gene_compare_with_weights: compares the mean expression between two groups of genes in a...

Description Usage Arguments Value Examples

Description

compares the mean expression between two groups of genes in a given expression matrix, optionally makes boxplots of the comparison, faceted by gene used in a pipeline, so checks are pretty rudimentary

Usage

1
2
3
4
5
6
7
8
9
gene_compare_with_weights(
  dat,
  camera2output,
  genesetidx,
  grp1name = "group1",
  grp2name = "group2",
  grp1idx,
  grp2idx
)

Arguments

dat

A matrix of gene expression

camera2output

output matrix from camera2()

genesetidx

a list of lists of gene indices

grp1name

String indicating name of the first group for the comparison

grp2name

String indicating name of the second group for the comparison

grp1idx

A vector of integer indices indicating columns belonging to group1

grp2idx

A vector of integer indices indicating columns belonging to group1

Value

table of genes' mean expression values for each group

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
## Not run: 
out <- gene_compare(expressionMatrix
                   ,c("LAP3","NCOA7","IFIT2","STAT1","TRIM21")
                   ,grp1name = "control"
                   ,grp2name = "experimental"
                   ,which(grepl("control", colnames(expressionMatrix)))
                   ,which(grepl("experimental", colnames(expressionMatrix)))
                                  )

## End(Not run)

ptvan/Pmisc documentation built on Nov. 19, 2020, 10:27 p.m.