intern.score {globalSeq}R Documentation

Internal function

Description

This function calculates the test statistic. It is called by the function omnibus.

Usage

intern.score(y, R, mu, phi)

Arguments

y

response variable: numeric vector of length n

R

numeric matrix of dimensions n*n (see example)

mu

mean parameters: numeric vector of length 1 or n

phi

dispersion parameter: non-negative real number

Value

The function returns a real number.

References

A Rauschenberger, MA Jonker, MA van de Wiel, and RX Menezes (2016). "Testing for association between RNA-Seq and high-dimensional data", BMC Bioinformatics. 17:118. html pdf (open access)

See Also

This is an internal function. The user functions are cursus, omnibus, and proprius.

Examples

# simulate high-dimensional data
n <- 30
p <- 100
set.seed(1)
y <- rnbinom(n,mu=10,size=1/0.25)
X <- matrix(rnorm(n*p),nrow=n,ncol=p)

# calculate test statistic
R <- X %*% t(X) / ncol(X)
mu <- mean(y)
phi <- (var(y)-mu)/mu^2
intern.score(y,R,mu,phi)


[Package globalSeq version 1.16.1 Index]