TestStatSP: Calculation of test statistic

In jkruppa/RepeatedHighDim: Global tests for expression data of high-dimensional sets of molecular features.

Description

Calculation of test statistic

Usage

TestStatSP(Y1, Y2)

Arguments

Y1	Matrix of expression levels in first group. Rows represent features (e.g. genes, proteins,...), columns represent samples.
Y2	Matrix of expression levels in second group. Rows represent features (e.g. genes, proteins,...), columns represent samples.

Details

Calculates the test statistic for RepeatedHighDim in the case of unpaired samples.

Value

A list containing the following items:

k

Indicates whether the paired or unpaired case was tested.

d

Number of features.

n1

Number of samples in group 1.

n2

Number of samples in group 2.

Fn

Test statistic.

f

First degree of freedoms.

f2

Second degree of freedom.

p

p-value.

Author(s)

Klaus Jung

References

• Brunner, E (2009) Repeated measures under non-sphericity. Proceedings of the 6th St. Petersburg Workshop on Simulation, 605-609.

• Jung K, Becker B, Brunner B and Beissbarth T (2011) Comparison of Global Tests for Functional Gene Sets in Two-Group Designs and Selection of Potentially Effect-causing Genes. Bioinformatics, 27: 1377-1383.

Examples

###### Global comparison of a set of 100 genes between two experimental groups.
X1 = matrix(rnorm(1000, 0, 1), 10, 100)
X2 = matrix(rnorm(1000, 0.1, 1), 10, 100)
RHD = RepeatedHighDim(X1, X2, paired=FALSE)
summary(RHD)

jkruppa/RepeatedHighDim documentation built on May 19, 2019, 12:45 p.m.