Renata Diaz 2021-05-09
To show the 95% interval, we need to load the distribution of shape metric values from the samples from the feasible set for a few communities. See rov_metric.md.
## Joining, by = c("sim", "source", "dat", "site", "singletons", "s0", "n0", "nparts")
## Joining, by = c("dat", "site", "s0", "n0", "nparts")
## Joining, by = c("sim", "source", "dat", "site", "singletons", "s0", "n0", "nparts")
## Joining, by = c("dat", "site", "s0", "n0", "nparts")
## Joining, by = c("sim", "source", "dat", "site", "singletons", "s0", "n0", "nparts")
## Joining, by = c("dat", "site", "s0", "n0", "nparts")
## Warning: Removed 1 rows containing non-finite values (stat_bin).
## Warning: Removed 3 row(s) containing missing values (geom_path).
## Warning: Removed 1 rows containing non-finite values (stat_bin).
## Warning: Removed 1 row(s) containing missing values (geom_path).
## Warning: Removed 1 rows containing non-finite values (stat_bin).
## [1] 6.6e+70
## [1] 1.001277e+12
## [1] 297
## TableGrob (3 x 4) "arrange": 7 grobs
## z cells name grob
## 1 1 (1-1,2-2) arrange gtable[layout]
## 2 2 (1-1,3-3) arrange gtable[layout]
## 3 3 (1-1,4-4) arrange gtable[layout]
## 4 4 (2-2,2-2) arrange gtable[layout]
## 5 5 (2-2,3-3) arrange gtable[layout]
## 6 6 (3-3,2-4) arrange text[GRID.text.1674]
## 7 7 (1-3,1-1) arrange text[GRID.text.1675]
## Joining, by = "Dataset"
## Joining, by = "Dataset"
## Note: Using an external vector in selections is ambiguous.
## ℹ Use `all_of(cols1)` instead of `cols1` to silence this message.
## ℹ See <https://tidyselect.r-lib.org/reference/faq-external-vector.html>.
## This message is displayed once per session.
## Note: Using an external vector in selections is ambiguous.
## ℹ Use `all_of(cols2)` instead of `cols2` to silence this message.
## ℹ See <https://tidyselect.r-lib.org/reference/faq-external-vector.html>.
## This message is displayed once per session.
## Warning: Removed 8 rows containing missing values (geom_point).
## `summarise()` has grouped output by 'Dataset'. You can override using the `.groups` argument.
## `summarise()` has grouped output by 'Dataset'. You can override using the `.groups` argument.
## `summarise()` has grouped output by 'Dataset'. You can override using the `.groups` argument.
## Joining, by = c("Dataset", "Number of elements")
## Joining, by = c("Dataset", "Number of elements")
## Warning: Removed 14 rows containing non-finite values (stat_bin).
## TableGrob (2 x 1) "arrange": 2 grobs
## z cells name grob
## 1 1 (1-1,1-1) arrange gtable[layout]
## 2 2 (2-2,1-1) arrange gtable[layout]
## Joining, by = c("sim", "source", "dat", "site", "singletons", "s0", "n0", "nparts")
## Joining, by = c("sim", "source", "dat", "site", "singletons", "s0", "n0", "nparts")
## Warning: Unknown or uninitialised column: `njks_skew`.
## `summarise()` has grouped output by 'Dataset'. You can override using the `.groups` argument.
## `summarise()` has grouped output by 'Dataset'. You can override using the `.groups` argument.
## Joining, by = c("Dataset", "resampling")
## Joining, by = "ometric2"
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.