Displaying 1 result from an estimated 1 matches for "eblp_glimmix_quad".
2010 Apr 11
1
Matrix is not symmetric under lme4
...(list = ls())
beta0<-2.5
beta1<--0.3
sigs2<-0.5
I<-4
#numberpatients<-c(40,100,160,200,400,600)
numberpatients<-c(1000)
#numberpatients<-3
times<-1
MSE_eBLP_quad<-NULL
for (i3 in 1:length(numberpatients)) {
n<-numberpatients[i3]
b0<-NULL
b1<-NULL
sigma2<-NULL
eBLP_glimmix_quad<-NULL
for (kk in 1:times) {
s<-rnorm(I,0,sqrt(sigs2))
y<-NULL
u<-NULL
center<-NULL
for (i in 1:I) {
for (k in 1:n) {
u<-c(u,0)
center<-c(center,i)
lam<-exp(beta0+s[i])
obs<-rpois(1,lam)
y<-c(y,obs)
}
for (k in 1:n) {
u<-c(u,1)
center<-c(cente...