[R] Programmation pour MLE
Dirk Eddelbuettel
edd at debian.org
Thu Jul 22 16:26:34 CEST 2004
On Thu, Jul 22, 2004 at 06:48:58AM +0200, Ahmed Elhabti wrote:
> Je veux cherch? l?estimateur de vraisemblance maximal (MLE)d?une fonction ? 3 param?tre inconue ?tant donn? une ?chantillon de taille 50 (les observations des valeurs de x) alors comment je peux proc?d?
>
> La fonction de densit? est d?finie par :
>
>
>
> f(x)= 1/3(g(a1)+g(a2)+g(a3))
>
>
>
> avec g(ai)=(exp(ai)*ai^x)/x! pour i=1,2,3.
The list language is English, so please pardon my reply for using it.
This should be reasonably straightforward with optim() and a simple function
that sums the terms of the log-likelihood, see help("optim").
There is also the mle function / class in stats4 which may help you as it
wraps around optim(), see help("mle", package="stats4").
> Je vous remercie beaucoup.
De rien.
Dirk
--
Those are my principles, and if you don't like them... well, I have others.
-- Groucho Marx
More information about the R-help
mailing list