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。