利用modeler和python进行蛋白质预测

Question

我正在尝试运行一个python脚本，它是作为modeller软件的教程由sailab编写的。本教程中有一个包含以下代码的python文件：

import pylab
import modeller

def r_enumerate(seq):
    """Enumerate a sequence in reverse order"""
    # Note that we don't use reversed() since Python 2.3 doesn't have it
    num = len(seq) - 1
    while num >= 0:
        yield num, seq[num]
        num -= 1

def get_profile(profile_file, seq):
    """Read `profile_file` into a Python array, and add gaps corresponding to
       the alignment sequence `seq`."""
    # Read all non-comment and non-blank lines from the file:
    f = open(profile_file)
    vals = []
    for line in f:
        if not line.startswith('#') and len(line) > 10:
            spl = line.split()
            vals.append(float(spl[-1]))
    # Insert gaps into the profile corresponding to those in seq:
    for n, res in r_enumerate(seq.residues):
        for gap in range(res.get_leading_gaps()):
            vals.insert(n, None)
    # Add a gap at position '0', so that we effectively count from 1:
    vals.insert(0, None)
    return vals

e = modeller.environ()
a = modeller.alignment(e, file='TvLDH-1bdmA.ali')

template = get_profile('1bdmA.profile', a['1bdmA'])
model = get_profile('TvLDH.profile', a['TvLDH'])

# Plot the template and model profiles in the same plot for comparison:
pylab.figure(1, figsize=(10,6))
pylab.xlabel('Alignment position')
pylab.ylabel('DOPE per-residue score')
pylab.plot(model, color='red', linewidth=2, label='Model')
pylab.plot(template, color='green', linewidth=2, label='Template')
pylab.legend()
pylab.savefig('dope_profile.png', dpi=65)

当我尝试在cmd中运行程序时，出现了一个错误： 没有名为pylab的模块 但是我已经安装了matplotlib、numpy和scipy，而且当我在visualstudio中运行它时，也没有问题。有人能帮我理解问题是什么吗？你知道吗

Answer 1

您的python路径对于命令行来说可能是一团糟。您很可能安装了多个版本，而命令行的默认版本没有安装pylab。你知道吗

请看……
How to add to the PYTHONPATH in Windows, so it finds my modules/packages?