脚本 /TJPROJ6/RNA_SH/personal_dir/zhangxin/jiaoben/intron_ratio/intron_ratio.py

该分析首先是对gene进行定量，去掉fpkm小于1的gene，然后对gene中intron进行定量，计算fpkm。intron保留率的计算根据intron的fpkm比上该intron所在gene全部exon的reads计算的fpkm；在计算intron保留率的差异时，选择表达量差异在2倍以上，使用t检验，阈值一选择矫正后的p值小于0.05。

原始数据或bam文件