[R] Left censored responses in mixed effects models

Giovanni Parrinello parrinel at med.unibs.it
Tue May 13 09:37:36 CEST 2008


Dear Bert,
a solution is the 'package'
censre3 by  Hughes JP
Reference: Hughes JP: Mixed effects models with censored data with 
application to HIV RNA levels.  Biometrics, 55:625-629, 1999.

Giovanni

Don MacQueen ha scritto:
> I assume you've looked at the NADA package(?) While I don't believe it 
> goes as far as dealing the mixed effects models, it might give you a 
> starting point, and possibly some additional references.
>
> -Don
>
> At 9:08 AM -0700 5/12/08, Bert Gunter wrote:
>> Dear R Fellow-Travellers:
>>
>> What is your recommended way of dealing with a left-censored response
>> (non-detects) in (linear Gaussian) mixed effects models?
>>
>> Specifics: Response is a numeric positive measurement (of volume, 
>> actually);
>> but when it falls below some unknown and slightly random value 
>> (depending on
>> how the sample is prepared and measured), it cannot be measured and is
>> recorded as 0.
>>
>> There is some statistical literature on this, but I was unable to find
>> anything that appeared to me to implement a strategy in any R 
>> package. If it
>> matters, I am less interested in inference than in removing possible 
>> bias in
>> estimation.
>>
>> Feel free to respond off-list if you feel that this would not be of 
>> general
>> interest.
>>
>> Cheers,
>>
>> Bert Gunter
>> Genentech
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide 
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>
>


-- 
dr. Giovanni Parrinello
Department of Biotecnologies
Medical Statistics Unit
University of Brescia
Viale Europa, 11 25123 Brescia
email: parrinel at med.unibs.it
Phone: +390303717528
Fax: +390303717488



More information about the R-help mailing list