функция экспоненциального роста в R - подгонке данных

Question

У меня есть следующие данные:

> x = 1:40
> y =
 5    6   28   30   31   34   39   48   63   70   82   91  107  112  127  146  171  198  258  334  403  497  571  657  730  883 1024 1139
1329 1635 2059 2545 3105 3684 4289 4778 5351 5916 6729 7600

Как определить экспоненциальную функцию (включая коэффициенты и константы) этого график в R?

Answer 1

В вашем посте у вас есть только 38 значений для y, поэтому я в основном предполагаю, что x будет 1:38. Если y = exp(ax+b), вы можете изменить его на log(y) = ax + b и установить линейную модель. Ниже будут работать с правильными значениями:

x = 1:38
y = c(5, 6, 28, 30, 31, 34, 39, 48, 63, 70, 82, 91, 107, 112, 127, 
146, 171, 198, 258, 334, 403, 497, 571, 657, 730, 883, 1024, 
1635, 2059, 2545, 3105, 3684, 4289, 4778, 5351, 5916, 6729, 7600
)

fit = lm(log(y) ~ x)
plot(x,y)
lines(x,exp(fitted(fit)),col="blue")