Simulates gene intensities from the two components model of Durbin and Rocke.
simdurbin2(mu, alpha, seta, seps)
the means where the gene intensities are simulated from
background mean of gene intensities
standard deviation of high-level gene intensities
standard deviation of low-level gene intensities
simdurbin this function simulates gene
intensities from the Durbin and Rocke two component model (see
Durbin and Rocke, 2001). Unlike
simdurbin this function
generates each intensity with a mean specified by an entry in
length(mu) intensities are generated with
i having mean
Guy Nason <email@example.com>
Rocke, D.M. and Durbin, B.P. (2001), A model for measurement error for gene expression arrays, Journal of Computational Biology, 8, 557-569
# # Generate data with mu equal to 80 and mu equal to 80000 with 10 replicates each # # v <- simdurbin2( c(rep(80, 10), rep(80000,10)), alpha=24800, seta=0.227, seps=4800) # # Let's look at the mean of the first 10 and the last 10 # mean(v[1:10]) # 27361.95 # mean(v[11:20]) # 92455.66
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.