# Fitting APM response, Height, MPG, Age as predictors
APM.gam <- gam(APM~s(Height)+s(MPG)+s(Age),data=baskball.df)

# Plotting all predictors
par(mfrow=c(2,2))
plot(APM.gam,se=TRUE)