Displaying 1 result from an estimated 1 matches for "numberpati".
Did you mean:
numberpart
2010 Apr 11
1
Matrix is not symmetric under lme4
Dear all,
My code is presented as the following.
library(MASS)
library(rmutil)
library(repeated)
library(lme4)
library(arm)
#install.packages("Zelig", repos = "http://gking.harvard.edu")
library(Zelig)
rm(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)...