Displaying 1 result from an estimated 1 matches for "estimate_sigma2".
2010 Apr 11
1
Matrix is not symmetric under lme4
...(i2 in 1:I){
yi1<-c(yi1,sum(y[u==0 & center==i2]))
yi2<-c(yi2,sum(y[u==1 & center==i2]))
}
##method2##############
fit2 <- glmer(y ~ u + (1|center), family=poisson, data=df,nAGQ=3,verbose =
TRUE, REML = FALSE )
#fit2 <- glmer(y ~ u + (1|center), family=poisson, data=df,nAGQ=3)
estimate_sigma2<-sigma.hat(fit2)$sigma$center
b_estimate<-fit2 at fixef
alpha1=b_estimate[1]
alpha2=b_estimate[1]+b_estimate[2]
T1=exp(alpha1+estimate_sigma2/2)
T2=exp(alpha2+estimate_sigma2/2)
T3=exp(estimate_sigma2)-1
T4=exp(0.5*estimate_sigma2)
w=exp(s)
eBLP_glimmix_quad[kk]=mean((T4+n*T3*T4*(yi1/n+yi2...