Description Usage Arguments Value Author(s) Examples
This function automatically estimates weights
for the "weight"
and "dperm.weight"
options of metaseqR for combining pvalues from
multiple statistical tests. It creates simulated
dataset based on real data and then performs
statistical analysis with metaseqR several times
in order to derive False Discovery Curves. Then,
the average areas under the false discovery curves
are used to construct weights for each algorithm,
according to its performance when using simulated
data.
1 2 3 4 5 6 
counts 
the real raw counts table from which the simulation parameters will be estimated. It must not be normalized and must contain only integer counts, without any other annotation elements and unique gene identifiers as the rownames attribute. 
normalization 
same as 
statistics 
same as 
nsim 
the number of simulations to perform to estimate the weights. It default to 10. 
N 
the number of genes to produce.
See 
samples 
a vector with 2 integers, which are the number of samples for each condition (two conditions currently supported). 
ndeg 
a vector with 2 integers, which are the number of differentially expressed genes to be produced. The first element is the number of upregulated genes while the second is the number of downregulated genes. 
fc.basis 
the minimum foldchange for deregulation. 
top 
the top 
model.org 
the organism from which the
data are derived. It must be one of

seed 
a list of seed for reproducible
simulations. Defaults to 
draw.fpc 
draw the averaged false
discovery curves? Default to 
multic 
whether to run in parallel
(if package 
... 
Further arguments to be passed to

A vector of weights to be used in
metaseqr
with the
weights
option.
Panagiotis Moulos
1 2 3 4 5 6 7 8 9  data("mm9.gene.data",package="metaseqR")
multic < check.parallel(0.8)
weights < estimate.aufc.weights(
counts=as.matrix(mm9.gene.counts[,9:12]),
normalization="edaseq",
statistics=c("deseq","edger"),
nsim=3,N=100,ndeg=c(10,10),top=10,model.org="mm9",
seed=10,multic=multic,libsize.gt=1e+5
)

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.