...t; ABC = lapply (tem,function(x)read.table(x,sep="\t"))
ABC 是一个列表数据,可以通过下面的方式获得每个文件内的内容,例如第一个文件:
> ABC[[1]] V1 V21 A 1>
把文件的名字(前缀)赋给每个读取到的数据
> for (i in 1:length(ABC)){+ assign(m...
从TCGA数据中提取lncRNA的表达量时,需要知道lncRNA的编号和对应的名称。这些信息可以从GTF文件中提取。提取的话,可以采用如下的代码实现。
#!/usr/bin/perl -w
use strict;
my $biotype_file = shift @ARGV;
my $gtf = shift @ARGV;
my $biotype = shift @A...