求助？用R语言比较数字大小出错了

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

想挑出一个range内，存在的位点，编程如下，但是不知道为什么有的结果正确有的结果错的。求指出！！
cgi <- read.csv("HumanMethylation450_15017482_v1-2.csv",sep=",", header = F, stringsAsFactors = F)
gene<- read.csv("sample_Q.csv", sep=",", header = F, stringsAsFactors = F)
name <- c(paste0("V",c(1:37)))
info <- data.frame(sapply(name, function(x) x = character(0)), stringsAsFactors = F)
as.numeric(cgi_tmp$V13)
for(i in c(1:22,"X","Y")){
  gene_tmp <- gene[gene$V2 == paste0(i),]
  cgi_tmp <- cgi[cgi$V12 == paste0(i),]
   fun <- function(x){
    logic1 <- x <= as.numeric(gene_tmp$V4)
    logic2 <- x >= as.numeric(gene_tmp$V3)
    idx1 <- which(logic1)
    idx2 <- which(logic2)
    c1 <- intersect(idx1, idx2)
      if(length(c1) != 0){
      info[(nrow(info)+1):(nrow(info)+length(c1)),] <<- c(cgi_tmp[cgi_tmp$V13 == x,1:33], gene_tmp[c1, 1:4])
      return(info)
    }
  }
  mapply(fun,cgi_tmp$V13,SIMPLIFY = T)
}

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：R语言 range 存在的 不知道 RAN