我想计算GAM模型中每个预测变量的相对重要性,因此我尝试了relaimpo
软件包。下面是一个示例:
library(relaimpo)
library(mgcv)
gam1 <- gam(mpg ~ s(drat) + s(wt) + s(qsec), data = mtcars, method = "REML")
summary(gam1) # Here is one R^2
calc.relimp(gam1) # Here is the other R^2
很奇怪,summary()
中的R ^ 2为0.867,而calc.relimp()
中的R ^ 2为0.837(83.7%)。为什么会发生这种分歧?