Re: [R] Maximum likelihood estimation in R with censored Data

From: Vincent Goulet <vincent.goulet_at_act.ulaval.ca>
Date: Fri, 13 Jun 2008 14:26:46 -0400

Le ven. 13 juin à 13:55, Ben Bolker a écrit :

> Bluder Olivia <olivia.bluder <at> k-ai.at> writes:
>
>>
>> Hello,
>>
>> I'm trying to calculate the Maximum likelihood estimators for a
>> dataset
>> which contains censored data.
>>
>> I started by using the function "nlm", but isn't there a separate
>> method
>> for doing this for e.g. the "weibull" and the "log-normal"
>> distribution?
>>
>> Thanks,
>>
>> Olivia
>
> This is not *quite* enough detail about what you
> want to do. Can you (as the posting guide suggests!)
> give us a small example of what you want to do? You may be able
> to do this via the survreg() command in the survival
> package, or you may want to do it yourself by constructing
> a log-likelihood function with dweibull() for uncensored
> data and pweibull() for censored data [or dlnorm/plnorm].

If you want to go the second route, function coverage() in package actuar will build the censored density function for you. You can then feed this function to fitdistr() just like for "usual" ML estimation.

HTH Vincent

>
>
> Ben Bolker
>
> ______________________________________________
> 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.



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. Received on Fri 13 Jun 2008 - 18:35:47 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 13 Jun 2008 - 19:30:39 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive