Displaying 1 result from an estimated 1 matches for "b_estim".
Did you mean:
estim
2010 Apr 11
1
Matrix is not symmetric under lme4
...r==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/n-T1-T2)/(1+T3*n*(T1+T2))-w)^2)
}
MSE...