test/helpfile_correct_test.r

library(eforensics);
## 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.