-
Notifications
You must be signed in to change notification settings - Fork 4
/
Fine_Mapping_PD_results.Rmd
104 lines (84 loc) · 3.04 KB
/
Fine_Mapping_PD_results.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
---
title: "<center><h1>Fine Mapping:</h1>Parkinson's Disease</h1></center>"
# author:
# "<div class='container'>
# <h3>Brian M. Schilder, Bioinformatician II<br>
# Raj Lab<br>
# Department of Neuroscience<br>
# Icahn School of Medicine at Mount Sinai<br>
# NYC, New York<br>
# </h3>
# <a href='https://github.com/RajLabMSSM/Fine_Mapping' target='_blank'><img src='./echolocatoR/images/echo_logo_sm.png'></a>
# <a href='https://github.com/RajLabMSSM' target='_blank'><img src='./web/images/github.png'></a>
# <a class='item' href='https://rajlabmssm.github.io/RajLab_website/' target='_blank'>
# <img src='./web/images/brain-icon.png'>
# <span class='caption'>RAJ LAB</span>
# <a href='https://icahn.mssm.edu/' target='_blank'><img src='./web/images/sinai.png'></a>
# </div>"
# date: "<br>Most Recent Update:<br> `r Sys.Date()`"
# output: rmarkdown::html_vignette
output:
html_document:
theme: cerulean
highlight: zenburn
code_folding: show
toc: false
toc_float: true
smooth_scroll: true
number_sections: false
self_contained: true
css: ./web/css/style.css
editor_options:
chunk_output_type: inline
---
```{r setup, message=F, warning=F, class.output = "pre"}
# knitr::opts_chunk$set(error = T) # Knit even with errors
# source("echolocatoR/R/MAIN.R")
library(dplyr)
library(stringr); library(bsselectR) # devtools::install_github("walkerke/bsselectR")
```
# Nalls et. al. (2019) w/ 23andMe
Summarize Fine-mapping Results
## Show all plots
```{r Gather plots, results='asis'}
# .tabset-dropdown
loci <- list.dirs("./Data/GWAS/Nalls23andMe_2019", recursive = F) %>% basename()
loci <- loci[loci!="_genome_wide"]
plot.paths <- file.path("./Data/GWAS/Nalls23andMe_2019",
loci, "Multi-finemap",
paste0(loci,"_ggbio.png"))
# Add link to plot
# baseURL <- "https://github.com/RajLabMSSM/Fine_Mapping/blob/master/Data/GWAS/Nalls23andMe_2019"
# links <- data.frame(Link=file.path(baseURL, top_loci_sort$Locus,"Multi-finemap",
# paste0(top_loci_sort$Locus,"_ggbio.png")))
# top_loci_sort <- cbind(links, top_loci_sort)
# file.path("https://rajlabmssm.github.io/Fine_Mapping/Data/Nalls23andMe_2019",
# genes,"Multi-finemap", paste0(genes,"_ggbio.png"))
# Dropdown
# genes <- basename(dirname(dirname(plot.paths)))
# plot.paths <- file.path("https://rajlabmssm.github.io/Fine_Mapping/Data/GWAS/Nalls23andMe_2019",genes,"Multi-finemap",paste0(genes,"_ggbio.png"))
# names(plot.paths) <- genes
#
#
# bsselect(vector = plot.paths,
# type = "img",
# selected = "LRRK2",
# frame_height = "1000",
# live_search = T,
# show_tick = T)
for (p in plot.paths){
if(file.exists(p)){
try({
gene <- basename(dirname(dirname(p)))
cat(' \n###', gene, ' \n')
cat(paste0('![',gene,'](',p,')'))
cat(' \n')
})
}
}
```
# Session Information
```{r Session Information}
sessionInfo()
print(paste("susieR ", packageVersion("susieR")))
```