我有一個關於R. std中student-t分佈的std和t的定義的問題,它帶有rugamo包,而t來自stats包。繪圖時std和t分佈之間的差異
qqplot(qstd(c(1:1000)/1001, nu=5),qt(c(1:1000)/1001,df=5))
abline(0,1)
很明顯看到兩個定義是不同的。任何人都可以告訴我爲什麼會有這種差異,哪個功能給出了正確的值?
我有一個關於R. std中student-t分佈的std和t的定義的問題,它帶有rugamo包,而t來自stats包。繪圖時std和t分佈之間的差異
qqplot(qstd(c(1:1000)/1001, nu=5),qt(c(1:1000)/1001,df=5))
abline(0,1)
很明顯看到兩個定義是不同的。任何人都可以告訴我爲什麼會有這種差異,哪個功能給出了正確的值?
如果你打開qstd
功能(這實際上是從fGarch),你會看到,它的修改stats::qt
:
> qstd
function (p, mean = 0, sd = 1, nu = 5)
{
s = sqrt(nu/(nu - 2))
result = qt(p = p, df = nu) * sd/s + mean
result
}
<environment: namespace:fGarch>
那麼,是什麼它給你的是一個"non-standardized Student's t-distribution"作爲反對標準化分佈可從統計。如果s=1
和自由度達到無窮大,它們將產生相同的結果。
難道不是簡單的抽樣問題嗎?如果你將你的樣本增加到100萬個點,看起來如何? – 2015-02-24 08:52:27
謝謝你的評論。不幸的是,它不是一個抽樣問題,qstd和qt函數僅僅爲t分佈返回不同的分位數值。 – Rosenlund 2015-02-24 09:05:46
'qstd'是包'fGarch'的函數。我現在還沒有意義,但'qstd(c(1:1000)/ 1001,nu = 5),qt(c(1:1000)/ 1001,df = 5)'是一個常數。 – 2015-02-24 09:22:33