library(gam);
# Load the Pima Indians Diabetes dataset
data(PimaIndiansDiabetes2, package="mlbench");
# Calculate smooth odds ratios using flexOR
mod1 <- flexOR(
data=PimaIndiansDiabetes2,
response="diabetes",
formula= ~ s(age) + s(mass) + s(pedigree) + pressure + glucose
);
# Predict the probabilities using predict.OR
predict(mod1, predictor="age", ref.value=40)
Run the code above in your browser using DataLab