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

Input data from file "183-phylip/chr6_8786.phylip" (29 taxa, 372 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	-1233.4531	2610.4432	0.4879	0.4879
HKYI	60	-1233.8461	2611.2291	0.3294	0.8173
HKYIG	61	-1233.4627	2613.3254	0.1155	0.9327
GTRG	64	-1230.3368	2615.7746	0.0339	0.9667
GTRI	64	-1230.8826	2616.8662	0.0197	0.9863
GTRIG	65	-1230.3262	2618.6916	0.0079	0.9942
HKY	59	-1239.5090	2619.7103	0.0047	0.9990
F81G	59	-1241.8100	2624.3123	0.0005	0.9994
F81I	59	-1242.1888	2625.0699	0.0003	0.9998
GTR	63	-1237.4330	2627.0479	0.0001	0.9999
F81IG	60	-1241.8221	2627.1811	0.0001	1.0000
F81	58	-1247.5463	2632.9584	0.0000	1.0000
K2PG	57	-1255.6148	2646.2869	0.0000	1.0000
K2PI	57	-1255.9388	2646.9350	0.0000	1.0000
K2PIG	58	-1255.6189	2649.1037	0.0000	1.0000
K2P	56	-1260.9789	2654.2244	0.0000	1.0000
JC69G	56	-1261.4331	2655.1328	0.0000	1.0000
JC69I	56	-1261.7728	2655.8123	0.0000	1.0000
SYMG	61	-1255.3749	2657.1498	0.0000	1.0000
SYMI	61	-1255.6941	2657.7882	0.0000	1.0000
JC69IG	57	-1261.4348	2657.9269	0.0000	1.0000
SYMIG	62	-1255.3836	2660.0488	0.0000	1.0000
JC69	55	-1266.7861	2663.0660	0.0000	1.0000
SYM	60	-1260.6408	2664.8186	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1233.4531	2586.9062	0.3503	0.3503
HKYI	60	-1233.8461	2587.6922	0.2365	0.5868
GTRG	64	-1230.3368	2588.6736	0.1448	0.7316
HKYIG	61	-1233.4627	2588.9254	0.1277	0.8593
GTRI	64	-1230.8826	2589.7652	0.0839	0.9432
GTRIG	65	-1230.3262	2590.6524	0.0538	0.9970
HKY	59	-1239.5090	2597.0180	0.0022	0.9992
GTR	63	-1237.4330	2600.8661	0.0003	0.9995
F81G	59	-1241.8100	2601.6200	0.0002	0.9998
F81I	59	-1242.1888	2602.3776	0.0002	0.9999
F81IG	60	-1241.8221	2603.6441	0.0001	1.0000
F81	58	-1247.5463	2611.0926	0.0000	1.0000
K2PG	57	-1255.6148	2625.2296	0.0000	1.0000
K2PI	57	-1255.9388	2625.8776	0.0000	1.0000
K2PIG	58	-1255.6189	2627.2378	0.0000	1.0000
SYMG	61	-1255.3749	2632.7498	0.0000	1.0000
SYMI	61	-1255.6941	2633.3882	0.0000	1.0000
K2P	56	-1260.9789	2633.9578	0.0000	1.0000
SYMIG	62	-1255.3836	2634.7672	0.0000	1.0000
JC69G	56	-1261.4331	2634.8662	0.0000	1.0000
JC69I	56	-1261.7728	2635.5456	0.0000	1.0000
JC69IG	57	-1261.4348	2636.8696	0.0000	1.0000
SYM	60	-1260.6408	2641.2816	0.0000	1.0000
JC69	55	-1266.7861	2643.5723	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1233.4531	2822.0398	0.5618	0.5618
HKYI	60	-1233.8461	2822.8258	0.3792	0.9410
HKYIG	61	-1233.4627	2827.9779	0.0288	0.9699
HKY	59	-1239.5090	2828.2328	0.0254	0.9953
F81G	59	-1241.8100	2832.8347	0.0025	0.9978
F81I	59	-1242.1888	2833.5923	0.0017	0.9996
F81	58	-1247.5463	2838.3884	0.0002	0.9997
F81IG	60	-1241.8221	2838.7778	0.0001	0.9998
GTRG	64	-1230.3368	2839.4828	0.0001	0.9999
GTRI	64	-1230.8826	2840.5744	0.0001	1.0000
GTRIG	65	-1230.3262	2845.3805	0.0000	1.0000
GTR	63	-1237.4330	2847.7564	0.0000	1.0000
K2PG	57	-1255.6148	2848.6065	0.0000	1.0000
K2PI	57	-1255.9388	2849.2546	0.0000	1.0000
K2P	56	-1260.9789	2853.4158	0.0000	1.0000
JC69G	56	-1261.4331	2854.3242	0.0000	1.0000
K2PIG	58	-1255.6189	2854.5337	0.0000	1.0000
JC69I	56	-1261.7728	2855.0037	0.0000	1.0000
JC69	55	-1266.7861	2859.1115	0.0000	1.0000
JC69IG	57	-1261.4348	2860.2465	0.0000	1.0000
SYMG	61	-1255.3749	2871.8023	0.0000	1.0000
SYMI	61	-1255.6941	2872.4407	0.0000	1.0000
SYM	60	-1260.6408	2876.4152	0.0000	1.0000
SYMIG	62	-1255.3836	2877.7387	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_8786.phylip = 1 - 372;
 Partition Dummy = 1:183-phylip/chr6_8786.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_8786.phylip = 1 - 372;
 Partition Dummy = 1:183-phylip/chr6_8786.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_8786.phylip = 1 - 372;
 Partition Dummy = 1:183-phylip/chr6_8786.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

