plotConvergence {birta}R Documentation

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.

Author(s)

Benedikt Zacher zacher@lmb.uni-muenchen.de

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

[Package birta version 1.26.0 Index]