Xmcupo.effectsize: Effect Size for Xmcupo Statistic

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

Description

This function computes the Cramer's Phi and Modified Cramer's Phi Criterion for the test statistic Xmcupo.sevsample.

Usage

Xmcupo.effectsize(group.data)

Arguments

group.data

A list where each element is a matrix of taxonomic counts(columns) for each sample(rows).

Details

Note: The matrices in group.data must contain the same taxa, in the same order.

Value

A vector containing the Chi-Squared statistic value, the Cramer's Phi Criterion, and the modified Cramer's Phi Criterion.

Examples

	data(saliva)
	data(throat)

	### Combine the data sets into a single list
	group.data <- list(saliva, throat)
	
	effect <- Xmcupo.effectsize(group.data)
	effect

Example output

Loading required package: dirmult

Attaching package: 'HMP'

The following object is masked from 'package:dirmult':

    weirMoM

        Chi-Squared          Cramer Phi Modified-Cramer Phi 
       98.397385359         0.009882751         0.173844356 

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