R二项分布检验与FDR校正

二项分布是重复n次的实验，且每次实验都是独立的，只有两种结果，并且相互对立的，生活中最常见的是投硬币~~~在生物领域内也有很多符合此类分布的，如二倍体动物等位基因，来源于父本和母本的重组等。具体公式什么的博主就不写了，写个关于ASE的例子吧。

# cat binom.r | R --slave --args <file>
args <- commandArgs()
fa <- read.table(args[4], header=FALSE, sep="\t")
n1 = fa$V8
n2 = fa$V12
len = length(n1)
pv = numeric(len)
for(i in 1:len){
pv[i] = 0
if(n1[i] > n2[i]) pv[i]=binom.test(n1[i], n1[i]+n2[i], p=1/2, alternative="greater")$p.value else pv[i]=binom.test(n2[i], n1[i]+n2[i], p=1/2, alternative="greater")$p.value
}
qv <- p.adjust(pv, method="fdr")#fdr校正
fa$pv = formatC(pv, digits=4)
fa$fdr = formatC(qv, digits=4)
write.table(fa, file=paste(args[4],".out",sep=""), sep="\t", col.names=FALSE, row.names=FALSE, quote=FALSE)

其实也不算是什么例子&#