问题描述
我下载了.bed文件,并使用此https://genome.ucsc.edu/cgi-bin/hgLiftOver
对其进行了转换
现在我想使用此代码将转换后的文件"datahg38.bed"
加载到R中
library(rtracklayer)
dataset<-import.bed("datahg38.bed")
我收到此错误
error : $ operator is invalid for atomic vectors
有什么解决办法吗?
所以我发现输出格式是300,000 x 1的列表,应该是300,000的3个变量的表
格式是这样的
chr1:10142-10351
chr1:10453-10563
chr1:13044-13104
.
.
.
所以基本上我需要使用read.table
将其读取为表格,然后将其转换为300,000 of 3的矩阵。字符:
和-
应该用于分割数据。
我正在寻找这样的输出
chr1 10142 10351
chr1 10453 10563
chr1 13044 13104
...
那么有人可以建议一种有效的方法来分离这些数据吗?
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)