exactWrapper: Wrapper function to easily perform an exactTest analysis on...

Description Usage Arguments Details Examples

View source: R/RNASeqSuite.R

Description

Performs an exactTest and stores the results in a DGEExact object.

Usage

1
exactWrapper(data, group, htsfilter=TRUE, cfilter=0, cutoff=0, adjust.method="BH", sort.by="FDR", decreasing=FALSE)

Arguments

data

A data frame of the raw counts

group

The factor created from grpSelection containing the different group names

htsfilter

A boolean determining whether or not to use HTSFilter to filter count matrix

cfilter

A numeric value controlling the use of cfilter, a custom outlier detector

cutoff

Filters based on a hard cutoff of read counts

adjust.method

The method used to adjust p-values. Valid options are "holm", "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr", and "none"

sort.by

What criteria to sort the DGEExact object. Valid options are "logFC", "logCPM", "F", "PValue", "FDR", and "none"

decreasing

A boolean to tell sort to ascend or descend across values upon sort

Details

htsfilter calls the HTSFilter library if set to TRUE. This will remove genes based on a similarity ranking via a Jaccard index. Genes with constant, low counts will be removed.

cfilter is a custom outlier detector. cfilter normalizes the count vectors and then takes the difference between the two groups. The standard deviation and mean of the difference vector is stored and the genes are subsequently filtered based on the value of cfilter. Gene count values that lie outside of the standard deviation times the cfilter will be filtered. The lower the value of cfilter, the more stringent the filter will be. The default value of 0 disables this filter.

cutoff is a hard count cutoff. The average of every group is calculated and if the average for all groups falls below the specified cutoff value, the gene is removed. The default value of 0 disables this filter.

Examples

1
y <- exactWrapper(data, group, htsfilter=TRUE, cfilter=2, cutoff=100)

Dudemanguy911/RNASeqSuite documentation built on Dec. 9, 2019, 8:43 a.m.