R/ZT.statistics.Hnull.Ha.R

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

Defines functions ZT.statistics.Hnull.Ha

Documented in ZT.statistics.Hnull.Ha

ZT.statistics.Hnull.Ha <-
function(Nrs, fit, type){
	if(tolower(type) == "hnull"){
		genData <- Multinomial(Nrs, fit$pi)
	}else{
		genData <- Dirichlet.multinomial(Nrs, fit$gamma)
	}
	
	ZT <- c(Z.statistics(genData), T.statistics(genData))
	
	return(ZT)
}

Try the HMP package in your browser