plotConvergence {birta} | R Documentation |
Plots the log likelihood along MCMC sampling.
plotConvergence(res, nburnin=NULL, title="")
res |
The result from birta.run (a list). |
nburnin |
Number of iterations used for the burn in. |
title |
Optional title of the plot. |
Benedikt Zacher zacher@lmb.uni-muenchen.de
data(humanSim) data(humanSim) design = model.matrix(~0+factor(c(rep("control", 5), rep("treated", 5)))) colnames(design) = c("control", "treated") contrasts = "treated - control" limmamRNA = limmaAnalysis(sim$dat.mRNA, design, contrasts) limmamiRNA = limmaAnalysis(sim$dat.miRNA, design, contrasts) sim_result = birta(sim$dat.mRNA, sim$dat.miRNA, limmamRNA=limmamRNA, limmamiRNA=limmamiRNA, nrep=c(5,5,5,5), genesets=genesets, model="all-plug-in", niter=50000, nburnin=10000, sample.weights=FALSE, potential_swaps=potential_swaps) plotConvergence(sim_result, nburnin=10000, title="simulation")