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
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.
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”) |
/project/mstephens/test_rss/data/load2013/load2013_sumstat.mat
/scratch/PI/whwong/zduren/share/PECA_human/PECA2/Liver_network.txt
h=0.3; piva=10.^(-(0.25:0.25:8))
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 |
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