# #——风格,回声= FALSE,结果= '飞机 '------------------------------- BiocStyle:减价()suppressPackageStartupMessages(库(grasp2db )) ## ---- loadup2 ------------------------------------------------------------- 库(grasp2db) # #——grasp2 -------------------------------------------------------------- grasp2 < - grasp2 () grasp2 # #——lkp ----------------------------------------------------------------- 变异< -台(grasp2“变种”)q1 % > % =(变体选择(Pvalue,NegativeLog10PBin) % > %过滤器(NegativeLog10PBin > 8) % > %总结(maxp = max (Pvalue), n = n ())) q1 # #——lkex ---------------------------------------------------------------- 解释(q1) # #——lkp2 ---------------------------------------------------------------- 研究< -台(grasp2,“研究”)large_effect < -变体% > %选择(PMID, SNPid_dbSNP134 NegativeLog10PBin) % > %过滤器(NegativeLog10PBin > 5)表型< - left_join (large_effect,研究% > %选择(PMID,PaperPhenotypeDescription)表型# #——doasth -------------------------------------------------------------- lkaw < - semi_join(变体% > %过滤器(NegativeLog10PBin < = 4) % > %选择(PMID、chr_hg19 SNPid_dbSNP134, PolyPhen2),研究% > %过滤器(PaperPhenotypeDescription = =“哮喘 ")) ## ---- dogre --------------------------------------------------------------- lkaw % > %过滤器(" % amaging PolyPhen2 %和% %") ## ---- lkbasic ------------------------------------------------------------- grasp2 # #——lkcon --------------------------------------------------------------- gcon = grasp2反对美元库(RSQLite) gcon dbListTables (gcon) # #——getgw --------------------------------------------------------------- 库(gwascat)Data (gwrngs19) # hg19地址;NHGRI船只hg38 gwrngs19 # #——lkanti -------------------------------------------------------------- gr22 % > % =变体过滤器(chr_hg19 = =“22”)abs22 = checkAnti (22) 1 - (abs22 % > % nrow ()) / (gr22 % > % % > %收集计数 %>% `[[`(" n ")) ## ---- lkabs --------------------------------------------------------------- abs22