線形モデルから F 統計の p 値にアクセスするために R で使用している便利な関数を見つけました
lmp<-function(modelobject){
if(class(modelobject !="lm" stop ("Not an object of class 'lm ")
f<-(modelobject)$fstatistic
p<-pf(f[1],f[2],f[3],lower.tail=F)
attributes(p)<-NULL
return(p)
}
ただし、この関数を使用するために必要なデータは、クラス lm のオブジェクトではなく、12 個の要素を持つリストです。リストの各要素は、lm から生成された要約です。リストデータの各要約には、次を使用してアクセスできます
data[[1]]$fstatistic
リスト内の各要素をループして F 統計の p 値を返すように関数を変更したいので、関数を変更して for ループ lmpmod<-function(modelobject){
for (i in 1:12) {
f<-(modelobject)[[i]]$fstatistic
p<-pf(f[1], f[2], f[3], lower.tail=F)
attributes(p)<-NULL
return(p)
}
}
しかし、出力はそれ自体を上書きし続けます。リスト内の各要素にアクセスして F 統計の p 値を計算できるように、この関数を変更する方法を誰か教えてもらえますか?