par中的参数srt和crt都是对应text显示的。
而我看了par中的参数只有las类似于这样的功能，但是也仅限于平行或者垂直与坐标轴。
显然不符合我的要求
遂用text来做labels，下面是例子
par(srt=45,mai=c(2,1,1,1),xpd=T,adj=1)
barplot(rbind(Aco.comp1[1:10,3],Aco.comp1[1:10,4]),beside=T,horiz=F)
text(seq(2,29,by=3),rep(-0.01,10),labels=Aco.comp1[1:10,2],cex=1)
至于text的坐标选择，就是自己画几次看看，然后确定一下最终的选择值。比如我就是花了几次发现-0.01比较合适
数据分析中常碰见多重检验问题(multiple
testing).Benjamini于1995年提出一种方法,通过控制FDR(False Discovery
Rate)来决定P值的域值.
假设你挑选了R个差异表达的基因，其中有S个是真正有差异表达的，另外有V个其实是没有差异表达的，是假阳性的.实践中希望错误比例Q＝V/R平均而言不能超过某个预先设定的值（比如0.05），在统计学上，这也就等价于控制FDR不能超过5％.
根据Benjamini在他的文章中所证明的定理，控制fdr的步骤实际上非常简单。
设总共有m个候选基因，每个基因对应的p值从小到大排列分别是p(1),p(2),...,p(m),则若想控制fdr不能超过q，则只需找到最大的正整数i，使得
p(i)<=
(i*q)/m.然后，挑选对应p(1),p(2),...,p(i)的基因做为差异表达基因，这样就能从统计学上保证fdr不超过q。
在网上找了半天，没有找到好的答案，只好看manual
不过太长了， 也不理解
最后搞懂了以下两条希望对别人能有所帮助：
BAM＝＝>SAM
samtools view -h -o out.sam out.bam
SAM＝＝>BAM
samtools view -bS out.sam >out.bam
-b 意思使输出使BAM format
-S 意思使输入使SAM，如果@SQ 缺剩， 要写-t
所以如果没有@SQ
samtools faidx ref.fa
samtools view -bt ref.fa.fai out.sam > out.bam
加载中…
