Я пытаюсь согласовать данные о мировых продажах BEV и PHEV с функцией логистического роста.
Например, я использую R-код для модели логистического роста населения Австралии вR (Элиз Гулд)
Но с помощью этого метода я получаю ошибку:
Ошибка в qr.default (.swts * gr): NA / NaN/ Inf в вызове внешней функции (arg 1) Дополнительно: предупреждающее сообщение: In Ops.factor (xmid, input): '-' не имеет значения для факторов
Мои данные:
EV_temp<-structure(list(Country = structure(c(21L, 21L, 21L, 21L, 21L,
21L, 21L, 21L, 21L, 21L, 21L, 21L, 21L), .Label = c("Australia",
"Brazil", "Canada", "Chile", "China", "Finland", "France", "Germany",
"India", "Japan", "Korea", "Mexico", "Netherlands", "New Zealand",
"Norway", "Others", "Portugal", "South Africa", "Sweden", "Thailand",
"Total", "United Kingdom", "United States"), class = "factor"),
Year = structure(1:13, .Label = c("2005", "2006", "2007",
"2008", "2009", "2010", "2011", "2012", "2013", "2014", "2015",
"2016", "2017"), class = "factor"), TotalSales = c(1.89,
2.23, 2.69, 5.15, 7.48, 14.26, 61.33, 179.03, 381.3, 703.65,
1239.45, 1982.04, 3109.05)), class = "data.frame", row.names = c(NA,
-13L), .Names = c("Country", "Year", "TotalSales"))
Код:
EV.ss <- nls(TotalSales ~ SSlogis(Year, phi1, phi2, phi3), data = EV_temp)
Если бы кто-то мог указать мне правильное направление в решении этого, это было бы очень признательно.