diff --git a/docs/data-preparation.html b/docs/data-preparation.html
index 7fdc09f..68504bd 100644
--- a/docs/data-preparation.html
+++ b/docs/data-preparation.html
@@ -358,8 +358,8 @@
diff --git a/docs/introduction.html b/docs/introduction.html
index 87fbcd8..8ead5bd 100644
--- a/docs/introduction.html
+++ b/docs/introduction.html
@@ -261,7 +261,7 @@ Figure 6.3: Illustration to the sensitivity of estimated process rates to a change in prior bulk density.
-
-
-
-
- |
-
-
-
- |
-
-Dependent variable:
- |
-
-
-
- |
-
- |
-
-
-
- |
-
-Nitrification
- |
-
-Denitrification
- |
-
-Reduction
- |
-
-
-
- |
-
-(1)
- |
-
-(2)
- |
-
-(3)
- |
-
-
-
- |
-
-
-
-BD
- |
-
--20.907***
- |
-
--4.672***
- |
-
-31.267***
- |
-
-
-
- |
-
-(0.993)
- |
-
-(0.619)
- |
-
-(0.926)
- |
-
-
-
- |
-
- |
-
- |
-
- |
-
-
-
-Constant
- |
-
-47.404***
- |
-
-26.665***
- |
-
--27.614***
- |
-
-
-
- |
-
-(1.678)
- |
-
-(1.045)
- |
-
-(1.565)
- |
-
-
-
- |
-
- |
-
- |
-
- |
-
-
-
- |
-
-
-
-Observations
- |
-
-7
- |
-
-7
- |
-
-7
- |
-
-
-
-R2
- |
-
-0.989
- |
-
-0.919
- |
-
-0.996
- |
+
+
+
+
+ |
+Nitrification |
+Denitrification |
+Reduction |
-
-
-Adjusted R2
- |
-
-0.987
- |
-
-0.903
- |
-
-0.995
- |
+
+Predictors |
+Estimates |
+p |
+Estimates |
+p |
+Estimates |
+p |
-
-
-Residual Std. Error (df = 5)
- |
-
-0.263
- |
-
-0.164
- |
-
-0.245
- |
+
+(Intercept) |
+47.4 |
+<0.001 |
+26.7 |
+<0.001 |
+-27.6 |
+<0.001 |
-
-
-F Statistic (df = 1; 5)
- |
-
-442.914***
- |
-
-57.032***
- |
-
-1,139.130***
- |
+
+BD |
+-20.9 |
+<0.001 |
+-4.7 |
+0.001 |
+31.3 |
+<0.001 |
-
-
- |
+
+Observations |
+7 |
+7 |
+7 |
-
-
-Note:
- |
-
-p<0.1; p<0.05; p<0.01
- |
+
+R2 / R2 adjusted |
+0.989 / 0.987 |
+0.919 / 0.903 |
+0.996 / 0.995 |
+
+
-
+
+
diff --git a/scripts/sensitivity-analysis/04-bulk-density.R b/scripts/sensitivity-analysis/04-bulk-density.R
index 609d44d..e138b27 100644
--- a/scripts/sensitivity-analysis/04-bulk-density.R
+++ b/scripts/sensitivity-analysis/04-bulk-density.R
@@ -1,4 +1,7 @@
+# attach packages to search path
+library(PRE)
+
# determine a range of bulk density to be searched
BD <- getParameters()$BD + seq(-0.15, 0.15, 0.05)
@@ -9,7 +12,8 @@ f <- function(x) {
getN2ON(parameters = parameters) |>
getMissing() |>
calculateFluxes(parameters = parameters)
- return(longPRE(data, 1, 7.5, n = 15)$processes)
+ result <- t(sapply(getParameters()$depth, function(d) longPRE(data, 1, d, n = 3)$processes))
+ return(data.frame(depth = getParameters()$depth, result))
}
# apply said function to all bulk densities
diff --git a/sensitivity-analysis.qmd b/sensitivity-analysis.qmd
index 43fdd9d..b35f783 100644
--- a/sensitivity-analysis.qmd
+++ b/sensitivity-analysis.qmd
@@ -156,5 +156,6 @@ data <- read.csv(file.path("scripts","sensitivity-analysis","output","bulk-densi
m1 <- lm(Nitrification ~ BD, data)
m2 <- lm(Denitrification ~ BD, data)
m3 <- lm(Reduction ~ BD, data)
-stargazer::stargazer(m1, m2, m3, type = "html")
-```
+#table <- stargazer::stargazer(m1, m2, m3, type = "html")
+sjPlot::tab_model(m1, m2, m3, show.p = TRUE, digits = 1, show.ci = FALSE)
+```
\ No newline at end of file