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

Description Usage Arguments Details Examples

View source: R/RNASeqSuite.R


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


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



A data frame of the raw counts


The factor created from grpSelection containing the different group names


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


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


Filters based on a hard cutoff of read counts


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


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


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


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.


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

Dudemanguy/RNASeqSuite documentation built on Dec. 5, 2019, 12:45 a.m.