Compute score test statistics.
compute_stat(Y, K_int, y_fixed, K0, sigma2_hat, tau_hat)
(matrix, n*1) The vector of response variable.
(matrix, n*n) The kernel matrix to be tested.
(vector of length n) Estimated fixed effect of the response.
(numeric) The estimated noise of the fixed effect.
(numeric) The estimated noise of the kernel effect.
(matrix, n*n) Estimated ensemble kernel matrix.
(numeric) The computed test statistic.
The test statistic is distributed as a scaled Chi-squared distribution.
Arnab Maity and Xihong Lin. Powerful tests for detecting a gene effect in the presence of possible gene-gene interactions using garrote kernel machines. December 2011.