[R] ..lsmeans and coxme..

[Matthias Kuhn]

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]]