私は多くの回帰を実行していますが、特定の変数の係数と p 値への影響のみに関心があります。したがって、私のスクリプトでは、glm サマリーから p 値を抽出できるようにしたいと考えています (係数自体を取得するのは簡単です)。私が知っている p 値を表示する唯一の方法は、summary(myReg) を使用することです。他の方法はありますか？
例えば：
fit <- glm(y ~ x1 + x2, myData)
x1Coeff <- fit$coefficients[2] # only returns coefficient, of course
x1pValue <- ???
行列として扱ってみましfit$coefficientsたが、p 値を単純に抽出することはまだできません。
これを行うことは可能ですか？
ありがとう！
番号の代わりに名前を直接入力できます
coef(summary(fit))[,'Pr(>|z|)']
係数の要約から利用可能な他のもの：
Estimate Std. Error z value Pr(>|z|)
summary私は過去にこの手法を使用して、適合モデル オブジェクトから、またはモデル オブジェクトから予測データを引き出しました。
coef(summary(m))[grepl("var_i_want$",row.names(coef(summary(m)))), 4]
これにより、データを取得する変数を簡単に編集できます。
または、指摘されているように @Ben、使用matchor %in%、よりややクリーンgrepl：
coef(summary(m))[row.names(coef(summary(m))) %in% "var_i_want" , 4]