Talk:Fdrsim.R

source("fdrsim.R") source("plotfun.R")

nseq <- 201 fdrsim("pi0-1000-sam", Estpi0 = TRUE, iiter=1000, seqpi=seq(0,1, length=nseq)) fdrsim("pi0-10000-sam", Estpi0 = TRUE, iiter=10000, seqpi=seq(0,1, length=nseq)) fdrsim("tmp2-1000-mwd", Estpi0 = TRUE, iiter=1000, use.siggenes=F, seqpi=seq(0,1, length=nseq)) fdrsim("tmp2-10000-mwd", Estpi0 = TRUE, iiter=10000, use.siggenes=F, seqpi=seq(0,1, length=nseq))

load(file.path("Data", "pi0-1000-sam")) summary(FDlist$p01) plotfun(FDlist, what="p0", binNo=100) plotfun(FDlist, what="p0NW", binNo=100) X11() load(file.path("Data", "pi0-10000-sam")) summary(FDlist$p01) plotfun(FDlist, what="p0", binNo=100)

X11() load(file.path("Data", "tmp2-1000-mwd")) summary(FDlist$p01) plotfun(FDlist, what="p0")

fdrsim("tmp", m=100, Estpi0 = FALSE, iiter=100, seqpi=seq(0,1, length=100)) load(file.path("Data", "tmp")) plotfun(FDlist)

load(file.path("Data", "tmp2-10000-mwd"))