Displaying 1 result from an estimated 1 matches for "mat3_1".
Did you mean:
mat1_1
2008 Dec 28
1
Random coefficients model with a covariate: coxme function
...s?
ugroup<-paste(rep(1:32, each=3), rep(c(0:1,"z"), 32), sep='/') #unique groups
mat1_1<-bdsmatrix(rep(c(1,1,1,1,1,1), 32),
blocksize=rep(3,32), dimnames=list(ugroup, ugroup))
mat2_1<-bdsmatrix(rep(c(0,0,0,1,0,0), 32),
blocksize=rep(3,32), dimnames=list(ugroup, ugroup))
mat3_1<-bdsmatrix(rep(c(0,0,0,0,0,1), 32),
blocksize=rep(3,32), dimnames=list(ugroup, ugroup))
group<-paste(data1$centers, data1$treat, data1$covar, sep="/")
fit3<- coxme(Surv(y, uncens)~ treat + covar,
data1, random=~1|group1, varlist=list(mat1_1,
mat2_1, mat3_1), rescale=F, pdcheck...