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 (rossini@blindglobe.net)
Date: Thu 12 Jun 2003 - 11:53:39 EST


Message-id: <8765nchxi4.fsf@jeeves.blindglobe.net>

Thomas Lumley <tlumley@u.washington.edu> 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.

best,
-tony

-- 
A.J. Rossini  /  rossini@u.washington.edu  /  rossini@scharp.org
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}}

______________________________________________ R-help@stat.math.ethz.ch mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help


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