...dlist for changing###############
open(IN,"$opts{idlist}")|| die "open ID file failed\n";
while(<IN>){
chomp;
my($newID,$gene)=split(/\s+/,$_);
if($gff=~/$gene/){
$gff=~ s/$gene/$newID/g;
}
}
close(IN);
#####################print out new gff#######################
ope...
...codecsimport uuidfrom PIL import ImageBin=os.path.split(os.path.realpath(__file__))[0]
##########################################################################
dom = xml.dom.minidom.parse(sys.argv[1])root = dom.documentElementff=open(sys.argv[2],"w")if root.nodeName=='pathway': for i i...
...rd
在癌症研究中：
hazard ratio > 1 is called bad prognostic factorhazard ratio < 1 is called good prognostic factor
z(-3.176)值代表Wald统计量，其值等于回归系数coef除以其标准误se(coef)，即z = coef/se(coef)；有统计量必有其对应的假设检验的显著性...
...9.chromosome.22.gtf -p gene_length
Traceback (most recent call last):
File "/home/zpp/work/my_rnaseq/scripts/get_gene_length_from_gtf.py", line 49, in <module>
if kvs['gene_id'] in geneL and kvs['transcript_id'] in geneL[kvs['gene_id']]:
KeyError: 'gene_id'
....validations.output这个文件
sh: alignment.validations.output: No such file or directoryError, cmd: update_alignment_status.dbi -M 'pasa.sqlite' < alignment.validations.output > pasa_run.log.dir/alignment.validation_loading.output died with ret 256 No such file or directory at PASAp...