Monte Carlo based contingency table tests considering the compositional approach to contingency tables.
rSDev.test(x, R = 999, method = "multinom")A list with class “htest” containing the following components:
statisticthe value of the relative simplicial deviance (test statistic).
methoda character string indicating what type of rSDev.test was performed.
p.valuethe p-value for the test.
matrix, data.frame or table
an integer specifying the number of replicates used in the Monte Carlo test.
either “rmultinom” (default) or “permutation”.
Matthias Templ, Karel Hron
Method “rmultinom” generate multinomially distributed samples
from the independent probability table, which is estimated from x using geometric mean marginals.
The relative simplicial deviance of the original data are then compared to the generated ones.
Method “permutation” permutes the entries of x and compares the relative simplicial deviance estimated from
the original data to the ones of the permuted data (the independence table is unchanged and originates on x).
Method “rmultinom” should be preferred, while method “permutation” can be used for comparisons.
Egozcue, J.J., Pawlowsky-Glahn, V., Templ, M., Hron, K. (2015) Independence in contingency tables using simplicial geometry. Communications in Statistics - Theory and Methods, 44 (18), 3978--3996.
rSDev
data(precipitation)
rSDev.test(precipitation)
Run the code above in your browser using DataLab