Re: [R] MLE for lambda of Poisson distribution using fitdistr

2009-10-28 Thread Terry Therneau
In general Poisson data consists of a pair of numbers (y,n), where y is the event count for the unit and n is the size of the unit. The Poisson MLE is sum(y)/sum(n). A general example is county level data where y is the number of events (rare cancer) and n is the county size. Two special cases a

Re: [R] MLE for lambda of Poisson distribution using fitdistr

2009-10-27 Thread Peter Ehlers
_ From: Peter Ehlers To: Kjetil Halvorsen Cc: r-h...@stat.math.ethz.ch; ankush...@yahoo.com Sent: Tue, October 27, 2009 10:15:31 AM Subject: Re: [R] MLE for lambda of Poisson distribution using fitdistr Kjetil Halvorsen wrote: What is wrong with using mean(x) to get the MLE of the poisson lam

Re: [R] MLE for lambda of Poisson distribution using fitdistr

2009-10-27 Thread Ankush Bhargava
he poisson lambda? Thanks again! Ankush From: Peter Ehlers To: Kjetil Halvorsen Sent: Tue, October 27, 2009 10:15:31 AM Subject: Re: [R] MLE for lambda of Poisson distribution using fitdistr Kjetil Halvorsen wrote: > What is wrong with using > mean(x)

Re: [R] MLE for lambda of Poisson distribution using fitdistr

2009-10-27 Thread Peter Ehlers
Kjetil Halvorsen wrote: What is wrong with using mean(x) to get the MLE of the poisson lambda? and mean(x)/length(x) to get its estimated variance. -Peter Ehlers Kjetil On Tue, Oct 27, 2009 at 9:17 AM, David Winsemius wrote: On Oct 26, 2009, at 11:25 PM, ankush...@yahoo.com wrote:

Re: [R] MLE for lambda of Poisson distribution using fitdistr

2009-10-27 Thread Kjetil Halvorsen
What is wrong with using mean(x) to get the MLE of the poisson lambda? Kjetil On Tue, Oct 27, 2009 at 9:17 AM, David Winsemius wrote: > > On Oct 26, 2009, at 11:25 PM, ankush...@yahoo.com wrote: > >> Hi, >> >> I am using the fitdistr of MASS to get the MLE for the lambda of a Poisson >> distribu

Re: [R] MLE for lambda of Poisson distribution using fitdistr

2009-10-27 Thread David Winsemius
On Oct 26, 2009, at 11:25 PM, ankush...@yahoo.com wrote: Hi, I am using the fitdistr of MASS to get the MLE for the lambda of a Poisson distribution. When i run the fitdistr command, i get an output that looks like - lambda 3.75 (0.03343) Couple of questions - 1. is the ML