1

Rヘルプhelp(quantile)で、あなたは見ることができます

Type 7 m = 1-p.  p[k] = (k - 1) / (n - 1).  In this case, p[k] =mode[F(x[k])]. 
This is used by S.

今、私は例を持っています：

w<-c(75,64,47.4,66.9,62.2,62.2,58.7,63.5,66.6,64,57,69,56.9,50,72.0)  
sort(w)  
[1] 47.4 50.0 56.9 57.0 58.7 62.2 62.2 63.5 64.0 64.0 66.6 66.9 69.0 72.0 75.0  
quantile(w)  
   0%   25%   50%   75%  100%   
47.40 57.85 63.50 66.75 75.00

タイプ7の数式を使用して結果を取得するにはどうすればよいですか？

4

1 に答える 1

2

答えが正しいかどうかを判断するのに問題があります。

> quantile(w, type=7)
   0%   25%   50%   75%  100% 
47.40 57.85 63.50 66.75 75.00

私の問題は、分位数のデフォルトがtype = 7であり、すでにその結果が得られていることです。Quantile.defaultのコードを見ると、type=7のセクションがあります。

        index <- 1 + (n - 1) * probs
        lo <- floor(index)
        hi <- ceiling(index)
        x <- sort(x, partial = unique(c(lo, hi)))
        qs <- x[lo]
        i <- which(index > lo)
        h <- (index - lo)[i]
        qs[i] <- (1 - h) * qs[i] + h * x[hi[i]]
于 2013-03-11T00:26:16.323 に答える