This function computes the Generalized Wald-type test statistic (Wilson and Koehler, 1984) and corresponding p-value to assess whether the sample RAD probability-means from multiple populations are the same or different. The statistics assumes that a common RAD probability-mean vector for comparison under the null hypothesis is known.
A list where each element is a matrix of taxonomic counts(columns) for each sample(rows).
The RAD-probability mean vector.
Note: The matrices in
group.data must contain the same taxa, in the same order.
A list containing the Generalized Wald-type statistics and p-value.
Wilson, J. R., and Koehler, K. J. (1984). Testing of equality of vectors of proportions for several cluster samples. Proceedings of Joint Statistical Association Meetings. Survey Research Methods.
1 2 3 4 5 6 7 8 9 10 11 12
Loading required package: dirmult Attaching package: 'HMP' The following object is masked from 'package:dirmult': weirMoM Iteration 1: Log-likelihood value: -1426219.55743915 Iteration 2: Log-likelihood value: -1426159.69392069 Iteration 3: Log-likelihood value: -1426137.57581644 Iteration 4: Log-likelihood value: -1426134.36208055 Iteration 5: Log-likelihood value: -1426134.28420196 Iteration 6: Log-likelihood value: -1426134.28414898 $`Xmc statistics`  295.529 $`p value`  0
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.