Subject: [modeller_usage] how to evaluate 100 profiles?
From: hori koshii <>
Date: Fri, 20 Oct 2006 12:58:19 -0700 (PDT)
Hi, does anyone know python and modeller very well and give us some help to run a model evaluation script for 100 homology models.
(how to set 'EGPh.B99990001.pdb' as a variable and make it like a loop so the number can be increased by 1 in the next run), or another other ways you know?
It will be so appreciated and thanks.
The items need to be updated in the loop are bolded in the following script.
---------------------------------------------------------------------------------------------------------------------- from modeller import * from modeller.automodel import * # Load the automodel class
log.verbose() # request verbose output env = environ() # create a new MODELLER environment to build this model in env.libs.topology.read(file='$(LIB)/top_heav.lib') # read
topology env.libs.parameters.read(file='$(LIB)/par.lib') # read parameters
# directories for input atom files env.io.atom_files_directory = './:../atom_files'