R ppoints 概率图的坐标

说明

生成概率点序列 (1:m - a)/(m + (1-a)-a)，其中 m 是 n (如果是 length(n)==1 )或 length(n) 。

用法

ppoints(n, a = if(n <= 10) 3/8 else 1/2)

参数

n	生成的点数或观察向量。
a	要使用的偏移分数；通常在 (0,1) 中。

细节

如果 0 < a < 1 ，则结果值在 (0,1) 范围内(不包括边界)。在任何情况下，生成的序列在 [0,1] 中都是对称的，即 p + rev(p) == 1 。

ppoints() 在 qqplot 和 qqnorm 中使用来生成用于评估逆分布的概率集。

a 的选择遵循 Becker 等人 (1988) 中同名函数的文档，并且似乎受到 Blom (1958) 关于期望正态顺序统计的近似值的结果的启发(另请参阅 quantile )。

将a分别取为1/2、0、1、1/3和3/8，可以获得连续样本分位数类型5至9(参见quantile)的概率点。

例子

ppoints(4) # the same as  ppoints(1:4)
ppoints(10)
ppoints(10, a = 1/2)

## Visualize including the fractions :
require(graphics)
p.ppoints <- function(n, ..., add = FALSE, col = par("col")) {
  pn <- ppoints(n, ...)
  if(add)
      points(pn, pn, col = col)
  else {
      tit <- match.call(); tit[[1]] <- quote(ppoints)
      plot(pn,pn, main = deparse(tit), col=col,
           xlim = 0:1, ylim = 0:1, xaxs = "i", yaxs = "i")
      abline(0, 1, col = adjustcolor(1, 1/4), lty = 3)
  }
  if(!add && requireNamespace("MASS", quietly = TRUE))
    text(pn, pn, as.character(MASS::fractions(pn)),
         adj = c(0,0)-1/4, cex = 3/4, xpd = NA, col=col)
  abline(h = pn, v = pn, col = adjustcolor(col, 1/2), lty = 2, lwd = 1/2)
}

p.ppoints(4)
p.ppoints(10)
p.ppoints(10, a = 1/2)
p.ppoints(21)
p.ppoints(8) ; p.ppoints(8, a = 1/2, add=TRUE, col="tomato")

参考

Becker, R. A., Chambers, J. M. and Wilks, A. R. (1988) The New S Language. Wadsworth & Brooks/Cole.

Blom, G. (1958) Statistical Estimates and Transformed Beta Variables. Wiley

也可以看看

qqplot、qqnorm。