• R语言
  • 请教:R中如何做P-P图?

在R中,做Q-Q图有相关的函数:qqnorm , qqplot ,qqline,但是做P-P图有什么函数呢?请高手指教,谢谢!

如果你知道什么是PP图,就不会问这个问题了,因为画它实在是太简单了。例:

x=rnorm(30)
plot(seq(0,1,length=length(x)),pnorm(sort(x)))

若无异议,本帖将正式改编为一道习题(<https://bookdown.org/xiangyun/msg/gallery.html>)。

谢谢老大指教!如果作为习题,能否弄一个较完善的函数可以完成画PP图和分布拟合线,更多一点可以绘制分布拟合的置信区间线,下面是利用Minitab完成的正态概率图:



[attachment=339]





数据是:



0.328

0.346

0.372

0.33

0.356

0.42

0.346

0.324

0.306

0.322

0.25

0.328

0.368

0.288

0.374

0.382

0.332

0.288

0.304

0.37

0.352

0.318

0.366

0.316

0.342

0.328

0.328

0.296

0.322

0.35
[quote]引用第2楼谢益辉于2010-02-22 14:54发表的  :

如果你知道什么是PP图,就不会问这个问题了,因为画它实在是太简单了。例:
x=rnorm(30)<br />
plot(seq(0,1,length=length(x)),pnorm(sort(x)))




若无异议,本帖将正式改编为一道习题(http://yihui.name/cn/publication/)。

[/quote]

弱问个题外话,谢老大这本画图的书好久可以出来,这样子是不是R Graphics这本书就可以在中国没市场了,嘿嘿
你这个概率图倒是把我画糊涂了,概率图难道不是理论概率对实际概率的散点图么?怎么变成概率对数据的散点图了?还有那直线和置信区间分别是什么意思?……
R Graphics在国内有卖么?貌似我都没看到……我这书写完怕是遥遥无期,慢慢磨吧,反正放网上了,觉得有用就瞅瞅,没用可以拍砖



我一看你说“好久”就晓得是四川/重庆人,哈哈
下面是Minitab中关于概率图的说明:



与轴为线性的经验累积分布函数图不同,概率图的 Y 轴(有时是 X 轴)将变换,以使拟合分布线(中心的蓝线)形成直线(请参见获取概率图点的方法)。



Minitab 还显示拟合分布的大约 95% 的置信区间(弯曲的蓝线)。这些置信区间是以逐点方式给出的,表示它们是对拟合分布上的每个点单独进行计算的,从而无法控制全族误差。因此,如果使用它们对每个样本估计多个参数,则发生类型 I 错误的几率比所选的 a 水平大。



通常,置信区间之外的点出现在尾部。在图的下半部分,置信带右侧的点指示左尾的数据比基于拟合分布所预期的要少。在上半部分,置信带右侧的点指示右尾的数据比所预期的多。上述示例中的数据为正(右)偏斜,意味着左尾的观测值较比所预期的少,右尾的观测值比所预期的多。



拟合分布线

用于确定样本数据服从特定分布的程度。Minitab 使用从样本导出的参数估计值或用户输入的历史值生成拟合分布线。一般情况下,这些分布线会与实际数据重叠,因此您可以将经验数据与假设分布直接进行比较。拟合分布线可以出现在直方图、概率图和经验累积分布函数图中。
在谢老师网站上无法留言,学生在此感谢谢老师。只是第七章和第八章还未写好吗?
呃,让我想起郭德纲常说的一句话,“这是人话吗?”我看了半天实在没看懂…… [s:12]