NLME сходится в R - PullRequest
Новая
       46

NLME сходится в R

0 голосов
/ 24 апреля 2020

Я боролся за то, чтобы последний этап моего NLME сходился, находя правильные начальные значения для констант a, b и c. Как и в большинстве кодов, последний шаг, вероятно, является наиболее важным и будет необходим.

Пожалуйста, помогите.

код: 

Type3 <- function(DBHcm, a, b, c){
y <- (a*DBHcm^b)*exp(c*DBHcm)
return(y)
}

fit_data<- groupedData(Prod~DBHcm|Machine_type, data=NLS_data)
Model3f <- nlme (Prod ~ Type3(DBHcm, a, b, c),
                 data = fit_data,
                 fixed = a +b +c ~ 1, 
                 random = a+ b + c ~1|Machine_type, 
                 weights = varPower(),
                 na.action=na.exclude,
                 start = c(a = c(0.000039), b = c(5.387), c = c(-0.132))) #this is where i am battling
...