[R-es] Ajuste con exponencial

[Jorge I Velez]

Hola Hector,

No soy experto, pero en

http://r.789695.n4.nabble.com/Fitting-weibull-exponential-and-lognormal-distributions-to-left-truncated-data-td869977.html
http://www.r-bloggers.com/r-help-follow-up-truncated-exponential/
http://www.jstatsoft.org/v16/c02/paper

hay algunas ideas.    Espero te sirvan.

Saludos,
Jorge.-



2015-01-28 22:27 GMT+11:00 "Hector Gómez Fuerte" <hector3 en gmx.es>:

> Saludos cordiales.
>
> Lamentablemente lo que dice Carlos no es correcto. Cuando la distribución
> exponencial la truncamos en un intervalo la constante de la función de
> densidad (para que integre 1) tiene una dependencia (complicada) del
> parámetro que multiplica al exponente, con lo cual la ecuación de
> verosimilitud no es nada sencila ni se puede resolver exactamente . La
> estimación maximo-verosímil requeriría de algoritmos numéricos y por tanto
> de software para su cálculo. Yo creo que esta no debe ser la mejor
> solución, y me sorprende no haber encontrado (puede que por mi torpeza)
> nada para ello en el R. De aquí la pregunta que hacia en este foro.
>
> Héctor Gómez
>
> *Enviar:* martes 27 de enero de 2015 a las 22:07
> *De:* "Carlos J. Gil Bellosta " <cgb en datanalytics.com>
> *Para:* "Hector Gómez Fuerte" <hector3 en gmx.es>
> *CC:* "Lista R" <r-help-es en r-project.org>
> *Asunto:* Re: [R-es] Ajuste con exponencial
> Hola, ¿qué tal?
>
> Creo que el ajuste (por máxima verosimilitud) de lambda es el inverso
> de la media de tus datos. Tu densidad en el intervalo de interés es
> como la de la exponencial (dividida por una constante de
> normalización). El logaritmo de la verosimitud es, por lo tanto, como
> el de la exponencial sin truncar más una constante.
>
> Luego la teoría habitual (de cómo el inverso de la media es el
> estimador por MV de lambda) aplica con cambios mínimos.
>
> Un saludo,
>
> Carlos J. Gil Bellosta
> http://www.datanalytics.com
>
> El día 27 de enero de 2015, 19:53, "Hector Gómez Fuerte"
> <hector3 en gmx.es> escribió:
> > Buenas tardes,
> > ¿cómo puedo con el R ajustar una distribución exponencial truancada (en
> el
> > intervalo [10,60]) a un vector de datos?
> > Muchas gracias.
> > Héctor Gómez
> >
> > _______________________________________________
> > R-help-es mailing list
> > R-help-es en r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-help-es
> >
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>

	[[alternative HTML version deleted]]