
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
01:03 02/17/2011

Input data from file "183-phylip/chrZ_6615.phylip" (29 taxa, 340 characters)

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


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
K2PG	57	-1182.5717	2502.5903	0.3143	0.3143
K2PI	57	-1182.7453	2502.9374	0.2642	0.5786
K2PIG	58	-1181.5937	2503.5432	0.1952	0.7737
HKYG	60	-1179.4217	2505.0800	0.0905	0.8643
HKYI	60	-1179.8159	2505.8684	0.0610	0.9253
HKYIG	61	-1178.4607	2506.1300	0.0535	0.9788
SYMG	61	-1180.6558	2510.5203	0.0060	0.9848
SYMI	61	-1180.9957	2511.2000	0.0042	0.9890
GTRG	64	-1176.5941	2511.4427	0.0038	0.9928
SYMIG	62	-1179.6637	2511.5296	0.0036	0.9964
GTRIG	65	-1175.5530	2512.4198	0.0023	0.9987
GTRI	64	-1177.6530	2513.5606	0.0013	1.0000
JC69G	56	-1202.2371	2539.0326	0.0000	1.0000
JC69I	56	-1202.6361	2539.8305	0.0000	1.0000
JC69IG	57	-1201.8474	2541.1415	0.0000	1.0000
F81G	59	-1199.6286	2542.5429	0.0000	1.0000
F81I	59	-1200.0548	2543.3954	0.0000	1.0000
F81IG	60	-1199.2349	2544.7063	0.0000	1.0000
K2P	56	-1220.5489	2575.6562	0.0000	1.0000
HKY	59	-1218.9307	2581.1472	0.0000	1.0000
SYM	60	-1218.6597	2583.5560	0.0000	1.0000
GTR	63	-1216.4589	2588.1351	0.0000	1.0000
JC69	55	-1243.4715	2618.6332	0.0000	1.0000
F81	58	-1240.4987	2621.3533	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1179.4217	2478.8434	0.1642	0.1642
HKYIG	61	-1178.4607	2478.9214	0.1579	0.3221
K2PG	57	-1182.5717	2479.1435	0.1413	0.4633
K2PIG	58	-1181.5937	2479.1873	0.1382	0.6016
K2PI	57	-1182.7453	2479.4906	0.1188	0.7204
HKYI	60	-1179.8159	2479.6319	0.1107	0.8310
GTRIG	65	-1175.5530	2481.1059	0.0530	0.8840
GTRG	64	-1176.5941	2481.1881	0.0508	0.9348
GTRI	64	-1177.6530	2483.3061	0.0176	0.9525
SYMG	61	-1180.6558	2483.3117	0.0176	0.9700
SYMIG	62	-1179.6637	2483.3274	0.0174	0.9875
SYMI	61	-1180.9957	2483.9914	0.0125	1.0000
JC69G	56	-1202.2371	2516.4743	0.0000	1.0000
F81G	59	-1199.6286	2517.2572	0.0000	1.0000
JC69I	56	-1202.6361	2517.2722	0.0000	1.0000
JC69IG	57	-1201.8474	2517.6947	0.0000	1.0000
F81I	59	-1200.0548	2518.1097	0.0000	1.0000
F81IG	60	-1199.2349	2518.4697	0.0000	1.0000
K2P	56	-1220.5489	2553.0979	0.0000	1.0000
HKY	59	-1218.9307	2555.8615	0.0000	1.0000
SYM	60	-1218.6597	2557.3194	0.0000	1.0000
GTR	63	-1216.4589	2558.9177	0.0000	1.0000
JC69	55	-1243.4715	2596.9431	0.0000	1.0000
F81	58	-1240.4987	2596.9974	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1182.5717	2697.3934	0.5021	0.5021
K2PI	57	-1182.7453	2697.7405	0.4221	0.9241
K2PIG	58	-1181.5937	2701.2662	0.0724	0.9966
HKYG	60	-1179.4217	2708.5802	0.0019	0.9984
HKYI	60	-1179.8159	2709.3686	0.0013	0.9997
HKYIG	61	-1178.4607	2712.4870	0.0003	0.9999
SYMG	61	-1180.6558	2716.8774	0.0000	1.0000
SYMI	61	-1180.9957	2717.5570	0.0000	1.0000
SYMIG	62	-1179.6637	2720.7220	0.0000	1.0000
GTRG	64	-1176.5941	2726.2406	0.0000	1.0000
GTRI	64	-1177.6530	2728.3586	0.0000	1.0000
GTRIG	65	-1175.5530	2729.9874	0.0000	1.0000
JC69G	56	-1202.2371	2730.8952	0.0000	1.0000
JC69I	56	-1202.6361	2731.6932	0.0000	1.0000
JC69IG	57	-1201.8474	2735.9446	0.0000	1.0000
F81G	59	-1199.6286	2743.1650	0.0000	1.0000
F81I	59	-1200.0548	2744.0175	0.0000	1.0000
F81IG	60	-1199.2349	2748.2065	0.0000	1.0000
K2P	56	-1220.5489	2767.5188	0.0000	1.0000
HKY	59	-1218.9307	2781.7693	0.0000	1.0000
SYM	60	-1218.6597	2787.0561	0.0000	1.0000
GTR	63	-1216.4589	2800.1413	0.0000	1.0000
JC69	55	-1243.4715	2807.5351	0.0000	1.0000
F81	58	-1240.4987	2819.0763	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 183-phylip/chrZ_6615.phylip = 1 - 340;
 Partition Dummy = 1:183-phylip/chrZ_6615.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 (K2PG)]

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_6615.phylip = 1 - 340;
 Partition Dummy = 1:183-phylip/chrZ_6615.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_6615.phylip = 1 - 340;
 Partition Dummy = 1:183-phylip/chrZ_6615.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

