[R] fitting mixed models to censored data?

Bert Gunter gunter.berton at gene.com
Mon Apr 23 20:16:31 CEST 2007


AFAIK, this is subject area of active current research. Diggle, Heagerty,
Liang, and Zeger , 2002, (ANALYSIS OF LONGITUDINAL DATA) say on p.316: "An
emerging consensus is that analysis of data with potentially informative
dropouts necessarily involves assumptions which are difficult, or even
impossible, to check from the observed data."  This was ca 1994, I believe,
so I don't know whether this view is still held among experts (which I am
not). But if it is, you may do well to be careful of whatever SAS does even
if you do have to go running off to it.


Bert Gunter
Genentech Nonclinical Statistics

-----Original Message-----
From: r-help-bounces at stat.math.ethz.ch
[mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Douglas Grove
Sent: Monday, April 23, 2007 10:58 AM
To: r-help at stat.math.ethz.ch
Subject: [R] fitting mixed models to censored data?


I'm trying to figure out if there are any packages allowing
one to fit mixed models (or non-linear mixed models) to data
that includes censoring.

I've done some searching already on CRAN and through the mailing
list archives, but haven't discovered anything.  Since I may well
have done a poor job searching I thought I'd ask here prior to
giving up.

I understand that SAS's proc nlmixed can accomodate censoring
(though proc mixed apparently can't), so if I can't find 
something available in R, I'll have to break down and use
that.  Please, save me from having to use SAS!

Thanks much,

R-help at stat.math.ethz.ch mailing list
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

More information about the R-help mailing list