MCtest: Monte Carlo permutation test

Description

This function performs Monte Carlo permutation test to asses the statistical significance of R and Q statistics.

Usage

MCtest(RAN, RQ, nper = 100)

Arguments

RAN

matrix of ranks to be permuted, columns refer to studies, rows refer to genes

observed values of R (average rank) and Q (heterogenity) - as produced by compute.RQ

nper

number of permutations

Value

Returns a matrix with four columns. First (Second) column represents significance level of high (low) avarage rank. Third (fourth) represents significance level of high (low) heterogenity.

References

Zintzaras, E., Ioannidis, J.P.A 2008 Meta-analysis for ranked discovery datasets: Theoretical framework and empirical demonstration for microarrays, Computational Biology and Chemistry 32, 39-47

Examples

Run this code

RANK<-cbind(sample(100), sample(100), sample(100))
RQ<-compute.RQ(RANK)
head(RQ)       
MCtest(RANK, RQ, nper=100)