[R-es] Cálculo - intervalo de confianza - modelo nls - predict

Jorge I Velez jorge|v@nve|ez @end|ng |rom gm@||@com
Mar Sep 1 16:34:20 CEST 2020


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 using 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 using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>

	[[alternative HTML version deleted]]



Más información sobre la lista de distribución R-help-es