在使用mgcv软件包对其中一个协变量使用micv平滑函数(来自scam软件包)构建的GAM模型进行预测时,我得到了以下错误:
Error in predict.gam(mod_gam, newdata = dfGAMDev) :
number of items to replace is not a multiple of replacement length代码:
library(dplyr)
library(ggplot2)
library(car)
library(mgcv)
library(scam)
mod_gam <- gam(Y ~ X1 + X2+ X3+ s(X4, bs='micv') + X5, data = dfGAMDev )
predDevGAM <- predict(mod_gam, newdata =dfGAMDev)我也尝试了以下几种方法:
predDevGAM <- predict.gam(mod_gam, newdata =dfGAMDev)所有的变量都是数字的
发布于 2015-12-17 18:49:50
我一直在尝试建立一个有形状约束的GAM模型,用无约束的GAM模型函数,用gam代替scam函数。以下代码起作用了:
mod_gam <- scam(Y ~ X1 + X2+ X3+ s(X4, bs='micv') + X5, data = dfGAMDev)很抱歉提出了这个问题。我误解了,scam包中的平滑函数可以在gam函数中使用。
https://stackoverflow.com/questions/34330110
复制相似问题