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

Input data from file "183-phylip/chr2_25941.phylip" (29 taxa, 317 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	-1142.0017	2432.5972	0.8304	0.8304
HKYIG	61	-1143.1447	2437.9521	0.0571	0.8875
GTRG	64	-1138.6021	2438.2202	0.0499	0.9374
HKYI	60	-1145.3117	2439.2171	0.0303	0.9678
GTRIG	65	-1137.7667	2439.7167	0.0236	0.9914
GTRI	64	-1140.3604	2441.7368	0.0086	1.0000
K2PG	57	-1160.3335	2460.1960	0.0000	1.0000
K2PIG	58	-1161.9238	2466.3748	0.0000	1.0000
SYMG	61	-1158.0334	2467.7295	0.0000	1.0000
K2PI	57	-1164.2592	2468.0474	0.0000	1.0000
SYMIG	62	-1157.1920	2469.1400	0.0000	1.0000
HKY	59	-1161.9226	2469.3939	0.0000	1.0000
SYMI	61	-1159.1391	2469.9409	0.0000	1.0000
GTR	63	-1157.2163	2472.3060	0.0000	1.0000
F81IG	60	-1166.2776	2481.1489	0.0000	1.0000
F81I	59	-1167.8573	2481.2633	0.0000	1.0000
F81G	59	-1168.8514	2483.2515	0.0000	1.0000
K2P	56	-1181.5495	2499.6529	0.0000	1.0000
JC69G	56	-1183.8836	2504.3210	0.0000	1.0000
SYM	60	-1179.0023	2506.5984	0.0000	1.0000
JC69IG	57	-1186.0248	2511.5785	0.0000	1.0000
JC69I	56	-1188.3890	2513.3318	0.0000	1.0000
F81	58	-1189.0282	2520.5835	0.0000	1.0000
JC69	55	-1204.8229	2543.2472	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1142.0017	2404.0034	0.4420	0.4420
GTRG	64	-1138.6021	2405.2043	0.2425	0.6845
GTRIG	65	-1137.7667	2405.5335	0.2057	0.8902
HKYIG	61	-1143.1447	2408.2893	0.0519	0.9421
GTRI	64	-1140.3604	2408.7209	0.0418	0.9839
HKYI	60	-1145.3117	2410.6233	0.0161	1.0000
K2PG	57	-1160.3335	2434.6670	0.0000	1.0000
SYMG	61	-1158.0334	2438.0668	0.0000	1.0000
SYMIG	62	-1157.1920	2438.3841	0.0000	1.0000
K2PIG	58	-1161.9238	2439.8476	0.0000	1.0000
SYMI	61	-1159.1391	2440.2782	0.0000	1.0000
GTR	63	-1157.2163	2440.4325	0.0000	1.0000
HKY	59	-1161.9226	2441.8453	0.0000	1.0000
K2PI	57	-1164.2592	2442.5185	0.0000	1.0000
F81IG	60	-1166.2776	2452.5552	0.0000	1.0000
F81I	59	-1167.8573	2453.7147	0.0000	1.0000
F81G	59	-1168.8514	2455.7029	0.0000	1.0000
K2P	56	-1181.5495	2475.0991	0.0000	1.0000
SYM	60	-1179.0023	2478.0047	0.0000	1.0000
JC69G	56	-1183.8836	2479.7672	0.0000	1.0000
JC69IG	57	-1186.0248	2486.0496	0.0000	1.0000
JC69I	56	-1188.3890	2488.7780	0.0000	1.0000
F81	58	-1189.0282	2494.0564	0.0000	1.0000
JC69	55	-1204.8229	2519.6457	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1142.0017	2629.5375	0.9480	0.9480
HKYI	60	-1145.3117	2636.1574	0.0346	0.9826
HKYIG	61	-1143.1447	2637.5823	0.0170	0.9996
GTRG	64	-1138.6021	2645.7740	0.0003	0.9999
K2PG	57	-1160.3335	2648.9244	0.0001	0.9999
GTRI	64	-1140.3604	2649.2906	0.0000	1.0000
GTRIG	65	-1137.7667	2649.8621	0.0000	1.0000
K2PI	57	-1164.2592	2656.7759	0.0000	1.0000
K2PIG	58	-1161.9238	2657.8639	0.0000	1.0000
HKY	59	-1161.9226	2663.6205	0.0000	1.0000
SYMG	61	-1158.0334	2667.3598	0.0000	1.0000
SYMI	61	-1159.1391	2669.5712	0.0000	1.0000
SYMIG	62	-1157.1920	2671.4360	0.0000	1.0000
F81I	59	-1167.8573	2675.4899	0.0000	1.0000
GTR	63	-1157.2163	2677.2433	0.0000	1.0000
F81G	59	-1168.8514	2677.4781	0.0000	1.0000
F81IG	60	-1166.2776	2678.0893	0.0000	1.0000
K2P	56	-1181.5495	2685.5976	0.0000	1.0000
JC69G	56	-1183.8836	2690.2657	0.0000	1.0000
JC69I	56	-1188.3890	2699.2765	0.0000	1.0000
JC69IG	57	-1186.0248	2700.3070	0.0000	1.0000
SYM	60	-1179.0023	2703.5388	0.0000	1.0000
F81	58	-1189.0282	2712.0727	0.0000	1.0000
JC69	55	-1204.8229	2726.3853	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/chr2_25941.phylip = 1 - 317;
 Partition Dummy = 1:183-phylip/chr2_25941.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/chr2_25941.phylip = 1 - 317;
 Partition Dummy = 1:183-phylip/chr2_25941.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/chr2_25941.phylip = 1 - 317;
 Partition Dummy = 1:183-phylip/chr2_25941.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

