qlfWrapper: Wrapper function to easily perform a glmQLFTest analysis on...

In Dudemanguy911/RNASeqSuite: Convenient R wrapper of popular RNAseq pipelines designed to simplify data analysis

Description

Performs a glmQLFTest and stores the results in a DGELRT object.

Usage

qlfWrapper(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

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

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