У меня недавно возникла проблема с запуском модели GAM из ранее работающего кода. Я считаю, что это связано с обновленной версией R и обновленной версией пакета mgcv. Так что было бы здорово узнать, есть ли у кого-то одна проблема или есть решение.Запуск функции GAM из пакета mgcv с R 3.2.2
Я в настоящее время работает: R version 3.2.2 (2015-08-14) - «Пожарная безопасность» на Windows. И я использую пакет mgcv 1.8-7.
Ниже приведен пример кода, который повторно выводит сообщение об ошибке при запуске на моем компьютере.
###Load package
library(mgcv)
Это mgcv 1.8-7.
###Simulate some example data
set.seed(2) ## simulate some data...
dat <- gamSim(1,n=400,dist="normal",scale=2)
###Run normal model
b <- gam(y~s(x0)+s(x1)+s(x2)+s(x3),data=dat, family=gaussian())
Это работает.
###change the smoothness selection method to REML
b0 <- gam(y~s(x0)+s(x1)+s(x2)+s(x3),data=dat,method="REML")
дает следующее сообщение об ошибке:
Ошибка в .C (C_gdi1, X = as.double (х [хороший,]), Е = as.double (Sr), Eb = а. двойной (Eb),: Неверное количество аргументов (48), ожидая, что 47 для «gdi1»
Спасибо за вашу помощь
FWIW, я не получаю ошибку wh ru Я запускаю ваш код в R 3.2.1 с mgcv 1.8-7. Я также не получил предупреждения при загрузке mgcv о пакете, который был построен в R 3.2.2. Таким образом, это похоже на проблему с версией. – ulfelder