[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]]

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,

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