#fit a model
deathm1 <- mixed_model(data_2w_Tdeath,
"PI", c("Genotype", "Time"),
"Experiment")
#plot model
plot_lm_predict(deathm1,
Genotype, PI, Time)
#fit zooplankton data
z1 <- ga_model(data = data_zooplankton,
Y_value = "log(density_adj)",
Fixed_Factor = "taxon",
Smooth_Factor = "day")
#plot fitted data
plot_lm_predict(Model = z1,
xcol = day,
ycol = `log(density_adj)`,
ByFactor = taxon)
Run the code above in your browser using DataLab