GPseq (version 0.5)

perform_permutation_de: Perform a Permutation Test to fit the null distribution of the test statistic to a Gamma Distribution

Description

This function calculates the shape and scale parameters for a Gamma Distribution which is used to fit the null distribution of the test statistic calculated for the Differential Expression of two genes or exons

Usage

perform_permutation_de(x,y,num_permute)

Arguments

x

Vector for counts of gene/exon in sample 1

y

Vector for counts of gene/exon in sample 2

num_permute

Number of times the permutation is done

Value

shape

The shape parameter of the Gamma distribution

scale

The scale parameter of the Gamma distribution

References

Consul, P. C. (1989) Generalized Poisson Distributions: Properties and Applications. New York: Marcel Dekker. Sudeep Srivastava, Liang Chen A two-parameter generalized Poisson model to improve the analysis of RNA-Seq data Nucleic Acids Research Advance Access published July 29,2010 doi : 10.1093/nar/gkq670

See Also

generalized_poisson_likelihood

Examples

x = rpois(100,10);
y = rpois(100,5);

gamma = perform_permutation_de(x,y,10);

cat("Shape parameter = ",gamma$shape,"Scale parameter = ",gamma$scale);