Total C
TCnull <- lmer(log(TotalC_percent+1) ~ 1 + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
TCT <- lmer(log(TotalC_percent+1) ~ TSF + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
TCE <- lmer(log(TotalC_percent+1) ~ ESD + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
TCTE <- lmer(log(TotalC_percent+1) ~ TSF + ESD + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
anova(TCnull, TCT, TCE, TCTE)
## Data: SoilJuly181920
## Models:
## TCnull: log(TotalC_percent + 1) ~ 1 + (1 | Location/Year)
## TCT: log(TotalC_percent + 1) ~ TSF + (1 | Location/Year)
## TCE: log(TotalC_percent + 1) ~ ESD + (1 | Location/Year)
## TCTE: log(TotalC_percent + 1) ~ TSF + ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## TCnull 4 -68.582 -53.930 38.291 -76.582
## TCT 8 -71.060 -41.756 43.530 -87.060 10.478 4 0.03311 *
## TCE 8 -137.872 -108.568 76.936 -153.872 66.812 0
## TCTE 12 -141.023 -97.068 82.512 -165.023 11.152 4 0.02491 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
anova(TCnull, TCT)
## Data: SoilJuly181920
## Models:
## TCnull: log(TotalC_percent + 1) ~ 1 + (1 | Location/Year)
## TCT: log(TotalC_percent + 1) ~ TSF + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## TCnull 4 -68.582 -53.930 38.291 -76.582
## TCT 8 -71.060 -41.756 43.530 -87.060 10.478 4 0.03311 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
anova(TCnull, TCE)
## Data: SoilJuly181920
## Models:
## TCnull: log(TotalC_percent + 1) ~ 1 + (1 | Location/Year)
## TCE: log(TotalC_percent + 1) ~ ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## TCnull 4 -68.582 -53.93 38.291 -76.582
## TCE 8 -137.872 -108.57 76.936 -153.872 77.29 4 6.53e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
anova(TCnull, TCTE)
## Data: SoilJuly181920
## Models:
## TCnull: log(TotalC_percent + 1) ~ 1 + (1 | Location/Year)
## TCTE: log(TotalC_percent + 1) ~ TSF + ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## TCnull 4 -68.582 -53.930 38.291 -76.582
## TCTE 12 -141.023 -97.068 82.512 -165.023 88.441 8 9.632e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Mult_TC_TSF <- glht(TCTE, linfct=mcp(TSF="Tukey"))
summary(Mult_TC_TSF)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(TotalC_percent + 1) ~ TSF + ESD + (1 | Location/Year),
## data = SoilJuly181920, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1yr - 2yr - 1yr < == 0 0.008654 0.029044 0.298 0.9982
## 2yr - 3yr - 1yr < == 0 -0.032951 0.032971 -0.999 0.8520
## 3yr - 4yr - 1yr < == 0 -0.025971 0.042914 -0.605 0.9735
## Unburned - 1yr < == 0 -0.082160 0.029841 -2.753 0.0449 *
## 2yr - 3yr - 1yr - 2yr == 0 -0.041604 0.032924 -1.264 0.7075
## 3yr - 4yr - 1yr - 2yr == 0 -0.034625 0.042760 -0.810 0.9256
## Unburned - 1yr - 2yr == 0 -0.090814 0.029846 -3.043 0.0190 *
## 3yr - 4yr - 2yr - 3yr == 0 0.006979 0.044776 0.156 0.9999
## Unburned - 2yr - 3yr == 0 -0.049210 0.034687 -1.419 0.6084
## Unburned - 3yr - 4yr == 0 -0.056189 0.044731 -1.256 0.7120
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
Mult_TC_ESD <- glht(TCTE, linfct=mcp(ESD="Tukey"))
summary(Mult_TC_ESD)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(TotalC_percent + 1) ~ TSF + ESD + (1 | Location/Year),
## data = SoilJuly181920, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 0.077199 0.042648 1.810 0.33889
## SaLo - Clayey == 0 0.123282 0.033014 3.734 0.00142 **
## Sandy - Clayey == 0 -0.112471 0.030461 -3.692 0.00174 **
## ThCl - Clayey == 0 0.070941 0.081428 0.871 0.89656
## SaLo - Loamy == 0 0.046083 0.043149 1.068 0.80477
## Sandy - Loamy == 0 -0.189670 0.038479 -4.929 < 0.001 ***
## ThCl - Loamy == 0 -0.006258 0.085231 -0.073 0.99999
## Sandy - SaLo == 0 -0.235753 0.026717 -8.824 < 0.001 ***
## ThCl - SaLo == 0 -0.052340 0.077089 -0.679 0.95592
## ThCl - Sandy == 0 0.183412 0.077709 2.360 0.11139
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
#2018
TCTE18 <- lmer(log(TotalC_percent+1) ~ TSF + ESD + (1|Location), data=SoilJuly18, REML = FALSE)
Mult_TC18_TSF <- glht(TCTE18, linfct=mcp(TSF="Tukey"))
summary(Mult_TC18_TSF)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(TotalC_percent + 1) ~ TSF + ESD + (1 | Location),
## data = SoilJuly18, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1yr - 2yr - 1yr < == 0 -0.04531 0.04632 -0.978 0.5881
## Unburned - 1yr < == 0 -0.08429 0.03964 -2.126 0.0835 .
## Unburned - 1yr - 2yr == 0 -0.03898 0.03836 -1.016 0.5641
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
Mult_TC18_ESD <- glht(TCTE18, linfct=mcp(ESD="Tukey"))
summary(Mult_TC18_ESD)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(TotalC_percent + 1) ~ TSF + ESD + (1 | Location),
## data = SoilJuly18, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 0.05952 0.06433 0.925 0.8744
## SaLo - Clayey == 0 0.10514 0.04991 2.106 0.1953
## Sandy - Clayey == 0 -0.11992 0.04573 -2.623 0.0575 .
## ThCl - Clayey == 0 0.03281 0.12643 0.260 0.9989
## SaLo - Loamy == 0 0.04562 0.06432 0.709 0.9486
## Sandy - Loamy == 0 -0.17944 0.05821 -3.083 0.0149 *
## ThCl - Loamy == 0 -0.02671 0.12993 -0.206 0.9995
## Sandy - SaLo == 0 -0.22506 0.04078 -5.519 <0.001 ***
## ThCl - SaLo == 0 -0.07232 0.11865 -0.610 0.9701
## ThCl - Sandy == 0 0.15274 0.12105 1.262 0.6900
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
#2019
TCTE19 <- lmer(log(TotalC_percent+1) ~ TSF + ESD + (1|Location), data=SoilJuly19, REML = FALSE)
Mult_TC19_TSF <- glht(TCTE19, linfct=mcp(TSF="Tukey"))
summary(Mult_TC19_TSF)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(TotalC_percent + 1) ~ TSF + ESD + (1 | Location),
## data = SoilJuly19, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1yr - 2yr - 1yr < == 0 -0.03972 0.04986 -0.797 0.8559
## 2yr - 3yr - 1yr < == 0 -0.09165 0.04865 -1.884 0.2350
## Unburned - 1yr < == 0 -0.14811 0.04811 -3.079 0.0114 *
## 2yr - 3yr - 1yr - 2yr == 0 -0.05193 0.05102 -1.018 0.7389
## Unburned - 1yr - 2yr == 0 -0.10840 0.04986 -2.174 0.1303
## Unburned - 2yr - 3yr == 0 -0.05646 0.04865 -1.161 0.6517
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
Mult_TC19_ESD <- glht(TCTE19, linfct=mcp(ESD="Tukey"))
summary(Mult_TC19_ESD)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(TotalC_percent + 1) ~ TSF + ESD + (1 | Location),
## data = SoilJuly19, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 0.04495 0.07122 0.631 0.9660
## SaLo - Clayey == 0 0.08829 0.05553 1.590 0.4735
## Sandy - Clayey == 0 -0.13501 0.05040 -2.679 0.0493 *
## ThCl - Clayey == 0 -0.09567 0.13913 -0.688 0.9539
## SaLo - Loamy == 0 0.04334 0.07065 0.613 0.9694
## Sandy - Loamy == 0 -0.17996 0.06383 -2.819 0.0333 *
## ThCl - Loamy == 0 -0.14062 0.14284 -0.984 0.8474
## Sandy - SaLo == 0 -0.22330 0.04456 -5.011 <0.001 ***
## ThCl - SaLo == 0 -0.18396 0.13057 -1.409 0.5934
## ThCl - Sandy == 0 0.03934 0.13318 0.295 0.9981
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
#2020
TCTE20 <- lmer(log(TotalC_percent+1) ~ TSF + ESD + (1|Location), data=SoilJuly20, REML = FALSE)
Mult_TC20_TSF <- glht(TCTE20, linfct=mcp(TSF="Tukey"))
summary(Mult_TC20_TSF)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(TotalC_percent + 1) ~ TSF + ESD + (1 | Location),
## data = SoilJuly20, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1yr - 2yr - 1yr < == 0 0.117888 0.055123 2.139 0.141
## 2yr - 3yr - 1yr < == 0 0.031811 0.057011 0.558 0.944
## 3yr - 4yr - 1yr < == 0 0.009232 0.055625 0.166 0.998
## 2yr - 3yr - 1yr - 2yr == 0 -0.086077 0.057006 -1.510 0.431
## 3yr - 4yr - 1yr - 2yr == 0 -0.108656 0.055650 -1.952 0.206
## 3yr - 4yr - 2yr - 3yr == 0 -0.022579 0.058243 -0.388 0.980
## (Adjusted p values reported -- single-step method)
Mult_TC20_ESD <- glht(TCTE20, linfct=mcp(ESD="Tukey"))
summary(Mult_TC20_ESD)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(TotalC_percent + 1) ~ TSF + ESD + (1 | Location),
## data = SoilJuly20, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 0.14122 0.07896 1.788 0.35284
## SaLo - Clayey == 0 0.12107 0.06079 1.992 0.24619
## Sandy - Clayey == 0 -0.09306 0.05653 -1.646 0.43915
## ThCl - Clayey == 0 0.19990 0.15432 1.295 0.66995
## SaLo - Loamy == 0 -0.02016 0.07620 -0.265 0.99878
## Sandy - Loamy == 0 -0.23428 0.07140 -3.282 0.00757 **
## ThCl - Loamy == 0 0.05867 0.15783 0.372 0.99540
## Sandy - SaLo == 0 -0.21412 0.05002 -4.281 < 0.001 ***
## ThCl - SaLo == 0 0.07883 0.14696 0.536 0.98147
## ThCl - Sandy == 0 0.29296 0.14894 1.967 0.25804
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
Total N
#Overall
TNnull <- lmer(TotalN_percent ~ 1 + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
TNT <- lmer(TotalN_percent ~ TSF + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
TNE <- lmer(TotalN_percent ~ ESD + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
TNTE <- lmer(TotalN_percent ~ TSF + ESD + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
anova(TNnull, TNT, TNE, TNTE)
## Data: SoilJuly181920
## Models:
## TNnull: TotalN_percent ~ 1 + (1 | Location/Year)
## TNT: TotalN_percent ~ TSF + (1 | Location/Year)
## TNE: TotalN_percent ~ ESD + (1 | Location/Year)
## TNTE: TotalN_percent ~ TSF + ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## TNnull 4 -859.54 -844.89 433.77 -867.54
## TNT 8 -858.41 -829.11 437.21 -874.41 6.8724 4 0.1428
## TNE 8 -920.65 -891.35 468.32 -936.65 62.2386 0
## TNTE 12 -919.64 -875.68 471.82 -943.64 6.9901 4 0.1364
anova(TNnull, TNT)
## Data: SoilJuly181920
## Models:
## TNnull: TotalN_percent ~ 1 + (1 | Location/Year)
## TNT: TotalN_percent ~ TSF + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## TNnull 4 -859.54 -844.89 433.77 -867.54
## TNT 8 -858.41 -829.11 437.21 -874.41 6.8724 4 0.1428
anova(TNnull, TNE)
## Data: SoilJuly181920
## Models:
## TNnull: TotalN_percent ~ 1 + (1 | Location/Year)
## TNE: TotalN_percent ~ ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## TNnull 4 -859.54 -844.89 433.77 -867.54
## TNE 8 -920.65 -891.35 468.32 -936.65 69.111 4 3.497e-14 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
anova(TNnull, TNTE)
## Data: SoilJuly181920
## Models:
## TNnull: TotalN_percent ~ 1 + (1 | Location/Year)
## TNTE: TotalN_percent ~ TSF + ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## TNnull 4 -859.54 -844.89 433.77 -867.54
## TNTE 12 -919.64 -875.68 471.82 -943.64 76.101 8 2.968e-13 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Mult_TN_TSF <- glht(TNTE, linfct=mcp(TSF="Tukey"))
summary(Mult_TN_TSF)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = TotalN_percent ~ TSF + ESD + (1 | Location/Year),
## data = SoilJuly181920, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1yr - 2yr - 1yr < == 0 -0.0020246 0.0074100 -0.273 0.999
## 2yr - 3yr - 1yr < == 0 -0.0042672 0.0084912 -0.503 0.987
## 3yr - 4yr - 1yr < == 0 -0.0032920 0.0112100 -0.294 0.998
## Unburned - 1yr < == 0 -0.0192555 0.0077430 -2.487 0.090 .
## 2yr - 3yr - 1yr - 2yr == 0 -0.0022426 0.0084780 -0.265 0.999
## 3yr - 4yr - 1yr - 2yr == 0 -0.0012675 0.0111689 -0.113 1.000
## Unburned - 1yr - 2yr == 0 -0.0172310 0.0077466 -2.224 0.165
## 3yr - 4yr - 2yr - 3yr == 0 0.0009752 0.0115767 0.084 1.000
## Unburned - 2yr - 3yr == 0 -0.0149883 0.0091991 -1.629 0.470
## Unburned - 3yr - 4yr == 0 -0.0159635 0.0120079 -1.329 0.665
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
Mult_TN_ESD <- glht(TNTE, linfct=mcp(ESD="Tukey"))
summary(Mult_TN_ESD)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = TotalN_percent ~ TSF + ESD + (1 | Location/Year),
## data = SoilJuly181920, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 0.004694 0.010877 0.432 0.9918
## SaLo - Clayey == 0 0.022302 0.008423 2.648 0.0539 .
## Sandy - Clayey == 0 -0.035008 0.007771 -4.505 <0.001 ***
## ThCl - Clayey == 0 0.002296 0.020767 0.111 1.0000
## SaLo - Loamy == 0 0.017608 0.011001 1.601 0.4666
## Sandy - Loamy == 0 -0.039702 0.009813 -4.046 <0.001 ***
## ThCl - Loamy == 0 -0.002398 0.021734 -0.110 1.0000
## Sandy - SaLo == 0 -0.057310 0.006817 -8.407 <0.001 ***
## ThCl - SaLo == 0 -0.020006 0.019669 -1.017 0.8313
## ThCl - Sandy == 0 0.037304 0.019822 1.882 0.2996
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
#2018
TNTE18 <- lmer(TotalN_percent ~ TSF + ESD + (1|Location), data=SoilJuly18, REML = FALSE)
Mult_TN18_TSF <- glht(TNTE18, linfct=mcp(TSF="Tukey"))
summary(Mult_TN18_TSF)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = TotalN_percent ~ TSF + ESD + (1 | Location), data = SoilJuly18,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1yr - 2yr - 1yr < == 0 -0.022372 0.011936 -1.874 0.1446
## Unburned - 1yr < == 0 -0.024943 0.010217 -2.441 0.0384 *
## Unburned - 1yr - 2yr == 0 -0.002571 0.009887 -0.260 0.9631
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
Mult_TN18_ESD <- glht(TNTE18, linfct=mcp(ESD="Tukey"))
summary(Mult_TN18_ESD)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = TotalN_percent ~ TSF + ESD + (1 | Location), data = SoilJuly18,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 0.010252 0.016541 0.620 0.9683
## SaLo - Clayey == 0 0.015956 0.012823 1.244 0.7012
## Sandy - Clayey == 0 -0.029151 0.011765 -2.478 0.0836 .
## ThCl - Clayey == 0 -0.016056 0.032492 -0.494 0.9863
## SaLo - Loamy == 0 0.005705 0.016484 0.346 0.9965
## Sandy - Loamy == 0 -0.039403 0.014965 -2.633 0.0560 .
## ThCl - Loamy == 0 -0.026308 0.033374 -0.788 0.9261
## Sandy - SaLo == 0 -0.045108 0.010485 -4.302 <0.001 ***
## ThCl - SaLo == 0 -0.032012 0.030535 -1.048 0.8155
## ThCl - Sandy == 0 0.013095 0.031131 0.421 0.9925
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
#2019
TNTE19 <- lmer(TotalN_percent ~ TSF + ESD + (1|Location), data=SoilJuly19, REML = FALSE)
Mult_TN19_TSF <- glht(TNTE19, linfct=mcp(TSF="Tukey"))
summary(Mult_TN19_TSF)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = TotalN_percent ~ TSF + ESD + (1 | Location), data = SoilJuly19,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1yr - 2yr - 1yr < == 0 -0.005244 0.012972 -0.404 0.978
## 2yr - 3yr - 1yr < == 0 -0.017343 0.012660 -1.370 0.518
## Unburned - 1yr < == 0 -0.023750 0.012520 -1.897 0.229
## 2yr - 3yr - 1yr - 2yr == 0 -0.012099 0.013272 -0.912 0.799
## Unburned - 1yr - 2yr == 0 -0.018506 0.012972 -1.427 0.482
## Unburned - 2yr - 3yr == 0 -0.006407 0.012660 -0.506 0.958
## (Adjusted p values reported -- single-step method)
Mult_TN19_ESD <- glht(TNTE19, linfct=mcp(ESD="Tukey"))
summary(Mult_TN19_ESD)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = TotalN_percent ~ TSF + ESD + (1 | Location), data = SoilJuly19,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 -0.005051 0.018455 -0.274 0.99859
## SaLo - Clayey == 0 0.003574 0.014374 0.249 0.99903
## Sandy - Clayey == 0 -0.049732 0.013085 -3.801 0.00113 **
## ThCl - Clayey == 0 -0.048271 0.036016 -1.340 0.63920
## SaLo - Loamy == 0 0.008625 0.018194 0.474 0.98824
## Sandy - Loamy == 0 -0.044681 0.016534 -2.702 0.04642 *
## ThCl - Loamy == 0 -0.043221 0.036941 -1.170 0.74686
## Sandy - SaLo == 0 -0.053306 0.011541 -4.619 < 0.001 ***
## ThCl - SaLo == 0 -0.051846 0.033880 -1.530 0.51280
## ThCl - Sandy == 0 0.001460 0.034517 0.042 1.00000
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
#2020
TNTE20 <- lmer(TotalN_percent ~ TSF + ESD + (1|Location), data=SoilJuly20, REML = FALSE)
Mult_TN20_TSF <- glht(TNTE20, linfct=mcp(TSF="Tukey"))
summary(Mult_TN20_TSF)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = TotalN_percent ~ TSF + ESD + (1 | Location), data = SoilJuly20,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1yr - 2yr - 1yr < == 0 0.020729 0.013775 1.505 0.434
## 2yr - 3yr - 1yr < == 0 0.008671 0.014255 0.608 0.929
## 3yr - 4yr - 1yr < == 0 0.003166 0.013903 0.228 0.996
## 2yr - 3yr - 1yr - 2yr == 0 -0.012058 0.014252 -0.846 0.832
## 3yr - 4yr - 1yr - 2yr == 0 -0.017563 0.013910 -1.263 0.587
## 3yr - 4yr - 2yr - 3yr == 0 -0.005505 0.014571 -0.378 0.982
## (Adjusted p values reported -- single-step method)
Mult_TN20_ESD <- glht(TNTE20, linfct=mcp(ESD="Tukey"))
summary(Mult_TN20_ESD)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = TotalN_percent ~ TSF + ESD + (1 | Location), data = SoilJuly20,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 0.01356 0.01996 0.679 0.956
## SaLo - Clayey == 0 0.03088 0.01541 2.003 0.240
## Sandy - Clayey == 0 -0.02749 0.01425 -1.929 0.276
## ThCl - Clayey == 0 0.04571 0.03910 1.169 0.748
## SaLo - Loamy == 0 0.01731 0.01955 0.886 0.891
## Sandy - Loamy == 0 -0.04105 0.01805 -2.275 0.136
## ThCl - Loamy == 0 0.03215 0.04006 0.803 0.922
## Sandy - SaLo == 0 -0.05837 0.01264 -4.616 <0.001 ***
## ThCl - SaLo == 0 0.01483 0.03701 0.401 0.994
## ThCl - Sandy == 0 0.07320 0.03760 1.947 0.267
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
P
#Overall
Pnull <- lmer(log(P_ppm+1) ~ 1 + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
PT <- lmer(log(P_ppm+1) ~ TSF + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
PE <- lmer(log(P_ppm+1) ~ ESD + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
PTE <- lmer(log(P_ppm+1) ~ TSF + ESD + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
anova(Pnull, PT, PE, PTE)
## Data: SoilJuly181920
## Models:
## Pnull: log(P_ppm + 1) ~ 1 + (1 | Location/Year)
## PT: log(P_ppm + 1) ~ TSF + (1 | Location/Year)
## PE: log(P_ppm + 1) ~ ESD + (1 | Location/Year)
## PTE: log(P_ppm + 1) ~ TSF + ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Pnull 4 371.12 385.77 -181.56 363.12
## PT 8 372.68 401.98 -178.34 356.68 6.4405 4 0.1686
## PE 8 287.15 316.46 -135.58 271.15 85.5280 0
## PTE 12 290.15 334.11 -133.08 266.15 4.9994 4 0.2874
anova(Pnull, PT)
## Data: SoilJuly181920
## Models:
## Pnull: log(P_ppm + 1) ~ 1 + (1 | Location/Year)
## PT: log(P_ppm + 1) ~ TSF + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Pnull 4 371.12 385.77 -181.56 363.12
## PT 8 372.68 401.98 -178.34 356.68 6.4405 4 0.1686
anova(Pnull, PE)
## Data: SoilJuly181920
## Models:
## Pnull: log(P_ppm + 1) ~ 1 + (1 | Location/Year)
## PE: log(P_ppm + 1) ~ ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Pnull 4 371.12 385.77 -181.56 363.12
## PE 8 287.15 316.46 -135.58 271.15 91.969 4 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
anova(Pnull, PTE)
## Data: SoilJuly181920
## Models:
## Pnull: log(P_ppm + 1) ~ 1 + (1 | Location/Year)
## PTE: log(P_ppm + 1) ~ TSF + ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Pnull 4 371.12 385.77 -181.56 363.12
## PTE 12 290.15 334.11 -133.08 266.15 96.968 8 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Mult_P_TSF <- glht(PTE, linfct=mcp(TSF="Tukey"))
summary(Mult_P_TSF)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(P_ppm + 1) ~ TSF + ESD + (1 | Location/Year),
## data = SoilJuly181920, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1yr - 2yr - 1yr < == 0 0.06356 0.06054 1.050 0.827
## 2yr - 3yr - 1yr < == 0 0.01565 0.06939 0.225 0.999
## 3yr - 4yr - 1yr < == 0 -0.09733 0.09166 -1.062 0.821
## Unburned - 1yr < == 0 -0.04628 0.06330 -0.731 0.947
## 2yr - 3yr - 1yr - 2yr == 0 -0.04792 0.06929 -0.692 0.957
## 3yr - 4yr - 1yr - 2yr == 0 -0.16089 0.09133 -1.762 0.387
## Unburned - 1yr - 2yr == 0 -0.10984 0.06333 -1.734 0.404
## 3yr - 4yr - 2yr - 3yr == 0 -0.11297 0.09462 -1.194 0.748
## Unburned - 2yr - 3yr == 0 -0.06192 0.07526 -0.823 0.921
## Unburned - 3yr - 4yr == 0 0.05105 0.09829 0.519 0.985
## (Adjusted p values reported -- single-step method)
Mult_P_ESD <- glht(PTE, linfct=mcp(ESD="Tukey"))
summary(Mult_P_ESD)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(P_ppm + 1) ~ TSF + ESD + (1 | Location/Year),
## data = SoilJuly181920, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 -0.11414 0.08881 -1.285 0.6749
## SaLo - Clayey == 0 -0.03555 0.06877 -0.517 0.9837
## Sandy - Clayey == 0 -0.46652 0.06346 -7.351 <0.001 ***
## ThCl - Clayey == 0 0.44662 0.16956 2.634 0.0560 .
## SaLo - Loamy == 0 0.07859 0.08977 0.876 0.8949
## Sandy - Loamy == 0 -0.35238 0.08012 -4.398 <0.001 ***
## ThCl - Loamy == 0 0.56076 0.17743 3.160 0.0115 *
## Sandy - SaLo == 0 -0.43097 0.05566 -7.743 <0.001 ***
## ThCl - SaLo == 0 0.48217 0.16064 3.002 0.0193 *
## ThCl - Sandy == 0 0.91314 0.16187 5.641 <0.001 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
#2018
PTE18 <- lmer(log(P_ppm+1) ~ TSF + ESD + (1|Location), data=SoilJuly18, REML = FALSE)
Mult_P18_TSF <- glht(PTE18, linfct=mcp(TSF="Tukey"))
summary(Mult_P18_TSF)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(P_ppm + 1) ~ TSF + ESD + (1 | Location), data = SoilJuly18,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1yr - 2yr - 1yr < == 0 -0.06963 0.10030 -0.694 0.7652
## Unburned - 1yr < == 0 -0.20194 0.08592 -2.350 0.0485 *
## Unburned - 1yr - 2yr == 0 -0.13232 0.08320 -1.590 0.2476
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
Mult_P18_ESD <- glht(PTE18, linfct=mcp(ESD="Tukey"))
summary(Mult_P18_ESD)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(P_ppm + 1) ~ TSF + ESD + (1 | Location), data = SoilJuly18,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 -0.30770 0.13657 -2.253 0.1431
## SaLo - Clayey == 0 -0.11112 0.10526 -1.056 0.8125
## Sandy - Clayey == 0 -0.57677 0.09760 -5.910 <0.001 ***
## ThCl - Clayey == 0 0.20348 0.26721 0.762 0.9347
## SaLo - Loamy == 0 0.19658 0.13267 1.482 0.5463
## Sandy - Loamy == 0 -0.26907 0.12349 -2.179 0.1685
## ThCl - Loamy == 0 0.51118 0.27350 1.869 0.3077
## Sandy - SaLo == 0 -0.46565 0.08652 -5.382 <0.001 ***
## ThCl - SaLo == 0 0.31460 0.25381 1.240 0.7053
## ThCl - Sandy == 0 0.78025 0.25752 3.030 0.0177 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
#2018
PTE18 <- lmer(log(P_ppm+1) ~ TSF + ESD + (1|Location), data=SoilJuly18, REML = FALSE)
Mult_P18_TSF <- glht(PTE18, linfct=mcp(TSF="Tukey"))
summary(Mult_P18_TSF)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(P_ppm + 1) ~ TSF + ESD + (1 | Location), data = SoilJuly18,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1yr - 2yr - 1yr < == 0 -0.06963 0.10030 -0.694 0.7652
## Unburned - 1yr < == 0 -0.20194 0.08592 -2.350 0.0485 *
## Unburned - 1yr - 2yr == 0 -0.13232 0.08320 -1.590 0.2475
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
Mult_P18_ESD <- glht(PTE18, linfct=mcp(ESD="Tukey"))
summary(Mult_P18_ESD)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(P_ppm + 1) ~ TSF + ESD + (1 | Location), data = SoilJuly18,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 -0.30770 0.13657 -2.253 0.1432
## SaLo - Clayey == 0 -0.11112 0.10526 -1.056 0.8125
## Sandy - Clayey == 0 -0.57677 0.09760 -5.910 <0.001 ***
## ThCl - Clayey == 0 0.20348 0.26721 0.762 0.9347
## SaLo - Loamy == 0 0.19658 0.13267 1.482 0.5463
## Sandy - Loamy == 0 -0.26907 0.12349 -2.179 0.1685
## ThCl - Loamy == 0 0.51118 0.27350 1.869 0.3078
## Sandy - SaLo == 0 -0.46565 0.08652 -5.382 <0.001 ***
## ThCl - SaLo == 0 0.31460 0.25381 1.240 0.7053
## ThCl - Sandy == 0 0.78025 0.25752 3.030 0.0176 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
K
Knull <- lmer(log(K_ppm+1) ~ 1 + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
## boundary (singular) fit: see ?isSingular
KT <- lmer(log(K_ppm+1) ~ TSF + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
## boundary (singular) fit: see ?isSingular
KE <- lmer(log(K_ppm+1) ~ ESD + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
## boundary (singular) fit: see ?isSingular
KTE <- lmer(log(K_ppm+1) ~ TSF + ESD + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
## boundary (singular) fit: see ?isSingular
anova(Knull, KT, KE, KTE)
## Data: SoilJuly181920
## Models:
## Knull: log(K_ppm + 1) ~ 1 + (1 | Location/Year)
## KT: log(K_ppm + 1) ~ TSF + (1 | Location/Year)
## KE: log(K_ppm + 1) ~ ESD + (1 | Location/Year)
## KTE: log(K_ppm + 1) ~ TSF + ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Knull 4 204.54 219.20 -98.272 196.54
## KT 8 200.31 229.61 -92.155 184.31 12.235 4 0.015684 *
## KE 8 145.69 174.99 -64.843 129.69 54.623 0
## KTE 12 138.73 182.69 -57.366 114.73 14.954 4 0.004797 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
anova(Knull, KT)
## Data: SoilJuly181920
## Models:
## Knull: log(K_ppm + 1) ~ 1 + (1 | Location/Year)
## KT: log(K_ppm + 1) ~ TSF + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Knull 4 204.54 219.20 -98.272 196.54
## KT 8 200.31 229.61 -92.155 184.31 12.235 4 0.01568 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
anova(Knull, KE)
## Data: SoilJuly181920
## Models:
## Knull: log(K_ppm + 1) ~ 1 + (1 | Location/Year)
## KE: log(K_ppm + 1) ~ ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Knull 4 204.54 219.20 -98.272 196.54
## KE 8 145.69 174.99 -64.843 129.69 66.858 4 1.044e-13 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
anova(Knull, KTE)
## Data: SoilJuly181920
## Models:
## Knull: log(K_ppm + 1) ~ 1 + (1 | Location/Year)
## KTE: log(K_ppm + 1) ~ TSF + ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Knull 4 204.54 219.20 -98.272 196.54
## KTE 12 138.73 182.69 -57.366 114.73 81.812 8 2.109e-14 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Mult_K_TSF <- glht(KTE, linfct=mcp(TSF="Tukey"))
summary(Mult_K_TSF)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(K_ppm + 1) ~ TSF + ESD + (1 | Location/Year),
## data = SoilJuly181920, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1yr - 2yr - 1yr < == 0 0.0009085 0.0474769 0.019 1.00000
## 2yr - 3yr - 1yr < == 0 -0.0327548 0.0542782 -0.603 0.97367
## 3yr - 4yr - 1yr < == 0 0.0435284 0.0713872 0.610 0.97265
## Unburned - 1yr < == 0 -0.1705685 0.0494040 -3.453 0.00479 **
## 2yr - 3yr - 1yr - 2yr == 0 -0.0336633 0.0541995 -0.621 0.97075
## 3yr - 4yr - 1yr - 2yr == 0 0.0426199 0.0711383 0.599 0.97436
## Unburned - 1yr - 2yr == 0 -0.1714771 0.0494218 -3.470 0.00448 **
## 3yr - 4yr - 2yr - 3yr == 0 0.0762832 0.0739279 1.032 0.83607
## Unburned - 2yr - 3yr == 0 -0.1378138 0.0583869 -2.360 0.12180
## Unburned - 3yr - 4yr == 0 -0.2140970 0.0759859 -2.818 0.03722 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
Mult_K_ESD <- glht(KTE, linfct=mcp(ESD="Tukey"))
summary(Mult_K_ESD)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(K_ppm + 1) ~ TSF + ESD + (1 | Location/Year),
## data = SoilJuly181920, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 0.046007 0.068343 0.673 0.958
## SaLo - Clayey == 0 0.050670 0.052700 0.961 0.859
## Sandy - Clayey == 0 -0.274544 0.049149 -5.586 <0.001 ***
## ThCl - Clayey == 0 -0.225801 0.130225 -1.734 0.385
## SaLo - Loamy == 0 0.004663 0.067374 0.069 1.000
## Sandy - Loamy == 0 -0.320551 0.061550 -5.208 <0.001 ***
## ThCl - Loamy == 0 -0.271808 0.135696 -2.003 0.240
## Sandy - SaLo == 0 -0.325214 0.042818 -7.595 <0.001 ***
## ThCl - SaLo == 0 -0.276471 0.124584 -2.219 0.154
## ThCl - Sandy == 0 0.048743 0.124974 0.390 0.994
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
Ca
#Overall
Canull <- lmer(log(Ca_ppm+1) ~ 1 + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
CaT <- lmer(log(Ca_ppm+1) ~ TSF + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
CaE <- lmer(log(Ca_ppm+1) ~ ESD + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
CaTE <- lmer(log(Ca_ppm+1) ~ TSF + ESD + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
anova(Canull, CaT, CaE, CaTE)
## Data: SoilJuly181920
## Models:
## Canull: log(Ca_ppm + 1) ~ 1 + (1 | Location/Year)
## CaT: log(Ca_ppm + 1) ~ TSF + (1 | Location/Year)
## CaE: log(Ca_ppm + 1) ~ ESD + (1 | Location/Year)
## CaTE: log(Ca_ppm + 1) ~ TSF + ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Canull 4 432.47 447.13 -212.24 424.47
## CaT 8 438.20 467.50 -211.10 422.20 2.2766 4 0.685
## CaE 8 417.34 446.64 -200.67 401.34 20.8607 0
## CaTE 12 421.71 465.67 -198.85 397.71 3.6260 4 0.459
anova(Canull, CaT)
## Data: SoilJuly181920
## Models:
## Canull: log(Ca_ppm + 1) ~ 1 + (1 | Location/Year)
## CaT: log(Ca_ppm + 1) ~ TSF + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Canull 4 432.47 447.13 -212.24 424.47
## CaT 8 438.20 467.50 -211.10 422.20 2.2766 4 0.685
anova(Canull, CaE)
## Data: SoilJuly181920
## Models:
## Canull: log(Ca_ppm + 1) ~ 1 + (1 | Location/Year)
## CaE: log(Ca_ppm + 1) ~ ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Canull 4 432.47 447.13 -212.24 424.47
## CaE 8 417.34 446.64 -200.67 401.34 23.137 4 0.0001189 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
anova(Canull, CaTE)
## Data: SoilJuly181920
## Models:
## Canull: log(Ca_ppm + 1) ~ 1 + (1 | Location/Year)
## CaTE: log(Ca_ppm + 1) ~ TSF + ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Canull 4 432.47 447.13 -212.24 424.47
## CaTE 12 421.71 465.67 -198.85 397.71 26.763 8 0.0007767 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Mult_Ca_TSF <- glht(CaTE, linfct=mcp(TSF="Tukey"))
summary(Mult_Ca_TSF)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(Ca_ppm + 1) ~ TSF + ESD + (1 | Location/Year),
## data = SoilJuly181920, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1yr - 2yr - 1yr < == 0 0.007164 0.076702 0.093 1.000
## 2yr - 3yr - 1yr < == 0 0.023771 0.087342 0.272 0.999
## 3yr - 4yr - 1yr < == 0 -0.029954 0.114214 -0.262 0.999
## Unburned - 1yr < == 0 -0.120975 0.079239 -1.527 0.537
## 2yr - 3yr - 1yr - 2yr == 0 0.016607 0.087212 0.190 1.000
## 3yr - 4yr - 1yr - 2yr == 0 -0.037118 0.113799 -0.326 0.997
## Unburned - 1yr - 2yr == 0 -0.128139 0.079261 -1.617 0.478
## 3yr - 4yr - 2yr - 3yr == 0 -0.053725 0.118780 -0.452 0.991
## Unburned - 2yr - 3yr == 0 -0.144746 0.092786 -1.560 0.515
## Unburned - 3yr - 4yr == 0 -0.091021 0.120144 -0.758 0.941
## (Adjusted p values reported -- single-step method)
Mult_Ca_ESD <- glht(CaTE, linfct=mcp(ESD="Tukey"))
summary(Mult_Ca_ESD)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(Ca_ppm + 1) ~ TSF + ESD + (1 | Location/Year),
## data = SoilJuly181920, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 0.34230 0.11273 3.037 0.01728 *
## SaLo - Clayey == 0 0.22611 0.08729 2.590 0.06264 .
## Sandy - Clayey == 0 0.25682 0.08049 3.191 0.01038 *
## ThCl - Clayey == 0 -0.46675 0.21525 -2.168 0.17155
## SaLo - Loamy == 0 -0.11619 0.11419 -1.017 0.83105
## Sandy - Loamy == 0 -0.08548 0.10172 -0.840 0.90814
## ThCl - Loamy == 0 -0.80905 0.22535 -3.590 0.00263 **
## Sandy - SaLo == 0 0.03071 0.07063 0.435 0.99151
## ThCl - SaLo == 0 -0.69286 0.20371 -3.401 0.00517 **
## ThCl - Sandy == 0 -0.72357 0.20538 -3.523 0.00316 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
Mg
#Overall
Mgnull <- lmer(log(Mg_ppm+1) ~ 1 + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
MgT <- lmer(log(Mg_ppm+1) ~ TSF + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
MgE <- lmer(log(Mg_ppm+1) ~ ESD + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
MgTE <- lmer(log(Mg_ppm+1) ~ TSF + ESD + (1|Location/Year), data=SoilJuly181920, REML = FALSE)
anova(Mgnull, MgT, MgE, MgTE)
## Data: SoilJuly181920
## Models:
## Mgnull: log(Mg_ppm + 1) ~ 1 + (1 | Location/Year)
## MgT: log(Mg_ppm + 1) ~ TSF + (1 | Location/Year)
## MgE: log(Mg_ppm + 1) ~ ESD + (1 | Location/Year)
## MgTE: log(Mg_ppm + 1) ~ TSF + ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Mgnull 4 240.69 255.34 -116.35 232.69
## MgT 8 243.52 272.83 -113.76 227.52 5.1702 4 0.2703
## MgE 8 226.90 256.20 -105.45 210.90 16.6245 0
## MgTE 12 230.91 274.86 -103.45 206.91 3.9926 4 0.4070
anova(Mgnull, MgT)
## Data: SoilJuly181920
## Models:
## Mgnull: log(Mg_ppm + 1) ~ 1 + (1 | Location/Year)
## MgT: log(Mg_ppm + 1) ~ TSF + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Mgnull 4 240.69 255.34 -116.35 232.69
## MgT 8 243.52 272.83 -113.76 227.52 5.1702 4 0.2703
anova(Mgnull, MgE)
## Data: SoilJuly181920
## Models:
## Mgnull: log(Mg_ppm + 1) ~ 1 + (1 | Location/Year)
## MgE: log(Mg_ppm + 1) ~ ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Mgnull 4 240.69 255.34 -116.35 232.69
## MgE 8 226.90 256.20 -105.45 210.90 21.795 4 0.0002202 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
anova(Mgnull, MgTE)
## Data: SoilJuly181920
## Models:
## Mgnull: log(Mg_ppm + 1) ~ 1 + (1 | Location/Year)
## MgTE: log(Mg_ppm + 1) ~ TSF + ESD + (1 | Location/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## Mgnull 4 240.69 255.34 -116.35 232.69
## MgTE 12 230.91 274.86 -103.45 206.91 25.787 8 0.001142 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Mult_Mg_TSF <- glht(MgTE, linfct=mcp(TSF="Tukey"))
summary(Mult_Mg_TSF)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(Mg_ppm + 1) ~ TSF + ESD + (1 | Location/Year),
## data = SoilJuly181920, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1yr - 2yr - 1yr < == 0 0.024008 0.055207 0.435 0.992
## 2yr - 3yr - 1yr < == 0 0.019006 0.063173 0.301 0.998
## 3yr - 4yr - 1yr < == 0 -0.039770 0.083214 -0.478 0.989
## Unburned - 1yr < == 0 -0.079791 0.057541 -1.387 0.628
## 2yr - 3yr - 1yr - 2yr == 0 -0.005002 0.063077 -0.079 1.000
## 3yr - 4yr - 1yr - 2yr == 0 -0.063778 0.082915 -0.769 0.937
## Unburned - 1yr - 2yr == 0 -0.103799 0.057564 -1.803 0.363
## 3yr - 4yr - 2yr - 3yr == 0 -0.058776 0.086079 -0.683 0.959
## Unburned - 2yr - 3yr == 0 -0.098797 0.068144 -1.450 0.587
## Unburned - 3yr - 4yr == 0 -0.040021 0.088790 -0.451 0.991
## (Adjusted p values reported -- single-step method)
Mult_Mg_ESD <- glht(MgTE, linfct=mcp(ESD="Tukey"))
summary(Mult_Mg_ESD)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(Mg_ppm + 1) ~ TSF + ESD + (1 | Location/Year),
## data = SoilJuly181920, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 -0.01398 0.08049 -0.174 0.9998
## SaLo - Clayey == 0 0.08347 0.06223 1.341 0.6388
## Sandy - Clayey == 0 -0.12214 0.05764 -2.119 0.1906
## ThCl - Clayey == 0 0.28875 0.15358 1.880 0.3009
## SaLo - Loamy == 0 0.09745 0.08069 1.208 0.7242
## Sandy - Loamy == 0 -0.10816 0.07257 -1.491 0.5394
## ThCl - Loamy == 0 0.30273 0.16047 1.886 0.2976
## Sandy - SaLo == 0 -0.20562 0.05043 -4.077 <0.001 ***
## ThCl - SaLo == 0 0.20528 0.14597 1.406 0.5956
## ThCl - Sandy == 0 0.41089 0.14686 2.798 0.0356 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
Ammonium
#Overall
NH4null <- lmer(log(NH4_ppm+1) ~ 1 + (1|Location/Month/Year), data=Soil181920, REML = FALSE)
## boundary (singular) fit: see ?isSingular
NH4T <- lmer(log(NH4_ppm+1) ~ TSF + (1|Location/Month/Year), data=Soil181920, REML = FALSE)
## boundary (singular) fit: see ?isSingular
NH4E <- lmer(log(NH4_ppm+1) ~ ESD + (1|Location/Month/Year), data=Soil181920, REML = FALSE)
## boundary (singular) fit: see ?isSingular
NH4TE <- lmer(log(NH4_ppm+1) ~ TSF + ESD + (1|Location/Month/Year), data=Soil181920, REML = FALSE)
## boundary (singular) fit: see ?isSingular
anova(NH4null, NH4T, NH4E, NH4TE)
## Data: Soil181920
## Models:
## NH4null: log(NH4_ppm + 1) ~ 1 + (1 | Location/Month/Year)
## NH4T: log(NH4_ppm + 1) ~ TSF + (1 | Location/Month/Year)
## NH4E: log(NH4_ppm + 1) ~ ESD + (1 | Location/Month/Year)
## NH4TE: log(NH4_ppm + 1) ~ TSF + ESD + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## NH4null 5 1757.2 1782.4 -873.60 1747.2
## NH4T 9 1753.4 1798.8 -867.70 1735.4 11.796 4 0.01893 *
## NH4E 9 1756.6 1802.1 -869.32 1738.6 0.000 0
## NH4TE 13 1752.0 1817.6 -862.99 1726.0 12.668 4 0.01302 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
anova(NH4null, NH4T)
## Data: Soil181920
## Models:
## NH4null: log(NH4_ppm + 1) ~ 1 + (1 | Location/Month/Year)
## NH4T: log(NH4_ppm + 1) ~ TSF + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## NH4null 5 1757.2 1782.4 -873.6 1747.2
## NH4T 9 1753.4 1798.8 -867.7 1735.4 11.796 4 0.01893 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
anova(NH4null, NH4E)
## Data: Soil181920
## Models:
## NH4null: log(NH4_ppm + 1) ~ 1 + (1 | Location/Month/Year)
## NH4E: log(NH4_ppm + 1) ~ ESD + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## NH4null 5 1757.2 1782.4 -873.60 1747.2
## NH4E 9 1756.6 1802.1 -869.32 1738.6 8.5564 4 0.0732 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
anova(NH4null, NH4TE)
## Data: Soil181920
## Models:
## NH4null: log(NH4_ppm + 1) ~ 1 + (1 | Location/Month/Year)
## NH4TE: log(NH4_ppm + 1) ~ TSF + ESD + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## NH4null 5 1757.2 1782.4 -873.60 1747.2
## NH4TE 13 1752.0 1817.6 -862.99 1726.0 21.224 8 0.006574 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Mult_NH4_TSF <- glht(NH4TE, linfct=mcp(TSF="Tukey"))
summary(Mult_NH4_TSF)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(NH4_ppm + 1) ~ TSF + ESD + (1 | Location/Month/Year),
## data = Soil181920, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1yr - 2yr - 1yr < == 0 -0.10003 0.03768 -2.655 0.0583 .
## 2yr - 3yr - 1yr < == 0 -0.07894 0.04361 -1.810 0.3576
## 3yr - 4yr - 1yr < == 0 -0.17222 0.05831 -2.954 0.0248 *
## Unburned - 1yr < == 0 -0.02274 0.04007 -0.567 0.9788
## 2yr - 3yr - 1yr - 2yr == 0 0.02108 0.04349 0.485 0.9882
## 3yr - 4yr - 1yr - 2yr == 0 -0.07220 0.05801 -1.245 0.7169
## Unburned - 1yr - 2yr == 0 0.07729 0.04010 1.928 0.2933
## 3yr - 4yr - 2yr - 3yr == 0 -0.09328 0.05943 -1.570 0.5067
## Unburned - 2yr - 3yr == 0 0.05620 0.04863 1.156 0.7693
## Unburned - 3yr - 4yr == 0 0.14949 0.06410 2.332 0.1294
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
Mult_NH4_ESD <- glht(NH4TE, linfct=mcp(ESD="Tukey"))
summary(Mult_NH4_ESD)
##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(NH4_ppm + 1) ~ TSF + ESD + (1 | Location/Month/Year),
## data = Soil181920, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 0.074733 0.055540 1.346 0.6351
## SaLo - Clayey == 0 -0.041609 0.043120 -0.965 0.8564
## Sandy - Clayey == 0 -0.005689 0.039594 -0.144 0.9999
## ThCl - Clayey == 0 0.201403 0.106051 1.899 0.2903
## SaLo - Loamy == 0 -0.116342 0.056489 -2.060 0.2145
## Sandy - Loamy == 0 -0.080422 0.050086 -1.606 0.4630
## ThCl - Loamy == 0 0.126670 0.111078 1.140 0.7639
## Sandy - SaLo == 0 0.035920 0.034785 1.033 0.8232
## ThCl - SaLo == 0 0.243012 0.100206 2.425 0.0954 .
## ThCl - Sandy == 0 0.207092 0.101078 2.049 0.2190
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)
Nitrate
Soil Moisture
Decomposition
Filler