Plot.Pi: Plot the Pi Vector

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

Description

Plots the taxa proportions for every group.

Usage

	Plot.PI(estPi, errorBars = TRUE, logScale = FALSE, 
		main = "PI Vector", ylab = "Fractional Abundance")

Arguments

estPi	The results for either MLE or MOM from the function 'Est.Pi'.
errorBars	A boolean to display the error bars or not.
logScale	A boolean to log the y scale or not.
main	A string to be used as the plots title.
ylab	A string to be used as the plots y-axis title.

Value

A plot of the pi vectors for every group.

Examples

	## Not run: 
		data(saliva)
		data(throat)
		data(tonsils)
		
		### Combine the data sets into a single list
		group.data <- list(saliva, throat, tonsils)
		
		### Get PI using MLE with CI
		mle <- Est.PI(group.data)$MLE
		
		### Plot with Error Bars
		Plot.PI(mle)
		
		### Plot without Error Bars
		Plot.PI(mle, FALSE)
		
		### Plot with Error Bars and scaling
		Plot.PI(mle, TRUE, TRUE)
	
## End(Not run)

Example output

Loading required package: dirmult

Attaching package: ‘HMP’

The following object is masked from ‘package:dirmult’:

    weirMoM

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