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

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)
#> [,1] #> [1,] -4.197374