Dear all, I am using "coxme" function in Kinship library to fit random treatment effect nested within centre. I got 3 treatments (0,1,2) and 3 centres. I used following commands, but got an error.
> ugroup=paste(rep(1:3,each=3),rep(0:2,3),sep='/') > mat1=bdsmatrix(rep(c(1,1,1,1,1,1,1,1,1),3),blocksize=rep(3,3),dimnames=list(ugroup,ugroup)) > mat2=bdsmatrix(rep(c(0,0,0,0,0,0,0,0,1),3),blocksize=rep(3,3),dimnames=list(ugroup,ugroup)) > group=paste(dat1$centre,dat1$treat,sep='/') > coxme(Surv(time,status) ~ as.factor(treat), data=dat1,random= > ~1|group,varlist=list(mat1,mat2),rescale=F,pdcheck=FALSE) Error in coxme.fit(X, Y, strats, offset, init, control, weights = weights, : Random effects variance is not spd Could anyone help me correcting this error? Many thanks in advance. Ruwanthi ____________________________________________________________________________________ Looking for last minute shopping deals? ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.