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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2352.3649	4736.7297	0.5052	0.5052
HKYIG	17	-2352.0423	4738.0846	0.2566	0.7618
HKYI	16	-2353.6756	4739.3512	0.1362	0.8980
GTRG	20	-2350.5326	4741.0651	0.0578	0.9558
GTRIG	21	-2350.2342	4742.4685	0.0287	0.9845
GTRI	20	-2351.8484	4743.6969	0.0155	1.0000
HKY	15	-2393.4878	4816.9757	0.0000	1.0000
GTR	19	-2391.9392	4821.8785	0.0000	1.0000
F81G	15	-2404.6685	4839.3371	0.0000	1.0000
F81I	15	-2405.4859	4840.9719	0.0000	1.0000
F81IG	16	-2404.4970	4840.9940	0.0000	1.0000
K2PG	13	-2422.1867	4870.3733	0.0000	1.0000
K2PIG	14	-2421.8964	4871.7928	0.0000	1.0000
K2PI	13	-2423.1919	4872.3839	0.0000	1.0000
SYMG	17	-2419.3454	4872.6909	0.0000	1.0000
SYMIG	18	-2418.9802	4873.9604	0.0000	1.0000
SYMI	17	-2420.5269	4875.0539	0.0000	1.0000
F81	14	-2442.5394	4913.0787	0.0000	1.0000
K2P	12	-2461.8594	4947.7187	0.0000	1.0000
JC69G	12	-2462.9270	4949.8540	0.0000	1.0000
SYM	16	-2459.0194	4950.0387	0.0000	1.0000
JC69I	12	-2463.7230	4951.4461	0.0000	1.0000
JC69IG	13	-2462.7410	4951.4820	0.0000	1.0000
JC69	11	-2500.6193	5023.2386	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2352.3649	4812.5212	0.7592	0.7592
HKYI	16	-2353.6756	4815.1427	0.2047	0.9639
HKYIG	17	-2352.0423	4818.6130	0.0361	1.0000
GTRG	20	-2350.5326	4835.8045	0.0000	1.0000
GTRI	20	-2351.8484	4838.4362	0.0000	1.0000
GTRIG	21	-2350.2342	4841.9448	0.0000	1.0000
HKY	15	-2393.4878	4888.0302	0.0000	1.0000
F81G	15	-2404.6685	4910.3916	0.0000	1.0000
GTR	19	-2391.9392	4911.8809	0.0000	1.0000
F81I	15	-2405.4859	4912.0264	0.0000	1.0000
F81IG	16	-2404.4970	4916.7855	0.0000	1.0000
K2PG	13	-2422.1867	4931.9539	0.0000	1.0000
K2PI	13	-2423.1919	4933.9644	0.0000	1.0000
K2PIG	14	-2421.8964	4938.1103	0.0000	1.0000
SYMG	17	-2419.3454	4953.2193	0.0000	1.0000
SYMI	17	-2420.5269	4955.5823	0.0000	1.0000
SYMIG	18	-2418.9802	4959.2258	0.0000	1.0000
F81	14	-2442.5394	4979.3963	0.0000	1.0000
K2P	12	-2461.8594	5004.5623	0.0000	1.0000
JC69G	12	-2462.9270	5006.6976	0.0000	1.0000
JC69I	12	-2463.7230	5008.2897	0.0000	1.0000
JC69IG	13	-2462.7410	5013.0626	0.0000	1.0000
SYM	16	-2459.0194	5025.8302	0.0000	1.0000
JC69	11	-2500.6193	5075.3452	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_5279.phylip = 1 - 843;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5279.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_5279.phylip = 1 - 843;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5279.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_5279.phylip = 1 - 843;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5279.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

