> ExomeCount[1:3,c('22S03466993','22S03467008','22S03467008')]
22S03466993 22S03467008 22S03467008.1
1 435 596 596
2 656 829 829
3 667 939 939
> aa = colnames(ExomeCount[1:3,c('22S03466993','22S03467008','22S03467008')])
> aa
[1] "22S03466993" "22S03467008" "22S03467008.1"
> ExomeCount[1:3,aa]
Error in `[.data.frame`(ExomeCount, 1:3, aa) : undefined columns selected
ExomeCount 的colnames里面本来只有一个22S03467008
我在调用的时候加入了两个 “22S03467008”, 那第二个“22S03467008”就成了“22S03467008.1”,
但再把“22S03467008.1” 拿去到ExomeCount中提取数据的时候就会报错。
我的代码运行的报错如下:
Optimization of the choice of aggregate reference set, this process can take some time
Number of selected bins: 1932
Error in `[.data.frame`(ExomeCount, , my.choice$reference.choice, drop = FALSE) :
undefined columns selected
Calls: as.matrix -> [ -> [.data.frame
In addition: There were 50 or more warnings (use warnings() to see the first 50)
Execution halted
这个坑太费时间了。