分散の信頼区間を 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
コメント