In DiogoFerrari/eforensics:

Setting up the library

library(eforensics);

I then run the help file just to see what happens.

## model
## -----
model    = 'bl'

## simulating data and parameters
## ------------------------------
help(ef_simulateData)
sim_data = ef_simulateData(n=7000, nCov=1, model=model)
data     = sim_data$data

## mcmc parameters
## ---------------
mcmc    = list(burnin=1, n.adapt=100, n.iter=100, n.chains=2)

## samples
## -------
help(eforensics)
samples    = eforensics(w ~.-a, a ~ .-w,  data=data, model=model, mcmc=mcmc)

summary(samples)
summary(samples, join.chains=T)

Notice that this throws an error: Error in check_mcmc(mcmc)'' This is from the fact thatburnin'' should probably be ``burn.in''. Now, if I correct this, this is what happens.

## model
## -----
model    = 'bl'

## simulating data and parameters
## ------------------------------
help(ef_simulateData)
sim_data = ef_simulateData(n=7000, nCov=1, model=model)
data     = sim_data$data

## mcmc parameters
## ---------------
mcmc = list(burn.in=1, n.adapt=100, n.iter=100, n.chains=2)

## samples
## -------
help(eforensics)
samples    = eforensics(w ~.-a, a ~ .-w,  data=data, model=model, mcmc=mcmc)

summary(samples)
summary(samples, join.chains=T)

DiogoFerrari/eforensics documentation built on Aug. 31, 2019, 3:44 p.m.