2017-08-06 27 views
0

私はベータ回帰モデル(パッケージ'betareg'を使用)とプロットを持っていますが、結果を報告するためにはR-平方とベータが必要です。私はlm方程式からBetaを見つけるためのlm.beta funtionとlm方程式からr-squaredを見つけるための概要(lm(DV~IV, data=mydata)$r.squaredを知っています。ベータ回帰モデルのこれらの値はどのようにして求められますか?RのBetaregモデルから、R-squaredとBetaの値はどうやって見つかりますか?

+0

Rsquaredを定義してください。ヒント:線形回帰以外は単純ではありません。 –

答えて

2

クラスbetaregのオブジェクトには、さまざまな抽出機能があります。表1を参照して​​を参照してください。簡単な例として

ReadingSkillsケーススタディ(5.1節)を検討:そのよう疑似などの特定の部分を抽出するには

summary(m) 
## Call: 
## betareg(formula = accuracy ~ iq * dyslexia | iq + dyslexia, data = ReadingSkills) 
## 
## Standardized weighted residuals 2: 
##  Min  1Q Median  3Q  Max 
## -2.3900 -0.6416 0.1572 0.8524 1.6446 
## 
## Coefficients (mean model with logit link): 
##    Estimate Std. Error z value Pr(>|z|)  
## (Intercept) 1.1232  0.1428 7.864 3.73e-15 *** 
## iq   0.4864  0.1331 3.653 0.000259 *** 
## dyslexia  -0.7416  0.1428 -5.195 2.04e-07 *** 
## iq:dyslexia -0.5813  0.1327 -4.381 1.18e-05 *** 
## 
## Phi coefficients (precision model with log link): 
##    Estimate Std. Error z value Pr(>|z|)  
## (Intercept) 3.3044  0.2227 14.835 < 2e-16 *** 
## iq   1.2291  0.2672 4.600 4.23e-06 *** 
## dyslexia  1.7466  0.2623 6.658 2.77e-11 *** 
## --- 
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
## 
## Type of estimator: ML (maximum likelihood) 
## Log-likelihood: 65.9 on 7 Df 
## Pseudo R-squared: 0.5756 
## Number of iterations: 25 (BFGS) + 1 (Fisher scoring) 

library("betareg") 
data("ReadingSkills", package = "betareg") 
m <- betareg(accuracy ~ iq * dyslexia | iq + dyslexia, data = ReadingSkills) 

を通常の概要は、あなたが探し情報を持っていますR-平方で、summary()の要素にアクセスできます。

summary(m)$pseudo.r.squared 
## 0.5756258 

または専用の方法があります。

coef(m) 
##  (Intercept)    iq   dyslexia  iq:dyslexia 
##   1.1232250   0.4863696  -0.7416450  -0.5812569 
## (phi)_(Intercept)   (phi)_iq (phi)_dyslexia 
##   3.3044312   1.2290731   1.7465642 
coef(m, model = "mean") 
## (Intercept)   iq dyslexia iq:dyslexia 
## 1.1232250 0.4863696 -0.7416450 -0.5812569 
coef(m, model = "precision") 
## (Intercept)   iq dyslexia 
## 3.304431 1.229073 1.746564 
関連する問題