R Under development (unstable) (2023-12-07 r85661 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2023 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library('bayescount')
> 
> data <- data.frame(Count=rpois(80,rep(c(10,10,2,2), 20)),
+ Subject=rep(1:20, each=4), Time=rep(rep(1:2,each=2),40),
+ Sample=1:2, Control=rep(c(0,1), each=40))
> # Compile the model - a paired model is required because
> # there are replicate samples within an individual:
> model <- fecrt.model(data, paired.model=TRUE)
module runjags loaded
Returning compiled model
> # Update the model - requires runjags:
> library('runjags')
> results <- extend.jags(model, burnin=5000)
Compiling rjags model...
Calling the simulation using the rjags method...
Adapting the model for 1000 iterations...
Burning in the model for 5000 iterations...
Running the model for 10000 iterations...
Simulation complete
Calculating summary statistics...
Calculating the Gelman-Rubin statistic for 6 variables....
Finished running the simulation
> 
> proc.time()
   user  system elapsed 
  33.09    0.35   35.43