degPatterns: Make groups of genes using expression profile

Description Usage Arguments Details Value Examples

View source: R/clustering.R

Description

Make groups of genes using expression profile

Usage

1
2
3
degPatterns(ma, metadata, minc = 15, summarize = "group", time = "time",
  col = "condition", reduce = FALSE, cutoff = 0.7, scale = TRUE,
  plot = TRUE, fixy = NULL)

Arguments

ma

log2 normalized count matrix

metadata

data frame with sample information. Rownames should match ma column names row number should be the same length than p-values vector.

minc

integer minimum number of genes in a group that will be return

summarize

character column name in metadata that will be used to gorup replicates. For instance, a merge between summarize and time parameters: control_point0 ... etc

time

character column name in metadata that will be used as variable that changes, normally a time variable.

col

character column name in metadata to separate samples. Normally control/mutant

reduce

boolean reduce number of clusters using correlation values between them.

cutoff

integer threshold for correlation expression to merge clusters (0 - 1)

scale

boolean scale the ma values by row

plot

boolean plot the clusters found

fixy

vector integers used as ylim in plot

Details

It would be used diana function to detect a value to cut the expression based clustering at certain height. It can work with one or more groups with 2 or more several time points. The different patterns can be merged to get similar ones into only one pattern. The expression correlation of the patterns will be used to decide whether some need to be merged or not.

Value

list wiht two items. df is a data.frame with two columns. The first one with genes, the second with the clusters they belong. pass_to_plot is a vector of the clusters that pass the minc cutoff.

Examples

1
2
3
4
5
data(humanSexDEedgeR)
ma <- humanSexDEedgeR$counts[1:100,]
des <- data.frame(row.names=colnames(ma), 
group=as.factor(humanSexDEedgeR$samples$group))
res <- degPatterns(ma, des, time="group", col=NULL)

DEGreport documentation built on May 20, 2017, 10:17 p.m.

Search within the DEGreport package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs in the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.