Я немного растерялся. Учитывая, что у меня есть несколько начальных случаев, как мне спрогнозировать остальную часть кривой, если я предполагаю, что она примет форму функции ошибки s-образная кривая .
Данные:
dat <- data.frame(
time = 1:10,
number = c(0.5, 1, 2, 4, 8, 16, 31, 60, 110, 200)
)
Поскольку кривая уже становится более плоской, должна быть возможность предсказать, когда она будет насыщенной.
Я думал, что это сработает, но я получаю ошибку "единый градиент":
model <-nls(
formula = time ~ 2 * pnorm(number * sqrt(2)) - 1,
data = dat,
start = list(number = 0)
)
predict(model, newdata = data.frame(time = 11:50, number = NA))
Я благодарен за любые подсказки!