An example gene intensity simulator with given mean values mu0
which then returns the variance, skew and kurtosis of the
DDHF transformed intensities. The gene intensity
simulation is performed using simdurbin2
dhhrcomp(nsims = 1024, nmu = 4,
mu0 = c(0, 5, 10, 15, 20, 25, 30, 40, 50, 60,
65, 70, 80, 100, 200, 300, 500, 1000) * 1000)
A vector of length nmu of the mean intensities considered
A vector of nmu variances, one for each mean in
mu0.
The variance is the variance of the DDHFm transformed simulated
intensity data.
A vector of nmu skewnesses, one for each mean in
mu0.
The skewness is the skewness of the DDHFm transformed simulated
intensity data.
A vector of nmu kurtoses, one for each mean in
mu0.
The kurtosis is the kurtosis of the DDHFm transformed simulated
intensity data.
Number of replicates for each mu
Selects first nmu genes from mu for use
The possible gene intensity means
Guy Nason <g.p.nason@bris.ac.uk.
This function simulates some gene intensities. A list of
possible intensity means are supplied in mu0.
The first nmu of these are selected. Then for each
of the selected means nsims gene intensities are
generated.
The intensities are them subjected to DDHF transformation.
Durbin, B.P., Hardin, J.S., Hawkins, D.M. and Rocke, D.M. (2002), A variance-stabilizing transformation for gene expression microarray data, Bioinformatics, 18, S105-S110
simdurbin2