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

Input data from file "261-aggressive-phylip/7217_hg18.phylip" (29 taxa, 512 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	-1712.3789	3571.3708	0.9725	0.9725
GTRIG	65	-1715.0784	3579.3944	0.0176	0.9901
HKYG	60	-1722.4096	3581.0498	0.0077	0.9978
GTRI	64	-1718.7839	3584.1808	0.0016	0.9994
HKYI	60	-1725.2907	3586.8120	0.0004	0.9999
SYMG	61	-1725.4489	3589.7067	0.0001	1.0000
HKYIG	61	-1726.8866	3592.5822	0.0000	1.0000
K2PG	57	-1733.8131	3596.1900	0.0000	1.0000
SYMIG	62	-1729.5192	3600.4371	0.0000	1.0000
K2PI	57	-1736.2804	3601.1247	0.0000	1.0000
SYMI	61	-1732.5646	3603.9381	0.0000	1.0000
K2PIG	58	-1738.4393	3607.9868	0.0000	1.0000
GTR	63	-1737.8480	3619.6960	0.0000	1.0000
HKY	59	-1750.1274	3633.9185	0.0000	1.0000
SYM	60	-1751.1697	3638.5701	0.0000	1.0000
K2P	56	-1761.1816	3648.3940	0.0000	1.0000
F81G	59	-1779.4627	3692.5890	0.0000	1.0000
F81IG	60	-1778.8681	3693.9668	0.0000	1.0000
F81I	59	-1781.9783	3697.6204	0.0000	1.0000
JC69G	56	-1787.6965	3701.4237	0.0000	1.0000
JC69IG	57	-1787.1303	3702.8244	0.0000	1.0000
JC69I	56	-1790.0878	3706.2064	0.0000	1.0000
F81	58	-1805.3808	3741.8697	0.0000	1.0000
JC69	55	-1813.5469	3750.6027	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1712.3789	3552.7578	0.9718	0.9718
GTRIG	65	-1715.0784	3560.1567	0.0240	0.9959
HKYG	60	-1722.4096	3564.8192	0.0023	0.9982
GTRI	64	-1718.7839	3565.5678	0.0016	0.9998
HKYI	60	-1725.2907	3570.5814	0.0001	0.9999
SYMG	61	-1725.4489	3572.8978	0.0000	1.0000
HKYIG	61	-1726.8866	3575.7733	0.0000	1.0000
K2PG	57	-1733.8131	3581.6262	0.0000	1.0000
SYMIG	62	-1729.5192	3583.0385	0.0000	1.0000
K2PI	57	-1736.2804	3586.5608	0.0000	1.0000
SYMI	61	-1732.5646	3587.1293	0.0000	1.0000
K2PIG	58	-1738.4393	3592.8786	0.0000	1.0000
GTR	63	-1737.8480	3601.6960	0.0000	1.0000
HKY	59	-1750.1274	3618.2548	0.0000	1.0000
SYM	60	-1751.1697	3622.3395	0.0000	1.0000
K2P	56	-1761.1816	3634.3632	0.0000	1.0000
F81G	59	-1779.4627	3676.9253	0.0000	1.0000
F81IG	60	-1778.8681	3677.7362	0.0000	1.0000
F81I	59	-1781.9783	3681.9567	0.0000	1.0000
JC69G	56	-1787.6965	3687.3930	0.0000	1.0000
JC69IG	57	-1787.1303	3688.2606	0.0000	1.0000
JC69I	56	-1790.0878	3692.1756	0.0000	1.0000
F81	58	-1805.3808	3726.7616	0.0000	1.0000
JC69	55	-1813.5469	3737.0939	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1722.4096	3819.1187	0.7776	0.7776
K2PG	57	-1733.8131	3823.2107	0.1005	0.8781
GTRG	64	-1712.3789	3824.0106	0.0674	0.9455
HKYI	60	-1725.2907	3824.8808	0.0436	0.9891
K2PI	57	-1736.2804	3828.1453	0.0085	0.9976
SYMG	61	-1725.4489	3831.4356	0.0016	0.9993
HKYIG	61	-1726.8866	3834.3111	0.0004	0.9996
GTRIG	65	-1715.0784	3835.6478	0.0002	0.9998
GTRI	64	-1718.7839	3836.8206	0.0001	1.0000
K2PIG	58	-1738.4393	3838.7014	0.0000	1.0000
SYMI	61	-1732.5646	3845.6671	0.0000	1.0000
SYMIG	62	-1729.5192	3845.8146	0.0000	1.0000
HKY	59	-1750.1274	3868.3160	0.0000	1.0000
GTR	63	-1737.8480	3868.7105	0.0000	1.0000
K2P	56	-1761.1816	3871.7094	0.0000	1.0000
SYM	60	-1751.1697	3876.6390	0.0000	1.0000
JC69G	56	-1787.6965	3924.7391	0.0000	1.0000
F81G	59	-1779.4627	3926.9865	0.0000	1.0000
JC69I	56	-1790.0878	3929.5218	0.0000	1.0000
JC69IG	57	-1787.1303	3929.8451	0.0000	1.0000
F81I	59	-1781.9783	3932.0178	0.0000	1.0000
F81IG	60	-1778.8681	3932.0357	0.0000	1.0000
JC69	55	-1813.5469	3970.2018	0.0000	1.0000
F81	58	-1805.3808	3972.5844	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/7217_hg18.phylip = 1 - 512;
 Partition Dummy = 1:261-aggressive-phylip/7217_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/7217_hg18.phylip = 1 - 512;
 Partition Dummy = 1:261-aggressive-phylip/7217_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/7217_hg18.phylip = 1 - 512;
 Partition Dummy = 1:261-aggressive-phylip/7217_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

