interactiontable: Returns a list of interactions with associated statistics.

Description Usage Arguments Value Warning Author(s) See Also Examples

View source: R/interactiontable.R

Description

This is a extanded wrapper around the topTable function from the limma package, as an option the ordinary t statistics can be calculated as well.

Usage

1
interactiontable(ebfit, sort = "none", ord.t = FALSE, correction = "BH")

Arguments

ebfit

ebfit a MArrayLM object produced by the eBayes function

sort

character string specifying which statistic to rank genes by, possible arguments are none, ID,size, t,B,adj.P.val,P.Value, and if ord.t = TRUE: ord.t, ord.p and ord.p.adj.

ord.t

Logical, should ordinary t statistics be calculted? Default is FALSE.

correction

method used to adjust the p-values for multiple testing. Default is BH. See p.adjust for the complete list of options.

Value

Returns a dataframe where the rows are the interaction pairs and the columns the statistics:

ID: Interaction pair if

size: the average interaction size

t: the moderated t statistics

P.Value: p-value for the moderated t statistics

adj.P.Val: adjusted p-value

B: the b statistics

if the ord.t=TRUE, the ordinary t statistics (ord.t), with correspnding p-values (ord.p) and adjusted p-values (ord.p.adj)

Warning

usage of the ordinary t statistics is not recommended for data sets with few replicates.

Author(s)

Elin Axelsson

See Also

p.adjust,topTable

Examples

1
2
3
4
5
6
7
8
9
## simulated data
     y <- matrix(rnorm(50*4,sd=1),50,4)
     rownames(y) <- paste("Pair",1:50)
    
     # fit and eBayes
     fit <- lmFit(y)
     fit <- eBayes(fit)
     tt = interactiontable(fit,sort="size")
     head(tt)

coRNAi documentation built on Nov. 17, 2017, 11:14 a.m.