# NOT RUN {
# }
# NOT RUN {
library(evgam)
data(COprcp)
COprcp$year <- format(COprcp$date, "%Y")
COprcp_gev <- aggregate(prcp ~ year + meta_row, COprcp, max)
COprcp_gev <- cbind(COprcp_gev, COprcp_meta[COprcp_gev$meta_row,])
fmla_gev <- list(prcp ~ s(lon, lat, k=30) + s(elev, bs="cr"), ~ s(lon, lat, k=20), ~ 1)
m_gev <- evgam(fmla_gev, data=COprcp_gev, family="gev")
logLik(m_gev)
AIC(m_gev)
BIC(m_gev)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab