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 < [email protected]> 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 > [email protected] > https://stat.ethz.ch/mailman/listinfo/r-help-es > [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
