免费观看又色又爽又黄的小说免费_美女福利视频国产片_亚洲欧美精品_美国一级大黄大色毛片

R語言畫棒棒糖圖展示snp在基因上的位置是怎樣的

這篇文章將為大家詳細(xì)講解有關(guān)R語言畫棒棒糖圖展示snp在基因上的位置是怎樣的,文章內(nèi)容質(zhì)量較高,因此小編分享給大家做個(gè)參考,希望大家閱讀完這篇文章后對相關(guān)知識(shí)有一定的了解。

創(chuàng)新互聯(lián)服務(wù)項(xiàng)目包括揭陽網(wǎng)站建設(shè)、揭陽網(wǎng)站制作、揭陽網(wǎng)頁制作以及揭陽網(wǎng)絡(luò)營銷策劃等。多年來,我們專注于互聯(lián)網(wǎng)行業(yè),利用自身積累的技術(shù)優(yōu)勢、行業(yè)經(jīng)驗(yàn)、深度合作伙伴關(guān)系等,向廣大中小型企業(yè)、政府機(jī)構(gòu)等提供互聯(lián)網(wǎng)行業(yè)的解決方案,揭陽網(wǎng)站推廣取得了明顯的社會(huì)效益與經(jīng)濟(jì)效益。目前,我們服務(wù)的客戶以成都為中心已經(jīng)輻射到揭陽省份的部分城市,未來相信會(huì)繼續(xù)擴(kuò)大服務(wù)區(qū)域并繼續(xù)獲得客戶的支持與信任!

今天在一個(gè)微信群里看見有人問下面兩幅圖用R語言如何實(shí)現(xiàn)

R語言畫棒棒糖圖展示snp在基因上的位置是怎樣的  
image.png
R語言畫棒棒糖圖展示snp在基因上的位置是怎樣的  
image.png
 第一幅就是散點(diǎn)圖

兩列數(shù)據(jù)x,y用來控制散點(diǎn)的位置,一列變量用來映射顏色,F(xiàn)requency用來映射大小。這個(gè)實(shí)現(xiàn)起來相對比較容易。還有一列用來添加文字標(biāo)簽

下面我們構(gòu)造一份數(shù)據(jù)集用來畫圖

df<-data.frame(pvalue=sample(seq(0,1,by=0.001),20,replace=F),
              Effect_size=sample(seq(-2,2,by=0.05),20,replace = F),
              frequency=sample(seq(0,1,by=0.01),20,replace=F),
              color=sample(c("A","B"),20,replace = T),
              text=LETTERS[1:20])
df              
library(ggplot2)
library(ggrepel)
ggplot(df,aes(x=Effect_size,y=pvalue,size=frequency))+
 geom_point(aes(color=color),alpha=0.5)+
 theme_bw()+
 scale_size_continuous(range = c(5,10))+
 geom_text_repel(aes(label=text),size=5)
 
R語言畫棒棒糖圖展示snp在基因上的位置是怎樣的  
image.png
 第二幅圖如何實(shí)現(xiàn)我自己還真沒有思路。

印象里好像是看見過有公眾號(hào)推文介紹這個(gè)圖。我翻了翻之前轉(zhuǎn)發(fā)過的朋友圈還真找到了。有一個(gè)專門的包 trackViewer

這個(gè)包對應(yīng)的論文是

trackViewer: a Bioconductor package for interactive and integrative visualization of multi-omics data. Nature Methods https://www.nature.com/articles/s41592-019-0430-y

github主頁 https://github.com/jianhong/trackViewer

一份很詳細(xì)的幫助文檔 https://bioconductor.org/packages/release/bioc/vignettes/trackViewer/inst/doc/trackViewer.html

文檔中實(shí)現(xiàn)的圖片有

R語言畫棒棒糖圖展示snp在基因上的位置是怎樣的  
image.png
R語言畫棒棒糖圖展示snp在基因上的位置是怎樣的  
image.png

等等。

今天就不重復(fù)文檔中的每個(gè)例子了。知道有這個(gè)包可以實(shí)現(xiàn)這個(gè)功能就可以了。重復(fù)一個(gè)最基本的例子

library(trackViewer)
SNP <- c(10, 12, 1400, 1402)
sample.gr <- GRanges("chr1", IRanges(SNP, width=1, names=paste0("snp", SNP)))
features <- GRanges("chr1", IRanges(c(1, 501, 1001),
                                   width=c(120, 400, 405),
                                   names=paste0("block", 1:3)))
features$fill <- c("#FF8833", "#51C6E6", "#DFA32D")
sample.gr$color <- sample.int(6, length(SNP), replace=TRUE)
sample.gr$border <- sample(c("gray80", "gray30"), length(SNP), replace=TRUE)
lolliplot(sample.gr, features)
 
R語言畫棒棒糖圖展示snp在基因上的位置是怎樣的

關(guān)于R語言畫棒棒糖圖展示snp在基因上的位置是怎樣的就分享到這里了,希望以上內(nèi)容可以對大家有一定的幫助,可以學(xué)到更多知識(shí)。如果覺得文章不錯(cuò),可以把它分享出去讓更多的人看到。

分享文章:R語言畫棒棒糖圖展示snp在基因上的位置是怎樣的
網(wǎng)站路徑:http://m.newbst.com/article26/jocccg.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供App開發(fā)標(biāo)簽優(yōu)化云服務(wù)器網(wǎng)站收錄網(wǎng)站策劃建站公司

廣告

聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請盡快告知,我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場,如需處理請聯(lián)系客服。電話:028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時(shí)需注明來源: 創(chuàng)新互聯(lián)

搜索引擎優(yōu)化