annotate deseq/differential_expression_analysis_pipeline_for_rnaseq_data-a03838a6eb54/DiffExpAnal/loadHTSeqResultFiles.py @ 10:6e573fd3c41b draft

Uploaded
author yufei-luo
date Mon, 13 May 2013 10:06:30 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
10
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
1 #!/usr/bin/env python
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
2
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
3 """
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
4 Yufei LUO
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
5 """
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
6
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
7
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
8
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
9 import optparse, sys
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
10
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
11
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
12 def __main__():
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
13 #Parse Command Line
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
14 parser = optparse.OptionParser()
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
15 parser.add_option('-i', '--inputs', dest='inputFiles', default=None, help='several input files. (seperated by @ or @@' )
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
16 parser.add_option( '-o', '--output', dest='outputFile', default=None, help='The output list of HTSeq results files(.tabular) on txt format.' )
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
17 ( options, args ) = parser.parse_args()
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
18
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
19
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
20 out = open(options.outputFile, 'w')
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
21 out.write("label\tfiles\tgroup\n")
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
22 if options.inputFiles == None:
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
23 raise Exception, 'input file name is not defined!'
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
24
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
25 groupCount = 1
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
26 fileCount = 0
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
27
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
28 inputFiles = sys.argv[6:]
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
29 print '\n\nthe length of inputfiles is : %s \n' % len(inputFiles)
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
30 i = 0
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
31 while i < (len(inputFiles)-1):
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
32 if inputFiles[i] == "@":
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
33 i += 1
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
34 fileCount = 1
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
35 groupCount += 1
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
36 out.write("Group%s_%s\t%s\t%s\n" % (groupCount, fileCount, inputFiles[i], groupCount))
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
37 else:
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
38 fileCount += 1
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
39 out.write("Group%s_%s\t%s\t%s\n" % (groupCount, fileCount, inputFiles[i], groupCount))
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
40 i += 1
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
41
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
42 out.close()
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
43
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
44
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
45
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
46 if __name__=="__main__": __main__()