HMP (version 2.0.1)

Xsc.onesample: 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

# NOT RUN {
	data(saliva)
	data(throat)
	
	### Get pi from the dirichlet-multinomial parameters
	pi0 <- dirmult(saliva)$pi
	
	xsc <- Xsc.onesample(throat, pi0)
	xsc
# }