genesimulator: Gene means simulator
In DDHFm: Variance Stabilization by Data-Driven Haar-Fisz (for Microarrays)

genesimulator

R Documentation

Gene means simulator

Description

Simulates means of the gene intensities for "nps" genes, each replicated "nreps" times. The mean is drawn from a gamma distribution with shape parameter "shape" and scale parameter "scale"

Usage

genesimulator(nreps = 3, nps = 100, shape = 4, scale = 100)

Arguments

`nreps`	Number of replicates
`nps`	Number of genes
`shape`	Shape parameter
`scale`	Scale parameter

Details

For many problems a set of reasonable gene mean intensities is useful for testing algorithms that later draw actual gene intensities with distributions that possess a gene mean intensity.

Value

A matrix containing nreps*nps rows and 3 columns. The first col contains all the gene mean intensities. The second and third col contain the gene replicate number and gene number respectively.