aldex.ttest {ALDEx2}R Documentation

calculate Welch's t-test and Wilcoxon test statistics

Description

calculates expected values of the Welch's t-test and Wilcoxon rank test on the data returned by clr_function.r

Usage

aldex.ttest(clr, conditions, paired.test = FALSE, hist.plot=FALSE)

Arguments

clr

clr is the data output of the aldex.clr function

conditions

a description of the data structure to be used for testing

paired.test

whether the Welch's test should be paired or not

hist.plot

whether to plot a histogram of P values for an individual Dirichlet Monte-Carlo instance. Plot is output to the standard R plotting device.

Details

An explicit example for two conditions is shown in the ‘Examples’ below.

Value

Outputs a dataframe with the following information:

we.ep

a vector containing the expected P value of the Welch's t-test for each feature

we.eBH

a vector containing the expected value of the Benjamini Hochberg corrected P value for each feature

wi.ep

a vector containing the expected P value of the Wilcoxon test for each feature

wi.eBH

a vector containing the expected value of the Benjamini Hochberg corrected P value for each feature

Author(s)

Greg Gloor

References

Please use the citation given by citation(package="ALDEx").

See Also

aldex.clr, aldex.glm, aldex.effect, selex

Examples

    # x is the output of the \code{x <- aldex.clr(data, mc.samples)} function
    # conditions is a description of the data
    # for the selex dataset, conditions <- c(rep("N", 7), rep("S", 7))
    data(selex)
    #subset for efficiency
    selex <- selex[1201:1600,]
    conds <- c(rep("NS", 7), rep("S", 7))
    x <- aldex.clr(selex, conds, mc.samples=2, denom="all")
    ttest.test <- aldex.ttest(x, conds)


[Package ALDEx2 version 1.12.0 Index]