Tengo una ecuación para predecir el peso de los manatíes a partir de su edad, en días (dias, en portugués):
R <- function(a, b, c, dias) c + a*(1 - exp(-b*dias))
Lo modelé en R, usando nls (), y obtuve este gráfico:
Ahora quiero calcular el intervalo de confianza del 95% y trazarlo en el gráfico. He usado los límites inferior y superior para cada variable a, byc, así:
lower a = a - 1.96*(standard error of a)
higher a = a + 1.96*(standard error of a)
(the same for b and c)
luego trazo una línea inferior usando a, b, c y una línea más alta usando a, b, c. Pero no estoy seguro de si esa es la forma correcta de hacerlo. Me está dando este gráfico:
¿Es esta la forma de hacerlo, o lo estoy haciendo mal?