R で分散の信頼区間を計算する方法

分散の信頼区間を R で計算する方法の解説。

• 分散の信頼区間を計算するスクリプト
• まとめ
• 参考文献

分散の信頼区間を計算するスクリプト

Rスクリプトは以下の通り。

var.interval <- function(data, conf.level=0.95){
df <- length(data) - 1
chilower <- qchisq((1-conf.level)/2, df)
chiupper <- qchisq((1-conf.level)/2, df, lower.tail=FALSE)
v <- var(data)
c(df * v/chiupper, df * v/chilower)
}

var.interval()のカッコ内に、分散の信頼区間を計算したい変数名を書いて実行する。

95％信頼区間の場合は、それ以外何も指定しなくてOK。

> var.interval(Imatinib_Stop$Age)
[1] 202.3442 354.2129

90％信頼区間が計算したい場合は、conf.level=0.90をカッコ内に書き足す。

その際、変数名とは半角カンマで区切る。

> var.interval(Imatinib_Stop$Age, conf.level=0.90)
[1] 210.8776 337.2703

まとめ

分散の信頼区間を R で求めるスクリプトを解説した。

参考になれば。

参考文献

Using R to Find Confidence Intervals

http://pages.stat.wisc.edu/~yandell/st571/R/append7.pdf