
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:26 04/04/2011

Input data from file "591-tomodel-phylip/chr11_1600.phylip" (7 taxa, 1137 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: GTRG


(Output sorted by AIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	21	-3358.9353	6759.8706	0.5678	0.5678
GTRG	20	-3360.2211	6760.4422	0.4267	0.9945
GTRI	20	-3364.5695	6769.1390	0.0055	1.0000
SYMIG	18	-3372.6799	6781.3598	0.0000	1.0000
SYMG	17	-3373.8740	6781.7479	0.0000	1.0000
SYMI	17	-3377.6821	6789.3643	0.0000	1.0000
HKYG	16	-3383.5967	6799.1935	0.0000	1.0000
HKYIG	17	-3382.6530	6799.3059	0.0000	1.0000
HKYI	16	-3386.8042	6805.6085	0.0000	1.0000
K2PG	13	-3390.8688	6807.7377	0.0000	1.0000
K2PIG	14	-3389.9381	6807.8762	0.0000	1.0000
K2PI	13	-3393.9185	6813.8370	0.0000	1.0000
GTR	19	-3425.7958	6889.5917	0.0000	1.0000
F81G	15	-3435.9530	6901.9060	0.0000	1.0000
F81IG	16	-3435.2803	6902.5606	0.0000	1.0000
F81I	15	-3438.4827	6906.9653	0.0000	1.0000
SYM	16	-3437.8607	6907.7214	0.0000	1.0000
JC69G	12	-3441.8750	6907.7500	0.0000	1.0000
JC69IG	13	-3441.1948	6908.3897	0.0000	1.0000
JC69I	12	-3444.3722	6912.7445	0.0000	1.0000
HKY	15	-3447.6916	6925.3832	0.0000	1.0000
K2P	12	-3454.7525	6933.5049	0.0000	1.0000
F81	14	-3496.1596	7020.3192	0.0000	1.0000
JC69	11	-3502.0408	7026.0816	0.0000	1.0000



(Output sorted by BIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-3360.2211	6861.1651	0.8521	0.8521
GTRIG	21	-3358.9353	6865.6297	0.0914	0.9435
SYMG	17	-3373.8740	6867.3624	0.0384	0.9820
GTRI	20	-3364.5695	6869.8620	0.0110	0.9930
SYMIG	18	-3372.6799	6872.0105	0.0038	0.9967
K2PG	13	-3390.8688	6873.2076	0.0021	0.9988
SYMI	17	-3377.6821	6874.9788	0.0009	0.9997
K2PIG	14	-3389.9381	6878.3823	0.0002	0.9998
K2PI	13	-3393.9185	6879.3070	0.0001	0.9999
HKYG	16	-3383.5967	6879.7718	0.0001	1.0000
HKYIG	17	-3382.6530	6884.9204	0.0000	1.0000
HKYI	16	-3386.8042	6886.1869	0.0000	1.0000
JC69G	12	-3441.8750	6968.1838	0.0000	1.0000
JC69I	12	-3444.3722	6973.1783	0.0000	1.0000
JC69IG	13	-3441.1948	6973.8596	0.0000	1.0000
F81G	15	-3435.9530	6977.4482	0.0000	1.0000
F81I	15	-3438.4827	6982.5075	0.0000	1.0000
F81IG	16	-3435.2803	6983.1389	0.0000	1.0000
GTR	19	-3425.7958	6985.2785	0.0000	1.0000
SYM	16	-3437.8607	6988.2997	0.0000	1.0000
K2P	12	-3454.7525	6993.9387	0.0000	1.0000
HKY	15	-3447.6916	7000.9254	0.0000	1.0000
JC69	11	-3502.0408	7081.4792	0.0000	1.0000
F81	14	-3496.1596	7090.8253	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRIG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_1600.phylip = 1 - 1137;
 Partition Dummy = 1:591-tomodel-phylip/chr11_1600.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRIG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_1600.phylip = 1 - 1137;
 Partition Dummy = 1:591-tomodel-phylip/chr11_1600.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_1600.phylip = 1 - 1137;
 Partition Dummy = 1:591-tomodel-phylip/chr11_1600.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


-------------------------------------------------------------
End of Output

