Подгонка gompertz к набору данных? - PullRequest
0 голосов
/ 29 ноября 2018

Я пытаюсь приспособить модель роста Гомперца к приведенным ниже данным, но она генерирует ошибку.Моя цель - увидеть, каким будет объем продаж в будущем, если он будет следовать той же тенденции по всей функции.Верхний предел, вероятно, должен быть около 11 000.

 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 = c(2005, 2006, 2007, 2008, 2009, 2010, 2011, 2012, 
     2013, 2014, 2015, 2016, 2017), 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)), row.names = c(NA, -13L), .Names = c("Country", 
 "Year", "TotalSales"), class = "data.frame")

EV.ss <- nls(TotalSales ~ SSgompertz(Year, Asym, b2, b3), data = EV_temp)

Дает ошибку:

Error in qr.qty(QR.rhs, .swts * ddot(attr(rhs, "gradient"), lin)) : 
   NA/NaN/Inf in foreign function call (arg 5)
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...