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

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

Minimum AIC  model: GTR
Minimum AICc model: K2P
Minimum BIC  model: K2P


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
K2P	56	-658.1673	1466.5621	0.5239	0.5239
HKY	59	-653.7112	1468.5932	0.1898	0.7136
K2PI	57	-658.1294	1470.0901	0.0898	0.8034
K2PG	57	-658.1599	1470.1510	0.0871	0.8905
HKYI	60	-653.6309	1472.1698	0.0317	0.9222
HKYG	60	-653.6790	1472.2660	0.0302	0.9525
K2PIG	58	-658.1325	1473.7437	0.0144	0.9669
SYM	60	-654.6872	1474.2823	0.0110	0.9780
GTR	63	-649.0735	1474.5470	0.0097	0.9876
HKYIG	61	-653.6343	1475.9600	0.0048	0.9924
SYMI	61	-654.6023	1477.8959	0.0018	0.9942
GTRI	64	-648.7994	1477.9258	0.0018	0.9960
SYMG	61	-654.6516	1477.9946	0.0017	0.9977
GTRG	64	-648.8935	1478.1141	0.0016	0.9993
SYMIG	62	-654.6054	1481.7326	0.0003	0.9996
GTRIG	65	-648.8039	1481.9116	0.0002	0.9999
JC69	55	-668.6797	1484.0260	0.0001	0.9999
F81	58	-664.5462	1486.5712	0.0000	1.0000
JC69I	56	-668.6528	1487.5332	0.0000	1.0000
JC69G	56	-668.6784	1487.5843	0.0000	1.0000
F81I	59	-664.5029	1490.1766	0.0000	1.0000
F81G	59	-664.5392	1490.2491	0.0000	1.0000
JC69IG	57	-668.6569	1491.1450	0.0000	1.0000
F81IG	60	-664.5062	1493.9203	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTR	63	-649.0735	1424.1470	0.2863	0.2863
HKY	59	-653.7112	1425.4225	0.1513	0.4375
GTRI	64	-648.7994	1425.5988	0.1385	0.5761
GTRG	64	-648.8935	1425.7871	0.1261	0.7021
HKYI	60	-653.6309	1427.2618	0.0603	0.7624
HKYG	60	-653.6790	1427.3581	0.0575	0.8199
GTRIG	65	-648.8039	1427.6078	0.0507	0.8707
K2P	56	-658.1673	1428.3346	0.0353	0.9059
HKYIG	61	-653.6343	1429.2686	0.0221	0.9280
SYM	60	-654.6872	1429.3744	0.0210	0.9490
K2PI	57	-658.1294	1430.2587	0.0135	0.9625
K2PG	57	-658.1599	1430.3197	0.0131	0.9756
SYMI	61	-654.6023	1431.2046	0.0084	0.9840
SYMG	61	-654.6516	1431.3033	0.0080	0.9920
K2PIG	58	-658.1325	1432.2650	0.0049	0.9969
SYMIG	62	-654.6054	1433.2108	0.0031	1.0000
F81	58	-664.5462	1445.0924	0.0000	1.0000
F81I	59	-664.5029	1447.0058	0.0000	1.0000
F81G	59	-664.5392	1447.0783	0.0000	1.0000
JC69	55	-668.6797	1447.3594	0.0000	1.0000
F81IG	60	-664.5062	1449.0123	0.0000	1.0000
JC69I	56	-668.6528	1449.3057	0.0000	1.0000
JC69G	56	-668.6784	1449.3568	0.0000	1.0000
JC69IG	57	-668.6569	1451.3137	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	56	-658.1673	1619.3867	0.8532	0.8532
K2PI	57	-658.1294	1624.7226	0.0592	0.9124
K2PG	57	-658.1599	1624.7835	0.0574	0.9698
HKY	59	-653.7112	1626.7096	0.0219	0.9918
K2PIG	58	-658.1325	1630.1404	0.0039	0.9957
HKYI	60	-653.6309	1631.9606	0.0016	0.9973
HKYG	60	-653.6790	1632.0568	0.0015	0.9988
SYM	60	-654.6872	1634.0731	0.0006	0.9994
JC69	55	-668.6797	1634.9999	0.0003	0.9997
HKYIG	61	-653.6343	1637.3790	0.0001	0.9998
GTR	63	-649.0735	1639.0807	0.0000	0.9999
SYMI	61	-654.6023	1639.3150	0.0000	0.9999
SYMG	61	-654.6516	1639.4137	0.0000	0.9999
JC69I	56	-668.6528	1640.3579	0.0000	1.0000
JC69G	56	-668.6784	1640.4090	0.0000	1.0000
F81	58	-664.5462	1642.9679	0.0000	1.0000
GTRI	64	-648.7994	1643.9441	0.0000	1.0000
GTRG	64	-648.8935	1644.1324	0.0000	1.0000
SYMIG	62	-654.6054	1644.7329	0.0000	1.0000
JC69IG	57	-668.6569	1645.7775	0.0000	1.0000
F81I	59	-664.5029	1648.2929	0.0000	1.0000
F81G	59	-664.5392	1648.3655	0.0000	1.0000
GTRIG	65	-648.8039	1649.3648	0.0000	1.0000
F81IG	60	-664.5062	1653.7111	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTR)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1308_hg18.phylip = 1 - 224;
 Partition Dummy = 1:261-aggressive-phylip/1308_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=equal;
 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);
END;


[Mrbayes block for the best AICc model (K2P)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1308_hg18.phylip = 1 - 224;
 Partition Dummy = 1:261-aggressive-phylip/1308_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


[Mrbayes block for the best BIC model (K2P)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1308_hg18.phylip = 1 - 224;
 Partition Dummy = 1:261-aggressive-phylip/1308_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

