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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	-1786.6402	3720.7197	0.6418	0.6418
GTRI	64	-1787.8993	3723.2379	0.1822	0.8240
GTRIG	65	-1786.6079	3723.3095	0.1758	0.9998
HKYG	60	-1800.2779	3737.5002	0.0001	0.9999
HKYI	60	-1801.5225	3739.9894	0.0000	1.0000
HKYIG	61	-1800.2657	3740.0813	0.0000	1.0000
GTR	63	-1817.9200	3780.6372	0.0000	1.0000
SYMG	61	-1821.2646	3782.0791	0.0000	1.0000
K2PG	57	-1827.4464	3784.0928	0.0000	1.0000
SYMIG	62	-1821.2430	3784.6534	0.0000	1.0000
SYMI	61	-1822.5535	3784.6568	0.0000	1.0000
K2PIG	58	-1827.4188	3786.6071	0.0000	1.0000
K2PI	57	-1828.9049	3787.0099	0.0000	1.0000
HKY	59	-1829.8709	3794.0928	0.0000	1.0000
SYM	60	-1848.8509	3834.6462	0.0000	1.0000
K2P	56	-1855.3844	3837.4111	0.0000	1.0000
F81G	59	-1860.5614	3855.4738	0.0000	1.0000
F81IG	60	-1860.5473	3858.0390	0.0000	1.0000
F81I	59	-1861.8928	3858.1366	0.0000	1.0000
JC69G	56	-1878.1730	3882.9881	0.0000	1.0000
JC69IG	57	-1878.1538	3885.5076	0.0000	1.0000
JC69I	56	-1879.5735	3885.7893	0.0000	1.0000
F81	58	-1888.2772	3908.3239	0.0000	1.0000
JC69	55	-1905.7003	3935.4966	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1786.6402	3701.2805	0.6010	0.6010
GTRIG	65	-1786.6079	3703.2159	0.2283	0.8293
GTRI	64	-1787.8993	3703.7986	0.1706	0.9999
HKYG	60	-1800.2779	3720.5558	0.0000	1.0000
HKYIG	61	-1800.2657	3722.5314	0.0000	1.0000
HKYI	60	-1801.5225	3723.0450	0.0000	1.0000
GTR	63	-1817.9200	3761.8400	0.0000	1.0000
SYMG	61	-1821.2646	3764.5292	0.0000	1.0000
SYMIG	62	-1821.2430	3766.4860	0.0000	1.0000
SYMI	61	-1822.5535	3767.1069	0.0000	1.0000
K2PG	57	-1827.4464	3768.8928	0.0000	1.0000
K2PIG	58	-1827.4188	3770.8375	0.0000	1.0000
K2PI	57	-1828.9049	3771.8099	0.0000	1.0000
HKY	59	-1829.8709	3777.7418	0.0000	1.0000
SYM	60	-1848.8509	3817.7018	0.0000	1.0000
K2P	56	-1855.3844	3822.7689	0.0000	1.0000
F81G	59	-1860.5614	3839.1228	0.0000	1.0000
F81IG	60	-1860.5473	3841.0946	0.0000	1.0000
F81I	59	-1861.8928	3841.7856	0.0000	1.0000
JC69G	56	-1878.1730	3868.3459	0.0000	1.0000
JC69IG	57	-1878.1538	3870.3076	0.0000	1.0000
JC69I	56	-1879.5735	3871.1471	0.0000	1.0000
F81	58	-1888.2772	3892.5543	0.0000	1.0000
JC69	55	-1905.7003	3921.4005	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-1786.6402	3970.1130	0.5822	0.5822
HKYG	60	-1800.2779	3972.5863	0.1690	0.7512
GTRI	64	-1787.8993	3972.6312	0.1653	0.9165
HKYI	60	-1801.5225	3975.0755	0.0487	0.9652
GTRIG	65	-1786.6079	3976.2490	0.0271	0.9923
HKYIG	61	-1800.2657	3978.7625	0.0077	1.0000
K2PG	57	-1827.4464	4008.3218	0.0000	1.0000
K2PI	57	-1828.9049	4011.2389	0.0000	1.0000
K2PIG	58	-1827.4188	4014.4671	0.0000	1.0000
SYMG	61	-1821.2646	4020.7602	0.0000	1.0000
SYMI	61	-1822.5535	4023.3380	0.0000	1.0000
HKY	59	-1829.8709	4025.5718	0.0000	1.0000
GTR	63	-1817.9200	4026.4721	0.0000	1.0000
SYMIG	62	-1821.2430	4026.9175	0.0000	1.0000
K2P	56	-1855.3844	4057.9974	0.0000	1.0000
SYM	60	-1848.8509	4069.7323	0.0000	1.0000
F81G	59	-1860.5614	4086.9528	0.0000	1.0000
F81I	59	-1861.8928	4089.6156	0.0000	1.0000
F81IG	60	-1860.5473	4093.1251	0.0000	1.0000
JC69G	56	-1878.1730	4103.5745	0.0000	1.0000
JC69I	56	-1879.5735	4106.3756	0.0000	1.0000
JC69IG	57	-1878.1538	4109.7366	0.0000	1.0000
F81	58	-1888.2772	4136.1839	0.0000	1.0000
JC69	55	-1905.7003	4152.4285	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/12539_hg18.phylip = 1 - 493;
 Partition Dummy = 1:261-aggressive-phylip/12539_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/12539_hg18.phylip = 1 - 493;
 Partition Dummy = 1:261-aggressive-phylip/12539_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 (GTRG)]

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


-------------------------------------------------------------
End of Output

