2015-08-18 5 views
1

Люди, Я пытаюсь напечатать таблицу {texreg} моделей lmer() {nlme} и lme() {lme4}, включая отклонения. Однако дисперсии существенно различаются между двумя моделями (на несколько порядков). Кажется, что lme() дисперсии являются квадратным корнем из lmer(). Какие из них правильные?texreg с объектами lmer и lme - отклонения отличаются

library(plm) 
library(lme4) 
library(nlme) 
library(texreg) 
data("Grunfeld", package="plm") 
reML0 <-lmer(inv ~ value + capital + (1|firm), data=Grunfeld) 
reML1 <- lme(inv ~ value + capital, data=Grunfeld, random=~1|firm) 
screenreg(list(reML0, reML1), digits=3, include.variance=TRUE) 

======================================================== 
          Model 1  Model 2  
-------------------------------------------------------- 
(Intercept)     -57.864 *  -57.864  
           (29.378)  (29.378) 
value       0.110 ***  0.110 *** 
           (0.011)  (0.011) 
capital       0.308 ***  0.308 *** 
           (0.017)  (0.017) 
-------------------------------------------------------- 
AIC       2205.851  2205.851  
Num. obs.      200   200   
Num. groups: firm    10      
Variance: firm.(Intercept) 7366.992     
Variance: Residual   2781.426     
Num. groups         10   
sigma           52.739  
sigma. RE          85.831  
======================================================== 
*** p < 0.001, ** p < 0.01, * p < 0.05 

ответ

2

Это кажется довольно ясно, что это дисперсия для модели 1 и сигмы для модели 2. Дисперсии квадратной сигмы отчетности, так что ответ на ваш вопрос «как»

 Смежные вопросы

  • Нет связанных вопросов^_^