如何绘制qq图_python画图

如何绘制qq图_python画图Q-Q图主要可以用来回答这些问题:两组数据是否来自同一分布PS:当然也可以用KS检验,利用python中scipy.stats.ks_2samp函数可以获得差值KSstatistic和P值从而实现判断。两组数据的尺度范围是否一致两组数据是否有类似的分布形状前面两个问题可以用样本数据集在Q-Q图上的点与参考线的距离判断;而后者则是用点的拟合线的斜率判断。用Q-Q图来分析分布的好处都有啥?(谁说对了…

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

Q-Q图主要可以用来回答这些问题:

两组数据是否来自同一分布

PS:当然也可以用KS检验,利用python中scipy.stats.ks_2samp函数可以获得差值KS statistic和P值从而实现判断。

两组数据的尺度范围是否一致

两组数据是否有类似的分布形状

前面两个问题可以用样本数据集在Q-Q图上的点与参考线的距离判断;而后者则是用点的拟合线的斜率判断。

用Q-Q图来分析分布的好处都有啥?(谁说对了就给他)

两组数据集的大小可以不同

可以回答上面的后两个问题,这是更深入的数据分布层面的信息。

那么,Q-Q图要怎么画呢?

将其中一组数据作为参考,另一组数据作为样本。样本数据每个值在样本数据集中的百分位数(percentile)作为其在Q-Q图上的横坐标值,而该值放到参考数据集中时的百分位数作为其在Q-Q图上的纵坐标。一般我们会在Q-Q图上做一条45度的参考线。如果两组数据来自同一分布,那么样本数据集的点应该都落在参考线附近;反之如果距离越远,这说明这两组数据很可能来自不同的分布。

python中利用scipy.stats.percentileofscore函数可以轻松计算上诉所需的百分位数;而利用numpy.polyfit函数和sklearn.linear_model.LinearRegression类可以用来拟合样本点的回归曲线

from scipy.stats import percentileofscore

from sklearn.linear_model import LinearRegression

import pandas as pd

import matplotlib.pyplot as plt

# df_samp, df_clu are two dataframes with input data set

ref = np.asarray(df_clu)

samp = np.asarray(df_samp)

ref_id = df_clu.columns

samp_id = df_samp.columns

# theoretical quantiles

samp_pct_x = np.asarray([percentileofscore(ref, x) for x in samp])

# sample quantiles

samp_pct_y = np.asarray([percentileofscore(samp, x) for x in samp])

# estimated linear regression model

p = np.polyfit(samp_pct_x, samp_pct_y, 1)

regr = LinearRegression()

model_x = samp_pct_x.reshape(len(samp_pct_x), 1)

model_y = samp_pct_y.reshape(len(samp_pct_y), 1)

regr.fit(model_x, model_y)

r2 = regr.score(model_x, model_y)

# get fit regression line

if p[1] > 0:

p_function = “y= %s x + %s, r-square = %s” %(str(p[0]), str(p[1]), str(r2))

elif p[1] < 0:

p_function = “y= %s x – %s, r-square = %s” %(str(p[0]), str(-p[1]), str(r2))

else:

p_function = “y= %s x, r-square = %s” %(str(p[0]), str(r2))

print “The fitted linear regression model in Q-Q plot using data from enterprises %s and cluster %s is %s” %(str(samp_id), str(ref_id), p_function)

# plot q-q plot

x_ticks = np.arange(0, 100, 20)

y_ticks = np.arange(0, 100, 20)

plt.scatter(x=samp_pct_x, y=samp_pct_y, color=’blue’)

plt.xlim((0, 100))

plt.ylim((0, 100))

# add fit regression line

plt.plot(samp_pct_x, regr.predict(model_x), color=’red’, linewidth=2)

# add 45-degree reference line

plt.plot([0, 100], [0, 100], linewidth=2)

plt.text(10, 70, p_function)

plt.xticks(x_ticks, x_ticks)

plt.yticks(y_ticks, y_ticks)

plt.xlabel(‘cluster quantiles – id: %s’ %str(ref_id))

plt.ylabel(‘sample quantiles – id: %s’ %str(samp_id))

plt.title(‘%s VS %s Q-Q plot’ %(str(ref_id), str(samp_id)))

plt.show()

9dd1645bfc79b5316789a798fd38df51.png

效果如上图所示,在本例中所用的样本数据在左下稀疏,在右上集中,且整体往上偏移,说明其分布应该与参考数据是不一样的(分布形状不同),用KS检验得到ks-statistic: 0.171464; p_value: 0.000000也验证了这一点;但是其斜率在约为1,且整体上偏的幅度不大,说明这两组数据的尺度是接近的。

PS: 这里的方法适用于不知道数据分布的情况。如果想检验数据是否符合某种已知的分布,例如正态分布请出门左转用scipy.stats.probplot函数。

参考:

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/169330.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • ubuntu16.04安装cuda10.2_opencv cuda

    ubuntu16.04安装cuda10.2_opencv cuda由于实验需要,在实验室电脑上搭建深度学习Caffee框架。一共花了两天的时间,其中遇到了不少的问题,记录一下。Caffee在配置上相对来说比较麻烦,需要前期安装的东西比较多,逐一介绍。CUDA:NVIDIA系列显卡支持的GPU编程框架,其实如果本身电脑是AMD的显卡,不用装也可用Caffee,只是速度会比较慢。所以最好有一块像样的显卡,最后我就败在这个上面了。。。MKA或是OpenAtlas

    2022年10月30日
  • httprunner(2)下载安装「建议收藏」

    httprunner(2)下载安装「建议收藏」环境要求HttpRunner是一个基于Python开发的测试框架,可以运行在macOS、Linux、Windows系统平台上。这里使用macOS系统进行演示对于python版本要求:py

  • 通过pycharm使用git[图文详解][通俗易懂]

    通过pycharm使用git[图文详解][通俗易懂]#如果之前已经正确clone后,可以通过gitpull进行同步代码,可指定分支gitpull#是pull当前分支的最新代码gitpullorigindev#表示fetchorigin的dev分支到当前目录,并与当前分支合并,相当于gitfetchorigindev&&gitmergeorigindev,取下来再合并gitpullorigindev:master#表示从远端origindev取下代码,并与本地的mas…

  • 阿里云矢量图标库用法_阿里矢量图库图标

    阿里云矢量图标库用法_阿里矢量图库图标1.登陆http://www.iconfont.cn/例如我需要购物车的图标2.按下回车后会出现一堆的购物车图标3.喜欢那个一个就鼠标移动到图标上去,加入购物车4点击加入购物车后,就会出现在右

  • bootoption没有启动项_javacontinue的用法

    bootoption没有启动项_javacontinue的用法一、现象从fetch说起,用fetch构造一个POST请求。1fetch(‘http://127.0.0.1:8000/api/login’,{2method:”POST”,3headers:({4’Content-Type’:’application/x-www-form-urlencoded’5}),6body:”name=”+name…

    2022年10月30日
  • 让你的网页更精彩 – Javascript 调用MSAgent

    让你的网页更精彩 – Javascript 调用MSAgent<!–让你的网页更精彩–MicrosoftAgent在网页上的简单应用(本文假设您使用WindowsXP或Windows2000操作系统)不知在你漫游互联网时可曾在他开某个网页时看到一个小巫师,蓝色的袍子上满是金黄的星星和月亮十分可爱。他会向你问好,给你介绍这个网站。你一定奇怪,那个巫师是怎么做出来的。其实他并不是网页实现的而是微软的一…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号