Xsc.onesample: Generalized Wald-Type Statistics: One Sample RAD...

In HMP: Hypothesis Testing and Power Calculations for Comparing Metagenomic Samples from HMP

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

Example output

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
$`Xsc statistics`
[1] 166.2704

$`p value`
[1] 0

HMP documentation built on Aug. 31, 2019, 5:05 p.m.