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" <hect...@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 " <c...@datanalytics.com> > *Para:* "Hector G�mez Fuerte" <hect...@gmx.es> > *CC:* "Lista R" <r-help-es@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" > <hect...@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@r-project.org > > https://stat.ethz.ch/mailman/listinfo/r-help-es > > > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > > [[alternative HTML version deleted]]
_______________________________________________ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es
