Flores-de-Gracia, Eric
2012-May-22 12:51 UTC
[R-es] Plotear intervalos de confianza en nls model
Hola Disculpen nuevamente por esta nueva pregunta. Tengo un modelo no lineal fiteado a mis datos, me interesa presentar en el plot las líneas de intervalos de confianza del modelo. Se que con los errores estandar de cada parametro puedo calcular los límites superior e inferior pero no quiero hacer eso mecánicamente y preferiría si puedo implementar una función dentro del script que me permita realizar ese trabajo. Este es mi script: kv<-seq(0.55,7.02,0.1) lumv<-predict(kluxreg3,list(klux3=kv),interval="confidence",level=0.95) lines(kv,lumv,lty=19,col="black") Agradezco cualquier ayuda, Eric Flores De Gracia School of Biosciences University of Exeter, Cornwall Campus Penryn, Cornwall TR10 9EZ United Kingdom Mobilephone: +044 07578724705 http://biosciences.exeter.ac.uk/cec/staff/postgradresearch/index.php?web_id=Eric_Flores_de_Gracia&tab=profile
Eric, Quizas la respuesta de Gabor Grothendieck en [1] sea lo que necesitas. Saludos, Jorge.- [1] http://r.789695.n4.nabble.com/plotting-confidence-bands-from-predict-nls-td3505012.html 2012/5/22 Flores-de-Gracia, Eric <>> > Hola > > Disculpen nuevamente por esta nueva pregunta. > > Tengo un modelo no lineal fiteado a mis datos, me interesa presentar en el > plot las líneas de intervalos de confianza del modelo. Se que con los > errores estandar de cada parametro puedo calcular los límites superior e > inferior pero no quiero hacer eso mecánicamente y preferiría si puedo > implementar una función dentro del script que me permita realizar ese > trabajo. > > Este es mi script: > > kv<-seq(0.55,7.02,0.1) > lumv<-predict(kluxreg3,list(klux3=kv),interval="confidence",level=0.95) > lines(kv,lumv,lty=19,col="black") > > Agradezco cualquier ayuda, > > Eric Flores De Gracia > School of Biosciences > University of Exeter, Cornwall Campus > Penryn, Cornwall > TR10 9EZ > United Kingdom > Mobilephone: +044 07578724705 > > http://biosciences.exeter.ac.uk/cec/staff/postgradresearch/index.php?web_id=Eric_Flores_de_Gracia&tab=profile > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es >[[alternative HTML version deleted]]
Hola, El año pasado también salió esta pregunta en la lista y se aportó (Argel Gastelum) código de cómo hacerlo. El correo está aquí: https://stat.ethz.ch/pipermail/r-help-es/2011-January/001782.html Saludos, Carlos Ortega www.qualityexcellence.es El 22 de mayo de 2012 14:51, Flores-de-Gracia, Eric <eef201@exeter.ac.uk>escribió:> > Hola > > Disculpen nuevamente por esta nueva pregunta. > > Tengo un modelo no lineal fiteado a mis datos, me interesa presentar en el > plot las líneas de intervalos de confianza del modelo. Se que con los > errores estandar de cada parametro puedo calcular los límites superior e > inferior pero no quiero hacer eso mecánicamente y preferiría si puedo > implementar una función dentro del script que me permita realizar ese > trabajo. > > Este es mi script: > > kv<-seq(0.55,7.02,0.1) > lumv<-predict(kluxreg3,list(klux3=kv),interval="confidence",level=0.95) > lines(kv,lumv,lty=19,col="black") > > Agradezco cualquier ayuda, > > Eric Flores De Gracia > School of Biosciences > University of Exeter, Cornwall Campus > Penryn, Cornwall > TR10 9EZ > United Kingdom > Mobilephone: +044 07578724705 > > http://biosciences.exeter.ac.uk/cec/staff/postgradresearch/index.php?web_id=Eric_Flores_de_Gracia&tab=profile > > _______________________________________________ > 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]]