#Necessary packages
library(lme4) 
library(car)
library(mgcv)


#Soil chemistry, liquid chemistry and ICP data models


model = lmer(variable | log10(variable) ~  Day*Basalt_bin*Bacillus_bin + (1|Mesocosm_no),data=data,REML= TRUE)

summary(model)
Anova(model, type="II")

