Exemplo n.º 1
0
import os
from query import consistency
from query import getdata
# from query import queryrna
# queryrna()

orig_pdb  = int(os.system("ls data/Pdb/ | wc -l"))
rna_pdb   = int(os.system("ls data/OnlyNA/ | wc -l"))
inp_pdb   = int(os.system("ls data/Inp/ | wc -l"))
rel_date  = int(os.system("grep \"REVDAT\" data/Pdb/*.pdb | awk '{ if ($5 ~ /^0/) print $0}' | wc -l"))

consistency(orig_pdb,rna_pdb,inp_pdb)

# Uncomment to check that the release dates match
# consistency(orig_pdb,rna_pdb,rel_date)

getdata()
Exemplo n.º 2
0
================================================================================
"""
import subprocess
import query
# from query import analysis
# from query import consistency

query.queryrna()

orig_pdb = subprocess.check_output('ls data/Pdb/ | wc -l', shell=True)
print "total number of pdb files = %d" % int(orig_pdb)
rna_pdb = subprocess.check_output('ls data/OnlyNA/ | wc -l', shell=True)
print "rna pdbs = %d" % int(rna_pdb)
inp_pdb = subprocess.check_output('ls data/Inp/ | wc -l', shell=True)
print "inp pdbs = %d" % int(inp_pdb)
query.consistency(orig_pdb, rna_pdb, inp_pdb)

# query.makedirs()

# query.download()

query.analysis()

query.helices()

query.getdata()

#####################################################################
# -4-
# Plot number of bases vs. year
# This module automatically generates plots after subsetting the