Basic plots of BGmix parameters and data.

Share:

Description

Plots gene effects and variances versus their corresponding data sufficient statistics (to show the effect of smoothing and shrinkage). Also plots "volcano plots": posterior probabilities of being classified in each mixture component versus the log fold change parameters.

Usage

1
plotBasic(res, ybar, ss, q.mean = T, q.diff = T, q.sig = T, q.volcano = T)

Arguments

res

list object output from 'ccParams'

ybar

ybar data (see BGmix help for details)

ss

ss data (see BGmix help for details)

q.mean

logical. Include mean plot?

q.diff

logical. Include log fold change plot?

q.sig

logical. Include variance plot?

q.volcano

logical. Include volcana plot (posterior classification v. fold change)?

Details

Note this plotting function is designed for model output from the unpaired differential expression design.

Value

No value is returned to R. Results from BGmix model are output to files.

Author(s)

Alex Lewin

Examples

1
2
3
4
5
6
## Note this is a very short MCMC run!
## For good analysis need proper burn-in period.
data(ybar,ss)
outdir <- BGmix(ybar, ss, c(8,8), nburn=0, niter=100, nthin=1)
params <- ccParams(outdir)
plotBasic(params,ybar,ss)