data(ybar, ss)
nreps <- c(8,8)
## Note this is a very short MCMC run!
## For good analysis need proper burn-in period.
outdir <- BGmix(ybar, ss, nreps, jstar=-1, nburn=0, niter=100, nthin=1)
params <- ccParams(outdir)
res <- ccTrace(outdir)
tpp.res <- TailPP(res, nreps, params, plots = FALSE)
FDR.res = FDRforTailPP(tpp.res$tpp, a1 = params$maa[1],
a2 = params$maa[2], n.rep1=nreps[1], n.rep2=nreps[2], p.cut = 0.8)
Run the code above in your browser using DataCamp Workspace