The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.
Some geese isotope data is included with this package. Find where it is with:
Load into R with:
library(readxl)
path <- system.file("extdata", "geese_data.xls", package = "simmr")
geese_data <- lapply(excel_sheets(path), read_excel, path = path)
If you want to see what the original Excel sheet looks like you can
run system(paste('open',path))
.
We can now separate out the data into parts
targets <- geese_data[[1]]
sources <- geese_data[[2]]
TEFs <- geese_data[[3]]
concdep <- geese_data[[4]]
Note that if you don’t have TEFs or concentration dependence you can set these all to the value 0 or just leave them blank in the step below.
simmr
simmr
and check convergencegeese_simmr_out <- simmr_mcmc(geese_simmr)
summary(geese_simmr_out,
type = "diagnostics",
group = 1
)
Check that the model fitted well:
## Compiling model graph
## Resolving undeclared variables
## Allocating nodes
## Graph information:
## Observed stochastic nodes: 40
## Unobserved stochastic nodes: 46
## Total graph size: 198
##
## Initializing model
Look at the influence of the prior:
Look at the histogram of the dietary proportions:
## Most popular orderings are as follows:
## Probability
## Day 428 > Day 124 > Day 398 > Day 1 0.2011
## Day 428 > Day 398 > Day 124 > Day 1 0.1811
## Day 428 > Day 124 > Day 1 > Day 398 0.1742
## Day 428 > Day 398 > Day 1 > Day 124 0.0950
## Day 428 > Day 1 > Day 124 > Day 398 0.0772
## Day 428 > Day 1 > Day 398 > Day 124 0.0667
## Day 398 > Day 428 > Day 124 > Day 1 0.0467
## Day 124 > Day 428 > Day 398 > Day 1 0.0425
## Day 124 > Day 428 > Day 1 > Day 398 0.0286
## Day 398 > Day 428 > Day 1 > Day 124 0.0222
## Day 398 > Day 124 > Day 428 > Day 1 0.0131
## Day 124 > Day 398 > Day 428 > Day 1 0.0128
## Day 1 > Day 428 > Day 124 > Day 398 0.0083
## Day 1 > Day 428 > Day 398 > Day 124 0.0064
## Day 124 > Day 1 > Day 428 > Day 398 0.0058
## Day 1 > Day 124 > Day 428 > Day 398 0.0050
## Day 398 > Day 1 > Day 428 > Day 124 0.0044
## Day 1 > Day 398 > Day 428 > Day 124 0.0028
## Day 124 > Day 398 > Day 1 > Day 428 0.0017
## Day 398 > Day 124 > Day 1 > Day 428 0.0017
## Day 124 > Day 1 > Day 398 > Day 428 0.0011
## Day 398 > Day 1 > Day 124 > Day 428 0.0008
## Day 1 > Day 124 > Day 398 > Day 428 0.0006
## Day 1 > Day 398 > Day 124 > Day 428 0.0003
For the many more options available to run and analyse output, see
the main vignette via vignette('simmr')
These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.