# Load dataset
library(gam)
data(PimaIndiansDiabetes2, package="mlbench");
# Calculate degrees of freedom using AIC
df2 <- dfgam(
response="diabetes",
nl.predictors=c("age", "mass"),
other.predictors=c("pedigree"),
smoother="s",
method="AIC",
data=PimaIndiansDiabetes2
);
print(df2$df);
Run the code above in your browser using DataLab