Last updated: 2018-08-31

workflowr checks: (Click a bullet for more information)
  • R Markdown file: up-to-date

    Great! Since the R Markdown file has been committed to the Git repository, you know the exact version of the code that produced these results.

  • Environment: empty

    Great job! The global environment was empty. Objects defined in the global environment can affect the analysis in your R Markdown file in unknown ways. For reproduciblity it’s best to always run the code in an empty environment.

  • Seed: set.seed(20180719)

    The command set.seed(20180719) was run prior to running the code in the R Markdown file. Setting a seed ensures that any results that rely on randomness, e.g. subsampling or permutations, are reproducible.

  • Session information: recorded

    Great job! Recording the operating system, R version, and package versions is critical for reproducibility.

  • Repository version: 4af3002

    Great! You are using Git for version control. Tracking code development and connecting the code version to the results is critical for reproducibility. The version displayed above was the version of the Git repository at the time these results were generated.

    Note that you need to be careful to ensure that all relevant files for the analysis have been committed to Git prior to generating the results (you can use wflow_publish or wflow_git_commit). workflowr only checks the R Markdown file, but you know if there are other scripts or data files that it depends on. Below is the status of the Git repository when the results were generated:
    
    Ignored files:
        Ignored:    .Rhistory
        Ignored:    .Rproj.user/
    
    
    Note that any generated files, e.g. HTML, png, CSS, etc., are not included in this status report because it is ok for generated content to have uncommitted changes.
Expand here to see past versions:
    File Version Author Date Message
    Rmd 4af3002 Xiang Zhu 2018-08-31 wflow_publish(“analysis/load2013_liver.Rmd”)
    html 5a3c37b Xiang Zhu 2018-08-31 Build site.
    Rmd 86e52f5 Xiang Zhu 2018-08-31 wflow_publish(“analysis/load2013_liver.Rmd”)
    html 1734dec Xiang Zhu 2018-08-30 Build site.
    Rmd 6b3b280 Xiang Zhu 2018-08-30 wflow_publish(“analysis/load2013_liver.Rmd”)
    Rmd dbd04eb Xiang Zhu 2018-08-30 minor update chunk option
    html 8843641 Xiang Zhu 2018-08-30 Build site.
    Rmd e5fe6e0 Xiang Zhu 2018-08-30 wflow_publish(“analysis/load2013_liver.Rmd”)

Input data

  • Processed GWAS file: /project/mstephens/test_rss/data/load2013/load2013_sumstat.mat
  • PECA2 liver-specific network: /scratch/PI/whwong/zduren/share/PECA_human/PECA2/Liver_network.txt
  • \(L_0\): 40 kb; \(L_1\): 100 kb
  • SNP-to-gene window: 10 Mb
  • SNP-to-RE window: exact zero
  • Hyperparameter grid: h=0.3; piva=10.^(-(0.25:0.25:8))

Results summary

First look at the approximated log marginal likelihoods (elbo column below).

log10.piva sigb elbo time
-7.75 3.8688047 3.983163e+15 67.84746
-8.00 5.1591340 3.116692e+11 68.16840
-7.25 2.1755888 3.978008e+10 58.15310
-7.50 2.9011943 3.797890e+10 45.59751
-7.00 1.6314614 1.465382e+10 61.22707
-6.75 1.2234235 4.938467e+09 24.58350
-5.75 0.3868805 1.478323e+09 112.50603
-6.50 0.9174382 1.470216e+09 26.28262
-5.50 0.2901194 1.068584e+09 1111.35141
-5.25 0.2175589 8.074074e+08 94.58703
-6.25 0.6879816 5.243371e+08 1825.34849
-5.00 0.1631461 4.497749e+08 64.31936
-6.00 0.5159134 4.489631e+08 25.41265
-4.75 0.1223423 2.190502e+08 1027.00462
-4.00 0.0515913 1.895669e+08 1047.94328
-4.25 0.0687982 1.065638e+08 1664.43326
-3.75 0.0386880 8.842670e+07 957.18869
-4.50 0.0917438 6.921695e+07 1258.99004
-3.50 0.0290119 6.456807e+07 1302.29656
-3.25 0.0217559 3.595338e+07 1276.26215
-3.00 0.0163146 1.699945e+07 1443.66553
-2.75 0.0122342 6.239009e+06 1536.16943
-2.50 0.0091744 1.241460e+06 1233.59871
-1.00 0.0016315 -6.400220e+05 4096.34627
-0.50 0.0009174 -7.316516e+05 3479.09999
-0.75 0.0012234 -7.666871e+05 3141.62185
-0.25 0.0006880 -8.983778e+05 2900.37023
-2.25 0.0068798 -9.468763e+05 3797.00448
-1.25 0.0021756 -9.971599e+05 3615.59358
-1.50 0.0029012 -1.100972e+06 4003.91446
-2.00 0.0051591 -1.141247e+06 4012.43595
-1.75 0.0038688 -1.435637e+06 3133.88371

Next look at the gene-level posterior statistics when the marginal likelihood is maximized.

hgnc_symbol chromosome_name start_position end_position gene_nid vb_weight vb_mean vb_var
STH 17 44076616 44077060 15105 1 82230388 0.0907612
ARL17B 17 44352150 44439130 15107 1 75155535 0.4658901
LRRC37A 17 44370099 44415160 15108 1 70458266 0.7754152
POTEH 22 16256441 16287937 17913 1 62692235 14.9676118
FAM27C 9 44990314 44991483 8578 1 9903777 0.0945179
FOXD4L5 9 70175707 70178815 8590 1 4213858 14.9676299
NBPF11 1 146032647 146082765 1073 1 4206964 0.0505033
FAM27A 9 45727107 45728274 8579 1 3620112 14.9676502
HLA-A 6 29909037 29913661 6055 1 3154370 5.0805388
FOXD4L6 9 69199480 69202204 8587 1 2762005 14.9676502

Session information

R version 3.5.1 (2018-07-02)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: macOS High Sierra 10.13.6

Matrix products: default
BLAS: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRblas.0.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/3.5/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     

other attached packages:
[1] bindrcpp_0.2.2 DT_0.4         knitr_1.20     dplyr_0.7.5   
[5] R.matlab_3.6.1

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.17      bindr_0.1.1       whisker_0.3-2    
 [4] magrittr_1.5      workflowr_1.1.1   tidyselect_0.2.4 
 [7] R6_2.2.2          rlang_0.2.1       highr_0.7        
[10] stringr_1.3.1     tools_3.5.1       R.oo_1.22.0      
[13] git2r_0.21.0      htmltools_0.3.6   yaml_2.1.19      
[16] rprojroot_1.3-2   digest_0.6.15     assertthat_0.2.0 
[19] tibble_1.4.2      purrr_0.2.5       htmlwidgets_1.2  
[22] R.utils_2.6.0     glue_1.2.0        evaluate_0.10.1  
[25] rmarkdown_1.10    stringi_1.2.3     pillar_1.2.3     
[28] compiler_3.5.1    backports_1.1.2   R.methodsS3_1.7.1
[31] pkgconfig_2.0.1  

This reproducible R Markdown analysis was created with workflowr 1.1.1