Message-ID: <A4E5A0B016B8CB41A485FC629B633CED3FF1B25853@GOLD.corp.lgc-group.com>
Date: 2012-07-05T12:39:19Z
From: S Ellison
Subject: Maximum Likelihood Estimation Poisson distribution mle {stats4}
In-Reply-To: <1341478091451-4635464.post@n4.nabble.com>
> -----Original Message-----
> > sample.mean<- sum(x*y)/sum(y)
> > sample.mean
> [1] 3.5433
>
> *This is the contradiction!! *
> Here I am getting the estimate as 3.5433(which is reasonable
> as most of the values are clustered around 3), but mle code
> gives the estimate 11.545(which may not be correct as this is
> out side the range 0:10)
>
> Why this contradiction?
>
>
Ermm.. the sample mean is mean(y) which is 11.545
I deduce that you and the help page have different views on what the sample y represents.
i also note that x does not figure at all in the help page's log-likelihood, suggesting that y is a simple set of counts, whereas you have interpreted x as the number of instances of each y. That appears not to be the case.
S
*******************************************************************
This email and any attachments are confidential. Any use...{{dropped:8}}