cp.fun.chisq.test {FunChisq} | R Documentation |
Comparative functional chi-squared tests on two or more contingency tables.
cp.fun.chisq.test( x, method = c("fchisq", "nfchisq", "default", "normalized"), log.p = FALSE )
x |
a list of at least two matrices representing contingency tables of the same dimensionality. |
method |
a character string to specify the method to compute the functional chi-squared statistic and its p-value. The default is Note: |
log.p |
logical; if |
The comparative functional chi-squared test determines whether the patterns underlying the contingency tables are heterogeneous in a functional way. Specifically, it evaluates whether the column variable is a changed function of the row variable across the contingency tables.
Two methods are provided to compute the functional chi-squared statistic and its p-value. When method = "fchisq"
(or "default"
), the p-value is computed using the chi-squared distribution; when method =
"nfchisq"
(or "normalized"
) a normalized statistic is obtained by shifting and scaling the original statistic and a p-value is computed using the standard normal distribution (Box et al., 2005). The normalized test is more conservative on the degrees of freedom.
A list with class "htest
" containing the following components:
statistic |
functional heterogeneity statistic if |
parameter |
degrees of freedom. |
p.value |
p-value of the comparative functional chi-squared test. By default, it is computed by the chi-squared distribution. If |
Yang Zhang and Joe Song
Box, G. E., Hunter, J. S. and Hunter, W. G. (2005) Statistics for Experimenters: Design, Innovation and Discovery, 2nd Ed., New York: Wiley-Interscience.
Zhang, Y. (2014) Nonparametric Statistical Methods for Biological Network Inference. Unpublished doctoral dissertation, Department of Computer Science, New Mexico State University, Las Cruces, USA.
Zhang, Y. and Song, M. (2013) Deciphering interactions in causal networks without parametric assumptions. arXiv Molecular Networks, arXiv:1311.2707. https://arxiv.org/abs/1311.2707
For comparative chi-squared test that does not consider functional dependencies, cp.chisq.test
.
## Not run: x <- matrix(c(4,0,4,0,4,0,1,0,1), 3) y <- t(x) z <- matrix(c(1,0,1,4,0,4,0,4,0), 3) data <- list(x,y,z) cp.fun.chisq.test(data) cp.fun.chisq.test(data, method="nfchisq") ## End(Not run)