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

Input data from file "261-aggressive-phylip/7659_hg18.phylip" (29 taxa, 615 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	-1748.0647	3639.2567	0.5145	0.5145
GTRI	64	-1748.7939	3640.7151	0.2482	0.7627
GTRIG	65	-1747.6508	3640.9301	0.2229	0.9856
HKYG	60	-1757.2521	3647.7172	0.0075	0.9930
HKYIG	61	-1756.7408	3649.1597	0.0036	0.9967
HKYI	60	-1758.0665	3649.3461	0.0033	1.0000
SYMG	61	-1768.3446	3672.3672	0.0000	1.0000
SYMIG	62	-1767.9195	3673.9911	0.0000	1.0000
K2PG	57	-1774.4219	3674.7146	0.0000	1.0000
SYMI	61	-1769.7609	3675.2000	0.0000	1.0000
K2PIG	58	-1773.8400	3675.9893	0.0000	1.0000
K2PI	57	-1775.2734	3676.4175	0.0000	1.0000
GTR	63	-1770.9937	3682.6226	0.0000	1.0000
HKY	59	-1779.6498	3690.0563	0.0000	1.0000
SYM	60	-1791.0573	3715.3275	0.0000	1.0000
F81G	59	-1792.2929	3715.3426	0.0000	1.0000
F81IG	60	-1791.7803	3716.7736	0.0000	1.0000
F81I	59	-1793.1005	3716.9578	0.0000	1.0000
K2P	56	-1797.3608	3718.1625	0.0000	1.0000
JC69G	56	-1807.1624	3737.7657	0.0000	1.0000
JC69IG	57	-1806.5129	3738.8965	0.0000	1.0000
JC69I	56	-1808.0134	3739.4678	0.0000	1.0000
F81	58	-1814.6729	3757.6551	0.0000	1.0000
JC69	55	-1829.9950	3781.0096	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1748.0647	3624.1294	0.4877	0.4877
GTRIG	65	-1747.6508	3625.3017	0.2714	0.7592
GTRI	64	-1748.7939	3625.5878	0.2352	0.9944
HKYG	60	-1757.2521	3634.5042	0.0027	0.9971
HKYIG	61	-1756.7408	3635.4816	0.0017	0.9988
HKYI	60	-1758.0665	3636.1331	0.0012	1.0000
SYMG	61	-1768.3446	3658.6891	0.0000	1.0000
SYMIG	62	-1767.9195	3659.8389	0.0000	1.0000
SYMI	61	-1769.7609	3661.5218	0.0000	1.0000
K2PG	57	-1774.4219	3662.8439	0.0000	1.0000
K2PIG	58	-1773.8400	3663.6799	0.0000	1.0000
K2PI	57	-1775.2734	3664.5468	0.0000	1.0000
GTR	63	-1770.9937	3667.9874	0.0000	1.0000
HKY	59	-1779.6498	3677.2995	0.0000	1.0000
SYM	60	-1791.0573	3702.1145	0.0000	1.0000
F81G	59	-1792.2929	3702.5859	0.0000	1.0000
F81IG	60	-1791.7803	3703.5606	0.0000	1.0000
F81I	59	-1793.1005	3704.2010	0.0000	1.0000
K2P	56	-1797.3608	3706.7216	0.0000	1.0000
JC69G	56	-1807.1624	3726.3248	0.0000	1.0000
JC69IG	57	-1806.5129	3727.0258	0.0000	1.0000
JC69I	56	-1808.0134	3728.0269	0.0000	1.0000
F81	58	-1814.6729	3745.3458	0.0000	1.0000
JC69	55	-1829.9950	3769.9899	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1757.2521	3899.8016	0.6448	0.6448
HKYI	60	-1758.0665	3901.4304	0.2856	0.9304
HKYIG	61	-1756.7408	3905.2005	0.0434	0.9738
GTRG	64	-1748.0647	3907.1132	0.0167	0.9904
GTRI	64	-1748.7939	3908.5717	0.0080	0.9985
GTRIG	65	-1747.6508	3912.7071	0.0010	0.9995
K2PG	57	-1774.4219	3914.8764	0.0003	0.9998
K2PI	57	-1775.2734	3916.5792	0.0001	1.0000
K2PIG	58	-1773.8400	3920.1340	0.0000	1.0000
SYMG	61	-1768.3446	3928.4081	0.0000	1.0000
SYMI	61	-1769.7609	3931.2408	0.0000	1.0000
SYMIG	62	-1767.9195	3933.9795	0.0000	1.0000
HKY	59	-1779.6498	3938.1752	0.0000	1.0000
GTR	63	-1770.9937	3946.5496	0.0000	1.0000
K2P	56	-1797.3608	3954.3325	0.0000	1.0000
F81G	59	-1792.2929	3963.4616	0.0000	1.0000
F81I	59	-1793.1005	3965.0767	0.0000	1.0000
SYM	60	-1791.0573	3967.4119	0.0000	1.0000
F81IG	60	-1791.7803	3968.8579	0.0000	1.0000
JC69G	56	-1807.1624	3973.9357	0.0000	1.0000
JC69I	56	-1808.0134	3975.6377	0.0000	1.0000
JC69IG	57	-1806.5129	3979.0582	0.0000	1.0000
F81	58	-1814.6729	4001.7999	0.0000	1.0000
JC69	55	-1829.9950	4013.1791	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/7659_hg18.phylip = 1 - 615;
 Partition Dummy = 1:261-aggressive-phylip/7659_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/7659_hg18.phylip = 1 - 615;
 Partition Dummy = 1:261-aggressive-phylip/7659_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/7659_hg18.phylip = 1 - 615;
 Partition Dummy = 1:261-aggressive-phylip/7659_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

