María Ángeles Onieva
2020-Sep-01 10:57 UTC
[R-es] Cálculo - intervalo de confianza - modelo nls - predict
Buenas tardes, Quisiera obtener el intervalo de confianza (y también intervalos de predicción) para los valores predichos en un modelo nls. ¿Hay alguna manera que no sea por ggplot2 (me interesaría obtener el valor listado -además de en el gráfico-) o por bootstrap? Os copio el código del ajuste del modelo y predicción para los 3 días siguientes: *#Ajuste del modelo* model = nls(formula = N~K*exp(-log(K/N0)*exp(-a*(t-t0))), data = datos, start = list(K=300, a = 0.25)) *#Predicción para 3 días* new_juliano c(juliano,juliano[(length(juliano))]+1,juliano[(length(juliano))]+2,juliano[(length(juliano))]+3) casos_predichos = predict(model,data.frame(t = new_juliano),interval "conficende",level = 0.95) Teóricamente debería devolver los intervalos con esto último, sin embargo, no los obtengo. Muchas gracias de antemano. Un cordial saludo, -- María Ángeles [[alternative HTML version deleted]]
Carlos Ortega
2020-Sep-01 11:26 UTC
[R-es] Cálculo - intervalo de confianza - modelo nls - predict
Hola, A través del paquete "broom" vas a poder obtenerlos de una forma muy compacta y manejable. Mira cómo hacerlo en esta función de "broom": - tidy.nls {broom} Gracias, Carlos Ortega www.qualityexcellence.es El mar., 1 sept. 2020 a las 12:57, María Ángeles Onieva (< m.onieva.medica3 en gmail.com>) escribió:> Buenas tardes, > > Quisiera obtener el intervalo de confianza (y también intervalos de > predicción) para los valores predichos en un modelo nls. > ¿Hay alguna manera que no sea por ggplot2 (me interesaría obtener el valor > listado -además de en el gráfico-) o por bootstrap? > > Os copio el código del ajuste del modelo y predicción para los 3 días > siguientes: > > *#Ajuste del modelo* > > model = nls(formula = N~K*exp(-log(K/N0)*exp(-a*(t-t0))), > data = datos, > start = list(K=300, a = 0.25)) > > *#Predicción para 3 días* > > new_juliano > > c(juliano,juliano[(length(juliano))]+1,juliano[(length(juliano))]+2,juliano[(length(juliano))]+3) > casos_predichos = predict(model,data.frame(t = new_juliano),interval > "conficende",level = 0.95) > > Teóricamente debería devolver los intervalos con esto último, sin embargo, > no los obtengo. > > Muchas gracias de antemano. > Un cordial saludo, > > -- > > María Ángeles > > [[alternative HTML version deleted]] > > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es >-- Saludos, Carlos Ortega www.qualityexcellence.es [[alternative HTML version deleted]]
Marcelino de la Cruz Rot
2020-Sep-01 12:43 UTC
[R-es] Cálculo - intervalo de confianza - modelo nls - predict
Hola: Igual es porque has escrito: predict(model,data.frame(t = new_juliano),interval "conficende",level = 0.95) en vez de: predict(model,data.frame(t = new_juliano),interval "confidence",level = 0.95) Un saludo, Marcelino El 01/09/2020 a las 12:57, María Ángeles Onieva escribió:> Buenas tardes, > > Quisiera obtener el intervalo de confianza (y también intervalos de > predicción) para los valores predichos en un modelo nls. > ¿Hay alguna manera que no sea por ggplot2 (me interesaría obtener el valor > listado -además de en el gráfico-) o por bootstrap? > > Os copio el código del ajuste del modelo y predicción para los 3 días > siguientes: > > *#Ajuste del modelo* > > model = nls(formula = N~K*exp(-log(K/N0)*exp(-a*(t-t0))), > data = datos, > start = list(K=300, a = 0.25)) > > *#Predicción para 3 días* > > new_juliano > c(juliano,juliano[(length(juliano))]+1,juliano[(length(juliano))]+2,juliano[(length(juliano))]+3) > casos_predichos = predict(model,data.frame(t = new_juliano),interval > "conficende",level = 0.95) > > Teóricamente debería devolver los intervalos con esto último, sin embargo, > no los obtengo. > > Muchas gracias de antemano. > Un cordial saludo, >-- Marcelino de la Cruz Rot Depto. de Biología y Geología Física y Química Inorgánica Universidad Rey Juan Carlos Móstoles España
Jorge I Velez
2020-Sep-01 14:34 UTC
[R-es] Cálculo - intervalo de confianza - modelo nls - predict
Marcelino dio en el blanco! @María: Para más detalles sobre cómo realizar predicción en modelos nls desde base, dale una mirada a la ayuda de ?predict.nls Felicidades, Jorge.- On Tue, Sep 1, 2020 at 7:43 AM Marcelino de la Cruz Rot < marcelino.delacruz en urjc.es> wrote:> Hola: > > Igual es porque has escrito: > > predict(model,data.frame(t = new_juliano),interval > "conficende",level = 0.95) > > > en vez de: > > predict(model,data.frame(t = new_juliano),interval > "confidence",level = 0.95) > > > Un saludo, > > Marcelino > > El 01/09/2020 a las 12:57, María Ángeles Onieva escribió: > > Buenas tardes, > > > > Quisiera obtener el intervalo de confianza (y también intervalos de > > predicción) para los valores predichos en un modelo nls. > > ¿Hay alguna manera que no sea por ggplot2 (me interesaría obtener el > valor > > listado -además de en el gráfico-) o por bootstrap? > > > > Os copio el código del ajuste del modelo y predicción para los 3 días > > siguientes: > > > > *#Ajuste del modelo* > > > > model = nls(formula = N~K*exp(-log(K/N0)*exp(-a*(t-t0))), > > data = datos, > > start = list(K=300, a = 0.25)) > > > > *#Predicción para 3 días* > > > > new_juliano > > > c(juliano,juliano[(length(juliano))]+1,juliano[(length(juliano))]+2,juliano[(length(juliano))]+3) > > casos_predichos = predict(model,data.frame(t = new_juliano),interval > > "conficende",level = 0.95) > > > > Teóricamente debería devolver los intervalos con esto último, sin > embargo, > > no los obtengo. > > > > Muchas gracias de antemano. > > Un cordial saludo, > > > > -- > Marcelino de la Cruz Rot > Depto. de Biología y Geología > Física y Química Inorgánica > Universidad Rey Juan Carlos > Móstoles España > > _______________________________________________ > 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]]