check.dist: Check for Distribution

In pmartinezarbizu/dada2pp: Miscelaneous Functions for Post-Processing dada2 Results.

Check for Distribution

Description

Use MASS fitdistr function to find best fitting distribution

Usage

check.dist(x)

Arguments

x	vector of values

Details

If x contains integer values, fitdistr is used to check for "negative binomial" or "Poisson" distribution. If x is float following distributions are tested "lognormal", "normal", "exponential","geometric",'gamma'. Solution are ranked according to AIC information criterium. Lowest is best.

Value

Ranked solutions

Author(s)

Pedro Martinez Arbizu

See Also

fitdistr

Examples

check.dist(rnorm(30,20,1)) 

# negative binomial distribution from a collection of poisson samples 
d <- c()
for(i in 1:20){d=c(d,rpois(50,sample(10:1000,1)))}
check.dist(d)

pmartinezarbizu/dada2pp documentation built on Feb. 7, 2024, 7:01 a.m.