从表格数据到 Genomic Ranges
在视频中,您学习了创建 GRanges 对象的方法。您可以用区间的名称、起始位点和结束位点(seqnames、start 和 end)来定义一个 GRange。若数据是表格格式,也可以将其转换为 GRanges 对象。我们将使用一个名为 seq_intervals 的数据框,因为序列区间很可能就是以这种形式保存的。注意:如果您更习惯使用 tibble,也同样适用。
您将使用预先定义好的数据框 seq_intervals,通过 as() 函数将其转换为 GRanges 对象。as() 函数在上一个视频中介绍过——它接收一个对象以及要转换成的类名。
本练习是课程的一部分
R 中的 Bioconductor 入门
练习说明
- 加载
GenomicRanges。 - 打印
seq_intervals,查看其结构。 - 将
seq_intervals转换为GRanges对象,并将新对象命名为myGR。 - 打印
myGR。
交互式实操练习
通过完成这段示例代码来试试这个练习。
# Load GenomicRanges package
library(GenomicRanges)
# Print seq_intervals
___
# Create myGR
___ <- ___(___, "___")
# Print myGR
___