Dear list-eners, I run into the following problem when I want to get contrasts from a coxme model using the lsmeans package: A call to lsmeans on the coxme model throws the following error:
Error in if (adjustSigma && object$method == "ML") V = V * object$dims$N/(object$dims$N - : missing value where TRUE/FALSE needed I give an example: library(coxme) library(lsmeans) fm <- coxme(Surv(y, uncens) ~ trt + (trt | center) + strata(center), data=eortc) summary(fm) lsmeans(fm, ~ trt) traceback points to the internal function lsmeans:::lsm.basis.lme() Any ideas? Thanks in advance. -m [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.