> -----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}}
______________________________________________
R-help@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.
