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

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

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


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
GTRG	20	-1383.0574	2807.4785	0.4358	0.4358
GTRI	20	-1383.4272	2808.2180	0.3011	0.7369
GTRIG	21	-1382.9881	2809.4787	0.1603	0.8972
HKYG	16	-1389.4987	2811.8749	0.0484	0.9455
HKYIG	17	-1389.0326	2813.0539	0.0268	0.9724
HKYI	16	-1390.6707	2814.2189	0.0150	0.9873
GTR	19	-1387.7060	2814.6437	0.0121	0.9995
HKY	15	-1395.0520	2820.8770	0.0005	1.0000
SYMG	17	-1403.1863	2841.3613	0.0000	1.0000
SYMI	17	-1403.5016	2841.9918	0.0000	1.0000
SYMIG	18	-1403.1261	2843.3590	0.0000	1.0000
K2PG	13	-1409.2519	2845.0880	0.0000	1.0000
K2PI	13	-1409.4035	2845.3913	0.0000	1.0000
K2PIG	14	-1409.2061	2847.0875	0.0000	1.0000
SYM	16	-1407.6703	2848.2181	0.0000	1.0000
K2P	12	-1414.2118	2852.9237	0.0000	1.0000
F81G	15	-1417.0093	2864.7915	0.0000	1.0000
F81I	15	-1417.2395	2865.2519	0.0000	1.0000
F81IG	16	-1416.9727	2866.8228	0.0000	1.0000
F81	14	-1421.7578	2872.1908	0.0000	1.0000
JC69G	12	-1435.0124	2894.5248	0.0000	1.0000
JC69I	12	-1435.2385	2894.9771	0.0000	1.0000
JC69IG	13	-1434.9783	2896.5409	0.0000	1.0000
JC69	11	-1439.6862	2901.7948	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1383.0574	2806.1148	0.4391	0.4391
GTRI	20	-1383.4272	2806.8544	0.3034	0.7425
GTRIG	21	-1382.9881	2807.9762	0.1731	0.9157
HKYG	16	-1389.4987	2810.9975	0.0382	0.9539
HKYIG	17	-1389.0326	2812.0652	0.0224	0.9763
HKYI	16	-1390.6707	2813.3415	0.0118	0.9882
GTR	19	-1387.7060	2813.4119	0.0114	0.9996
HKY	15	-1395.0520	2820.1041	0.0004	1.0000
SYMG	17	-1403.1863	2840.3726	0.0000	1.0000
SYMI	17	-1403.5016	2841.0031	0.0000	1.0000
SYMIG	18	-1403.1261	2842.2522	0.0000	1.0000
K2PG	13	-1409.2519	2844.5038	0.0000	1.0000
K2PI	13	-1409.4035	2844.8070	0.0000	1.0000
K2PIG	14	-1409.2061	2846.4123	0.0000	1.0000
SYM	16	-1407.6703	2847.3406	0.0000	1.0000
K2P	12	-1414.2118	2852.4237	0.0000	1.0000
F81G	15	-1417.0093	2864.0186	0.0000	1.0000
F81I	15	-1417.2395	2864.4789	0.0000	1.0000
F81IG	16	-1416.9727	2865.9454	0.0000	1.0000
F81	14	-1421.7578	2871.5155	0.0000	1.0000
JC69G	12	-1435.0124	2894.0248	0.0000	1.0000
JC69I	12	-1435.2385	2894.4771	0.0000	1.0000
JC69IG	13	-1434.9783	2895.9566	0.0000	1.0000
JC69	11	-1439.6862	2901.3724	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1389.4987	2882.3058	0.6785	0.6785
HKYI	16	-1390.6707	2884.6498	0.2102	0.8887
HKY	15	-1395.0520	2886.9556	0.0664	0.9550
HKYIG	17	-1389.0326	2887.8303	0.0428	0.9979
GTRG	20	-1383.0574	2895.2502	0.0010	0.9989
GTRI	20	-1383.4272	2895.9898	0.0007	0.9996
GTR	19	-1387.7060	2898.0905	0.0003	0.9999
GTRIG	21	-1382.9881	2901.5684	0.0000	0.9999
K2PG	13	-1409.2519	2902.4418	0.0000	1.0000
K2PI	13	-1409.4035	2902.7450	0.0000	1.0000
K2P	12	-1414.2118	2905.9049	0.0000	1.0000
K2PIG	14	-1409.2061	2908.8071	0.0000	1.0000
SYMG	17	-1403.1863	2916.1377	0.0000	1.0000
SYMI	17	-1403.5016	2916.7682	0.0000	1.0000
SYM	16	-1407.6703	2918.6490	0.0000	1.0000
SYMIG	18	-1403.1261	2922.4741	0.0000	1.0000
F81G	15	-1417.0093	2930.8701	0.0000	1.0000
F81I	15	-1417.2395	2931.3305	0.0000	1.0000
F81	14	-1421.7578	2933.9103	0.0000	1.0000
F81IG	16	-1416.9727	2937.2537	0.0000	1.0000
JC69G	12	-1435.0124	2947.5060	0.0000	1.0000
JC69I	12	-1435.2385	2947.9583	0.0000	1.0000
JC69	11	-1439.6862	2950.3968	0.0000	1.0000
JC69IG	13	-1434.9783	2953.8946	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr17_3088.phylip = 1 - 637;
 Partition Dummy = 1:591-tomodel-phylip/chr17_3088.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr17_3088.phylip = 1 - 637;
 Partition Dummy = 1:591-tomodel-phylip/chr17_3088.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr17_3088.phylip = 1 - 637;
 Partition Dummy = 1:591-tomodel-phylip/chr17_3088.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

