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"
## precip = 472:
##  treatment lsmean   SE   df lower.CL upper.CL .group
##  NP         28.06 1.99 24.7   23.970    32.16  a    
##  EA          6.40 1.72 23.0    2.843     9.96   b   
##  AG          4.36 1.72 23.0    0.798     7.92   b   
## 
## precip = 579:
##  treatment lsmean   SE   df lower.CL upper.CL .group
##  NP          8.89 1.72 23.0    5.334    12.45  a    
##  EA          6.03 1.72 23.0    2.474     9.59  ab   
##  AG          1.84 1.72 23.0   -1.716     5.40   b   
## 
## precip = 850:
##  treatment lsmean   SE   df lower.CL upper.CL .group
##  NP          8.51 1.72 23.0    4.950    12.07  a    
##  AG          2.67 1.72 23.0   -0.892     6.23   b   
##  EA          2.54 1.72 23.0   -1.019     6.10   b   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   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
## 4        NP    579  8.892644 1.720460 23.01103  5.3336945 12.451593     a 
## 5        EA    579  6.032500 1.720460 23.01103  2.4735509  9.591449     ab
## 6        AG    579  1.843352 1.720460 23.01103 -1.7155968  5.402301      b
## 8        NP    850  8.509000 1.720460 23.01103  4.9500509 12.067949     a 
## 7        AG    850  2.667000 1.720460 23.01103 -0.8919491  6.225949      b
## 9        EA    850  2.540000 1.720460 23.01103 -1.0189491  6.098949      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
##  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  
## 
## precip = 850:
##  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   
## 
## Degrees-of-freedom method: kenward-roger 
## Confidence level used: 0.95 
## significance level used: alpha = 0.05
##   treatment precip   lsmean       SE df lower.CL upper.CL .group
## 2        NP    472 22.16082 1.010455 25 20.07975 24.24189     a 
## 1        AG    472 18.46192 1.010455 25 16.38085 20.54299      b
## 3        EA    472 17.99958 1.188525 25 15.55177 20.44740      b
## 4        EA    579 22.53130 1.010455 25 20.45023 24.61238    a  
## 6        NP    579 17.76469 1.010455 25 15.68362 19.84576     b 
## 5        AG    579 13.02804 1.010455 25 10.94697 15.10911      c
## 8        AG    850 21.73883 1.010455 25 19.65776 23.81990     a 
## 9        NP    850 20.48556 1.188525 25 18.03774 22.93337     a 
## 7        EA    850 14.00480 1.010455 25 11.92373 16.08587      b

4 Single Ring Plotting

5 Cornell Sprinkle Infiltrometer Plotting