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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-1225.2034	2593.2106	0.8534	0.8534
HKYI	60	-1227.6834	2598.1704	0.0715	0.9249
HKYIG	61	-1226.3028	2598.2431	0.0689	0.9938
GTRG	64	-1224.5947	2603.4354	0.0051	0.9990
GTRI	64	-1226.9172	2608.0805	0.0005	0.9995
GTRIG	65	-1225.4958	2608.1436	0.0005	1.0000
K2PG	57	-1240.3669	2615.1412	0.0000	1.0000
K2PI	57	-1240.9239	2616.2552	0.0000	1.0000
K2PIG	58	-1240.0011	2617.1910	0.0000	1.0000
SYMG	61	-1239.6358	2624.9091	0.0000	1.0000
SYMI	61	-1240.0941	2625.8257	0.0000	1.0000
SYMIG	62	-1239.2360	2626.9610	0.0000	1.0000
HKY	59	-1251.6207	2643.2289	0.0000	1.0000
GTR	63	-1250.9834	2653.3253	0.0000	1.0000
K2P	56	-1268.0857	2667.8145	0.0000	1.0000
F81G	59	-1264.5583	2669.1042	0.0000	1.0000
F81IG	60	-1265.2268	2673.2573	0.0000	1.0000
F81I	59	-1266.9772	2673.9420	0.0000	1.0000
SYM	60	-1267.7608	2678.3253	0.0000	1.0000
JC69G	56	-1277.3770	2686.3971	0.0000	1.0000
JC69IG	57	-1278.0636	2690.5347	0.0000	1.0000
JC69I	56	-1279.7525	2691.1481	0.0000	1.0000
F81	58	-1291.4330	2720.0549	0.0000	1.0000
JC69	55	-1304.8985	2738.6926	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1225.2034	2570.4068	0.8011	0.8011
HKYIG	61	-1226.3028	2574.6056	0.0982	0.8993
HKYI	60	-1227.6834	2575.3667	0.0671	0.9663
GTRG	64	-1224.5947	2577.1894	0.0270	0.9933
GTRIG	65	-1225.4958	2580.9917	0.0040	0.9973
GTRI	64	-1226.9172	2581.8344	0.0026	1.0000
K2PG	57	-1240.3669	2594.7338	0.0000	1.0000
K2PI	57	-1240.9239	2595.8478	0.0000	1.0000
K2PIG	58	-1240.0011	2596.0022	0.0000	1.0000
SYMG	61	-1239.6358	2601.2716	0.0000	1.0000
SYMI	61	-1240.0941	2602.1882	0.0000	1.0000
SYMIG	62	-1239.2360	2602.4720	0.0000	1.0000
HKY	59	-1251.6207	2621.2413	0.0000	1.0000
GTR	63	-1250.9834	2627.9668	0.0000	1.0000
F81G	59	-1264.5583	2647.1167	0.0000	1.0000
K2P	56	-1268.0857	2648.1715	0.0000	1.0000
F81IG	60	-1265.2268	2650.4535	0.0000	1.0000
F81I	59	-1266.9772	2651.9545	0.0000	1.0000
SYM	60	-1267.7608	2655.5215	0.0000	1.0000
JC69G	56	-1277.3770	2666.7540	0.0000	1.0000
JC69IG	57	-1278.0636	2670.1273	0.0000	1.0000
JC69I	56	-1279.7525	2671.5050	0.0000	1.0000
F81	58	-1291.4330	2698.8660	0.0000	1.0000
JC69	55	-1304.8985	2719.7969	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1225.2034	2807.1321	0.9058	0.9058
HKYI	60	-1227.6834	2812.0919	0.0759	0.9817
HKYIG	61	-1226.3028	2815.2762	0.0154	0.9971
K2PG	57	-1240.3669	2819.6228	0.0018	0.9989
K2PI	57	-1240.9239	2820.7368	0.0010	0.9999
K2PIG	58	-1240.0011	2824.8366	0.0001	1.0000
GTRG	64	-1224.5947	2829.6963	0.0000	1.0000
GTRI	64	-1226.9172	2834.3413	0.0000	1.0000
GTRIG	65	-1225.4958	2837.4440	0.0000	1.0000
SYMG	61	-1239.6358	2841.9423	0.0000	1.0000
SYMI	61	-1240.0941	2842.8589	0.0000	1.0000
SYMIG	62	-1239.2360	2847.0880	0.0000	1.0000
HKY	59	-1251.6207	2854.0212	0.0000	1.0000
K2P	56	-1268.0857	2869.1150	0.0000	1.0000
GTR	63	-1250.9834	2876.5283	0.0000	1.0000
F81G	59	-1264.5583	2879.8965	0.0000	1.0000
F81I	59	-1266.9772	2884.7343	0.0000	1.0000
F81IG	60	-1265.2268	2887.1788	0.0000	1.0000
JC69G	56	-1277.3770	2887.6976	0.0000	1.0000
SYM	60	-1267.7608	2892.2468	0.0000	1.0000
JC69I	56	-1279.7525	2892.4485	0.0000	1.0000
JC69IG	57	-1278.0636	2895.0162	0.0000	1.0000
F81	58	-1291.4330	2927.7004	0.0000	1.0000
JC69	55	-1304.8985	2936.7951	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/chr6_3228.phylip = 1 - 382;
 Partition Dummy = 1:183-phylip/chr6_3228.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 183-phylip/chr6_3228.phylip = 1 - 382;
 Partition Dummy = 1:183-phylip/chr6_3228.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 183-phylip/chr6_3228.phylip = 1 - 382;
 Partition Dummy = 1:183-phylip/chr6_3228.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

