elif figure == '3':
gf.fig3()
elif figure == '4':
gf.fig4()
elif figure == 'S1':
gf.figS1()
elif figure == 'S2':
gf.figS2()
elif figure == 'S3':
gf.figS3()
else:
print "Command not recognized"
else:
gd.HMP_OTU_to_sparese_SbyS()
gd.Match_NAP_to_sparse_SbyS()
gd.get_SADs_HMP()
gd.get_SADs_mgrast()
gd.merge_SADs_mgrast()
methods = ['geom','lognorm', 'mete','zipf']
gd.generate_obs_pred_data(['EMPclosed','HMP', 'MGRAST'], methods, remove_obs = 0)
gd.generate_obs_pred_data(['EMPclosed','HMP', 'MGRAST'], methods, remove_obs = 1)
gd.generate_obs_pred_data(['95','97', '99'], methods, remove_obs = 0)
gd.stratifyDataOnce(['EMPclosed','HMP', 'MGRAST'], remove_obs = 0)
gd.stratifyDataOnce(['EMPclosed','HMP', 'MGRAST'], remove_obs = 1)
gd.stratifyDataBootstrap(remove_obs = 0)
gd.stratifyDataBootstrap(remove_obs = 1)
gd.stratifyDataBootstrap(remove_obs = 0, seqSim = '95')
gd.stratifyDataBootstrap(remove_obs = 0, seqSim = '97')
gd.stratifyDataBootstrap(remove_obs = 0, seqSim = '99')
from __future__ import division
import mete as mete
import os, math
import numpy as np
import modelsAndMetrics as mo
import macroeco_distributions as md
from scipy import stats
import generateData as gd
mydir = os.path.expanduser("~/github/MicroMETE/")
import importData as im
gd.stratifyDataBootstrap(remove_obs = 0)
#gd.stratifyDataBootstrap(remove_obs = 1)
#gd.stratifyDataBootstrap(remove_obs = 0, seqSim = '95')
#gd.stratifyDataBootstrap(remove_obs = 0, seqSim = '97')
#gd.stratifyDataBootstrap(remove_obs = 0, seqSim = '99')
#gd.stratifyDataOnce(['EMPclosed', 'MGRAST', 'HMP'], remove_obs = 0)
elif figure == 'TS4':
gf.tableS4()
elif figure == 'TS5':
gf.tableS5()
elif figure == 'TS6':
gf.tableS6()
else:
print "Command not recognized"
else:
gd.HMP_OTU_to_sparese_SbyS()
gd.Match_NAP_to_sparse_SbyS()
gd.get_SADs_HMP()
gd.get_SADs_mgrast()
gd.merge_SADs_mgrast()
methods = ['geom', 'lognorm', 'mete', 'zipf']
gd.generate_obs_pred_data(['EMPclosed', 'HMP', 'MGRAST'],
methods,
remove_obs=0)
gd.generate_obs_pred_data(['EMPclosed', 'HMP', 'MGRAST'],
methods,
remove_obs=1)
gd.generate_obs_pred_data(['95', '97', '99'], methods, remove_obs=0)
gd.stratifyDataOnce(['EMPclosed', 'HMP', 'MGRAST'], remove_obs=0)
gd.stratifyDataOnce(['EMPclosed', 'HMP', 'MGRAST'], remove_obs=1)
gd.stratifyDataBootstrap(remove_obs=0)
gd.stratifyDataBootstrap(remove_obs=1)
gd.stratifyDataBootstrap(remove_obs=0, seqSim='95')
gd.stratifyDataBootstrap(remove_obs=0, seqSim='97')
gd.stratifyDataBootstrap(remove_obs=0, seqSim='99')