Description Usage Arguments Details Value Author(s) References See Also Examples
Generate plots for the results of simuComp
. The power or FPR of limma voom, limma trans, integrated likelihood and Quasi likelihood method are compared.
1 |
intcomp |
An object of class "intcomp" produced by |
ylim |
The ylim of the plot. Should be a numeric vector of length 2. |
text |
The text shown on the plot. |
pos |
The position of the text on the plot. Should be a numeric vector of length 2. |
... |
Additional arguments passed to plot. |
If the null is set to be true when simuComp is called, then the p-value thresholds are plotted against the empirical FPR, which is the ratio of rejections. Otherwise the empirical discover rate (average number of discoveries divided by the number of genes) is plotted against p-value thresholds. Four plots are shown corresponds to the small, medium, large difference genes and all genes.
A plot is created on the current graphics device.
Yilun Zhang, David Rocke
Our paper
1 2 3 4 5 6 7 8 9 10 11 12 13 | ## Not run:
data(count.data)
data(condition)
count=count.data[,1:10]
cond=rep(0:1,each=5)
res=intSEQ(count, cond,nneighbour=400,lambda1=10,lambda2=0.05,meanmeth="local.mean"
,"loess",normalize=TRUE, coeff =c(2))
nsimu <- simuComp(res, ntime = 2, null = TRUE)
asimu <- simuComp(res, ntime = 2, null = FALSE)
plotComp(nsimu, text = " 3 vs 3 null simulation")
plotComp(asimu, text = " 3 vs 3 power simulation")
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.