私は、独裁者のゲームで寄付を説明する回帰モデルに合っています。私は多くの変数を持っているので、 'for'ループを使ってプロセスを自動化したいと思います。今のところ私は単変量モデルから始まります。 フィッティングfit[1:24]
を印刷/要約すると、切片と係数だけが表示されます。 p値が格納されていないようですね?forループにはp値が格納されていませんか? (in R)
sumres[i]=summary(fit[i])
pf(sumres[i]$fstatistic[1L], sumres[i]$fstatistic[2L],sumres[i]$fstatistic[3L], lower.tail = FALSE)
と
pVal[i] <- (fit[i])$coefficients[,4]
いつもはなって終わる:
predictor<-0
fit<-0
dictatorgame<-mydata$dictatorgame
sumres<-0
pVal<-0
for(i in 1:24) #24 predictor variables stored in column 1-24 in mydata
{
predictor<-mydata[i]
unlist(predictor)
fit[i]<-lm(dictatorgame~unlist(predictor))
}
は、私はそれらの両方はオブジェクトがアトミックであることを考えるように見せかけ、私はSOにここで見つける二つの異なる解決策を試してみましたエラーメッセージ$ operator is invalid for atomic vectors
。
'predictor <-mydata [i]'が間違っているようです。 'mydata'がdata.frameなら' predictor <-mydata [、i] 'でなければなりません。 – Alex
ええ、それは私がやったことのよりクリーンな解決策です( 'unlist(predictor)')。私のコードでそれを変更した後、私はまだp値を見ません(インターセプトと係数は正しいと思われます)。 – YaeVo