Last updated: 2017-11-07

Code version: db998b3

library(ashr)
library(edgeR)
library(limma)
library(pROC)
source("../code/gdash.R")

Read in data

data = readRDS("../data/liver.rds")
ngene = 1e4
data = top_gene_selection(ngene, data)$data

Global Null

\[ g \equiv 0 \ . \]

SNR = 0

\[ g = 0.6\delta_0 + 0.3N\left(0, \sigma^2\right) + 0.1N\left(0, \left(2.65\sigma\right)^2\right) \ . \]

SNR = 2

\[ g = 0.6\delta_0 + 0.3N\left(0, \sigma^2\right) + 0.1N\left(0, \left(3.58\sigma\right)^2\right) \ . \]

SNR = 5

\[ g = 0.6\delta_0 + 0.3N\left(0, \left(2\sigma\right)^2\right) + 0.1N\left(0, \left(4.43\sigma\right)^2\right) \ . \]

SNR = 10

\[ g = 0.6\delta_0 + 0.3N\left(0, \left(3\sigma\right)^2\right) + 0.1N\left(0, \left(8.54\sigma\right)^2\right) \ . \]

Session information

sessionInfo()
R version 3.4.2 (2017-09-28)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: macOS Sierra 10.12.6

Matrix products: default
BLAS: /Library/Frameworks/R.framework/Versions/3.4/Resources/lib/libRblas.0.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/3.4/Resources/lib/libRlapack.dylib

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

loaded via a namespace (and not attached):
 [1] compiler_3.4.2  backports_1.1.1 magrittr_1.5    rprojroot_1.2  
 [5] tools_3.4.2     htmltools_0.3.6 yaml_2.1.14     Rcpp_0.12.13   
 [9] stringi_1.1.5   rmarkdown_1.6   knitr_1.17      git2r_0.19.0   
[13] stringr_1.2.0   digest_0.6.12   evaluate_0.10.1

This R Markdown site was created with workflowr