Displaying 1 result from an estimated 1 matches for "plkw".
Did you mean:
pkw
2002 Dec 18
1
problem with 'gnls'
...and I get the following Error:
>> Error in "coef<-.corARMA"(*tmp*, value = c(174.267493382104, 173.412740072763 :
>> Coefficient matrix not invertible
Here is my R-code:
data <- d.plabutsch.neu
# calculating the starting coeficients:
r.CO2.rlm <- rlm(EF.CO2 ~ pLKW.total + Sat + Sun + Fz.total,
data = data, method = "MM", na.action =
na.exclude)
START <- coef(r.CO2.rlm); names(START) <- letters[1:5]
# 'gnls':
r.CO2.gnls <- gnls(log(EF.CO2) ~ log(a + b * pLKW.total) + c * Sat +...