test/new/extreme_simulation.R

In DiogoFerrari/eforensics:

### Extreme Data ###

model = 'bl'
sim_data = ef_simulateData(n=10000, nCov=20, model=model)
data = sim_data$data
mcmc = list(burn.in=100, n.adapt=100, n.iter=1000, n.chains=2)
samples = eforensics(w ~.-a, a ~ .-w,  data=data, model=model, mcmc=mcmc)
print(summary(samples));