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
誰かがこれを彼らのために働かせましたか？