Forage Quality and Grazer Selection Dissertation and Models Scratchpad
Jonathan Spiess
April 14, 2021
KgHa
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular## Data: HRECSIOmit
## Models:
## KgHaNull: log(KgHa + 1) ~ 1 + (1 | Location/Month/Year)
## KgHaG: log(KgHa + 1) ~ Treatment + (1 | Location/Month/Year)
## KgHaT: log(KgHa + 1) ~ TSF + (1 | Location/Month/Year)
## KgHaE: log(KgHa + 1) ~ ESD + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## KgHaNull 5 4031.2 4059.2 -2010.6 4021.2
## KgHaG 6 4032.3 4065.9 -2010.2 4020.3 0.8545 1 0.3553
## KgHaT 9 3663.1 3713.5 -1822.6 3645.1 375.2253 3 <2e-16 ***
## KgHaE 9 3990.3 4040.7 -1986.2 3972.3 0.0000 0
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Data: HRECSIOmit
## Models:
## KgHaNull: log(KgHa + 1) ~ 1 + (1 | Location/Month/Year)
## KgHaT: log(KgHa + 1) ~ TSF + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## KgHaNull 5 4031.2 4059.2 -2010.6 4021.2
## KgHaT 9 3663.1 3713.5 -1822.6 3645.1 376.08 4 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Data: HRECSIOmit
## Models:
## KgHaNull: log(KgHa + 1) ~ 1 + (1 | Location/Month/Year)
## KgHaG: log(KgHa + 1) ~ Treatment + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## KgHaNull 5 4031.2 4059.2 -2010.6 4021.2
## KgHaG 6 4032.3 4065.9 -2010.2 4020.3 0.8545 1 0.3553## Data: HRECSIOmit
## Models:
## KgHaNull: log(KgHa + 1) ~ 1 + (1 | Location/Month/Year)
## KgHaE: log(KgHa + 1) ~ ESD + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## KgHaNull 5 4031.2 4059.2 -2010.6 4021.2
## KgHaE 9 3990.3 4040.7 -1986.2 3972.3 48.909 4 6.099e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(KgHa + 1) ~ TSF + (1 | Location/Month/Year),
## data = HRECSIOmit, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 0.54567 0.03934 13.871 < 0.001 ***
## 2to3 - RB == 0 0.70380 0.04565 15.418 < 0.001 ***
## 3plus - RB == 0 0.78886 0.05974 13.205 < 0.001 ***
## NYB - RB == 0 0.52298 0.03624 14.429 < 0.001 ***
## 2to3 - 1to2 == 0 0.15813 0.04682 3.378 0.00621 **
## 3plus - 1to2 == 0 0.24319 0.06053 4.018 < 0.001 ***
## NYB - 1to2 == 0 -0.02269 0.04078 -0.556 0.98029
## 3plus - 2to3 == 0 0.08507 0.06245 1.362 0.64221
## NYB - 2to3 == 0 -0.18081 0.04869 -3.713 0.00178 **
## NYB - 3plus == 0 -0.26588 0.06364 -4.178 < 0.001 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(KgHa + 1) ~ Treatment + (1 | Location/Month/Year),
## data = HRECSIOmit, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Sheep - Cattle == 0 -0.1088 0.1135 -0.959 0.338
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(KgHa + 1) ~ ESD + (1 | Location/Month/Year),
## data = HRECSIOmit, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 0.06305 0.06051 1.042 0.81677
## Saline Lowland - Clayey == 0 -0.09777 0.04466 -2.189 0.16200
## Sandy - Clayey == 0 -0.20672 0.04205 -4.916 < 0.001 ***
## Thin Claypan - Clayey == 0 -0.45591 0.12277 -3.714 0.00162 **
## Saline Lowland - Loamy == 0 -0.16081 0.06019 -2.672 0.04954 *
## Sandy - Loamy == 0 -0.26976 0.05316 -5.075 < 0.001 ***
## Thin Claypan - Loamy == 0 -0.51896 0.12826 -4.046 < 0.001 ***
## Sandy - Saline Lowland == 0 -0.10895 0.03656 -2.980 0.02023 *
## Thin Claypan - Saline Lowland == 0 -0.35815 0.11742 -3.050 0.01615 *
## Thin Claypan - Sandy == 0 -0.24920 0.11819 -2.109 0.19216
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(KgHa + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit17,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## NYB - RB == 0 0.67086 0.08263 8.119 4.44e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(KgHa + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit17,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## NYB - RB == 0 0.21208 0.08145 2.604 0.00922 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(KgHa + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit18,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -0.08556 0.08740 -0.979 0.58865
## NYB - RB == 0 0.23026 0.07529 3.058 0.00634 **
## NYB - 1to2 == 0 0.31582 0.07688 4.108 < 0.001 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(KgHa + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit18,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 0.56662 0.10863 5.216 <1e-05 ***
## NYB - RB == 0 0.58637 0.09293 6.310 <1e-05 ***
## NYB - 1to2 == 0 0.01975 0.09293 0.212 0.975
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(KgHa + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit19,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -0.01770 0.08806 -0.201 0.9971
## 2to3 - RB == 0 -0.13261 0.08947 -1.482 0.4483
## NYB - RB == 0 0.13377 0.08781 1.523 0.4234
## 2to3 - 1to2 == 0 -0.11491 0.08782 -1.308 0.5574
## NYB - 1to2 == 0 0.15147 0.08688 1.743 0.3011
## NYB - 2to3 == 0 0.26638 0.08827 3.018 0.0139 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(KgHa + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit19,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 0.16255 0.09310 1.746 0.2999
## 2to3 - RB == 0 0.37507 0.09310 4.029 <0.001 ***
## NYB - RB == 0 0.24638 0.09013 2.734 0.0315 *
## 2to3 - 1to2 == 0 0.21253 0.09432 2.253 0.1095
## NYB - 1to2 == 0 0.08384 0.09154 0.916 0.7964
## NYB - 2to3 == 0 -0.12869 0.09154 -1.406 0.4955
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(KgHa + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit20,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 0.98254 0.08417 11.673 <0.001 ***
## 2to3 - RB == 0 1.05309 0.08330 12.641 <0.001 ***
## 3plus - RB == 0 0.88961 0.08464 10.510 <0.001 ***
## 2to3 - 1to2 == 0 0.07054 0.08443 0.836 0.838
## 3plus - 1to2 == 0 -0.09294 0.08578 -1.083 0.700
## 3plus - 2to3 == 0 -0.16348 0.08418 -1.942 0.210
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(KgHa + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit20,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 1.48437 0.09490 15.641 <1e-05 ***
## 2to3 - RB == 0 1.55129 0.09639 16.095 <1e-05 ***
## 3plus - RB == 0 1.51967 0.09639 15.767 <1e-05 ***
## 2to3 - 1to2 == 0 0.06693 0.09803 0.683 0.904
## 3plus - 1to2 == 0 0.03531 0.09803 0.360 0.984
## 3plus - 2to3 == 0 -0.03162 0.09931 -0.318 0.989
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)ADF
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular## Data: HRECSIOmit
## Models:
## ADFNull: ADF ~ 1 + (1 | Location/Month/Year)
## ADFG: ADF ~ Treatment + (1 | Location/Month/Year)
## ADFT: ADF ~ TSF + (1 | Location/Month/Year)
## ADFE: ADF ~ ESD + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## ADFNull 5 9757.3 9785.3 -4873.7 9747.3
## ADFG 6 9757.5 9791.1 -4872.8 9745.5 1.8206 1 0.1772
## ADFT 9 9402.6 9452.9 -4692.3 9384.6 360.9582 3 <2e-16 ***
## ADFE 9 9749.8 9800.2 -4865.9 9731.8 0.0000 0
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Data: HRECSIOmit
## Models:
## ADFNull: ADF ~ 1 + (1 | Location/Month/Year)
## ADFT: ADF ~ TSF + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## ADFNull 5 9757.3 9785.3 -4873.7 9747.3
## ADFT 9 9402.6 9452.9 -4692.3 9384.6 362.78 4 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Data: HRECSIOmit
## Models:
## ADFNull: ADF ~ 1 + (1 | Location/Month/Year)
## ADFG: ADF ~ Treatment + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## ADFNull 5 9757.3 9785.3 -4873.7 9747.3
## ADFG 6 9757.5 9791.1 -4872.8 9745.5 1.8206 1 0.1772## Data: HRECSIOmit
## Models:
## ADFNull: ADF ~ 1 + (1 | Location/Month/Year)
## ADFE: ADF ~ ESD + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## ADFNull 5 9757.3 9785.3 -4873.7 9747.3
## ADFE 9 9749.8 9800.2 -4865.9 9731.8 15.508 4 0.003755 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADF ~ TSF + (1 | Location/Month/Year), data = HRECSIOmit,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 2.4017 0.1653 14.532 <0.001 ***
## 2to3 - RB == 0 3.0023 0.1917 15.659 <0.001 ***
## 3plus - RB == 0 3.2299 0.2509 12.872 <0.001 ***
## NYB - RB == 0 1.9565 0.1523 12.850 <0.001 ***
## 2to3 - 1to2 == 0 0.6006 0.1966 3.054 0.0182 *
## 3plus - 1to2 == 0 0.8281 0.2542 3.257 0.0096 **
## NYB - 1to2 == 0 -0.4452 0.1713 -2.599 0.0672 .
## 3plus - 2to3 == 0 0.2276 0.2623 0.868 0.9051
## NYB - 2to3 == 0 -1.0458 0.2045 -5.114 <0.001 ***
## NYB - 3plus == 0 -1.2733 0.2673 -4.764 <0.001 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADF ~ Treatment + (1 | Location/Month/Year), data = HRECSIOmit,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Sheep - Cattle == 0 0.7423 0.5400 1.375 0.169
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADF ~ ESD + (1 | Location/Month/Year), data = HRECSIOmit,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 -0.1705 0.2555 -0.668 0.9580
## Saline Lowland - Clayey == 0 -0.3549 0.1885 -1.882 0.2967
## Sandy - Clayey == 0 -0.4575 0.1776 -2.577 0.0641 .
## Thin Claypan - Clayey == 0 1.0163 0.5183 1.961 0.2572
## Saline Lowland - Loamy == 0 -0.1843 0.2541 -0.725 0.9437
## Sandy - Loamy == 0 -0.2870 0.2244 -1.279 0.6763
## Thin Claypan - Loamy == 0 1.1868 0.5414 2.192 0.1611
## Sandy - Saline Lowland == 0 -0.1027 0.1543 -0.665 0.9585
## Thin Claypan - Saline Lowland == 0 1.3712 0.4957 2.766 0.0381 *
## Thin Claypan - Sandy == 0 1.4738 0.4989 2.954 0.0219 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit17,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## NYB - RB == 0 2.9514 0.3791 7.786 6.88e-15 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit17,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## NYB - RB == 0 0.7642 0.3877 1.971 0.0487 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit18,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -0.1657 0.3635 -0.456 0.8913
## NYB - RB == 0 0.7339 0.3137 2.339 0.0502 .
## NYB - 1to2 == 0 0.8995 0.3205 2.807 0.0138 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit18,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 1.79676 0.27333 6.574 <1e-06 ***
## NYB - RB == 0 1.77203 0.23391 7.576 <1e-06 ***
## NYB - 1to2 == 0 -0.02473 0.23391 -0.106 0.994
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit19,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 1.06154 0.40897 2.596 0.0468 *
## 2to3 - RB == 0 0.89684 0.41553 2.158 0.1351
## NYB - RB == 0 0.85235 0.40776 2.090 0.1561
## 2to3 - 1to2 == 0 -0.16470 0.40779 -0.404 0.9777
## NYB - 1to2 == 0 -0.20919 0.40351 -0.518 0.9547
## NYB - 2to3 == 0 -0.04449 0.41000 -0.109 0.9995
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit19,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 0.48533 0.31271 1.552 0.406
## 2to3 - RB == 0 0.52966 0.31271 1.694 0.327
## NYB - RB == 0 0.50522 0.30269 1.669 0.340
## 2to3 - 1to2 == 0 0.04433 0.31671 0.140 0.999
## NYB - 1to2 == 0 0.01989 0.30739 0.065 1.000
## NYB - 2to3 == 0 -0.02444 0.30739 -0.080 1.000
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit20,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 5.18445 0.45988 11.273 <1e-04 ***
## 2to3 - RB == 0 5.75131 0.45512 12.637 <1e-04 ***
## 3plus - RB == 0 5.20893 0.46244 11.264 <1e-04 ***
## 2to3 - 1to2 == 0 0.56686 0.46126 1.229 0.608
## 3plus - 1to2 == 0 0.02448 0.46867 0.052 1.000
## 3plus - 2to3 == 0 -0.54238 0.45992 -1.179 0.640
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit20,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 5.2502 0.4361 12.040 <1e-04 ***
## 2to3 - RB == 0 5.3856 0.4429 12.160 <1e-04 ***
## 3plus - RB == 0 4.7556 0.4429 10.737 <1e-04 ***
## 2to3 - 1to2 == 0 0.1354 0.4504 0.301 0.991
## 3plus - 1to2 == 0 -0.4946 0.4504 -1.098 0.691
## 3plus - 2to3 == 0 -0.6300 0.4564 -1.380 0.512
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)NDF
## boundary (singular) fit: see ?isSingular## Data: HRECSIOmit
## Models:
## NDFNull: NDF ~ 1 + (1 | Location/Month/Year)
## NDFG: NDF ~ Treatment + (1 | Location/Month/Year)
## NDFT: NDF ~ TSF + (1 | Location/Month/Year)
## NDFE: NDF ~ ESD + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## NDFNull 5 11270 11298 -5629.9 11260
## NDFG 6 11265 11299 -5626.6 11253 6.5678 1 0.01038 *
## NDFT 9 11043 11094 -5512.5 11025 228.1075 3 < 2e-16 ***
## NDFE 9 11258 11308 -5619.9 11240 0.0000 0
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Data: HRECSIOmit
## Models:
## NDFNull: NDF ~ 1 + (1 | Location/Month/Year)
## NDFT: NDF ~ TSF + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## NDFNull 5 11270 11298 -5629.9 11260
## NDFT 9 11043 11094 -5512.5 11025 234.68 4 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Data: HRECSIOmit
## Models:
## NDFNull: NDF ~ 1 + (1 | Location/Month/Year)
## NDFG: NDF ~ Treatment + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## NDFNull 5 11270 11298 -5629.9 11260
## NDFG 6 11265 11299 -5626.6 11253 6.5678 1 0.01038 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Data: HRECSIOmit
## Models:
## NDFNull: NDF ~ 1 + (1 | Location/Month/Year)
## NDFE: NDF ~ ESD + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## NDFNull 5 11270 11298 -5629.9 11260
## NDFE 9 11258 11308 -5619.9 11240 19.887 4 0.0005258 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = NDF ~ TSF + (1 | Location/Month/Year), data = HRECSIOmit,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 2.95505 0.24886 11.874 < 0.001 ***
## 2to3 - RB == 0 3.65679 0.28907 12.650 < 0.001 ***
## 3plus - RB == 0 3.69400 0.37875 9.753 < 0.001 ***
## NYB - RB == 0 2.23734 0.22976 9.738 < 0.001 ***
## 2to3 - 1to2 == 0 0.70174 0.29622 2.369 0.11834
## 3plus - 1to2 == 0 0.73895 0.38349 1.927 0.29319
## NYB - 1to2 == 0 -0.71771 0.25889 -2.772 0.04200 *
## 3plus - 2to3 == 0 0.03721 0.39521 0.094 0.99998
## NYB - 2to3 == 0 -1.41945 0.30950 -4.586 < 0.001 ***
## NYB - 3plus == 0 -1.45666 0.40469 -3.599 0.00265 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = NDF ~ Treatment + (1 | Location/Month/Year), data = HRECSIOmit,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Sheep - Cattle == 0 2.1264 0.7348 2.894 0.00381 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = NDF ~ ESD + (1 | Location/Month/Year), data = HRECSIOmit,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 -0.87073 0.37278 -2.336 0.11641
## Saline Lowland - Clayey == 0 -0.60560 0.27480 -2.204 0.15695
## Sandy - Clayey == 0 -0.77698 0.25866 -3.004 0.01884 *
## Thin Claypan - Clayey == 0 1.59015 0.75526 2.105 0.19342
## Saline Lowland - Loamy == 0 0.26513 0.37146 0.714 0.94674
## Sandy - Loamy == 0 0.09376 0.32738 0.286 0.99829
## Thin Claypan - Loamy == 0 2.46089 0.78946 3.117 0.01308 *
## Sandy - Saline Lowland == 0 -0.17137 0.22498 -0.762 0.93331
## Thin Claypan - Saline Lowland == 0 2.19575 0.72188 3.042 0.01674 *
## Thin Claypan - Sandy == 0 2.36713 0.72682 3.257 0.00831 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = NDF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit17,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## NYB - RB == 0 3.7738 0.6527 5.782 7.37e-09 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = NDF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit17,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## NYB - RB == 0 1.4118 0.4477 3.153 0.00162 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = NDF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit18,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -1.5493 0.6695 -2.314 0.0535 .
## NYB - RB == 0 -0.1915 0.5777 -0.331 0.9409
## NYB - 1to2 == 0 1.3578 0.5900 2.301 0.0552 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = NDF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit18,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 2.1674 0.4274 5.071 < 1e-04 ***
## NYB - RB == 0 1.4920 0.3658 4.079 0.000147 ***
## NYB - 1to2 == 0 -0.6754 0.3658 -1.847 0.153260
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = NDF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit19,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 2.4446 0.6563 3.725 0.00101 **
## 2to3 - RB == 0 1.2629 0.6668 1.894 0.23047
## NYB - RB == 0 1.5855 0.6543 2.423 0.07272 .
## 2to3 - 1to2 == 0 -1.1817 0.6544 -1.806 0.27049
## NYB - 1to2 == 0 -0.8591 0.6475 -1.327 0.54575
## NYB - 2to3 == 0 0.3226 0.6579 0.490 0.96127
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = NDF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit19,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 1.2281 0.5222 2.352 0.0866 .
## 2to3 - RB == 0 0.9528 0.5222 1.824 0.2617
## NYB - RB == 0 0.3593 0.5055 0.711 0.8929
## 2to3 - 1to2 == 0 -0.2753 0.5289 -0.521 0.9542
## NYB - 1to2 == 0 -0.8688 0.5133 -1.693 0.3274
## NYB - 2to3 == 0 -0.5935 0.5133 -1.156 0.6545
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = NDF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit20,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 6.9917 0.7286 9.596 <1e-04 ***
## 2to3 - RB == 0 8.2861 0.7210 11.492 <1e-04 ***
## 3plus - RB == 0 7.3028 0.7326 9.968 <1e-04 ***
## 2to3 - 1to2 == 0 1.2944 0.7308 1.771 0.287
## 3plus - 1to2 == 0 0.3112 0.7425 0.419 0.975
## 3plus - 2to3 == 0 -0.9833 0.7286 -1.349 0.531
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = NDF ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit20,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 5.4207 0.5434 9.976 <0.001 ***
## 2to3 - RB == 0 5.7580 0.5519 10.433 <0.001 ***
## 3plus - RB == 0 4.5696 0.5519 8.280 <0.001 ***
## 2to3 - 1to2 == 0 0.3373 0.5613 0.601 0.932
## 3plus - 1to2 == 0 -0.8511 0.5613 -1.516 0.428
## 3plus - 2to3 == 0 -1.1884 0.5686 -2.090 0.156
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)ADL
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular## Data: HRECSIOmit
## Models:
## ADLNull: ADL ~ 1 + (1 | Location/Month/Year)
## ADLG: ADL ~ Treatment + (1 | Location/Month/Year)
## ADLT: ADL ~ TSF + (1 | Location/Month/Year)
## ADLE: ADL ~ ESD + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## ADLNull 5 4360.5 4388.5 -2175.2 4350.5
## ADLG 6 4361.9 4395.5 -2174.9 4349.9 0.5766 1 0.4476
## ADLT 9 4093.5 4143.9 -2037.8 4075.5 274.3681 3 <2e-16 ***
## ADLE 9 4363.9 4414.3 -2173.0 4345.9 0.0000 0
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Data: HRECSIOmit
## Models:
## ADLNull: ADL ~ 1 + (1 | Location/Month/Year)
## ADLT: ADL ~ TSF + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## ADLNull 5 4360.5 4388.5 -2175.2 4350.5
## ADLT 9 4093.5 4143.9 -2037.8 4075.5 274.94 4 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Data: HRECSIOmit
## Models:
## ADLNull: ADL ~ 1 + (1 | Location/Month/Year)
## ADLG: ADL ~ Treatment + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## ADLNull 5 4360.5 4388.5 -2175.2 4350.5
## ADLG 6 4361.9 4395.5 -2174.9 4349.9 0.5766 1 0.4476## Data: HRECSIOmit
## Models:
## ADLNull: ADL ~ 1 + (1 | Location/Month/Year)
## ADLE: ADL ~ ESD + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## ADLNull 5 4360.5 4388.5 -2175.2 4350.5
## ADLE 9 4363.9 4414.3 -2173.0 4345.9 4.5632 4 0.3351##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADL ~ TSF + (1 | Location/Month/Year), data = HRECSIOmit,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 0.47566 0.04206 11.309 <0.001 ***
## 2to3 - RB == 0 0.60803 0.04912 12.379 <0.001 ***
## 3plus - RB == 0 0.50297 0.06469 7.776 <0.001 ***
## NYB - RB == 0 0.55765 0.03920 14.225 <0.001 ***
## 2to3 - 1to2 == 0 0.13237 0.05015 2.640 0.060 .
## 3plus - 1to2 == 0 0.02731 0.06529 0.418 0.993
## NYB - 1to2 == 0 0.08200 0.04444 1.845 0.337
## 3plus - 2to3 == 0 -0.10506 0.06696 -1.569 0.506
## NYB - 2to3 == 0 -0.05037 0.05349 -0.942 0.875
## NYB - 3plus == 0 0.05469 0.07004 0.781 0.933
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADL ~ Treatment + (1 | Location/Month/Year), data = HRECSIOmit,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Sheep - Cattle == 0 -0.2245 0.2944 -0.763 0.446
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADL ~ ESD + (1 | Location/Month/Year), data = HRECSIOmit,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 -0.113315 0.063995 -1.771 0.358
## Saline Lowland - Clayey == 0 -0.040513 0.047116 -0.860 0.900
## Sandy - Clayey == 0 -0.044928 0.044326 -1.014 0.831
## Thin Claypan - Clayey == 0 -0.186613 0.129452 -1.442 0.568
## Saline Lowland - Loamy == 0 0.072803 0.063884 1.140 0.762
## Sandy - Loamy == 0 0.068388 0.056188 1.217 0.715
## Thin Claypan - Loamy == 0 -0.073297 0.135399 -0.541 0.980
## Sandy - Saline Lowland == 0 -0.004415 0.038579 -0.114 1.000
## Thin Claypan - Saline Lowland == 0 -0.146100 0.123647 -1.182 0.737
## Thin Claypan - Sandy == 0 -0.141685 0.124535 -1.138 0.763
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADL ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit17,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## NYB - RB == 0 0.4590 0.1168 3.93 8.51e-05 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADL ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit17,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## NYB - RB == 0 0.1179 0.1102 1.069 0.285
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADL ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit18,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 0.27097 0.09203 2.944 0.00908 **
## NYB - RB == 0 0.59388 0.07945 7.475 < 0.001 ***
## NYB - 1to2 == 0 0.32291 0.08115 3.979 < 0.001 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADL ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit18,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 0.46718 0.08526 5.480 <1e-04 ***
## NYB - RB == 0 0.58914 0.07296 8.074 <1e-04 ***
## NYB - 1to2 == 0 0.12196 0.07296 1.672 0.215
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADL ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit19,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 0.38437 0.13506 2.846 0.0228 *
## 2to3 - RB == 0 0.55543 0.13723 4.047 <0.001 ***
## NYB - RB == 0 0.54235 0.13465 4.028 <0.001 ***
## 2to3 - 1to2 == 0 0.17106 0.13466 1.270 0.5819
## NYB - 1to2 == 0 0.15797 0.13327 1.185 0.6361
## NYB - 2to3 == 0 -0.01308 0.13541 -0.097 0.9997
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADL ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit19,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 0.42182 0.11602 3.636 0.00145 **
## 2to3 - RB == 0 0.62798 0.11602 5.413 < 0.001 ***
## NYB - RB == 0 0.61186 0.11230 5.449 < 0.001 ***
## 2to3 - 1to2 == 0 0.20617 0.11749 1.755 0.29527
## NYB - 1to2 == 0 0.19004 0.11403 1.667 0.34144
## NYB - 2to3 == 0 -0.01612 0.11403 -0.141 0.99899
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADL ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit20,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 0.70337 0.10205 6.892 <1e-04 ***
## 2to3 - RB == 0 0.63793 0.10099 6.317 <1e-04 ***
## 3plus - RB == 0 0.57536 0.10262 5.607 <1e-04 ***
## 2to3 - 1to2 == 0 -0.06544 0.10236 -0.639 0.919
## 3plus - 1to2 == 0 -0.12801 0.10400 -1.231 0.607
## 3plus - 2to3 == 0 -0.06257 0.10206 -0.613 0.928
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = ADL ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit20,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 0.82340 0.07718 10.668 <0.001 ***
## 2to3 - RB == 0 0.86985 0.07838 11.098 <0.001 ***
## 3plus - RB == 0 0.68731 0.07838 8.769 <0.001 ***
## 2to3 - 1to2 == 0 0.04645 0.07974 0.582 0.937
## 3plus - 1to2 == 0 -0.13609 0.07974 -1.707 0.320
## 3plus - 2to3 == 0 -0.18253 0.08076 -2.260 0.107
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)CP
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular## Data: HRECSIOmit
## Models:
## CPNull: log(CP + 1) ~ 1 + (1 | Location/Month/Year)
## CPG: log(CP + 1) ~ Treatment + (1 | Location/Month/Year)
## CPT: log(CP + 1) ~ TSF + (1 | Location/Month/Year)
## CPE: log(CP + 1) ~ ESD + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## CPNull 5 -184.60 -156.62 97.301 -194.60
## CPG 6 -186.36 -152.78 99.180 -198.36 3.7589 1 0.05253 .
## CPT 9 -584.42 -534.05 301.211 -602.42 404.0611 3 < 2e-16 ***
## CPE 9 -194.74 -144.37 106.368 -212.74 0.0000 0
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Data: HRECSIOmit
## Models:
## CPNull: log(CP + 1) ~ 1 + (1 | Location/Month/Year)
## CPT: log(CP + 1) ~ TSF + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## CPNull 5 -184.60 -156.62 97.301 -194.60
## CPT 9 -584.42 -534.05 301.211 -602.42 407.82 4 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Data: HRECSIOmit
## Models:
## CPNull: log(CP + 1) ~ 1 + (1 | Location/Month/Year)
## CPG: log(CP + 1) ~ Treatment + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## CPNull 5 -184.60 -156.62 97.301 -194.60
## CPG 6 -186.36 -152.78 99.180 -198.36 3.7589 1 0.05253 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Data: HRECSIOmit
## Models:
## CPNull: log(CP + 1) ~ 1 + (1 | Location/Month/Year)
## CPE: log(CP + 1) ~ ESD + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## CPNull 5 -184.60 -156.62 97.301 -194.60
## CPE 9 -194.74 -144.37 106.368 -212.74 18.134 4 0.001162 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(CP + 1) ~ TSF + (1 | Location/Month/Year),
## data = HRECSIOmit, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -0.18625 0.01336 -13.937 <0.001 ***
## 2to3 - RB == 0 -0.26520 0.01553 -17.081 <0.001 ***
## 3plus - RB == 0 -0.32898 0.02035 -16.168 <0.001 ***
## NYB - RB == 0 -0.07181 0.01234 -5.818 <0.001 ***
## 2to3 - 1to2 == 0 -0.07895 0.01591 -4.963 <0.001 ***
## 3plus - 1to2 == 0 -0.14273 0.02060 -6.929 <0.001 ***
## NYB - 1to2 == 0 0.11444 0.01391 8.227 <0.001 ***
## 3plus - 2to3 == 0 -0.06378 0.02122 -3.005 0.021 *
## NYB - 2to3 == 0 0.19339 0.01664 11.626 <0.001 ***
## NYB - 3plus == 0 0.25718 0.02176 11.821 <0.001 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(CP + 1) ~ Treatment + (1 | Location/Month/Year),
## data = HRECSIOmit, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Sheep - Cattle == 0 -0.11318 0.05576 -2.03 0.0424 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(CP + 1) ~ ESD + (1 | Location/Month/Year),
## data = HRECSIOmit, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 0.003903 0.020932 0.186 1.000
## Saline Lowland - Clayey == 0 -0.030607 0.015429 -1.984 0.246
## Sandy - Clayey == 0 0.020272 0.014522 1.396 0.599
## Thin Claypan - Clayey == 0 0.045904 0.042402 1.083 0.795
## Saline Lowland - Loamy == 0 -0.034511 0.020861 -1.654 0.429
## Sandy - Loamy == 0 0.016369 0.018382 0.890 0.888
## Thin Claypan - Loamy == 0 0.042000 0.044324 0.948 0.863
## Sandy - Saline Lowland == 0 0.050879 0.012632 4.028 <0.001 ***
## Thin Claypan - Saline Lowland == 0 0.076511 0.040525 1.888 0.294
## Thin Claypan - Sandy == 0 0.025631 0.040803 0.628 0.966
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(CP + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit17,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## NYB - RB == 0 0.00219 0.03879 0.056 0.955
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(CP + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit17,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## NYB - RB == 0 -0.06899 0.02938 -2.349 0.0188 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(CP + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit18,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -0.009242 0.036006 -0.257 0.964150
## NYB - RB == 0 0.121182 0.031078 3.899 0.000279 ***
## NYB - 1to2 == 0 0.130423 0.031744 4.109 0.000112 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(CP + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit18,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -0.163612 0.030911 -5.293 <1e-06 ***
## NYB - RB == 0 -0.158746 0.026452 -6.001 <1e-06 ***
## NYB - 1to2 == 0 0.004866 0.026452 0.184 0.981
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(CP + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit19,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -0.18151 0.03068 -5.917 <1e-06 ***
## 2to3 - RB == 0 -0.18415 0.03117 -5.908 <1e-06 ***
## NYB - RB == 0 0.01119 0.03059 0.366 0.983
## 2to3 - 1to2 == 0 -0.00264 0.03059 -0.086 1.000
## NYB - 1to2 == 0 0.19270 0.03027 6.367 <1e-06 ***
## NYB - 2to3 == 0 0.19534 0.03075 6.352 <1e-06 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(CP + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit19,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -0.05752 0.02596 -2.215 0.11889
## 2to3 - RB == 0 -0.09008 0.02596 -3.469 0.00277 **
## NYB - RB == 0 -0.07626 0.02513 -3.034 0.01267 *
## 2to3 - 1to2 == 0 -0.03256 0.02629 -1.238 0.60242
## NYB - 1to2 == 0 -0.01874 0.02552 -0.734 0.88332
## NYB - 2to3 == 0 0.01382 0.02552 0.541 0.94883
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(CP + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit20,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -0.35627 0.03249 -10.966 <1e-04 ***
## 2to3 - RB == 0 -0.50797 0.03215 -15.799 <1e-04 ***
## 3plus - RB == 0 -0.54404 0.03267 -16.653 <1e-04 ***
## 2to3 - 1to2 == 0 -0.15171 0.03259 -4.656 <1e-04 ***
## 3plus - 1to2 == 0 -0.18777 0.03311 -5.671 <1e-04 ***
## 3plus - 2to3 == 0 -0.03606 0.03249 -1.110 0.683
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = log(CP + 1) ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit20,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -0.352450 0.028250 -12.476 <1e-05 ***
## 2to3 - RB == 0 -0.374421 0.028689 -13.051 <1e-05 ***
## 3plus - RB == 0 -0.361363 0.028689 -12.596 <1e-05 ***
## 2to3 - 1to2 == 0 -0.021972 0.029185 -0.753 0.875
## 3plus - 1to2 == 0 -0.008913 0.029185 -0.305 0.990
## 3plus - 2to3 == 0 0.013059 0.029559 0.442 0.971
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)Moisture
Higher in burn patches later in the season
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular
## boundary (singular) fit: see ?isSingular## Data: HRECSIOmit
## Models:
## MoistureNull: Moisture ~ 1 + (1 | Location/Month/Year)
## MoistureG: Moisture ~ Treatment + (1 | Location/Month/Year)
## MoistureT: Moisture ~ TSF + (1 | Location/Month/Year)
## MoistureE: Moisture ~ ESD + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## MoistureNull 5 14244 14272 -7116.8 14234
## MoistureG 6 14244 14278 -7116.1 14232 1.3522 1 0.2449
## MoistureT 9 13868 13918 -6924.7 13850 382.8398 3 <2e-16 ***
## MoistureE 9 14239 14289 -7110.4 14221 0.0000 0
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Data: HRECSIOmit
## Models:
## MoistureNull: Moisture ~ 1 + (1 | Location/Month/Year)
## MoistureT: Moisture ~ TSF + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## MoistureNull 5 14244 14272 -7116.8 14234
## MoistureT 9 13868 13918 -6924.7 13850 384.19 4 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Data: HRECSIOmit
## Models:
## MoistureNull: Moisture ~ 1 + (1 | Location/Month/Year)
## MoistureG: Moisture ~ Treatment + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## MoistureNull 5 14244 14272 -7116.8 14234
## MoistureG 6 14244 14278 -7116.1 14232 1.3522 1 0.2449## Data: HRECSIOmit
## Models:
## MoistureNull: Moisture ~ 1 + (1 | Location/Month/Year)
## MoistureE: Moisture ~ ESD + (1 | Location/Month/Year)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## MoistureNull 5 14244 14272 -7116.8 14234
## MoistureE 9 14239 14289 -7110.4 14221 12.894 4 0.01181 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = Moisture ~ TSF + (1 | Location/Month/Year), data = HRECSIOmit,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -6.5321 0.4891 -13.355 < 0.001 ***
## 2to3 - RB == 0 -8.6338 0.5712 -15.114 < 0.001 ***
## 3plus - RB == 0 -10.2354 0.7524 -13.604 < 0.001 ***
## NYB - RB == 0 -6.7200 0.4559 -14.738 < 0.001 ***
## 2to3 - 1to2 == 0 -2.1017 0.5832 -3.604 0.00265 **
## 3plus - 1to2 == 0 -3.7033 0.7593 -4.877 < 0.001 ***
## NYB - 1to2 == 0 -0.1879 0.5170 -0.363 0.99608
## 3plus - 2to3 == 0 -1.6016 0.7787 -2.057 0.23016
## NYB - 2to3 == 0 1.9139 0.6222 3.076 0.01667 *
## NYB - 3plus == 0 3.5154 0.8148 4.314 < 0.001 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = Moisture ~ Treatment + (1 | Location/Month/Year),
## data = HRECSIOmit, REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Sheep - Cattle == 0 -3.977 3.380 -1.177 0.239
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = Moisture ~ ESD + (1 | Location/Month/Year), data = HRECSIOmit,
## REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## Loamy - Clayey == 0 1.9473 0.7643 2.548 0.0690 .
## Saline Lowland - Clayey == 0 0.9303 0.5627 1.653 0.4297
## Sandy - Clayey == 0 1.4924 0.5294 2.819 0.0328 *
## Thin Claypan - Clayey == 0 -1.4697 1.5461 -0.951 0.8615
## Saline Lowland - Loamy == 0 -1.0170 0.7629 -1.333 0.6407
## Sandy - Loamy == 0 -0.4549 0.6710 -0.678 0.9556
## Thin Claypan - Loamy == 0 -3.4170 1.6171 -2.113 0.1903
## Sandy - Saline Lowland == 0 0.5621 0.4608 1.220 0.7135
## Thin Claypan - Saline Lowland == 0 -2.4001 1.4768 -1.625 0.4474
## Thin Claypan - Sandy == 0 -2.9621 1.4874 -1.992 0.2424
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = Moisture ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit17,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## NYB - RB == 0 -10.175 1.325 -7.678 1.62e-14 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = Moisture ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit17,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## NYB - RB == 0 -5.370 1.267 -4.237 2.27e-05 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = Moisture ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit18,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -0.8260 1.1098 -0.744 0.73598
## NYB - RB == 0 -3.1900 0.9582 -3.329 0.00267 **
## NYB - 1to2 == 0 -2.3640 0.9788 -2.415 0.04125 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = Moisture ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit18,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -4.0994 1.7897 -2.291 0.05650 .
## NYB - RB == 0 -4.7905 1.5316 -3.128 0.00503 **
## NYB - 1to2 == 0 -0.6911 1.5316 -0.451 0.89318
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = Moisture ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit19,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -5.3644 0.9056 -5.924 <0.001 ***
## 2to3 - RB == 0 -5.1979 0.9201 -5.649 <0.001 ***
## NYB - RB == 0 -3.4714 0.9028 -3.845 <0.001 ***
## 2to3 - 1to2 == 0 0.1665 0.9029 0.184 0.998
## NYB - 1to2 == 0 1.8930 0.8936 2.118 0.147
## NYB - 2to3 == 0 1.7266 0.9080 1.902 0.227
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = Moisture ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit19,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -2.7675 0.9072 -3.051 0.0122 *
## 2to3 - RB == 0 -4.1885 0.9072 -4.617 <0.001 ***
## NYB - RB == 0 -5.1625 0.8781 -5.879 <0.001 ***
## 2to3 - 1to2 == 0 -1.4210 0.9187 -1.547 0.4094
## NYB - 1to2 == 0 -2.3950 0.8916 -2.686 0.0359 *
## NYB - 2to3 == 0 -0.9740 0.8916 -1.092 0.6941
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = Moisture ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit20,
## Treatment == "Cattle"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -10.7298 0.9923 -10.813 <0.001 ***
## 2to3 - RB == 0 -13.7527 0.9821 -14.003 <0.001 ***
## 3plus - RB == 0 -13.0935 0.9979 -13.121 <0.001 ***
## 2to3 - 1to2 == 0 -3.0230 0.9953 -3.037 0.0125 *
## 3plus - 1to2 == 0 -2.3637 1.0113 -2.337 0.0898 .
## 3plus - 2to3 == 0 0.6592 0.9924 0.664 0.9105
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)## boundary (singular) fit: see ?isSingular##
## Simultaneous Tests for General Linear Hypotheses
##
## Multiple Comparisons of Means: Tukey Contrasts
##
##
## Fit: lmer(formula = Moisture ~ TSF + (1 | Location/Month), data = subset(HRECSIOmit20,
## Treatment == "Sheep"), REML = FALSE)
##
## Linear Hypotheses:
## Estimate Std. Error z value Pr(>|z|)
## 1to2 - RB == 0 -12.4905 1.1335 -11.019 <1e-05 ***
## 2to3 - RB == 0 -13.4682 1.1511 -11.700 <1e-05 ***
## 3plus - RB == 0 -13.8731 1.1511 -12.052 <1e-05 ***
## 2to3 - 1to2 == 0 -0.9776 1.1710 -0.835 0.838
## 3plus - 1to2 == 0 -1.3826 1.1710 -1.181 0.639
## 3plus - 2to3 == 0 -0.4050 1.1860 -0.341 0.986
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## (Adjusted p values reported -- single-step method)Selection Index
The selection index data is very skewed, going to come back for the stats later
