Plot a barplot of AIC weights for each model and each gene group

Share:

Description

This function will plot a barplot where each bar corresponds to a gene group and the proportion of each bar that is filled with a certain color corresponds to the AIC weight for a given model. Bars are sorted according to which model has the highest weight.

Usage

1
2
plot_wAICbarplot(dat, names, col = 2:(length(names) + 1), title = "",
  cex = 1.4, border = par("fg"), space = NULL, names.arg = 1:nrow(dat))

Arguments

dat

a wAIC matrix, rows are be gene groups and columns are models

names

a vector of names for each model

col

a vector of colors, one color corresponding to each model

title

name for the barplot

cex

character expansion factor

border

type of border around each bar in the barplots

space

amount of space between bars in the barplot

names.arg

names of bars in the barplot

Value

invisibly returns the ordering of gene groups in the barplot.

Examples

1
2
3
4
5
6
7
8
## Not run: 
data(yeast)
GO.groups.pruned = good.groups(colnames(yeast.homozygote),GO.groups,30)
test_groups = GO.groups[GO.groups.pruned[1:100]]
yeast.test = test.groups(yeast.tree,yeast.homozygote,test_groups,print_names=T)
plot_wAICbarplot(mytest$wAIC,1:7)

## End(Not run)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.