nbinomTestForMatricesSH: Exact test under Negative Binomial Test with Shrinkage...

In sSeq: Shrinkage estimation of dispersion in Negative Binomial models for RNA-seq experiments with small sample size

Description

This is an internal function used by nbTestSH. It calculates the exact per-gene probabilities for p-values, and tests the null hypothesis that the expected expression of a gene under two conditions are not different.

Usage

nbinomTestForMatricesSH(countsA, countsB, sizeFactorsA, sizeFactorsB, 
    numPart=1, SHonly=FALSE, propForSigma=c(0, 1), shrinkTarget=NULL, 
    shrinkQuantile=NULL, cLA, cLB, contrast=NULL, 
    keepLevelsConsistant=TRUE, 
    useMMdisp=FALSE, shrinkVariance=FALSE, pairedDesign=FALSE, 
    pairedDesign.dispMethod="per-pair", useFisher=FALSE, Dispersions=NULL,
    eSlope=NULL, plotASD=FALSE, lwd_ASD=4.5, cex_ASD=1.2)

Arguments

countsA	A counts table under condition “condA".
countsB	A counts table under condition “condB".
sizeFactorsA	A vector of size factors under condition “condA".
sizeFactorsB	A vector of size factors under condition “condB".
numPart	An integer indicating the number of targets for the shrinkage dispersion estimates. “numPart=1" is the default value. It assumes that all the genes share one common target, and then the method of moment estimates are shrunk toward one single target. When it is assumed that the genes share multiple targets, the value for “numPart" is the number of targets and the grouped shrinkage estimates for dispersion are calculated.
SHonly	If ‘SHonly’ is TRUE, then the function outputs the shrinkage estimates for dispersion without testing the differentiation between conditions. If FALSE, then the function outputs a data frame including the per-gene p-values of tests.
propForSigma	A range vector between 0 and 1 that is used to select a subset of data. It helps users to make a flexible choice on the subset of data when they believe only part of data should be used to estimate the variation among per-gene dispersion. An input “propForSigma=c(0.1, 0.9)" means that the genes having method of moment estimates for dispersion greater than the 10th quantile and less than the 90th quantile are used to estimate the dispersion variation. The default input “propForSigma=c(0, 1)" is recommended. It means that we want to use all the data to estimate the dispersion variation.
shrinkTarget	A value for the shrinkage target of dispersion estimates. If “shrinkTarget=NULL" and “shrinkQuantile" is a value instead of NULL, then the quantile value for “shrinkQuantile" is converted into the scale of dispersion estimates and used as the target. If both of them are NULL, then a value that is small and minimizes the average squared difference is automatically used as the target value. If both of them are not NULL, then the value of “shrinkTarget" is used as the target.
shrinkQuantile	A quantile value for the shrinkage target of dispersion estimates. If “shrinkTarget=NULL" and “shrinkQuantile" is a value instead of NULL, then the quantile value for “shrinkQuantile" is converted into the scale of dispersion estimates and used as the target. If both of them are NULL, then a value that is small and minimizes the average squared difference is automatically used as the target value. If both of them are not NULL, then the value of “shrinkTarget" is used as the target.
cLA	A data.frame indicating the levels or extra factors under condition “condA".
cLB	A data.frame indicating the levels or extra factors under condition “condB".
contrast	A contrast vector for testing in complex experiments. The length of this vector equals to the number of columns in countsTable.
keepLevelsConsistant	A logic TRUE/FALSE value. When “coLevels" is used to indicate a paired design experiment, “keepLevelsConsistant=TRUE" silences the genes that have different changing directions (i.e. positive and negative test statistics) among individual samples by setting their p-values as 1.
useMMdisp	A logic value. When “useMMdisp=TRUE" the method of moment (MM) estimates for dispersion without any shrinkage approach are used for testing the differentiation of genes between two conditions.
shrinkVariance	A logic value. When “shrinkVariance=TRUE", the testing is based on the shrinkage estimates for variance instead of dispersion.
pairedDesign	A logic value. When pairedDesign=TRUE is specified, the tests are performed specifically for the paired design experiment. The Null hypotheses ∑_l(μ_{gA,l}-μ_{gB,l})=0 will be tested.
pairedDesign.dispMethod	A character specifying the method of selecting data used for the paired design experiment. When the input is “per-pair" (the default input), the dispersion estimates are shrunk within each pair of samples. The shrinkage target is different in different pair of samples. When the input is "pooled", firstly method of moment estimates for dispersion are obtained within each pair of samples, and then the average estimates across all pairs of samples are shrunk toward a common targets among genes.
useFisher	A logic value specifying whether Fisher's method of combining multiple p-values for a gene is used in the paired design experiment. In detail the formula of calculating the Fisher's combined p-values is pval_g=χ^2_{df=2k}(X>x) where k is the number of pairs and x=-2*∑^k_{l=1}log_{e}(p_{l}). The default input is FALSE and the formulae pval_{g}=exp(∑^k_{l=1}log_{e}(p_{l})) is used.
Dispersions	If it is not null, then the input is a vector of known dispersion values. The length of the vector equals to the number of genes in the counts table. The default value is “NULL".
eSlope	A positive value near to zero. When selecting the shrinkage target that is small and minimizing the average squared difference (ASD), the value of “elope" is a threshold to stop the selection steps if the absolute value of a local slope for the ASD is less than the threshold. The default value is 0.05.
plotASD	A logic value. If plotASD=TRUE, then the plot of average squared difference (ASD) versus target points is produced. The shrinkage (SH) estimates are obtained by shrinking the method of moment (MM) estimates toward a target. In the figure, the vertical axis are ASD values when the shrinkage target (represented by the horizontal axis) varies within the range of dispersion estimates. The selected target is a small value minimizing ASD.
lwd_ASD	A value specifying the width of the curve shown in the plot for the average squared difference when “plotASD=TRUE". The default value is 4.5.
cex_ASD	A value specifying the size of label text shown in the plot for the average squared difference when “plotASD=TRUE". The default value is 1.2.

See Also

nbTestSH.

sSeq documentation built on Nov. 8, 2020, 5:52 p.m.