全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
11696 2
2016-11-21
现在有两个数据表,一个是染色体上的基因片段所在的位置,一个是结合在染色体上的RNA所在的位置。人体内23条染色体和两条性染色体的数据是乱序的,基因片段的数目和RNA的数目也不相等,如图:
吃啊.PNG
我需要先判断基因和RNA是否在同一条染色体上,再对它们的位置进行比较,位置有交集的,则输出RNA的名字,我的代码是:

lncRNA<-read.csv("E:/THU_research_intern/lncRNA/data/2012data/lncRNA.csv",header=TRUE,stringsAsFactors = FALSE)
ChP<-read.csv("E:/THU_research_intern/lncRNA/data/2012data/GSM1436260_ENCFF002ENI.csv",header=TRUE,stringsAsFactors = FALSE)
ChP<-ChP[order(ChP[,1]),]
for(x in 1:2708)
  for(i in 1:118484)
    if(lncRNA[i,1]==ChP[x,1]&&lncRNA[x,3]<=ChP[i,2]||lncRNA[x,1]==ChP[i,1]&&lncRNA[x,2]>=ChP[i,6]) lncRNA=lncRNA[-c(i),]
lncRNA[,4]

错误信息是,
Error in if (lncRNA[i, 1] == ChP[x, 1] && lncRNA[x, 3] <= ChP[i, 2] ||  :
  missing value where TRUE/FALSE needed

应该怎么写呢?
附件列表
然.PNG

原图尺寸 15.99 KB

然.PNG

吃啊.PNG

原图尺寸 17 KB

吃啊.PNG

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2016-11-23 10:44:52
不知道楼主的错误是什么问题,我用了另一种方法(没有原始数据,自己编了几个,楼主参考一下吧):
复制代码
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2016-11-24 16:08:04
能把2个数据表贴出来看看嘛? 不清楚你得数据结构哦
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群