Dear Bert,
you can also translate in nlme, as I'm trying to do, the approach of
Thiébaut and *Gadda( *Mixed models for longitudinal left-censored
repeated measures. Computer Methods and Programs in Biomedicine 74
<
http://www.informatik.uni-trier.de/%7Eley/db/journals/cmpb/cmpb74.html#ThiebautJ04>(3):
<javascript:void(0)>(2004)) written in nlmixed(SAS)
Giovanni
Bert Gunter ha scritto:
> 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@... 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@...
Phone: +390303717528
Fax: +390303717488
______________________________________________
R-help@... mailing list
https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide
http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
