Exemplo n.º 1
0
    # Run MCMC
    result.sample(model,
                  progress_bar=progress_bar,
                  **sample_params)


def cl(PROJECTS=sorted(PB_PROJECTS.keys()),
       MODELS=('gpa333',),
       FREQS=(4.0, 8.0, 0.5, 1.0, 2.0, 16.0, 32.0, 40.0),  # angular frequencies (rad/s)
       log_dir='~'):
    """Generates command lines to launch the script with different parameters."""
    for project in PROJECTS:
        pbproject = PB_PROJECTS[project]
        for model, af, genotype in product(MODELS, FREQS, pbproject.genotypes()):
            expid = '%s__%s__%s__%g' % (project, model, genotype, af)
            print('python2 -u thebas/sinefitting/samplers.py sample '
                  '--pbproject %s '
                  '--freq %g '
                  '--genotype-id %s '
                  '--model-id %s '
                  '&>%s/%s.log' %
                  (project, af, genotype, model, log_dir, expid))


if __name__ == '__main__':
    import argh
    parser = argh.ArghParser()
    parser.add_commands([cl, sample])
    parser.dispatch()
Exemplo n.º 2
0
def main():
    parser = argh.ArghParser()
    parser.prog = 'samlkeygen'
    parser.add_commands(
        [authenticate, list_profiles, select_profile, run_command, version])
    parser.dispatch()