У меня проблемы с запуском функции gls()
в nlme
(пакет R). Он продолжает выдавать одно и то же сообщение об ошибке:
"Error: object of type 'closure' is not subsettable".
Основная строка кода проблемы:
temp.gls <- gls(SS.Temp.Age[,2] ~ SS.Temp.Age[,1], data = SS.Temp.Age, correlation = corAR1, subset = c(1:5956), method = "ML")
Но позвольте мне просмотреть предыдущий код, включая соответствующие пакеты R:
library(RPANDA)
library(nlme)
data(InfTemp) # The data is temperature through geologic time
InfTempSmoothed <- smooth.spline(InfTemp[,1], InfTemp[,2])
# I use a smooth spline function on the data.
SS.Temp.Age <- cbind(InfTempSmoothed$x,InfTempSmoothed$y)
# And then the line where I try to run gls()
temp.gls <- gls(SS.Temp.Age[,2] ~ SS.Temp.Age[,1], data = SS.Temp.Age, correlation = corAR1, subset = c(1:5956), method = "ML")
Любая помощь будет принята с благодарностью! Пожалуйста, дайте мне знать, если есть какие-либо разъяснения.