1 Set up

1.1 1) Load packages

1.2 2) Import data

1.3 3) Data Wrangling

1.4 4) Theme James

2 Infiltration

2.1 Infiltration function

2.2 Infiltration

2.3 Infiltration statistics

## [1] "Single Ring"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## treatment        951.71  475.85     2 23.179  46.197 8.191e-09 ***
## precip           464.18  232.09     2 23.178  22.532 3.673e-06 ***
## treatment:precip 476.75  119.19     4 23.172  11.571 2.612e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## [1] "Cornell Sprinkle Infiltrometer"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## treatment         36.15  18.077     2    25  4.4263   0.02261 *  
## precip            17.86   8.928     2    25  2.1862   0.13336    
## treatment:precip 319.70  79.924     4    25 19.5697 2.086e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

3 Extracting model means and pairwise comparisons

## [1] "Single Ring"
##   treatment precip    lsmean       SE       df   lower.CL  upper.CL .group
## 1        NP    472 28.062914 1.986313 24.71643 23.9696447 32.156183     a 
## 2        EA    472  6.401788 1.720460 23.01103  2.8428386  9.960737      b
## 3        AG    472  4.356709 1.720460 23.01103  0.7977601  7.915658      b
## 5        NP    579  8.509000 1.720460 23.01103  4.9500509 12.067949     a 
## 4        AG    579  2.667000 1.720460 23.01103 -0.8919491  6.225949      b
## 6        EA    579  2.540000 1.720460 23.01103 -1.0189491  6.098949      b
## 7        NP    850  8.892644 1.720460 23.01103  5.3336945 12.451593     a 
## 8        EA    850  6.032500 1.720460 23.01103  2.4735509  9.591449     ab
## 9        AG    850  1.843352 1.720460 23.01103 -1.7155968  5.402301      b
## [1] "Cornell Sprinkle Infiltrometer"
## precip = 472:
##  treatment lsmean   SE df lower.CL upper.CL .group
##  NP          22.2 1.01 25     20.1     24.2  a    
##  AG          18.5 1.01 25     16.4     20.5   b   
##  EA          18.0 1.19 25     15.6     20.4   b   
## 
## precip = 579:
##  treatment lsmean   SE df lower.CL upper.CL .group
##  AG          21.7 1.01 25     19.7     23.8  a    
##  NP          20.5 1.19 25     18.0     22.9  a    
##  EA          14.0 1.01 25     11.9     16.1   b   
## 
## precip = 850:
##  treatment lsmean   SE df lower.CL upper.CL .group
##  EA          22.5 1.01 25     20.5     24.6  a    
##  NP          17.8 1.01 25     15.7     19.8   b   
##  AG          13.0 1.01 25     10.9     15.1    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05 
## NOTE: Compact letter displays can be misleading
##       because they show NON-findings rather than findings.
##       Consider using 'pairs()', 'pwpp()', or 'pwpm()' instead.

4 Single Ring Plotting

5 Cornell Sprinkle Infiltrometer Plotting

5.1 Combined graph

5.2 Final graph

6 Right post hoc test

6.1 Single Ring Infiltrometer

## [1] "Single Ring"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq Mean Sq NumDF  DenDF F value    Pr(>F)    
## treatment        951.71  475.85     2 23.179  46.197 8.191e-09 ***
## precip           464.18  232.09     2 23.178  22.532 3.673e-06 ***
## treatment:precip 476.75  119.19     4 23.172  11.571 2.612e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  treatment precip lsmean   SE   df lower.CL upper.CL .group
##  NP        472     28.06 1.99 24.7   23.970    32.16  a    
##  NP        850      8.89 1.72 23.0    5.334    12.45   b   
##  NP        579      8.51 1.72 23.0    4.950    12.07   b   
##  EA        472      6.40 1.72 23.0    2.843     9.96   bc  
##  EA        850      6.03 1.72 23.0    2.474     9.59   bc  
##  AG        472      4.36 1.72 23.0    0.798     7.92   bc  
##  AG        579      2.67 1.72 23.0   -0.892     6.23    c  
##  EA        579      2.54 1.72 23.0   -1.019     6.10    c  
##  AG        850      1.84 1.72 23.0   -1.716     5.40    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05 
## NOTE: Compact letter displays can be misleading
##       because they show NON-findings rather than findings.
##       Consider using 'pairs()', 'pwpp()', or 'pwpm()' instead.

6.2 Cornell Sprinkle Infiltrometer

## [1] "Cornell Sprinkle Infiltrometer"
## Type III Analysis of Variance Table with Satterthwaite's method
##                  Sum Sq Mean Sq NumDF DenDF F value    Pr(>F)    
## treatment         36.15  18.077     2    25  4.4263   0.02261 *  
## precip            17.86   8.928     2    25  2.1862   0.13336    
## treatment:precip 319.70  79.924     4    25 19.5697 2.086e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  treatment precip lsmean   SE df lower.CL upper.CL .group
##  EA        850      22.5 1.01 25     20.5     24.6  a    
##  NP        472      22.2 1.01 25     20.1     24.2  a    
##  AG        579      21.7 1.01 25     19.7     23.8  a    
##  NP        579      20.5 1.19 25     18.0     22.9  ab   
##  AG        472      18.5 1.01 25     16.4     20.5   b   
##  EA        472      18.0 1.19 25     15.6     20.4   b   
##  NP        850      17.8 1.01 25     15.7     19.8   b   
##  EA        579      14.0 1.01 25     11.9     16.1    c  
##  AG        850      13.0 1.01 25     10.9     15.1    c  
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05 
## NOTE: Compact letter displays can be misleading
##       because they show NON-findings rather than findings.
##       Consider using 'pairs()', 'pwpp()', or 'pwpm()' instead.