birta (version 1.16.0)

plotConvergence: Plotting the likelihood along MCMC sampling.

Description

Plots the log likelihood along MCMC sampling.

Usage

plotConvergence(res, nburnin=NULL, title="")

Arguments

res

The result from birta.run (a list).

nburnin

Number of iterations used for the burn in.

title

Optional title of the plot.

See Also

birta

Examples

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")