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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	65	-1079.3909	2322.6949	0.5553	0.5553
GTRG	64	-1081.2843	2323.3244	0.4054	0.9607
GTRI	64	-1084.1274	2329.0107	0.0236	0.9843
HKYIG	61	-1089.8395	2331.1109	0.0083	0.9926
HKYG	60	-1091.5660	2331.5041	0.0068	0.9994
HKYI	60	-1093.9448	2336.2616	0.0006	1.0000
SYMG	61	-1103.2098	2357.8516	0.0000	1.0000
K2PIG	58	-1108.0898	2358.5027	0.0000	1.0000
K2PG	57	-1109.6530	2358.6393	0.0000	1.0000
SYMIG	62	-1103.3259	2361.1674	0.0000	1.0000
K2PI	57	-1111.8218	2362.9768	0.0000	1.0000
F81G	59	-1110.0849	2365.5057	0.0000	1.0000
SYMI	61	-1107.2486	2365.9291	0.0000	1.0000
F81IG	60	-1109.3350	2367.0421	0.0000	1.0000
F81I	59	-1114.0598	2373.4556	0.0000	1.0000
JC69G	56	-1123.2515	2382.8693	0.0000	1.0000
GTR	63	-1113.1751	2383.9737	0.0000	1.0000
HKY	59	-1120.5287	2386.3934	0.0000	1.0000
JC69IG	57	-1124.8889	2389.1110	0.0000	1.0000
JC69I	56	-1129.3306	2395.0276	0.0000	1.0000
K2P	56	-1139.5556	2415.4776	0.0000	1.0000
SYM	60	-1134.6789	2417.7299	0.0000	1.0000
F81	58	-1139.8798	2422.0827	0.0000	1.0000
JC69	55	-1156.7682	2446.9585	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1079.3909	2288.7819	0.6966	0.6966
GTRG	64	-1081.2843	2290.5685	0.2851	0.9817
GTRI	64	-1084.1274	2296.2548	0.0166	0.9983
HKYIG	61	-1089.8395	2301.6790	0.0011	0.9994
HKYG	60	-1091.5660	2303.1320	0.0005	1.0000
HKYI	60	-1093.9448	2307.8895	0.0000	1.0000
SYMG	61	-1103.2098	2328.4197	0.0000	1.0000
SYMIG	62	-1103.3259	2330.6517	0.0000	1.0000
K2PIG	58	-1108.0898	2332.1797	0.0000	1.0000
K2PG	57	-1109.6530	2333.3060	0.0000	1.0000
SYMI	61	-1107.2486	2336.4972	0.0000	1.0000
K2PI	57	-1111.8218	2337.6435	0.0000	1.0000
F81G	59	-1110.0849	2338.1698	0.0000	1.0000
F81IG	60	-1109.3350	2338.6700	0.0000	1.0000
F81I	59	-1114.0598	2346.1197	0.0000	1.0000
GTR	63	-1113.1751	2352.3501	0.0000	1.0000
JC69G	56	-1123.2515	2358.5029	0.0000	1.0000
HKY	59	-1120.5287	2359.0575	0.0000	1.0000
JC69IG	57	-1124.8889	2363.7777	0.0000	1.0000
JC69I	56	-1129.3306	2370.6611	0.0000	1.0000
SYM	60	-1134.6789	2389.3578	0.0000	1.0000
K2P	56	-1139.5556	2391.1112	0.0000	1.0000
F81	58	-1139.8798	2395.7596	0.0000	1.0000
JC69	55	-1156.7682	2423.5364	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1091.5660	2529.0434	0.5501	0.5501
HKYIG	61	-1089.8395	2531.3556	0.1731	0.7233
GTRG	64	-1081.2843	2531.5408	0.1578	0.8811
GTRIG	65	-1079.3909	2533.5193	0.0587	0.9398
HKYI	60	-1093.9448	2533.8010	0.0510	0.9907
GTRI	64	-1084.1274	2537.2271	0.0092	0.9999
K2PG	57	-1109.6530	2547.9219	0.0000	1.0000
K2PIG	58	-1108.0898	2550.5607	0.0000	1.0000
K2PI	57	-1111.8218	2552.2594	0.0000	1.0000
SYMG	61	-1103.2098	2558.0963	0.0000	1.0000
F81G	59	-1110.0849	2560.3161	0.0000	1.0000
SYMIG	62	-1103.3259	2564.0936	0.0000	1.0000
F81IG	60	-1109.3350	2564.5814	0.0000	1.0000
SYMI	61	-1107.2486	2566.1738	0.0000	1.0000
F81I	59	-1114.0598	2568.2659	0.0000	1.0000
JC69G	56	-1123.2515	2569.3536	0.0000	1.0000
JC69IG	57	-1124.8889	2578.3936	0.0000	1.0000
HKY	59	-1120.5287	2581.2038	0.0000	1.0000
JC69I	56	-1129.3306	2581.5118	0.0000	1.0000
GTR	63	-1113.1751	2589.5572	0.0000	1.0000
K2P	56	-1139.5556	2601.9619	0.0000	1.0000
F81	58	-1139.8798	2614.1407	0.0000	1.0000
SYM	60	-1134.6789	2615.2692	0.0000	1.0000
JC69	55	-1156.7682	2630.6220	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRIG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr2_7954.phylip = 1 - 319;
 Partition Dummy = 1:183-phylip/chr2_7954.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRIG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr2_7954.phylip = 1 - 319;
 Partition Dummy = 1:183-phylip/chr2_7954.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

