请教：R中如何做P-P图？

thinkyou

在R中，做Q-Q图有相关的函数：qqnorm ， qqplot ，qqline，但是做P-P图有什么函数呢？请高手指教，谢谢！

tryshy

ppplot.das {StatDA}

http://finzi.psych.upenn.edu/R/library/StatDA/html/ppplot.das.html

yihui

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

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

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

thinkyou

谢谢老大指教！如果作为习题，能否弄一个较完善的函数可以完成画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

skyer_117

[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这本书就可以在中国没市场了，嘿嘿

yihui

你这个概率图倒是把我画糊涂了，概率图难道不是理论概率对实际概率的散点图么？怎么变成概率对数据的散点图了？还有那直线和置信区间分别是什么意思？……

yihui

R Graphics在国内有卖么？貌似我都没看到……我这书写完怕是遥遥无期，慢慢磨吧，反正放网上了，觉得有用就瞅瞅，没用可以拍砖

我一看你说“好久”就晓得是四川/重庆人，哈哈

thinkyou

下面是Minitab中关于概率图的说明：

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

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

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

拟合分布线

用于确定样本数据服从特定分布的程度。Minitab 使用从样本导出的参数估计值或用户输入的历史值生成拟合分布线。一般情况下，这些分布线会与实际数据重叠，因此您可以将经验数据与假设分布直接进行比较。拟合分布线可以出现在直方图、概率图和经验累积分布函数图中。

ryusukekenji

在谢老师网站上无法留言，学生在此感谢谢老师。只是第七章和第八章还未写好吗？

yihui

呃，让我想起郭德纲常说的一句话，“这是人话吗？”我看了半天实在没看懂…… [s:12]

yihui

可以留言啊。