We simulated a dataset and created a model.
set.seed(459) # seed mass n <- 1000 seed.mass <- round(rnorm(n, mean = 250, sd = 75),digits = 1) ## Settin