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

[Marcelino de la Cruz Rot]

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