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
