Mam równanie, aby przewidzieć wagę manatów na podstawie ich wieku, w dniach (dias, po portugalsku):
R <- function(a, b, c, dias) c + a*(1 - exp(-b*dias))
Modelowałem to w R, używając nls () i otrzymałem następującą grafikę:
Teraz chcę obliczyć 95% przedział ufności i wykreślić go na grafice. Użyłem dolnej i górnej granicy dla każdej zmiennej a, b i c, w ten sposób:
lower a = a - 1.96*(standard error of a)
higher a = a + 1.96*(standard error of a)
(the same for b and c)
następnie wykreślam dolną linię za pomocą niższej a, b, c, a wyższą linię za pomocą wyższej a, b, c. Ale nie jestem pewien, czy to właściwy sposób, aby to zrobić. Daje mi tę grafikę:
Czy to jest sposób, aby to zrobić, czy robię to źle?