Displaying 1 result from an estimated 1 matches for "fit_full".
Did you mean:
bt_full
2007 Jun 18
1
Question about lmer
...n(ff_count, fixed_factors, rf_count, random_factors, y_values)
{
SA<-matrix(as.array(c(fixed_factors, random_factors)), ncol=3)
data<-as.data.frame(SA)
y<-as.array(y_values)
dd<-data.frame(SA)
for(i in 1:(ff_count+rf_count)){
dd[,i]<-as.factor(data[,i])
}
fit_full=lmer(y~dd[,1]+dd[,2]+(1|dd[,3]),method="ML")
fit_full
}
A<-c(0,0,0,0,1,1,1,1,0,0,0,0,1,1,1,1,0,0,0,0,1,1,1,1,0,0,0,0,1,1,1,1)
B<-c(0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1)
C<-c(0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1)
Y<-c...