Re: [R-es] Ajuste con exponencial
Hola, Si tienes su expresión algebráica, podrías ajustarla con un ajuste no lineal nls(). nls() viene por defecto dentro del paquete stats y si no fuese posible hacerla converger, puedes utilizar otros paquetes (nlstools, nls2) que modifican nls() que utilizan otros algoritmos de convergencia. Saludos, Carlos Ortega www.qualityexcellence.es El 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 -- Saludos, Carlos Ortega www.qualityexcellence.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
Re: [R-es] Ajuste con exponencial
Hola Hector, buenos d�as:Hay un m�todo generalista para maximizar funciones en R que quiz�s te valga, prueba lo siguiente a ver qu� tal (por supuesto, si te fijas en la definici�n de la funci�n de verosimilitud, ves que la he cuadrado a mano al intervalo que trataslibrary(optimx)muestra - c(50, 20, 31, 40, 10)funcionExpTrun - function(l, x){n - length(x)f - rep(0, n) for (i in 1:n){ f[i]- -log(l*exp(-l*x[i])/(exp(-10*l)-exp(-60*l))) }sumaf - sum(f)return (sumaf)}resultado - optim(par = c(0.1), fn = funcionExpTrun, method = c(L-BFGS-B), lower = c(-Inf, 0), upper = c(Inf, Inf), x = muestra)El c�digo lo he sacado (con alguna adaptaci�n por mi parte) de:http://www.mat.uda.cl/jolivares/probabilidades/EMV.pdfEl par�metro al que converge ser�a: $par[1] 0.02356897Un saludo From: hect...@gmx.es To: c...@datanalytics.com; r-help-es@r-project.org Date: Wed, 28 Jan 2015 12:27:56 +0100 Subject: Re: [R-es] Ajuste con exponencial 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
Re: [R-es] Ajuste con exponencial
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
Re: [R-es] Ajuste con exponencial
Saludos cordiales. Lamentablemente lo que dice Carlos no es correcto. Cuando la distribucin exponencial la truncamos en un intervalo la constante de la funcin de densidad (para que integre 1) tiene una dependencia (complicada) del parmetro que multiplica al exponente, con lo cual la ecuacin de verosimilitud no es nada sencila ni se puede resolver exactamente . La estimacin maximo-verosmil requerira de algoritmos numricos y por tanto de software para su clculo. Yo creo que esta no debe ser la mejor solucin, 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. Hctor Gmez Enviar:martes 27 de enero de 2015 a las 22:07 De:Carlos J. Gil Bellosta c...@datanalytics.com Para:Hector Gmez 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 mxima verosimilitud) de lambda es el inverso de la media de tus datos. Tu densidad en el intervalo de inters es como la de la exponencial (dividida por una constante de normalizacin). El logaritmo de la verosimitud es, por lo tanto, como el de la exponencial sin truncar ms una constante. Luego la teora habitual (de cmo el inverso de la media es el estimador por MV de lambda) aplica con cambios mnimos. Un saludo, Carlos J. Gil Bellosta http://www.datanalytics.com El da 27 de enero de 2015, 19:53, Hector Gmez Fuerte hect...@gmx.es escribi: Buenas tardes, cmo puedo con el R ajustar una distribucin exponencial truancada (en el intervalo [10,60]) a un vector de datos? Muchas gracias. Hctor Gmez ___ 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
Re: [R-es] Ajuste con exponencial
Hola Hector, buenos d�as: Hay un m�todo generalista para maximizar funciones en R que quiz�s te valga, prueba lo siguiente a ver qu� tal (por supuesto, si te fijas en la definici�n de la funci�n de verosimilitud, ves que la he cuadrado a mano al intervalo que tratas library(optimx) muestra - c(50, 20, 31, 40, 10) funcionExpTrun - function(l, x){n - length(x)f - rep(0, n) for (i in 1:n){ f[i]- -log(l*exp(-l*x[i])/(exp(-10*l)-exp(-60*l))) }sumaf - sum(f)return (sumaf)} resultado - optim(par = c(0.1), fn = funcionExpTrun, method = c(L-BFGS-B), lower = c(-Inf, 0), upper = c(Inf, Inf), x = muestra) El c�digo lo he sacado (con alguna adaptaci�n por mi parte) de: http://www.mat.uda.cl/jolivares/probabilidades/EMV.pdf El par�metro al que converge ser�a: $par[1] 0.02356897 Un saludo From: hect...@gmx.es To: c...@datanalytics.com; r-help-es@r-project.org Date: Wed, 28 Jan 2015 12:27:56 +0100 Subject: Re: [R-es] Ajuste con exponencial 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
[R-es] Ajuste con exponencial
Buenas tardes, cmo puedo con el R ajustar una distribucin exponencial truancada (en el intervalo [10,60]) a un vector de datos? Muchas gracias. Hctor Gmez ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es
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