Xsc.onesample {HMP} | R Documentation |
Generalized Wald-Type Statistics: One Sample RAD Probability-Mean Test Comparison
Description
This routine provides the value of the Generalized Wald-type statistic to assess whether the RAD probability-mean observed in one group of samples is equal to a known RAD probability-mean.
Usage
Xsc.onesample(data, pi0)
Arguments
data |
A matrix of taxonomic counts(columns) for each sample(rows). |
pi0 |
The RAD-probability mean vector. |
Value
A list containing Generalized Wald-type statistics and p-value.
Examples
data(saliva)
data(throat)
### Get pi from the dirichlet-multinomial parameters
pi0 <- dirmult(saliva)$pi
xsc <- Xsc.onesample(throat, pi0)
xsc
[Package HMP version 2.0.1 Index]