Re: [R] Problem with coxme

From: Terry Therneau <>
Date: Thu 06 Jul 2006 - 01:11:23 EST

Date: Wed, 5 Jul 2006 09:09:14 -0500 (CDT) From: Terry Therneau <> Subject: RE: Problem with coxme
Cc:,, Content-MD5: BXKVsHtW/1I9mIUqrXBU0g==

  The original question involved a strange error message from coxme when trying to fit a random slopes model:

     coxme(Surv(time, status) ~ rx+x2, data=rats2, random=~ (1+x2)|litter)
   	a. coxme does not yet support random slopes models
  	b. it should recognize the request and print a sensible error, instead
 it just dies further down in the routine after it has gotton thoroughly  confused.

        c. the actual error message is interesting, as an example of just how  odd the error message from a confused routine can be; but has no other  use.  

   Thanks to Jing Zhao who forwarded the message on to me, as I prefer Splus (blasphemy!) and don't monitor the r-help list, and also to Spencer Graves and others who tried to further identify the bug.

		Terry Therneau
		Mayo Clinic
		(author of coxme)

  Random slopes are on my "to do" list, and in fact most of the infrastructure for them is already built into the underlying C code. Unfortunately it's been on the list for 3 years now and may stay there until a data set crosses my own desk that requires it (or there is sufficient clamor from without). In my defense, I was first sidetracked by trying to understand and support pdMat structures, something I've since abandoned all hope of. They are great examples of object-oriented programming: elegant, short, code resuse, and incomprehensible (see mailing list PLEASE do read the posting guide! Received on Thu Jul 06 01:16:12 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Thu 06 Jul 2006 - 02:18:31 EST.

Mailing list information is available at Please read the posting guide before posting to the list.