Re: [R] mixed-effects models for left-censored data?

About this list Date view Thread view Subject view Author view Attachment view

From: A.J. Rossini (
Date: Thu 12 Jun 2003 - 11:53:39 EST

Message-id: <>

Thomas Lumley <> writes:

> On Wed, 11 Jun 2003, Remko Duursma wrote:
> For a random intercept model you could use survreg() and frailty() in the
> survival package.
> In general the random effects tobit model will be quite hard to fit,
> involving a numerical integration whose dimension is the number of random
> effects. Some sort of EM algorithm might work.

One huge catch with that approach is heterscedasticity, which seems to
pop its head up all too often with limit-of-detection assay data.

> There is a paper by Pettit in Biometrics some time ago on censored linear
> mixed models -- I don't have the reference with me.

There is also a paper by a fellow named Jim Hughes, in Biometrics
(late 90s?), on this exact topic -- he used single imputation, whereas
he mentioned later (private communication) that a multiple imputation
approach would be better. The S-PLUS code (it isn't pretty) is
somewhere on his WWW page, buried deep in the U Washington
Biostatistcs WWW site.

At least it used to be.


A.J. Rossini  /  /
Biomedical/Health Informatics and Biostatistics, University of Washington.
Biostatistics, HVTN/SCHARP, Fred Hutchinson Cancer Research Center.
FHCRC: 206-667-7025 (fax=4812)|Voicemail is pretty sketchy/use Email 

CONFIDENTIALITY NOTICE: This e-mail message and any attachments ... {{dropped}}

______________________________________________ mailing list

About this list Date view Thread view Subject view Author view Attachment view

This archive was generated by hypermail 2.1.3 : Tue 01 Jul 2003 - 09:11:53 EST