1. 简单实现qq图
输入为一个vector,我们以a <- seq(1, 250, 1)做为示例数据
a <- seq(1, 250, 1)
利用qqnorm函数直接绘制出了如下正态检验qq图
qqnorm(a)
还可以进一步使用qqline命令在qq图上加上标准直线
qqline(a, col=2, lwd=2) # 设置为红色加粗
注:qqline的默认算法为向量a上四分位数和下四分位数对应两个点的连线
By default qqline draws a line through the first and third quartiles[1].
2. 了解基本原理,手动实现qq plot
Step 1: 首先我们算出vector中每一个数对应的百分位数
在向量a中,数字1对应的累积比例(即小于等于数字1的频率)为1/length(a) = 0.04,数字250对应的累积比例为250/length(a) = 100%