I am running the aaml program of PAML using a amino acid alignment of several species. When I run the aaml.ctl file with my tree and phylip file, I receive the following error:
Error: initial value alpha <= 0 for fix_alpha = 0
Any ideas on what may be causing this error?
Thanks for the help!
Here is part of my ctl file:
seqfile = /gene.phylip
treefile = /gene.nh
outfile = /gene.out
noisy = 9
verbose = 0
runmode = 0
seqtype = 2
aaRatefile = /dat/wag.dat
model = 3
Mgene = 0
fix_alpha = 0
alpha = 0.
Malpha = 0
ncatG = 2
clock = 0
getSE = 0
RateAncestor = 1
Small_Diff = .5e-6
cleandata = 0
method = 0
Error: initial value alpha <= 0 for fix_alpha = 0
Any ideas on what may be causing this error?
Thanks for the help!
Here is part of my ctl file:
seqfile = /gene.phylip
treefile = /gene.nh
outfile = /gene.out
noisy = 9
verbose = 0
runmode = 0
seqtype = 2
aaRatefile = /dat/wag.dat
model = 3
Mgene = 0
fix_alpha = 0
alpha = 0.
Malpha = 0
ncatG = 2
clock = 0
getSE = 0
RateAncestor = 1
Small_Diff = .5e-6
cleandata = 0
method = 0