Fits a linear model with interaction term and tests for significance of the interaction term.

LinearModelInteractionTest(dat, xname, yname, response = NULL)

Arguments

dat

dat is a gene expression matrix (rows genes columns samples)

xname

the name of the first gene

yname

name of the second gene

response

a vector of two class labels, if NULL then split into 50/50. odd number of samples will ignore the middle sample

Value

single value for p-value

Examples

data(STRING) data(SKCM) LinearModelInteractionTest(dat = SKCM,xname = "EIF3C",yname = "EIF5B")
#> [1] 0.5482295