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

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

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
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
HKYG	60	-1088.7371	2321.8741	0.3607	0.3607
SYMG	61	-1087.8236	2322.9449	0.2112	0.5719
K2PG	57	-1093.5763	2322.9743	0.2081	0.7800
GTRG	64	-1084.6269	2325.3618	0.0631	0.8431
HKYI	60	-1090.8198	2326.0397	0.0449	0.8880
SYMI	61	-1089.6519	2326.6015	0.0339	0.9220
K2PI	57	-1095.6790	2327.1799	0.0254	0.9474
HKYIG	61	-1090.5215	2328.3406	0.0142	0.9616
SYMIG	62	-1089.1130	2328.4408	0.0135	0.9751
GTRI	64	-1086.3171	2328.7423	0.0116	0.9868
K2PIG	58	-1095.2614	2329.1850	0.0093	0.9961
GTRIG	65	-1085.9216	2330.9280	0.0039	1.0000
JC69G	56	-1113.4159	2359.8319	0.0000	1.0000
JC69I	56	-1113.5085	2360.0169	0.0000	1.0000
F81G	59	-1109.6355	2360.7926	0.0000	1.0000
F81I	59	-1109.7128	2360.9471	0.0000	1.0000
JC69IG	57	-1113.0309	2361.8835	0.0000	1.0000
F81IG	60	-1109.2837	2362.9674	0.0000	1.0000
K2P	56	-1127.5304	2388.0608	0.0000	1.0000
HKY	59	-1123.7174	2388.9564	0.0000	1.0000
SYM	60	-1122.4898	2389.3797	0.0000	1.0000
GTR	63	-1119.9551	2393.0618	0.0000	1.0000
JC69	55	-1147.6892	2425.5751	0.0000	1.0000
F81	58	-1143.9036	2426.4694	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1084.6269	2297.2537	0.2814	0.2814
HKYG	60	-1088.7371	2297.4741	0.2520	0.5334
SYMG	61	-1087.8236	2297.6473	0.2311	0.7645
GTRI	64	-1086.3171	2300.6342	0.0519	0.8164
K2PG	57	-1093.5763	2301.1525	0.0401	0.8565
SYMI	61	-1089.6519	2301.3038	0.0371	0.8936
HKYI	60	-1090.8198	2301.6397	0.0314	0.9250
GTRIG	65	-1085.9216	2301.8432	0.0284	0.9534
SYMIG	62	-1089.1130	2302.2261	0.0234	0.9768
HKYIG	61	-1090.5215	2303.0429	0.0156	0.9924
K2PI	57	-1095.6790	2305.3581	0.0049	0.9973
K2PIG	58	-1095.2614	2306.5228	0.0027	1.0000
F81G	59	-1109.6355	2337.2710	0.0000	1.0000
F81I	59	-1109.7128	2337.4255	0.0000	1.0000
F81IG	60	-1109.2837	2338.5674	0.0000	1.0000
JC69G	56	-1113.4159	2338.8319	0.0000	1.0000
JC69I	56	-1113.5085	2339.0169	0.0000	1.0000
JC69IG	57	-1113.0309	2340.0617	0.0000	1.0000
SYM	60	-1122.4898	2364.9797	0.0000	1.0000
HKY	59	-1123.7174	2365.4348	0.0000	1.0000
GTR	63	-1119.9551	2365.9103	0.0000	1.0000
K2P	56	-1127.5304	2367.0608	0.0000	1.0000
F81	58	-1143.9036	2403.8071	0.0000	1.0000
JC69	55	-1147.6892	2405.3784	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1093.5763	2522.8186	0.8653	0.8653
K2PI	57	-1095.6790	2527.0241	0.1057	0.9710
HKYG	60	-1088.7371	2530.8068	0.0159	0.9870
K2PIG	58	-1095.2614	2532.0777	0.0084	0.9954
SYMG	61	-1087.8236	2534.8688	0.0021	0.9975
HKYI	60	-1090.8198	2534.9724	0.0020	0.9995
SYMI	61	-1089.6519	2538.5254	0.0003	0.9998
HKYIG	61	-1090.5215	2540.2645	0.0001	1.0000
SYMIG	62	-1089.1130	2543.3365	0.0000	1.0000
GTRG	64	-1084.6269	2546.1419	0.0000	1.0000
GTRI	64	-1086.3171	2549.5224	0.0000	1.0000
GTRIG	65	-1085.9216	2554.6203	0.0000	1.0000
JC69G	56	-1113.4159	2556.6090	0.0000	1.0000
JC69I	56	-1113.5085	2556.7941	0.0000	1.0000
JC69IG	57	-1113.0309	2561.7278	0.0000	1.0000
F81G	59	-1109.6355	2566.7148	0.0000	1.0000
F81I	59	-1109.7128	2566.8693	0.0000	1.0000
F81IG	60	-1109.2837	2571.9000	0.0000	1.0000
K2P	56	-1127.5304	2584.8379	0.0000	1.0000
HKY	59	-1123.7174	2594.8786	0.0000	1.0000
SYM	60	-1122.4898	2598.3124	0.0000	1.0000
GTR	63	-1119.9551	2610.9096	0.0000	1.0000
JC69	55	-1147.6892	2619.2667	0.0000	1.0000
F81	58	-1143.9036	2629.3620	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 183-phylip/chr2_22128.phylip = 1 - 361;
 Partition Dummy = 1:183-phylip/chr2_22128.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 (HKYG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr2_22128.phylip = 1 - 361;
 Partition Dummy = 1:183-phylip/chr2_22128.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 (K2PG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr2_22128.phylip = 1 - 361;
 Partition Dummy = 1:183-phylip/chr2_22128.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

