Skip to content

Instantly share code, notes, and snippets.

@shug3502
Created March 5, 2019 08:53
Show Gist options
  • Save shug3502/8bb72b66ac341e5859a9101d9a8b9024 to your computer and use it in GitHub Desktop.
Save shug3502/8bb72b66ac341e5859a9101d9a8b9024 to your computer and use it in GitHub Desktop.
Minimal example to recreate RBi error
#minimal example to recreate error
library(rbi)
set.seed(123)
model_file <- system.file(package="rbi", "SIR.bi")
SIRmodel <- bi_model(model_file)
SIRdata <- bi_generate_dataset(SIRmodel, end_time=16*7, noutputs=16)
bi <- libbi(SIRmodel)
bi_prior <- sample(bi, target="prior", nsamples=1000, end_time=16*7, noutputs=16)
bi <- sample(bi_prior, target="posterior", nparticles=32, obs=SIRdata, verbose=TRUE)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment