如何在R中输入数据使用GenomicRanges包进行分析?



我想让数据像这样

library(Gviz)
library(GenomicRanges)
#Load data : class = GRanges
data(cpgIslands)
cpgIslands
## GRanges with 10 ranges and 0 metadata columns:
##        seqnames               ranges strand
##                         
##    [1]     chr7 [26549019, 26550183]      *
##    [2]     chr7 [26564119, 26564500]      *
##    [3]     chr7 [26585667, 26586158]      *
##    [4]     chr7 [26591772, 26593309]      *
##    [5]     chr7 [26594192, 26594570]      *
##    [6]     chr7 [26623835, 26624150]      *
##    [7]     chr7 [26659284, 26660352]      *
##    [8]     chr7 [26721294, 26721717]      *
##    [9]     chr7 [26821518, 26823297]      *
##   [10]     chr7 [26991322, 26991841]      *
##   ---
##   seqlengths:
##    chr7
##      NA
#Annotation track, title ="CpG"
atrack <- AnnotationTrack(cpgIslands, name = "CpG")
plotTracks(atrack)

这是我的数据

id  chr start   end
wnt4    1A  14938532    14938648
wnt5    1A  14940394    14940633
wnt6    1A  14943914    14944217
wnt7    1A  14945438    14945867
wnt1    1A  14946238    14946546

我的问题是如何使我的数据看起来像上面的数据?

有人能帮我吗?我不太熟悉输入数据。

请帮忙!

您可以将R中的数据作为数据框导入,而不是像这里解释的那样应用makeGRangesFromDataFrame函数https://web.mit.edu/~r/current/arch/i386_linux26/lib/R/library/GenomicRanges/html/makeGRangesFromDataFrame.html

最新更新