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.medi...@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@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
--
Saludos,
Carlos Ortega
www.qualityexcellence.es
[[alternative HTML version deleted]]
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
