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

Input data from file "261-aggressive-phylip/4714_hg18.phylip" (29 taxa, 430 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	-1645.3169	3441.4283	0.9237	0.9237
GTRIG	65	-1646.5161	3446.6036	0.0695	0.9932
HKYG	60	-1655.7088	3451.2551	0.0068	1.0000
HKYI	60	-1662.4246	3464.6865	0.0000	1.0000
GTRI	64	-1657.1072	3465.0089	0.0000	1.0000
SYMG	61	-1663.0321	3468.6185	0.0000	1.0000
SYMIG	62	-1663.0732	3471.4326	0.0000	1.0000
HKYIG	61	-1664.8901	3472.3346	0.0000	1.0000
SYMI	61	-1669.4743	3481.5029	0.0000	1.0000
K2PG	57	-1683.1539	3498.0821	0.0000	1.0000
K2PI	57	-1687.0735	3505.9212	0.0000	1.0000
GTR	63	-1679.9422	3507.9171	0.0000	1.0000
K2PIG	58	-1692.6582	3519.7639	0.0000	1.0000
HKY	59	-1695.5936	3528.3223	0.0000	1.0000
SYM	60	-1698.8533	3537.5439	0.0000	1.0000
F81G	59	-1707.6073	3552.3497	0.0000	1.0000
F81IG	60	-1709.6772	3559.1919	0.0000	1.0000
F81I	59	-1711.7799	3560.6950	0.0000	1.0000
K2P	56	-1727.7970	3584.7092	0.0000	1.0000
JC69I	56	-1755.6170	3640.3493	0.0000	1.0000
F81	58	-1754.6021	3643.6517	0.0000	1.0000
JC69G	56	-1760.1751	3649.4655	0.0000	1.0000
JC69IG	57	-1760.0161	3651.8065	0.0000	1.0000
JC69	55	-1811.1175	3748.7056	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1645.3169	3418.6338	0.8988	0.8988
GTRIG	65	-1646.5161	3423.0321	0.0997	0.9985
HKYG	60	-1655.7088	3431.4177	0.0015	1.0000
GTRI	64	-1657.1072	3442.2144	0.0000	1.0000
HKYI	60	-1662.4246	3444.8491	0.0000	1.0000
SYMG	61	-1663.0321	3448.0641	0.0000	1.0000
SYMIG	62	-1663.0732	3450.1465	0.0000	1.0000
HKYIG	61	-1664.8901	3451.7803	0.0000	1.0000
SYMI	61	-1669.4743	3460.9486	0.0000	1.0000
K2PG	57	-1683.1539	3480.3079	0.0000	1.0000
GTR	63	-1679.9422	3485.8844	0.0000	1.0000
K2PI	57	-1687.0735	3488.1470	0.0000	1.0000
K2PIG	58	-1692.6582	3501.3164	0.0000	1.0000
HKY	59	-1695.5936	3509.1872	0.0000	1.0000
SYM	60	-1698.8533	3517.7065	0.0000	1.0000
F81G	59	-1707.6073	3533.2145	0.0000	1.0000
F81IG	60	-1709.6772	3539.3545	0.0000	1.0000
F81I	59	-1711.7799	3541.5598	0.0000	1.0000
K2P	56	-1727.7970	3567.5939	0.0000	1.0000
JC69I	56	-1755.6170	3623.2340	0.0000	1.0000
F81	58	-1754.6021	3625.2043	0.0000	1.0000
JC69G	56	-1760.1751	3632.3502	0.0000	1.0000
JC69IG	57	-1760.0161	3634.0323	0.0000	1.0000
JC69	55	-1811.1175	3732.2351	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1655.7088	3675.2448	0.8474	0.8474
GTRG	64	-1645.3169	3678.7160	0.1494	0.9968
GTRIG	65	-1646.5161	3687.1782	0.0022	0.9989
HKYI	60	-1662.4246	3688.6762	0.0010	1.0000
SYMG	61	-1663.0321	3695.9550	0.0000	1.0000
HKYIG	61	-1664.8901	3699.6712	0.0000	1.0000
SYMIG	62	-1663.0732	3702.1012	0.0000	1.0000
GTRI	64	-1657.1072	3702.2966	0.0000	1.0000
SYMI	61	-1669.4743	3708.8395	0.0000	1.0000
K2PG	57	-1683.1539	3711.9436	0.0000	1.0000
K2PI	57	-1687.0735	3719.7828	0.0000	1.0000
K2PIG	58	-1692.6582	3737.0160	0.0000	1.0000
GTR	63	-1679.9422	3741.9028	0.0000	1.0000
HKY	59	-1695.5936	3748.9505	0.0000	1.0000
SYM	60	-1698.8533	3761.5337	0.0000	1.0000
F81G	59	-1707.6073	3772.9779	0.0000	1.0000
F81I	59	-1711.7799	3781.3232	0.0000	1.0000
F81IG	60	-1709.6772	3783.1816	0.0000	1.0000
K2P	56	-1727.7970	3795.1659	0.0000	1.0000
JC69I	56	-1755.6170	3850.8060	0.0000	1.0000
JC69G	56	-1760.1751	3859.9222	0.0000	1.0000
F81	58	-1754.6021	3860.9038	0.0000	1.0000
JC69IG	57	-1760.0161	3865.6680	0.0000	1.0000
JC69	55	-1811.1175	3955.7432	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/4714_hg18.phylip = 1 - 430;
 Partition Dummy = 1:261-aggressive-phylip/4714_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/4714_hg18.phylip = 1 - 430;
 Partition Dummy = 1:261-aggressive-phylip/4714_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/4714_hg18.phylip = 1 - 430;
 Partition Dummy = 1:261-aggressive-phylip/4714_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

