[Bio] / FigTutorial / show_functional_coupling.py Repository:
ViewVC logotype

View of /FigTutorial/show_functional_coupling.py

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1 - (download) (as text) (annotate)
Thu Jun 10 14:51:56 2004 UTC (15 years, 6 months ago) by disz
Branch: MAIN
CVS Tags: merge-bodev_news-3, rast_rel_2014_0912, rast_rel_2008_06_18, rast_rel_2008_06_16, rast_rel_2008_07_21, rast_rel_2010_0928, rast_2008_0924, Root-bobdev_news, rast_rel_2008_09_30, caBIG-13Feb06-00, rast_rel_2010_0526, rast_rel_2014_0729, merge-trunktag-bobdev_news-1, rast_rel_2009_05_18, caBIG-05Apr06-00, rast_rel_2009_0925, rast_rel_2010_1206, rast_rel_2010_0118, caBIG-00-00-00, rast_rel_2009_02_05, rast_rel_2011_0119, rast_rel_2008_12_18, merge-trunktag-bodev_news-3, merge-bobdev_news-2, merge-bobdev_news-1, rast_rel_2008_10_09, rast_release_2008_09_29, rast_rel_2008_04_23, rast_rel_2008_08_07, rast_rel_2009_07_09, rast_rel_2010_0827, myrast_33, rast_rel_2011_0928, rast_rel_2008_09_29, rast_rel_2008_10_29, rast_rel_2009_03_26, merge-trunktag-bobdev_news-2, rast_rel_2008_11_24, HEAD
Branch point for: Branch-bobdev_news
Python version of tutorial code

from FigKernelPackages import FIG2
import sys, re

fig = FIG2.FIG()

def die(msg):
	print msg
	sys.exit(0)

if len(sys.argv) != 4:
	die("usage: fc Bound SimCutoff FcCutOff < PEGs > PEG1-Sc-PEGs")


bound = sys.argv[1]
sim_cutoff = sys.argv[2]
fc_cutoff = sys.argv[3]

for line in sys.stdin.readlines():	
    match =  re.search("^(fig\|\d+\.\d+\.peg\.\d+)", line)
    if match:
    	peg = match.group(1)
	fc  = fig.coupling_and_evidence(peg, bound, sim_cutoff, fc_cutoff)
	for tuple in fc:
	    print "%s\t%s\t%s" % ( peg, tuple[0], tuple[1])

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3