python 相关性fdr矫正_科学网—根据p-value列表计算FDR（Benjamini-Hochberg） - 谭月茜的博文...

#!/bin/bash

#===========================================================

#    File Name:  calfdr

#      Version:  1.0

#       Author:  tanyuexi

#      Contact:  tanyuexi@gmail.com

#         Date:  Thu Oct 18 14:09:48 CST 2012

#  Discription:  This script is used to calculate FDR using Benjamini-Hochberg method.

#===========================================================

bn=`basename $0`

if [ -z $2 ]

then

echo -e "Usage: $bn In_p-value_list OutFile"

exit 1

fi

infile=$1

outfile=$2

cat >$infile.R <<__eof__>

p.srt

p.rank

output

fdrlevels

for(fdr in fdrlevels){for(i in length(p.srt):1){tmp

write(output,file="$outfile",ncolumn=1)

quit(save="no")

__EOF__

R CMD BATCH $infile.R $infile.Rout && rm $infile.R $infile.Rout

#EOF

转载本文请联系原作者获取授权，同时请注明本文来自谭月茜科学网博客。

链接地址：http://blog.sciencenet.cn/blog-747197-767694.html

上一篇：定时按照原目录结构自动备份.pl和.sh等脚本文件

下一篇：更新系统及软件，清除旧内核