Ich möchte jede Komponente des Prädiktors aus einem GAM-Modell separat mit der Option type="terms"
bewerten. Als Plausibilitätsprüfung habe ich die Ergebnisse mit einer Auswertung der Gesamtvorhersage unter Verwendung der Option type="response"
verglichen.mgcv: predict.gam() gibt unterschiedliche Ergebnisse für type = "terms" und type = "response"
Es stellt sich heraus, dass die Ergebnisse unterschiedlich sind. Hier ist ein Beispiel:
library(mgcv)
n<-200
sig <- 2
dat <- gamSim(1,n=n,scale=sig)
b<-gam(y~x0+s(I(x1^2))+s(x2)+offset(x3),da=dat)
nd <- data.frame(x0=c(.25,.5),x1=c(.25,.5),x2=c(.25,.5),x3=c(.25,.5))
a1 <- predict.gam(b,newdata=nd,type="response")
a2 <- rowSums(predict.gam(b,newdata=nd,type="terms")) + b$coefficients[1]
a1 - a2 # Should be zero!
# 1 2
# 0.25 0.50
Kann mir jemand bei diesem Problem helfen? Vielen Dank für Ihre Hilfe!