The data

The results

Dry bulk density 15 cm

## Warning: Removed 12 rows containing non-finite values (`stat_boxplot()`).

## Warning: Removed 12 rows containing non-finite values (`stat_boxplot()`).

## Warning: Removed 12 rows containing non-finite values (`stat_boxplot()`).
## Warning: Removed 12 rows containing non-finite values (`stat_bin()`).

## Analysis of Variance Table
## 
## Response: bulkdensity2022[, var]
##              Df  Sum Sq   Mean Sq F value Pr(>F)
## treatment     2 0.00668 0.0033385  0.4594 0.6333
## blok          3 0.02414 0.0080459  1.1072 0.3510
## blok:subplot  8 0.04233 0.0052908  0.7281 0.6663
## Residuals    82 0.59589 0.0072670               
## 
##  Shapiro-Wilk normality test
## 
## data:  model$residuals
## W = 0.9881, p-value = 0.5457
## 
## $emmeans
##  contrast                 estimate     SE df t.ratio p.value
##  Fibrous roots - Mixture  0.017522 0.0213 82   0.822  0.6904
##  Fibrous roots - Taproots 0.017856 0.0213 82   0.838  0.6807
##  Mixture - Taproots       0.000334 0.0213 82   0.016  0.9999
## 
## Results are averaged over the levels of: subplot, blok 
## Note: contrasts are still on the [ scale 
## P value adjustment: tukey method for comparing a family of 3 estimates 
## 
## $contrasts
##  contrast                                                estimate     SE df
##  (Fibrous roots - Mixture) - (Fibrous roots - Taproots) -0.000334 0.0213 82
##  (Fibrous roots - Mixture) - (Mixture - Taproots)        0.017188 0.0369 82
##  (Fibrous roots - Taproots) - (Mixture - Taproots)       0.017522 0.0213 82
##  t.ratio p.value
##   -0.016  0.9999
##    0.466  0.8875
##    0.822  0.6904
## 
## Results are averaged over the levels of: subplot, blok 
## P value adjustment: tukey method for comparing a family of 3 estimates

Bulk density 30 cm

## Warning: Removed 12 rows containing non-finite values (`stat_boxplot()`).

## Warning: Removed 12 rows containing non-finite values (`stat_boxplot()`).

## Warning: Removed 12 rows containing non-finite values (`stat_boxplot()`).
## Warning: Removed 12 rows containing non-finite values (`stat_bin()`).

## Analysis of Variance Table
## 
## Response: bulkdensity2022[, var]
##              Df  Sum Sq   Mean Sq F value Pr(>F)
## treatment     2 0.00074 0.0003706  0.0426 0.9583
## blok          3 0.02209 0.0073640  0.8459 0.4727
## blok:subplot  8 0.05078 0.0063469  0.7291 0.6654
## Residuals    82 0.71384 0.0087054               
## 
##  Shapiro-Wilk normality test
## 
## data:  model$residuals
## W = 0.98298, p-value = 0.2491
## 
## $emmeans
##  contrast                 estimate     SE df t.ratio p.value
##  Fibrous roots - Mixture  -0.00211 0.0233 82  -0.091  0.9955
##  Fibrous roots - Taproots -0.00666 0.0233 82  -0.285  0.9561
##  Mixture - Taproots       -0.00455 0.0233 82  -0.195  0.9793
## 
## Results are averaged over the levels of: subplot, blok 
## Note: contrasts are still on the [ scale 
## P value adjustment: tukey method for comparing a family of 3 estimates 
## 
## $contrasts
##  contrast                                               estimate     SE df
##  (Fibrous roots - Mixture) - (Fibrous roots - Taproots)  0.00455 0.0233 82
##  (Fibrous roots - Mixture) - (Mixture - Taproots)        0.00243 0.0404 82
##  (Fibrous roots - Taproots) - (Mixture - Taproots)      -0.00211 0.0233 82
##  t.ratio p.value
##    0.195  0.9793
##    0.060  0.9980
##   -0.091  0.9955
## 
## Results are averaged over the levels of: subplot, blok 
## P value adjustment: tukey method for comparing a family of 3 estimates