2015-08-31 2 views
0

У меня недавно возникла проблема с запуском модели 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»

Спасибо за вашу помощь

+0

FWIW, я не получаю ошибку wh ru Я запускаю ваш код в R 3.2.1 с mgcv 1.8-7. Я также не получил предупреждения при загрузке mgcv о пакете, который был построен в R 3.2.2. Таким образом, это похоже на проблему с версией. – ulfelder

ответ

0

я заново установил R и пакет mgcv и, кажется, как будто это имеет разрешил проблему.