Larval density drives thermogenesis and affects microbiota and substrate properties in black soldier fly trials
Publication
Klammsteiner, T., Heussler, C.D., Insam, H., Schlick-Steiner, B.C.*, Steiner, F.M.* (2025). Larval density drives thermogenesis and affects microbiota and substrate properties in black soldier fly trials. iScience 28(7), 112794. doi.org/10.1016/j.isci.2025.112794
(*equal contribution)
Experimental design
Environmental Data Download
Biomass Input and Output
–> Download data here
Table 1. Raw data on larvae and substrate biomass input and output (in grams).
FM = fresh matter, DM = dry matter.
Performance Parameters
–> Download data here
Table 2. Measured and calculated performance parameters.
The peak temperature measurement of 2.5SORBN-2 was identified as an outlier and subsequently removed from the dataset. GR = growth rate, SGR = specific growth rate, FCR = feed conversion rate, WRI = waste reduction index, SR = substrate reduction, ECI = efficiency of conversion of ingested feed, ECD = efficiency of conversion of digested feed.
Temperature Profiles
–> Download data here
Figure 2. Animation of average (n = 4) substrate temperatures at four larval densities (0, 1.25, 2.5, 5 larvae/cm²) reared in food waste containing either no (sorb-) or 0.15% (w/w) potassium sorbate (sorb+).
Microbiome Data Download
All data was generated on the Illumina NovaSeq6000 platform, targetting the V4 region on the 16S rRNA gene using the primer pair 515F (5’-GTGCCAGCMGCCGCGGTAA-3’) and 806R (5’-GGACTACHVGGGTWTCTAAT-3’).
Larval Bacteriobiota
–> Download the filtered and rarefied data here)
–> Download the filtered and rarefied ASV table here
–> Download the filtered and rarefied taxonomy table here
–> Download the metadata table for larval gut samples here
The raw sequence data generated for this study can be found in the European Nucleotide Archive (ENA) under the bioproject accession number PRJEB67330.
Substrate Bacteriobiota
–> Download the filtered and rarefied data here
–> Download the filtered and rarefied ASV table here
–> Download the filtered and rarefied taxonomy table here
–> Download the metadata table for substrate samples here
The raw sequence data generated for this study can be found in the European Nucleotide Archive (ENA) under the bioproject accession number PRJEB67330.
Statistics
Environmental data
Show the code
# Load packages
library(tidyverse)
# Import data
data <- read.csv("../../data/20240712-summary-performance_parameters.csv", sep = ";")Boxplots
Show the code
# Plot grid
data %>%
mutate(GROUP = factor(GROUP,
levels = c('0SORBN', '1.25SORBN', '2.5SORBN', '5SORBN',
'0SORBP', '1.25SORBP', '2.5SORBP', '5SORBP'))) %>%
select(GROUP, 6:14) %>%
pivot_longer(cols = -GROUP, names_to = "PARAMETER", values_to = "VALUE") %>%
ggplot() +
geom_boxplot(aes(x = GROUP, y = VALUE), fill = "#d9d7bf") +
facet_wrap(.~PARAMETER, scales = "free", ncol = 5) +
labs(x = NULL, y = NULL) +
theme_classic() +
theme(axis.text.x = element_text(angle = 45, hjust = 0.95, vjust = 1),
strip.background = element_rect(colour = NA),
strip.text = element_text(face = "bold")) 
Density plots
Show the code
# Plot grid
data %>%
mutate(GROUP = factor(GROUP,
levels = c('0SORBN', '1.25SORBN', '2.5SORBN', '5SORBN',
'0SORBP', '1.25SORBP', '2.5SORBP', '5SORBP'))) %>%
select(GROUP, 6:14) %>%
pivot_longer(cols = -GROUP, names_to = "PARAMETER", values_to = "VALUE") %>%
ggplot() +
geom_density(aes(x = VALUE), fill = "#d9d7bf", colour = NA) +
facet_wrap(.~PARAMETER, scales = "free", ncol = 5) +
labs(x = NULL, y = NULL) +
theme_classic() +
theme(strip.background = element_rect(colour = NA),
strip.text = element_text(face = "bold")) 
Testing for Normality
Show the code
# Load packages
library(reactable)
# Import data
data <- read.csv("../../data/20240712-summary-performance_parameters.csv", sep = ";")
# Check for normality
normality_biomass <- shapiro.test(data$PEAK_BIOMASS)
normality_days <- shapiro.test(data$DAYS_PEAK)
normality_temp <- shapiro.test(data$PEAK_TEMP)
normality_gr <- shapiro.test(data$GR)
normality_sgr <- shapiro.test(data$SGR)
normality_fcr <- shapiro.test(data$FCR)
normality_bce <- shapiro.test(data$BCE)
normality_wri <- shapiro.test(data$WRI)
normality_sr <- shapiro.test(data$SR)
normality <- data.frame(
Parameter = c('Peak biomass', 'Days to peak', 'Peak temperature',
'GR', 'SGR', 'FCR', 'BCE', 'WRI', 'SR'),
pvalue = c(normality_biomass$p.value, normality_days$p.value,
normality_temp$p.value, normality_gr$p.value,
normality_sgr$p.value, normality_fcr$p.value,
normality_bce$p.value, normality_wri$p.value, normality_sr$p.value)
)The data for the parameters “Days to peak”, “Peak temperature”, “SGR”, “FCR”, “WRI”, and “SR” are not normally distributed.
Two-way ANOVA and Tukey’s HSD post-hoc test
For normally distributed data.
Peak biomass
Show the code
# ANOVA
summary(aov(PEAK_BIOMASS ~ SUBST * DENSITY, data = data)) Df Sum Sq Mean Sq F value Pr(>F)
SUBST 1 0.015050 0.015050 17.199 0.000499 ***
DENSITY 1 0.008571 0.008571 9.795 0.005277 **
SUBST:DENSITY 1 0.000469 0.000469 0.536 0.472566
Residuals 20 0.017501 0.000875
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
6 Beobachtungen als fehlend gelöschtShow the code
# Tukey's HSD post-hoc test
TukeyHSD(aov(PEAK_BIOMASS ~ GROUP, data = data)) Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = PEAK_BIOMASS ~ GROUP, data = data)
$GROUP
diff lwr upr p adj
1.25SORBP-1.25SORBN -0.04550 -0.106839756 0.01583976 0.2227613
2.5SORBN-1.25SORBN 0.05375 -0.007589756 0.11508976 0.1060864
2.5SORBP-1.25SORBN -0.01975 -0.081089756 0.04158976 0.9040283
5SORBN-1.25SORBN 0.04325 -0.018089756 0.10458976 0.2677750
5SORBP-1.25SORBN 0.01200 -0.049339756 0.07333976 0.9878957
2.5SORBN-1.25SORBP 0.09925 0.037910244 0.16058976 0.0008261
2.5SORBP-1.25SORBP 0.02575 -0.035589756 0.08708976 0.7631205
5SORBN-1.25SORBP 0.08875 0.027410244 0.15008976 0.0026001
5SORBP-1.25SORBP 0.05750 -0.003839756 0.11883976 0.0736458
2.5SORBP-2.5SORBN -0.07350 -0.134839756 -0.01216024 0.0138064
5SORBN-2.5SORBN -0.01050 -0.071839756 0.05083976 0.9934011
5SORBP-2.5SORBN -0.04175 -0.103089756 0.01958976 0.3011738
5SORBN-2.5SORBP 0.06300 0.001660244 0.12433976 0.0421427
5SORBP-2.5SORBP 0.03175 -0.029589756 0.09308976 0.5817312
5SORBP-5SORBN -0.03125 -0.092589756 0.03008976 0.5973153Growth rate
Show the code
# ANOVA
summary(aov(GR ~ SUBST * DENSITY, data = data)) Df Sum Sq Mean Sq F value Pr(>F)
SUBST 1 3.760 3.760 3.347 0.082291 .
DENSITY 1 21.757 21.757 19.363 0.000276 ***
SUBST:DENSITY 1 2.696 2.696 2.400 0.137033
Residuals 20 22.473 1.124
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
6 Beobachtungen als fehlend gelöschtShow the code
# Tukey's HSD post-hoc test
TukeyHSD(aov(GR ~ GROUP, data = data)) Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = GR ~ GROUP, data = data)
$GROUP
diff lwr upr p adj
1.25SORBP-1.25SORBN -0.150 -1.5420918 1.2420918 0.9992626
2.5SORBN-1.25SORBN 2.800 1.4079082 4.1920918 0.0000649
2.5SORBP-1.25SORBN 2.300 0.9079082 3.6920918 0.0006585
5SORBN-1.25SORBN 3.450 2.0579082 4.8420918 0.0000040
5SORBP-1.25SORBN 1.725 0.3329082 3.1170918 0.0105120
2.5SORBN-1.25SORBP 2.950 1.5579082 4.3420918 0.0000333
2.5SORBP-1.25SORBP 2.450 1.0579082 3.8420918 0.0003241
5SORBN-1.25SORBP 3.600 2.2079082 4.9920918 0.0000022
5SORBP-1.25SORBP 1.875 0.4829082 3.2670918 0.0050990
2.5SORBP-2.5SORBN -0.500 -1.8920918 0.8920918 0.8576766
5SORBN-2.5SORBN 0.650 -0.7420918 2.0420918 0.6780902
5SORBP-2.5SORBN -1.075 -2.4670918 0.3170918 0.1899579
5SORBN-2.5SORBP 1.150 -0.2420918 2.5420918 0.1414509
5SORBP-2.5SORBP -0.575 -1.9670918 0.8170918 0.7746145
5SORBP-5SORBN -1.725 -3.1170918 -0.3329082 0.0105120Scheirer-Ray-Hare test and Dunn’s post-hoc test
For not-normally distributed data.
Days to peak
Show the code
# Scheirer-Ray-Hare Test
scheirerRayHare(DAYS_PEAK ~ SUBST * DENSITY, data = data)
DV: DAYS_PEAK
Observations: 24
D: 0.8347826
MS total: 50 Df Sum Sq H p.value
SUBST 1 170.67 4.0889 0.043166
DENSITY 2 688.00 16.4833 0.000263
SUBST:DENSITY 2 101.33 2.4278 0.297040
Residuals 18 0.00 Show the code
# Dunn's post-hoc test
dunnTest(DAYS_PEAK ~ GROUP, data = data, method = 'bh') Comparison Z P.unadj P.adj
1 1.25SORBN - 1.25SORBP -2.188988 0.0285977421 0.047662903
2 1.25SORBN - 2.5SORBN 0.000000 1.0000000000 1.000000000
3 1.25SORBP - 2.5SORBN 2.188988 0.0285977421 0.053620766
4 1.25SORBN - 2.5SORBP 0.000000 1.0000000000 1.000000000
5 1.25SORBP - 2.5SORBP 2.188988 0.0285977421 0.061280876
6 2.5SORBN - 2.5SORBP 0.000000 1.0000000000 1.000000000
7 1.25SORBN - 5SORBN -2.188988 0.0285977421 0.071494355
8 1.25SORBP - 5SORBN 0.000000 1.0000000000 1.000000000
9 2.5SORBN - 5SORBN -2.188988 0.0285977421 0.085793226
10 2.5SORBP - 5SORBN -2.188988 0.0285977421 0.107241533
11 1.25SORBN - 5SORBP -3.502380 0.0004611212 0.002305606
12 1.25SORBP - 5SORBP -1.313393 0.1890507054 0.257796416
13 2.5SORBN - 5SORBP -3.502380 0.0004611212 0.003458409
14 2.5SORBP - 5SORBP -3.502380 0.0004611212 0.006916818
15 5SORBN - 5SORBP -1.313393 0.1890507054 0.283576058SGR
Show the code
# Scheirer-Ray-Hare Test
scheirerRayHare(SGR ~ SUBST * DENSITY, data = data)
DV: SGR
Observations: 24
D: 0.9578261
MS total: 50 Df Sum Sq H p.value
SUBST 1 121.50 2.5370 0.11121
DENSITY 2 769.00 16.0572 0.00033
SUBST:DENSITY 2 54.25 1.1328 0.56757
Residuals 18 156.75 Show the code
# Dunn's post-hoc test
dunnTest(SGR ~ GROUP, data = data, method = 'bh') Comparison Z P.unadj P.adj
1 1.25SORBN - 1.25SORBP 0.2554446 0.7983797491 0.79837975
2 1.25SORBN - 2.5SORBN -2.7588014 0.0058013784 0.02175517
3 1.25SORBP - 2.5SORBN -3.0142459 0.0025761887 0.01288094
4 1.25SORBN - 2.5SORBP -1.9924677 0.0463197702 0.11579943
5 1.25SORBP - 2.5SORBP -2.2479122 0.0245817871 0.07374536
6 2.5SORBN - 2.5SORBP 0.7663337 0.4434777623 0.55434720
7 1.25SORBN - 5SORBN -3.1419682 0.0016781624 0.01258622
8 1.25SORBP - 5SORBN -3.3974128 0.0006802626 0.01020394
9 2.5SORBN - 5SORBN -0.3831669 0.7015960501 0.75171005
10 2.5SORBP - 5SORBN -1.1495006 0.2503496308 0.34138586
11 1.25SORBN - 5SORBP -1.4049451 0.1600375910 0.26672932
12 1.25SORBP - 5SORBP -1.6603897 0.0968360777 0.18156765
13 2.5SORBN - 5SORBP 1.3538562 0.1757822502 0.26367338
14 2.5SORBP - 5SORBP 0.5875225 0.5568528338 0.64252250
15 5SORBN - 5SORBP 1.7370231 0.0823830958 0.17653521FCR
Show the code
# Scheirer-Ray-Hare Test
scheirerRayHare(FCR ~ SUBST * DENSITY, data = data)
DV: FCR
Observations: 24
D: 0.9904348
MS total: 50 Df Sum Sq H p.value
SUBST 1 30.37 0.6134 0.43352
DENSITY 2 768.25 15.5134 0.00043
SUBST:DENSITY 2 279.75 5.6490 0.05934
Residuals 18 60.62 Show the code
# Dunn's post-hoc test
dunnTest(FCR ~ GROUP, data = data, method = 'bh') Comparison Z P.unadj P.adj
1 1.25SORBN - 1.25SORBP 0.8038538 4.214814e-01 0.4515872573
2 1.25SORBN - 2.5SORBN 4.0192688 5.837904e-05 0.0008756856
3 1.25SORBP - 2.5SORBN 3.2154150 1.302561e-03 0.0097692040
4 1.25SORBN - 2.5SORBP 1.6579484 9.732789e-02 0.1824897894
5 1.25SORBP - 2.5SORBP 0.8540946 3.930526e-01 0.4535221929
6 2.5SORBN - 2.5SORBP -2.3613204 1.820999e-02 0.0546299692
7 1.25SORBN - 5SORBN 2.6878860 7.190594e-03 0.0269647282
8 1.25SORBP - 5SORBN 1.8840322 5.956061e-02 0.1276298886
9 2.5SORBN - 5SORBN -1.3313828 1.830631e-01 0.3051051528
10 2.5SORBP - 5SORBN 1.0299376 3.030393e-01 0.3787991092
11 1.25SORBN - 5SORBP 2.8888494 3.866542e-03 0.0193327081
12 1.25SORBP - 5SORBP 2.0849957 3.706969e-02 0.0926742320
13 2.5SORBN - 5SORBP -1.1304193 2.582996e-01 0.3522266822
14 2.5SORBP - 5SORBP 1.2309011 2.183599e-01 0.3275398048
15 5SORBN - 5SORBP 0.2009634 8.407272e-01 0.8407271623WRI
Show the code
# Scheirer-Ray-Hare Test
scheirerRayHare(WRI ~ SUBST * DENSITY, data = data)
DV: WRI
Observations: 24
D: 0.993913
MS total: 50 Df Sum Sq H p.value
SUBST 1 459.37 9.2438 0.002363
DENSITY 2 307.56 6.1889 0.045299
SUBST:DENSITY 2 189.44 3.8120 0.148677
Residuals 18 186.62 Show the code
# Dunn's post-hoc test
dunnTest(WRI ~ GROUP, data = data, method = 'bh') Comparison Z P.unadj P.adj
1 1.25SORBN - 1.25SORBP 1.9810385 4.758696e-02 0.1189674044
2 1.25SORBN - 2.5SORBN -2.0311913 4.223559e-02 0.1267067709
3 1.25SORBP - 2.5SORBN -4.0122298 6.014791e-05 0.0009022187
4 1.25SORBN - 2.5SORBP 0.9779810 3.280837e-01 0.4101046472
5 1.25SORBP - 2.5SORBP -1.0030574 3.158331e-01 0.4306815631
6 2.5SORBN - 2.5SORBP 3.0091723 2.619605e-03 0.0196470361
7 1.25SORBN - 5SORBN -0.6770638 4.983655e-01 0.5750371062
8 1.25SORBP - 5SORBN -2.6581022 7.858205e-03 0.0392910267
9 2.5SORBN - 5SORBN 1.3541275 1.756957e-01 0.2635435329
10 2.5SORBP - 5SORBN -1.6550448 9.791541e-02 0.2098187453
11 1.25SORBN - 5SORBP -0.4012230 6.882560e-01 0.7374171056
12 1.25SORBP - 5SORBP -2.3822614 1.720668e-02 0.0645250396
13 2.5SORBN - 5SORBP 1.6299683 1.031082e-01 0.1933278513
14 2.5SORBP - 5SORBP -1.3792040 1.678319e-01 0.2797197832
15 5SORBN - 5SORBP 0.2758408 7.826703e-01 0.7826703474SR
Show the code
# Scheirer-Ray-Hare Test
scheirerRayHare(SR ~ SUBST * DENSITY, data = data)
DV: SR
Observations: 24
D: 0.993913
MS total: 50 Df Sum Sq H p.value
SUBST 1 459.37 9.2438 0.002363
DENSITY 2 307.56 6.1889 0.045299
SUBST:DENSITY 2 189.44 3.8120 0.148677
Residuals 18 186.62 Show the code
# Dunn's post-hoc test
dunnTest(SR ~ GROUP, data = data, method = 'bh') Comparison Z P.unadj P.adj
1 1.25SORBN - 1.25SORBP 1.9810385 4.758696e-02 0.1189674044
2 1.25SORBN - 2.5SORBN -2.0311913 4.223559e-02 0.1267067709
3 1.25SORBP - 2.5SORBN -4.0122298 6.014791e-05 0.0009022187
4 1.25SORBN - 2.5SORBP 0.9779810 3.280837e-01 0.4101046472
5 1.25SORBP - 2.5SORBP -1.0030574 3.158331e-01 0.4306815631
6 2.5SORBN - 2.5SORBP 3.0091723 2.619605e-03 0.0196470361
7 1.25SORBN - 5SORBN -0.6770638 4.983655e-01 0.5750371062
8 1.25SORBP - 5SORBN -2.6581022 7.858205e-03 0.0392910267
9 2.5SORBN - 5SORBN 1.3541275 1.756957e-01 0.2635435329
10 2.5SORBP - 5SORBN -1.6550448 9.791541e-02 0.2098187453
11 1.25SORBN - 5SORBP -0.4012230 6.882560e-01 0.7374171056
12 1.25SORBP - 5SORBP -2.3822614 1.720668e-02 0.0645250396
13 2.5SORBN - 5SORBP 1.6299683 1.031082e-01 0.1933278513
14 2.5SORBP - 5SORBP -1.3792040 1.678319e-01 0.2797197832
15 5SORBN - 5SORBP 0.2758408 7.826703e-01 0.7826703474Peak temperature
Show the code
# Scheirer-Ray-Hare Test
scheirerRayHare(PEAK_TEMP ~ SUBST * DENSITY, data = data)
DV: PEAK_TEMP
Observations: 29
D: 0.9972906
MS total: 72.5 Df Sum Sq H p.value
SUBST 1 294.67 4.0755 0.04351
DENSITY 3 1546.77 21.3927 0.00009
SUBST:DENSITY 3 43.57 0.6025 0.89585
Residuals 21 148.27 Show the code
# Dunn's post-hoc test
dunnTest(PEAK_TEMP ~ GROUP, data = data %>% filter(DENSITY != 0 & ID != '2.5SORBN-2'), method = 'bh') Comparison Z P.unadj P.adj
1 1.25SORBN - 1.25SORBP 1.3848333 0.1661034583 0.276839097
2 1.25SORBN - 2.5SORBN -1.0240725 0.3058010507 0.417001433
3 1.25SORBP - 2.5SORBN -2.3061789 0.0211006365 0.063301910
4 1.25SORBN - 2.5SORBP 0.4180629 0.6759011505 0.675901151
5 1.25SORBP - 2.5SORBP -0.9667704 0.3336588302 0.417073538
6 2.5SORBN - 2.5SORBP 1.4111235 0.1582082108 0.296640395
7 1.25SORBN - 5SORBN -2.3516036 0.0186926810 0.070097554
8 1.25SORBP - 5SORBN -3.7364369 0.0001866463 0.002799694
9 2.5SORBN - 5SORBN -1.1530895 0.2488736657 0.373310499
10 2.5SORBP - 5SORBN -2.7696665 0.0056113712 0.028056856
11 1.25SORBN - 5SORBP -1.7767672 0.0756065503 0.162014036
12 1.25SORBP - 5SORBP -3.1616005 0.0015690468 0.011767851
13 2.5SORBN - 5SORBP -0.6208943 0.5346691542 0.616925947
14 2.5SORBP - 5SORBP -2.1948301 0.0281757894 0.070439474
15 5SORBN - 5SORBP 0.5748364 0.5654019144 0.605787765BCE
Show the code
# Scheirer-Ray-Hare Test
scheirerRayHare(BCE ~ SUBST * DENSITY, data = data)
DV: BCE
Observations: 24
D: 1
MS total: 50 Df Sum Sq H p.value
SUBST 1 160.17 3.2033 0.073488
DENSITY 2 768.25 15.3650 0.000461
SUBST:DENSITY 2 148.08 2.9617 0.227448
Residuals 18 73.50 Show the code
# Dunn's post-hoc test
dunnTest(BCE ~ GROUP, data = data, method = 'bh') Comparison Z P.unadj P.adj
1 1.25SORBN - 1.25SORBP 0.80 4.237108e-01 0.4539758541
2 1.25SORBN - 2.5SORBN -3.20 1.374276e-03 0.0103070691
3 1.25SORBP - 2.5SORBN -4.00 6.334248e-05 0.0009501373
4 1.25SORBN - 2.5SORBP -0.85 3.953251e-01 0.4561443303
5 1.25SORBP - 2.5SORBP -1.65 9.894294e-02 0.1855180051
6 2.5SORBN - 2.5SORBP 2.35 1.877341e-02 0.0563202332
7 1.25SORBN - 5SORBN -1.95 5.117612e-02 0.1096631122
8 1.25SORBP - 5SORBN -2.75 5.959526e-03 0.0223482243
9 2.5SORBN - 5SORBN 1.25 2.112995e-01 0.3521659122
10 2.5SORBP - 5SORBN -1.10 2.713321e-01 0.3391651524
11 1.25SORBN - 5SORBP -2.00 4.550026e-02 0.1137506597
12 1.25SORBP - 5SORBP -2.80 5.110261e-03 0.0255513033
13 2.5SORBN - 5SORBP 1.20 2.301393e-01 0.3452090107
14 2.5SORBP - 5SORBP -1.15 2.501439e-01 0.3411052790
15 5SORBN - 5SORBP -0.05 9.601224e-01 0.9601223883Physicochemical parameters
Show the code
# Load packages
library(tidyverse)
# Import data
physchem <- read.csv("../../data/20240712-summary-physchem.csv", sep = ";")Boxplots
Show the code
# Plot grid
physchem %>%
mutate(GROUP = factor(GROUP,
levels = c('FSORBN', 'FSORBP', '0SORBN', '1.25SORBN', '2.5SORBN',
'5SORBN', '0SORBP', '1.25SORBP', '2.5SORBP', '5SORBP'))) %>%
select(GROUP, 6:11) %>%
pivot_longer(cols = -GROUP, names_to = "PARAMETER", values_to = "VALUE") %>%
ggplot() +
geom_boxplot(aes(x = GROUP, y = VALUE), fill = "#d9d7bf") +
facet_wrap(.~PARAMETER, scales = "free", ncol = 5) +
labs(x = NULL, y = NULL) +
theme_classic() +
theme(axis.text.x = element_text(angle = 45, hjust = 0.95, vjust = 1),
strip.background = element_rect(colour = NA),
strip.text = element_text(face = "bold")) 
Density plots
Show the code
# Plot grid
physchem %>%
mutate(GROUP = factor(GROUP,
levels = c('FSORBN', 'FSORBP', '0SORBN', '1.25SORBN', '2.5SORBN',
'5SORBN', '0SORBP', '1.25SORBP', '2.5SORBP', '5SORBP'))) %>%
select(GROUP, 6:11) %>%
pivot_longer(cols = -GROUP, names_to = "PARAMETER", values_to = "VALUE") %>%
ggplot() +
geom_density(aes(x = VALUE), fill = "#d9d7bf", colour = NA) +
facet_wrap(.~PARAMETER, scales = "free", ncol = 5) +
labs(x = NULL, y = NULL) +
theme_classic() +
theme(strip.background = element_rect(colour = NA),
strip.text = element_text(face = "bold")) 
Testing for Normality
Show the code
# Load packages
library(reactable)
# Import data
physchem <- read.csv("../../data/20240712-summary-physchem.csv", sep = ";")
# Check for normality
normality_ph <- shapiro.test(physchem$PH)
normality_h2o_reduct <- shapiro.test(physchem$H2O_REDUCT_REL)
normality_vs_reduct <- shapiro.test(physchem$VS_REDUCT_REL)
normality_C <- shapiro.test(physchem$C_REL)
normality_N <- shapiro.test(physchem$N_REL)
normality_CN_ratio <- shapiro.test(physchem$CN_ratio)
normality <- data.frame(
Parameter = c('pH', 'H2O loss [%]', 'VS reduction [%]',
'C [%]', 'N [%]', 'C:N ratio'),
pvalue = c(normality_ph$p.value, normality_h2o_reduct$p.value,
normality_vs_reduct$p.value, normality_C$p.value,
normality_N$p.value, normality_CN_ratio$p.value)
)Relative Nitrogen content is normally distributed, all other tested paramteres are not normally distributed.
Two-way ANOVA and Tukey’s HSD post-hoc test
For normally distributed data.
Nitrogen content
Show the code
# ANOVA
summary(aov(N_REL ~ SUBST * DENSITY, data = physchem)) Df Sum Sq Mean Sq F value Pr(>F)
SUBST 1 0.033 0.0333 0.178 0.6766
DENSITY 4 7.625 1.9062 10.179 3.23e-05 ***
SUBST:DENSITY 4 2.859 0.7147 3.817 0.0134 *
Residuals 28 5.243 0.1873
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1Show the code
# Tukey's HSD post-hoc test
TukeyHSD(aov(N_REL ~ GROUP, data = physchem)) Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = N_REL ~ GROUP, data = physchem)
$GROUP
diff lwr upr p adj
0SORBP-0SORBN -4.781333e-01 -1.6898503 0.73358367 0.9313275
1.25SORBN-0SORBN -9.040667e-01 -2.0375241 0.22939081 0.2079130
1.25SORBP-0SORBN -7.336667e-02 -1.2068241 1.06009081 1.0000000
2.5SORBN-0SORBN -1.547792e+00 -2.6812491 -0.41433419 0.0022669
2.5SORBP-0SORBN -1.536542e+00 -2.6699991 -0.40308419 0.0024767
5SORBN-0SORBN -6.436667e-01 -1.7771241 0.48979081 0.6392249
5SORBP-0SORBN -1.408092e+00 -2.5415491 -0.27463419 0.0067237
FSORBN-0SORBN -6.577917e-01 -1.7912491 0.47566581 0.6121097
FSORBP-0SORBN -6.577917e-01 -1.7912491 0.47566581 0.6121097
1.25SORBN-0SORBP -4.259333e-01 -1.5593908 0.70752414 0.9481188
1.25SORBP-0SORBP 4.047667e-01 -0.7286908 1.53822414 0.9617652
2.5SORBN-0SORBP -1.069658e+00 -2.2031158 0.06379914 0.0765719
2.5SORBP-0SORBP -1.058408e+00 -2.1918658 0.07504914 0.0823756
5SORBN-0SORBP -1.655333e-01 -1.2989908 0.96792414 0.9999546
5SORBP-0SORBP -9.299583e-01 -2.0634158 0.20349914 0.1799850
FSORBN-0SORBP -1.796583e-01 -1.3131158 0.95379914 0.9999099
FSORBP-0SORBP -1.796583e-01 -1.3131158 0.95379914 0.9999099
1.25SORBP-1.25SORBN 8.307000e-01 -0.2186777 1.88007771 0.2157709
2.5SORBN-1.25SORBN -6.437250e-01 -1.6931027 0.40565271 0.5397922
2.5SORBP-1.25SORBN -6.324750e-01 -1.6818527 0.41690271 0.5631718
5SORBN-1.25SORBN 2.604000e-01 -0.7889777 1.30977771 0.9968266
5SORBP-1.25SORBN -5.040250e-01 -1.5534027 0.54535271 0.8140109
FSORBN-1.25SORBN 2.462750e-01 -0.8031027 1.29565271 0.9979144
FSORBP-1.25SORBN 2.462750e-01 -0.8031027 1.29565271 0.9979144
2.5SORBN-1.25SORBP -1.474425e+00 -2.5238027 -0.42504729 0.0015916
2.5SORBP-1.25SORBP -1.463175e+00 -2.5125527 -0.41379729 0.0017523
5SORBN-1.25SORBP -5.703000e-01 -1.6196777 0.47907771 0.6914078
5SORBP-1.25SORBP -1.334725e+00 -2.3841027 -0.28534729 0.0051919
FSORBN-1.25SORBP -5.844250e-01 -1.6338027 0.46495271 0.6628334
FSORBP-1.25SORBP -5.844250e-01 -1.6338027 0.46495271 0.6628334
2.5SORBP-2.5SORBN 1.125000e-02 -1.0381277 1.06062771 1.0000000
5SORBN-2.5SORBN 9.041250e-01 -0.1452527 1.95350271 0.1371789
5SORBP-2.5SORBN 1.397000e-01 -0.9096777 1.18907771 0.9999793
FSORBN-2.5SORBN 8.900000e-01 -0.1593777 1.93937771 0.1501757
FSORBP-2.5SORBN 8.900000e-01 -0.1593777 1.93937771 0.1501757
5SORBN-2.5SORBP 8.928750e-01 -0.1565027 1.94225271 0.1474525
5SORBP-2.5SORBP 1.284500e-01 -0.9209277 1.17782771 0.9999899
FSORBN-2.5SORBP 8.787500e-01 -0.1706277 1.92812771 0.1612234
FSORBP-2.5SORBP 8.787500e-01 -0.1706277 1.92812771 0.1612234
5SORBP-5SORBN -7.644250e-01 -1.8138027 0.28495271 0.3116212
FSORBN-5SORBN -1.412500e-02 -1.0635027 1.03525271 1.0000000
FSORBP-5SORBN -1.412500e-02 -1.0635027 1.03525271 1.0000000
FSORBN-5SORBP 7.503000e-01 -0.2990777 1.79967771 0.3351184
FSORBP-5SORBP 7.503000e-01 -0.2990777 1.79967771 0.3351184
FSORBP-FSORBN 4.440892e-16 -1.0493777 1.04937771 1.0000000Scheirer-Ray-Hare test and Dunn’s post-hoc test
pH
Show the code
# Scheirer-Ray-Hare Test
scheirerRayHare(PH ~ SUBST * DENSITY, data = physchem)
DV: PH
Observations: 38
D: 0.9986869
MS total: 123.5 Df Sum Sq H p.value
SUBST 1 28.7 0.232 0.62978
DENSITY 4 4338.2 35.174 0.00000
SUBST:DENSITY 4 21.1 0.171 0.99655
Residuals 28 175.5 Show the code
# Dunn's post-hoc test
dunnTest(PH ~ GROUP, data = physchem, method = 'bh') Comparison Z P.unadj P.adj
1 0SORBN - 0SORBP 0.33084054 0.7407649384 0.793676720
2 0SORBN - 1.25SORBN -1.82736430 0.0676450135 0.152201280
3 0SORBP - 1.25SORBN -2.18104771 0.0291798877 0.093792496
4 0SORBN - 1.25SORBP -1.35578642 0.1751671349 0.315300843
5 0SORBP - 1.25SORBP -1.70946983 0.0873639560 0.178699001
6 1.25SORBN - 1.25SORBP 0.50936233 0.6104982754 0.670059083
7 0SORBN - 2.5SORBN -2.62315198 0.0087120401 0.035640164
8 0SORBP - 2.5SORBN -2.97683539 0.0029124034 0.018722593
9 1.25SORBN - 2.5SORBN -0.85954893 0.3900377356 0.516226415
10 1.25SORBP - 2.5SORBN -1.36891126 0.1710270147 0.320675653
11 0SORBN - 2.5SORBP -2.53473112 0.0112533664 0.042200124
12 0SORBP - 2.5SORBP -2.88841454 0.0038718923 0.017423515
13 1.25SORBN - 2.5SORBP -0.76404350 0.4448413229 0.526785777
14 1.25SORBP - 2.5SORBP -1.27340583 0.2028740711 0.326047614
15 2.5SORBN - 2.5SORBP 0.09550544 0.9239133720 0.944911403
16 0SORBN - 5SORBN -3.50736051 0.0004525756 0.006788633
17 0SORBP - 5SORBN -3.86104392 0.0001129036 0.005080662
18 1.25SORBN - 5SORBN -1.81460330 0.0695848947 0.149110489
19 1.25SORBP - 5SORBN -2.32396563 0.0201273361 0.069671548
20 2.5SORBN - 5SORBN -0.95505437 0.3395502029 0.477492473
21 2.5SORBP - 5SORBN -1.05055981 0.2934608094 0.455370222
22 0SORBN - 5SORBP -3.35999242 0.0007794461 0.008768769
23 0SORBP - 5SORBP -3.71367583 0.0002042704 0.004596084
24 1.25SORBN - 5SORBP -1.65542757 0.0978377966 0.191421776
25 1.25SORBP - 5SORBP -2.16478990 0.0304037783 0.091211335
26 2.5SORBN - 5SORBP -0.79587864 0.4261025769 0.547846170
27 2.5SORBP - 5SORBP -0.89138408 0.3727231565 0.508258850
28 5SORBN - 5SORBP 0.15917573 0.8735304261 0.914159748
29 0SORBN - FSORBN -0.64841959 0.5167136033 0.581302804
30 0SORBP - FSORBN -1.00210300 0.3162938480 0.459136231
31 1.25SORBN - FSORBN 1.27340583 0.2028740711 0.338123452
32 1.25SORBP - FSORBN 0.76404350 0.4448413229 0.541023231
33 2.5SORBN - FSORBN 2.13295476 0.0329284377 0.087163512
34 2.5SORBP - FSORBN 2.03744932 0.0416050356 0.098538242
35 5SORBN - FSORBN 3.08800913 0.0020150228 0.015112671
36 5SORBP - FSORBN 2.92883340 0.0034023673 0.017011837
37 0SORBN - FSORBP -0.64841959 0.5167136033 0.596208004
38 0SORBP - FSORBP -1.00210300 0.3162938480 0.474440772
39 1.25SORBN - FSORBP 1.27340583 0.2028740711 0.351128200
40 1.25SORBP - FSORBP 0.76404350 0.4448413229 0.556051654
41 2.5SORBN - FSORBP 2.13295476 0.0329284377 0.092611231
42 2.5SORBP - FSORBP 2.03744932 0.0416050356 0.104012589
43 5SORBN - FSORBP 3.08800913 0.0020150228 0.018135205
44 5SORBP - FSORBP 2.92883340 0.0034023673 0.019138316
45 FSORBN - FSORBP 0.00000000 1.0000000000 1.000000000H2O loss
Show the code
# Scheirer-Ray-Hare Test
scheirerRayHare(H2O_REDUCT_REL ~ SUBST * DENSITY, data = physchem)
DV: H2O_REDUCT_REL
Observations: 30
D: 1
MS total: 77.5 Df Sum Sq H p.value
SUBST 1 73.63 0.9501 0.32969
DENSITY 3 2104.00 27.1484 0.00001
SUBST:DENSITY 3 8.37 0.1080 0.99087
Residuals 22 61.50 Show the code
# Dunn's post-hoc test
dunnTest(H2O_REDUCT_REL ~ GROUP, data = physchem, method = 'bh') Comparison Z P.unadj P.adj
1 0SORBN - 0SORBP 0.4173650 6.764114e-01 0.701463723
2 0SORBN - 1.25SORBN -3.4950920 4.738983e-04 0.006634577
3 0SORBP - 1.25SORBN -3.9412740 8.104999e-05 0.002269400
4 0SORBN - 1.25SORBP -2.9001827 3.729452e-03 0.020884929
5 0SORBP - 1.25SORBP -3.3463647 8.187863e-04 0.007642005
6 1.25SORBN - 1.25SORBP 0.6425755 5.204996e-01 0.560538040
7 0SORBN - 2.5SORBN -2.3052735 2.115126e-02 0.065803930
8 0SORBP - 2.5SORBN -2.7514554 5.933110e-03 0.023732438
9 1.25SORBN - 2.5SORBN 1.2851509 1.987395e-01 0.309150409
10 1.25SORBP - 2.5SORBN 0.6425755 5.204996e-01 0.582959562
11 0SORBN - 2.5SORBP -1.7103642 8.719855e-02 0.174397102
12 0SORBP - 2.5SORBP -2.1565461 3.104105e-02 0.086914926
13 1.25SORBN - 2.5SORBP 1.9277264 5.388917e-02 0.125741392
14 1.25SORBP - 2.5SORBP 1.2851509 1.987395e-01 0.327335727
15 2.5SORBN - 2.5SORBP 0.6425755 5.204996e-01 0.607249543
16 0SORBN - 5SORBN -0.9295457 3.526063e-01 0.470141770
17 0SORBP - 5SORBN -1.3757277 1.689060e-01 0.295585413
18 1.25SORBN - 5SORBN 2.7711067 5.586612e-03 0.026070856
19 1.25SORBP - 5SORBN 2.1285312 3.329306e-02 0.084745982
20 2.5SORBN - 5SORBN 1.4859558 1.372908e-01 0.256276204
21 2.5SORBP - 5SORBN 0.8433803 3.990158e-01 0.507838276
22 0SORBN - 5SORBP -0.7064548 4.799054e-01 0.584232626
23 0SORBP - 5SORBP -1.1526367 2.490595e-01 0.367035084
24 1.25SORBN - 5SORBP 3.0120725 2.594706e-03 0.018162944
25 1.25SORBP - 5SORBP 2.3694970 1.781230e-02 0.062343045
26 2.5SORBN - 5SORBP 1.7269216 8.418175e-02 0.181314543
27 2.5SORBP - 5SORBP 1.0843461 2.782114e-01 0.389495925
28 5SORBN - 5SORBP 0.2409658 8.095816e-01 0.809581624VS reduction
Show the code
# Scheirer-Ray-Hare Test
scheirerRayHare(VS_REDUCT_REL ~ SUBST * DENSITY, data = physchem)
DV: VS_REDUCT_REL
Observations: 30
D: 1
MS total: 77.5 Df Sum Sq H p.value
SUBST 1 381.63 4.9243 0.02648
DENSITY 3 1427.25 18.4161 0.00036
SUBST:DENSITY 3 248.78 3.2101 0.36035
Residuals 22 189.83 Show the code
# Dunn's post-hoc test
dunnTest(VS_REDUCT_REL ~ GROUP, data = physchem, method = 'bh') Comparison Z P.unadj P.adj
1 0SORBN - 0SORBP 0.3246172 7.454708e-01 0.773080817
2 0SORBN - 1.25SORBN -1.9458491 5.167287e-02 0.120570026
3 0SORBP - 1.25SORBN -2.2928795 2.185494e-02 0.067993155
4 0SORBN - 1.25SORBP -0.5701214 5.685954e-01 0.636826810
5 0SORBP - 1.25SORBP -0.9171518 3.590631e-01 0.478750796
6 1.25SORBN - 1.25SORBP 1.4859558 1.372908e-01 0.256276204
7 0SORBN - 2.5SORBN -3.5446678 3.931082e-04 0.005503515
8 0SORBP - 2.5SORBN -3.8916982 9.954502e-05 0.002787261
9 1.25SORBN - 2.5SORBN -1.7269216 8.418175e-02 0.181314543
10 1.25SORBP - 2.5SORBN -3.2128773 1.314124e-03 0.012265158
11 0SORBN - 2.5SORBP -1.3137580 1.889277e-01 0.311174967
12 0SORBP - 2.5SORBP -1.6607884 9.675595e-02 0.193511901
13 1.25SORBN - 2.5SORBP 0.6827364 4.947734e-01 0.602332842
14 1.25SORBP - 2.5SORBP -0.8032193 4.218480e-01 0.536897423
15 2.5SORBN - 2.5SORBP 2.4096580 1.596748e-02 0.055886184
16 0SORBN - 5SORBN -2.5407584 1.106123e-02 0.051619087
17 0SORBP - 5SORBN -2.8877888 3.879603e-03 0.027157220
18 1.25SORBN - 5SORBN -0.6425755 5.204996e-01 0.607249543
19 1.25SORBP - 5SORBN -2.1285312 3.329306e-02 0.093220580
20 2.5SORBN - 5SORBN 1.0843461 2.782114e-01 0.389495925
21 2.5SORBP - 5SORBN -1.3253119 1.850677e-01 0.323868521
22 0SORBN - 5SORBP -2.4292129 1.513164e-02 0.060526577
23 0SORBP - 5SORBP -2.7762433 5.499104e-03 0.030794982
24 1.25SORBN - 5SORBP -0.5220926 6.016059e-01 0.647883259
25 1.25SORBP - 5SORBP -2.0080483 4.463816e-02 0.113624400
26 2.5SORBN - 5SORBP 1.2048290 2.282693e-01 0.336396900
27 2.5SORBP - 5SORBP -1.2048290 2.282693e-01 0.355085616
28 5SORBN - 5SORBP 0.1204829 9.041006e-01 0.904100625Carbon content
Show the code
# Scheirer-Ray-Hare Test
scheirerRayHare(C_REL ~ SUBST * DENSITY, data = physchem)
DV: C_REL
Observations: 38
D: 0.9995623
MS total: 123.5 Df Sum Sq H p.value
SUBST 1 79.61 0.6449 0.42196
DENSITY 4 3059.67 24.7855 0.00006
SUBST:DENSITY 4 298.06 2.4145 0.66001
Residuals 28 1130.17 Show the code
# Dunn's post-hoc test
dunnTest(C_REL ~ GROUP, data = physchem, method = 'bh') Comparison Z P.unadj P.adj
1 0SORBN - 0SORBP 0.69813524 0.485092622 0.77961314
2 0SORBN - 1.25SORBN 2.41577815 0.015701625 0.07850813
3 0SORBP - 1.25SORBN 1.66944018 0.095030176 0.18592861
4 0SORBN - 1.25SORBP 1.06058553 0.288878303 0.49998168
5 0SORBP - 1.25SORBP 0.31424756 0.753333020 0.89210489
6 1.25SORBN - 1.25SORBP -1.46377533 0.143255346 0.25785962
7 0SORBN - 2.5SORBN 3.06391375 0.002184619 0.09830785
8 0SORBP - 2.5SORBN 2.31757578 0.020472388 0.07677145
9 1.25SORBN - 2.5SORBN 0.70006646 0.483885800 0.80647633
10 1.25SORBP - 2.5SORBN 2.16384179 0.030476494 0.08571514
11 0SORBN - 2.5SORBP 2.53362098 0.011289076 0.06350105
12 0SORBP - 2.5SORBP 1.78728302 0.073891757 0.15833948
13 1.25SORBN - 2.5SORBP 0.12728481 0.898714982 0.94051568
14 1.25SORBP - 2.5SORBP 1.59106014 0.111596043 0.20924258
15 2.5SORBN - 2.5SORBP -0.57278165 0.566792546 0.79705202
16 0SORBN - 5SORBN 2.97553162 0.002924812 0.06580827
17 0SORBP - 5SORBN 2.22919365 0.025801023 0.08931123
18 1.25SORBN - 5SORBN 0.60460285 0.545442911 0.79177197
19 1.25SORBP - 5SORBN 2.06837818 0.038604476 0.09143165
20 2.5SORBN - 5SORBN -0.09546361 0.923946595 0.94494538
21 2.5SORBP - 5SORBN 0.47731804 0.633135666 0.86336682
22 0SORBN - 5SORBP 2.68092453 0.007341907 0.04719798
23 0SORBP - 5SORBP 1.93458656 0.053041050 0.11934236
24 1.25SORBN - 5SORBP 0.28639082 0.774578801 0.89374477
25 1.25SORBP - 5SORBP 1.75016615 0.080089648 0.16381973
26 2.5SORBN - 5SORBP -0.41367564 0.679111673 0.84888959
27 2.5SORBP - 5SORBP 0.15910601 0.873585350 0.93598430
28 5SORBN - 5SORBP -0.31821203 0.750324110 0.91255635
29 0SORBN - FSORBN 0.47137135 0.637375575 0.81948288
30 0SORBP - FSORBN -0.27496662 0.783341885 0.85976548
31 1.25SORBN - FSORBN -2.10019938 0.035711306 0.08927826
32 1.25SORBP - FSORBN -0.63642405 0.524500066 0.78675010
33 2.5SORBN - FSORBN -2.80026584 0.005106054 0.05744310
34 2.5SORBP - FSORBN -2.22748419 0.025914931 0.07774479
35 5SORBN - FSORBN -2.70480223 0.006834507 0.05125880
36 5SORBP - FSORBN -2.38659021 0.017005435 0.06956769
37 0SORBN - FSORBP 0.47137135 0.637375575 0.84358532
38 0SORBP - FSORBP -0.27496662 0.783341885 0.88125962
39 1.25SORBN - FSORBP -2.10019938 0.035711306 0.09452993
40 1.25SORBP - FSORBP -0.63642405 0.524500066 0.81387941
41 2.5SORBN - FSORBP -2.80026584 0.005106054 0.07659081
42 2.5SORBP - FSORBP -2.22748419 0.025914931 0.08329799
43 5SORBN - FSORBP -2.70480223 0.006834507 0.06151056
44 5SORBP - FSORBP -2.38659021 0.017005435 0.07652446
45 FSORBN - FSORBP 0.00000000 1.000000000 1.00000000C:N ratio
Show the code
# Scheirer-Ray-Hare Test
scheirerRayHare(CN_ratio ~ SUBST * DENSITY, data = physchem)
DV: CN_ratio
Observations: 38
D: 0.9995623
MS total: 123.5 Df Sum Sq H p.value
SUBST 1 111.18 0.9007 0.34260
DENSITY 4 2198.92 17.8128 0.00134
SUBST:DENSITY 4 849.73 6.8834 0.14218
Residuals 28 1407.67 Show the code
# Dunn's post-hoc test
dunnTest(CN_ratio ~ GROUP, data = physchem, method = 'bh') Comparison Z P.unadj P.adj
1 0SORBN - 0SORBP -0.6981352 0.4850926223 0.64203435
2 0SORBN - 1.25SORBN -1.2078891 0.2270899333 0.40876188
3 0SORBP - 1.25SORBN -0.4615511 0.6444032620 0.76310913
4 0SORBN - 1.25SORBP 0.0000000 1.0000000000 1.00000000
5 0SORBP - 1.25SORBP 0.7463380 0.4554632838 0.62108630
6 1.25SORBN - 1.25SORBP 1.3046693 0.1920054754 0.36001027
7 0SORBN - 2.5SORBN -2.7103852 0.0067205101 0.05040383
8 0SORBP - 2.5SORBN -1.9640473 0.0495246100 0.18571729
9 1.25SORBN - 2.5SORBN -1.6228813 0.1046147733 0.26153693
10 1.25SORBP - 2.5SORBN -2.9275507 0.0034164341 0.03843488
11 0SORBN - 2.5SORBP -3.1228352 0.0017911803 0.04030156
12 0SORBP - 2.5SORBP -2.3764972 0.0174778940 0.08738947
13 1.25SORBN - 2.5SORBP -2.0683782 0.0386044756 0.15792740
14 1.25SORBP - 2.5SORBP -3.3730475 0.0007434113 0.03345351
15 2.5SORBN - 2.5SORBP -0.4454968 0.6559607441 0.75687778
16 0SORBN - 5SORBN -0.3829892 0.7017277583 0.77018900
17 0SORBP - 5SORBN 0.3633487 0.7163443842 0.76751184
18 1.25SORBN - 5SORBN 0.8909937 0.3729325626 0.55939884
19 1.25SORBP - 5SORBN -0.4136756 0.6791116728 0.76400063
20 2.5SORBN - 5SORBN 2.5138750 0.0119412744 0.06716967
21 2.5SORBP - 5SORBN 2.9593719 0.0030826688 0.04624003
22 0SORBN - 5SORBP -2.5336210 0.0112890763 0.07257263
23 0SORBP - 5SORBP -1.7872830 0.0738917573 0.22167527
24 1.25SORBN - 5SORBP -1.4319541 0.1521569512 0.29769838
25 1.25SORBP - 5SORBP -2.7366234 0.0062073303 0.05586597
26 2.5SORBN - 5SORBP 0.1909272 0.8485826168 0.88805158
27 2.5SORBP - 5SORBP 0.6364241 0.5245000665 0.67435723
28 5SORBN - 5SORBP -2.3229478 0.0201819568 0.09081881
29 0SORBN - FSORBN -1.7381818 0.0821787785 0.21753206
30 0SORBP - FSORBN -0.9918439 0.3212736922 0.49852814
31 1.25SORBN - FSORBN -0.5727816 0.5667925456 0.68934229
32 1.25SORBP - FSORBN -1.8774510 0.0604563124 0.19432386
33 2.5SORBN - FSORBN 1.0500997 0.2936722811 0.48945380
34 2.5SORBP - FSORBN 1.4955965 0.1347588283 0.30320736
35 5SORBN - FSORBN -1.4637753 0.1432553463 0.29302230
36 5SORBP - FSORBN 0.8591725 0.3902453674 0.54878255
37 0SORBN - FSORBP -1.7381818 0.0821787785 0.23112781
38 0SORBP - FSORBP -0.9918439 0.3212736922 0.51633272
39 1.25SORBN - FSORBP -0.5727816 0.5667925456 0.70849068
40 1.25SORBP - FSORBP -1.8774510 0.0604563124 0.20927185
41 2.5SORBN - FSORBP 1.0500997 0.2936722811 0.50827895
42 2.5SORBP - FSORBP 1.4955965 0.1347588283 0.31916565
43 5SORBN - FSORBP -1.4637753 0.1432553463 0.30697574
44 5SORBP - FSORBP 0.8591725 0.3902453674 0.56648521
45 FSORBN - FSORBP 0.0000000 1.0000000000 1.00000000Funding
This research was funded in whole or in part by the Austrian Science Fund (FWF) [grant DOI: 10.55776/P35401]. For open access purposes, the author has applied a CC BY public copyright license to any author accepted manuscript version arising from this submission.