GLMs and GAMs comparison in global study
layout(matrix(c(1,2,3,3), 2, 2, byrow=TRUE))
plot(GLMresults$aucs,xlab="number of covariates",ylab="AUC",
main="GLM",ylim=c(0.6,0.9))
lines(GLMresults$aucopt)
plot(GAMresults$aucs,xlab="number of covariates",ylab="AUC",
main="GAM",ylim=c(0.6,0.9))
lines(GAMresults$aucopt)
plot(GLMresults$aucopt,xlab="number of covariates",ylab="AUC",
ylim=c(0.77,0.90),type='b',pch=1,main="GLM vs. GAM",cex=1.24)
lines(GAMresults$aucopt,type='b',pch=2,cex=1.25)
legend("bottomleft",c("GLM","GAM"),lty=1,pch=c(1,2),
box.lty=0,cex=1.25)