Temperature*Flooding

TemperatureFlooding.model = lmer (Fluxes ~ Temperature*Flooding + (1|PID) + (1|Temperature) + (1|Flooding), data = BETADecompositioncsv )
    
Temperature.modelNULL = lmer (Fluxes ~ Temperature + (1|PID) + (1|Temperature) + (1|Flooding),data = BETADecompositioncsv, REML=FALSE)

TemperatureFlooding.modelFULL = lmer (Fluxes ~ Temperature*Flooding + (1|PID) + (1|Temperature) + (1|Flooding),
                        data = BETADecompositioncsv, REML=FALSE)

anova(Flooding.modelNULL,TemperatureFlooding.modelFULL)
    
plot(fitted(TemperatureFlooding.model),residuals(TemperatureFlooding.model))