2016-12-05 4 views
-2

下記のR出力から個々の統計を抽出したいと思います。実行する方法?GOFロジスティックモデルR

# R code 
library(pscl) 
full <- glm(RESPONSE ~ as.factor(SEX) + TOTCHOL + AGE + ..., 
     data=training, family=binomial(link="logit")) 
summary(full) 
# McFadden's GOF pseudo R-squared, 0.2 to 0.4 indicates good model fit 
McFadden=round(pR2(full),digits=4)['McFadden']; McFadden 

出力:

library(pscl) 
round(pR2(reduced),digits=4) 
     llh llhNull  G2 McFadden  r2ML r2CU 
-165.6085 -191.1523 51.0875 0.1336 0.0663 0.1651 

がエラー:

pR2$McFadden 
Error in pR2$McFadden : object of type 'closure' is not subsettable 

答えて

0

私はあなたが私たちのすべてのコードを与えるとreducedが回帰コールの値が割り当てられたことがなかったことを推測していますリストではなくベクトルであるように見えるので、この未テストコードのようなものを使用します(

)。
round(pR2(reduced),digits=4)['McFadden'] 
+0

[McFadden]というこの種の構文は一度も使用していません。どうもありがとうございました。それは働く... MM –

関連する問題