
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:24 08/01/2011

Input data from file "261-aggressive-phylip/4561_hg18.phylip" (29 taxa, 455 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-1280.8222	2710.9777	0.4205	0.4205
GTRIG	65	-1279.9659	2711.9884	0.2537	0.6741
GTRI	64	-1281.7022	2712.7378	0.1744	0.8485
HKYG	60	-1287.8947	2714.3681	0.0772	0.9257
HKYI	60	-1288.5498	2715.6782	0.0401	0.9658
HKYIG	61	-1287.3738	2715.9944	0.0342	1.0000
K2PG	57	-1302.0200	2734.6949	0.0000	1.0000
K2PI	57	-1302.6835	2736.0220	0.0000	1.0000
K2PIG	58	-1301.5405	2736.3638	0.0000	1.0000
SYMG	61	-1298.8526	2738.9521	0.0000	1.0000
SYMI	61	-1299.5557	2740.3582	0.0000	1.0000
SYMIG	62	-1298.3142	2740.5571	0.0000	1.0000
GTR	63	-1303.3943	2753.4126	0.0000	1.0000
HKY	59	-1308.9655	2753.8550	0.0000	1.0000
F81G	59	-1311.6226	2759.1693	0.0000	1.0000
F81I	59	-1312.2885	2760.5010	0.0000	1.0000
F81IG	60	-1311.1720	2760.9227	0.0000	1.0000
K2P	56	-1322.3454	2772.7309	0.0000	1.0000
JC69G	56	-1323.0362	2774.1126	0.0000	1.0000
JC69I	56	-1323.7170	2775.4743	0.0000	1.0000
JC69IG	57	-1322.5844	2775.8236	0.0000	1.0000
SYM	60	-1319.2409	2777.0605	0.0000	1.0000
F81	58	-1331.6885	2796.6599	0.0000	1.0000
JC69	55	-1342.9434	2811.3254	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1280.8222	2689.6443	0.4202	0.4202
GTRIG	65	-1279.9659	2689.9319	0.3639	0.7841
GTRI	64	-1281.7022	2691.4044	0.1743	0.9584
HKYG	60	-1287.8947	2695.7894	0.0195	0.9778
HKYIG	61	-1287.3738	2696.7476	0.0121	0.9899
HKYI	60	-1288.5498	2697.0996	0.0101	1.0000
K2PG	57	-1302.0200	2718.0400	0.0000	1.0000
K2PIG	58	-1301.5405	2719.0810	0.0000	1.0000
K2PI	57	-1302.6835	2719.3671	0.0000	1.0000
SYMG	61	-1298.8526	2719.7053	0.0000	1.0000
SYMIG	62	-1298.3142	2720.6285	0.0000	1.0000
SYMI	61	-1299.5557	2721.1113	0.0000	1.0000
GTR	63	-1303.3943	2732.7885	0.0000	1.0000
HKY	59	-1308.9655	2735.9309	0.0000	1.0000
F81G	59	-1311.6226	2741.2452	0.0000	1.0000
F81IG	60	-1311.1720	2742.3440	0.0000	1.0000
F81I	59	-1312.2885	2742.5769	0.0000	1.0000
K2P	56	-1322.3454	2756.6907	0.0000	1.0000
JC69G	56	-1323.0362	2758.0724	0.0000	1.0000
SYM	60	-1319.2409	2758.4818	0.0000	1.0000
JC69IG	57	-1322.5844	2759.1687	0.0000	1.0000
JC69I	56	-1323.7170	2759.4341	0.0000	1.0000
F81	58	-1331.6885	2779.3771	0.0000	1.0000
JC69	55	-1342.9434	2795.8868	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1287.8947	2943.0072	0.6179	0.6179
HKYI	60	-1288.5498	2944.3174	0.3209	0.9389
HKYIG	61	-1287.3738	2948.0857	0.0488	0.9876
K2PG	57	-1302.0200	2952.8969	0.0044	0.9920
GTRG	64	-1280.8222	2953.3434	0.0035	0.9956
K2PI	57	-1302.6835	2954.2241	0.0023	0.9978
GTRI	64	-1281.7022	2955.1035	0.0015	0.9993
GTRIG	65	-1279.9659	2957.7512	0.0004	0.9997
K2PIG	58	-1301.5405	2958.0583	0.0003	1.0000
SYMG	61	-1298.8526	2971.0434	0.0000	1.0000
SYMI	61	-1299.5557	2972.4495	0.0000	1.0000
SYMIG	62	-1298.3142	2976.0869	0.0000	1.0000
HKY	59	-1308.9655	2979.0285	0.0000	1.0000
F81G	59	-1311.6226	2984.3428	0.0000	1.0000
F81I	59	-1312.2885	2985.6744	0.0000	1.0000
K2P	56	-1322.3454	2987.4274	0.0000	1.0000
JC69G	56	-1323.0362	2988.8091	0.0000	1.0000
F81IG	60	-1311.1720	2989.5619	0.0000	1.0000
JC69I	56	-1323.7170	2990.1707	0.0000	1.0000
GTR	63	-1303.3943	2992.3673	0.0000	1.0000
JC69IG	57	-1322.5844	2994.0257	0.0000	1.0000
SYM	60	-1319.2409	3005.6996	0.0000	1.0000
F81	58	-1331.6885	3018.3543	0.0000	1.0000
JC69	55	-1342.9434	3022.5031	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 261-aggressive-phylip/4561_hg18.phylip = 1 - 455;
 Partition Dummy = 1:261-aggressive-phylip/4561_hg18.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 (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/4561_hg18.phylip = 1 - 455;
 Partition Dummy = 1:261-aggressive-phylip/4561_hg18.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 BIC model (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/4561_hg18.phylip = 1 - 455;
 Partition Dummy = 1:261-aggressive-phylip/4561_hg18.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

