如何获取基因的长度（非冗余外显子(EXON)长度之和）

  gle = lapply(split(k,k$V8),function(x){
    tmp=apply(x,1,function(y){
        y[1]:y[2]
    })
    length(unique(unlist(tmp)))
  })
  gle=data.frame(gene_name=names(gle),
               length=as.numeric(gle))
  save(gle,file = "gle.Rdata")

根据提供的代码段，可以看出这是R语言中的一段代码，使用了lapply函数和apply函数对数据进行处理。代码中使用了split函数对exon数据按照gene_name进行分组，并将每个组中的数据进行处理，最终保存到gle数据框中。

K

image.png