Python:多 QQ 绘图

2024-07-04

我是新人,通常来自 R。 我想创建一个包含多行的 QQ-Plot。

我有一个测试版分布式数据集 我想尝试不同的 beta 分布参数并在oneQQ-Plot 以便更好地进行比较。如果我尝试下面的代码,每个图都有相同的颜色,并且我得到 3 条 QQ 线。 有没有可能将这三个 QQ 情节合二为一?
我希望你能明白我的问题

import scipy.stats as stats
import numpy
x=numpy.random.beta(2, 3, size=100)
stats.probplot(x, dist=stats.beta, sparams=(2,3),plot=plt,fit=False)
stats.probplot(x, dist=stats.beta, sparams=(1,2),plot=plt,fit=False)
stats.probplot(x, dist=stats.beta, sparams=(1,4),plot=plt,fit=False)

亲切的问候


可以,然后呢stats.probplot让我有点困惑。该文档明确指出:

probplot生成概率图,不应将其与 Q-Q 或 P-P 图。

然而我能找到的所有资料都表明概率图指的是 Q-Q 图或 P-P 图。去搞清楚。

无论如何,就我而言,你所生成的isQ-Q 图。

在我看来,这个选项fit=False of stats.probplot被忽略,并且始终将回归线添加到数据中。

不管怎样,为了得到你想要的,我们可以显式创建一个 matplotlib 轴实例,并使用get_lines方法来删除不需要的回归线并更改标记颜色。

import scipy.stats as stats
import numpy as np
import matplotlib.pyplot as plt

plt.style.use('seaborn')

x = numpy.random.beta(2, 3, size=100)

fig, ax = plt.subplots(1, 1, figsize=(6, 4))
stats.probplot(x, dist=stats.beta, sparams=(2,3), plot=plt, fit=False)
stats.probplot(x, dist=stats.beta, sparams=(1,2), plot=plt, fit=False)
stats.probplot(x, dist=stats.beta, sparams=(1,4), plot=plt, fit=False)

# Remove the regression lines
ax.get_lines()[1].remove()
ax.get_lines()[2].remove()
ax.get_lines()[3].remove()

# Change colour of scatter
ax.get_lines()[0].set_markerfacecolor('C0')
ax.get_lines()[1].set_markerfacecolor('C1')
ax.get_lines()[2].set_markerfacecolor('C2')

# Add on y=x line
ax.plot([0, 1], [0, 1], c='C3')

这给了我以下内容,我认为这一次确实是您想要的:

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

Python:多 QQ 绘图 的相关文章

随机推荐

  • 如何在 Spring Boot 中使用 CommonsMultipartResolver

    我尝试在 Boot 中使用 CommonsMultipartResolver 将我的旧应用程序 WAR 转换为 Boot 现在它得到以下代码 Configuration public class TestConfig Bean public
  • SharePoint Web 部件与 ASP.NET Web 部件

    有人可以告诉我 SharePoint Web 部件和 ASP NET Web 部件之间的区别吗 我发现描述差异的信息很少 提前致谢 MSDN 在 ASP NET 2 0 Web 部件和 Windows SharePoint Services
  • MacOS 上的 Jenkins 发出“Docker:找不到命令”

    当使用 Pipeline 语法和 Docker 代理从 Jenkinsfile 运行作业时 管道失败并显示 Docker 找不到命令 我理解这意味着 1 Docker 未安装 2 Jenkins 未指向正确的 Docker 安装路径 我的情
  • 如何重试基于块的 URL 请求

    我正在使用 iOS7 的新 URL 请求方法获取数据 如下所示 NSMutableURLRequest request NSMutableURLRequest requestWithURL NSURL URLWithString self
  • 预期主表达式在“...”之前,C++ 编译错误

    SO 上有很多类似标题的帖子 但它们似乎是由各种语法错误触发的 而且我还没有看到一致的模式 using namespace std class A public A int a a a int a int main A x 3 A y 0
  • 预加载队列中的 mp3 文件,以避免播放队列中下一个文件时出现任何延迟

    我正在编写一个脚本 其中我正在播放多个 mp3 并且每个文件都在队列中 播放下一个 mp3 文件会出现轻微延迟 因为缓冲 加载文件需要时间 我如何缓冲队列中的下一个 mp3 文件 以便所有文件顺利运行而没有任何延迟 getData 1 fu
  • C++:你会选择 boost::date_time 还是 icu::date/time 库?

    我的应用程序需要自定义时间和日期设置功能 我检查了 ICU 和 boost date time 库 从完整性的角度来看 两者似乎都满足我的要求 我想知道两者之间是否有任何偏好 依据是什么 哪一个会在表现上得分 如果没有有关您的特定用例和环境
  • 如何在新窗口中获取dom元素?

    JavaScript 中的一个简单任务是打开一个新窗口并在其中写入 但我需要在一个dom元素中写入一个带有ID的div var novoForm window open somform html wFormx width 800 heigh
  • Windows 8 - 低延迟音频

    我正在考虑为即将推出的 Windows 8 开发一个应用程序 该应用程序需要低延迟音频录制和播放 我正在尝试找出操作系统是否支持这一点 而不是其他平台 http code google com p android issues detail
  • 替换 Apache POI XWPF 中的文本

    我刚刚发现 Apache POI 库对于使用 Java 编辑 Word 文件非常有用 具体来说 我想编辑一个DOCX使用 Apache POI 的 XWPF 类的文件 我发现没有合适的方法 文档可以执行此操作 有人可以分步骤解释一下如何替换
  • 上传大文件nginx + uwsgi

    堆栈 flask 0 10 uwsgi 1 4 5 nginx 1 2 3 我可以通过我的应用程序上传小文件 CONTENT LENGTH 无效 太大 跳过 nginx 日志没有显示任何有用的内容 我尝试了以下操作 但没有成功 nginx配
  • 如何在 Rust 中包装对使用 VarArgs 的 FFI 函数的调用?

    mexPrintf 就像printf 接受参数的可变参数列表 但我不知道用 Rust 包装它的最佳方法是什么 有一个可变泛型的 RFC https github com rust lang rfcs issues 376 但是今天我们能做什
  • 无法在单例对象中找到隐式值

    我有这个代码 trait Context implicit val e Encoder trait Encoder def write Unit println Test trait AsyncEncoders this Context g
  • getInstance() 不适用于实时数据库中除 us-central1 之外的其他位置

    我正在尝试将用户凭据保存到 Firebase 实时数据库 但是 当我执行该程序时 数据库没有更新 我已正确配置 Firebase 设置 因为身份验证和存储 均为 Firebase 正在运行 build gradle dependencies
  • 如何设置 UWP 应用的默认语言?

    我有2种语言 en US 和 ru Ru 当我将默认语言设置为 ru RU 时 出现错误 PRI257 0xdef00522 找到语言 ru ru en us 的资源 但未找到默认语言的资源 ru RU 改变 默认语言或使用默认语言限定资源
  • 如何在 .proto 文件中处理带有协议缓冲区的通用类型对象?

    我花了一些时间寻找一些替代方法来处理通用对象 我看到了与我类似的问题 但没有我想象的那么具体 协议缓冲区有多种我可以使用的标量类型 但它们大多是原始的 我希望我的消息是灵活的 并且能够有一个作为某种列表的字段 假设我的 proto 文件如下
  • XAML 页面中的 Bindable 值存在问题

    我有一个保龄球回顾示例代码 其中输入了 3 场比赛并总结在系列列中 我遇到的问题是系列列没有按照我的预期进行更新 显然 我做错了什么 但看不到我做错了什么 以下是支持该应用程序的代码
  • 如何覆盖 nltk 的 pos_tag 分配给文本的 POS 标签?

    我使用 nltk 中的 pos tag 来标记一组 未标记的 技术文档中的文本并获得良好的结果 但它总是将 authenticated 等单词标记为动词 而有时它可以用作形容词 换句话说 仅仅改变标签并不是每次都有效 是否有一个好方法来覆盖
  • 嵌入式阵列文档中的 Morphia Mongodb 更新失败

    我是 Morphia 的新手 正在尝试更新嵌入式文档中的字段 这是当前的结构 class A List b BList class B String field 所以我的结构在 MongoDb 中如下所示 id ObjectId 5bab8
  • Python:多 QQ 绘图

    我是新人 通常来自 R 我想创建一个包含多行的 QQ Plot 我有一个测试版分布式数据集 我想尝试不同的 beta 分布参数并在oneQQ Plot 以便更好地进行比较 如果我尝试下面的代码 每个图都有相同的颜色 并且我得到 3 条 QQ