我正在尝试使用 BioPython 中 codeml 的 aaml 部分在系统发育树中找到一堆氨基酸序列的 dN/ds 值。不幸的是，运行它时出现错误（结果文件为空）。
from Bio.Phylo.PAML import codeml
cml = codeml.Codeml(alignment = "family.phy", tree = "family.phy_phyml_tree.txt", out_file = "results.out")
cml.read_ctl_file("aaml.ctl")
cml.alignment = "family.phy"
cml.tree = "family.phy_phyml_tree.txt"
cml.out_file = "results.out"
results=cml.run()
我得到的错误是
Traceback (most recent call last):
File "<pyshell#57>", line 1, in <module>
cml.run()
File "C:\Python33\lib\site-packages\Bio\Phylo\PAML\codeml.py", line 189, in run
results = read(self.out_file)
File "C:\Python33\lib\site-packages\Bio\Phylo\PAML\codeml.py", line 210, in read
raise ValueError("Invalid results file")
ValueError: Invalid results file
有没有人得到这个为他们工作？