#load example data
data(vi)
#reduce size to speed-up example
vi <- vi[1:1000, ]
#integer counts response and continuous predictor
#to data frame without NAs
df <- data.frame(
y = vi[["vi_counts"]],
x = vi[["swi_max"]]
) |>
na.omit()
#GLM model with Poisson family
f_r2_glm_poisson(df = df)
#GLM model with second degree polynomials and Poisson family
f_r2_glm_poisson_poly2(df = df)
#GAM model with Poisson family
f_r2_gam_poisson(df = df)
Run the code above in your browser using DataLab