# Times of landfalling US hurricanes
data(USlandfall)
# convert dates to years, as a continuous variable
year <- as.integer(format(USlandfall$date, "%Y"))
day <- as.integer(format(USlandfall$date, "%j"))
USlandfall$year <- year + pmin(day / 365, 1)
hits <- subset(USlandfall, landfall == 1)
# this creates nodes in the default way
m1 <- ppgam( ~ s(year), hits)
print(m1)
Run the code above in your browser using DataLab