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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-1480.0800	3102.4770	0.7316	0.7316
HKYIG	61	-1480.1990	3105.5295	0.1590	0.8907
HKYI	60	-1482.0772	3106.4714	0.0993	0.9900
GTRG	64	-1479.0165	3111.7119	0.0072	0.9972
GTRIG	65	-1479.1273	3114.8181	0.0015	0.9987
GTRI	64	-1480.7428	3115.1645	0.0013	1.0000
K2PG	57	-1497.0016	3127.9790	0.0000	1.0000
K2PIG	58	-1497.6439	3132.0272	0.0000	1.0000
K2PI	57	-1499.5361	3133.0481	0.0000	1.0000
SYMG	61	-1496.2990	3137.7295	0.0000	1.0000
SYMIG	62	-1496.2986	3140.5605	0.0000	1.0000
SYMI	61	-1497.7303	3140.5922	0.0000	1.0000
HKY	59	-1511.4069	3162.3335	0.0000	1.0000
GTR	63	-1510.3880	3171.5883	0.0000	1.0000
F81G	59	-1525.2525	3190.0248	0.0000	1.0000
K2P	56	-1529.8839	3190.9968	0.0000	1.0000
F81IG	60	-1525.3237	3192.9645	0.0000	1.0000
F81I	59	-1527.4027	3194.3253	0.0000	1.0000
SYM	60	-1528.4873	3199.2917	0.0000	1.0000
JC69G	56	-1539.8097	3210.8484	0.0000	1.0000
JC69IG	57	-1539.8094	3213.5947	0.0000	1.0000
JC69I	56	-1541.7924	3214.8137	0.0000	1.0000
F81	58	-1556.5967	3249.9328	0.0000	1.0000
JC69	55	-1571.5260	3271.5506	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1480.0800	3080.1599	0.6484	0.6484
HKYIG	61	-1480.1990	3082.3980	0.2118	0.8602
HKYI	60	-1482.0772	3084.1544	0.0880	0.9482
GTRG	64	-1479.0165	3086.0329	0.0344	0.9826
GTRIG	65	-1479.1273	3088.2546	0.0113	0.9939
GTRI	64	-1480.7428	3089.4855	0.0061	1.0000
K2PG	57	-1497.0016	3108.0032	0.0000	1.0000
K2PIG	58	-1497.6439	3111.2878	0.0000	1.0000
K2PI	57	-1499.5361	3113.0722	0.0000	1.0000
SYMG	61	-1496.2990	3114.5980	0.0000	1.0000
SYMIG	62	-1496.2986	3116.5973	0.0000	1.0000
SYMI	61	-1497.7303	3117.4607	0.0000	1.0000
HKY	59	-1511.4069	3140.8138	0.0000	1.0000
GTR	63	-1510.3880	3146.7760	0.0000	1.0000
F81G	59	-1525.2525	3168.5051	0.0000	1.0000
F81IG	60	-1525.3237	3170.6474	0.0000	1.0000
K2P	56	-1529.8839	3171.7678	0.0000	1.0000
F81I	59	-1527.4027	3172.8055	0.0000	1.0000
SYM	60	-1528.4873	3176.9746	0.0000	1.0000
JC69G	56	-1539.8097	3191.6195	0.0000	1.0000
JC69IG	57	-1539.8094	3193.6189	0.0000	1.0000
JC69I	56	-1541.7924	3195.5848	0.0000	1.0000
F81	58	-1556.5967	3229.1934	0.0000	1.0000
JC69	55	-1571.5260	3253.0521	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1480.0800	3317.9747	0.8466	0.8466
HKYI	60	-1482.0772	3321.9691	0.1149	0.9615
HKYIG	61	-1480.1990	3324.1763	0.0381	0.9997
K2PG	57	-1497.0016	3333.9272	0.0003	0.9999
K2PI	57	-1499.5361	3338.9963	0.0000	1.0000
GTRG	64	-1479.0165	3339.7020	0.0000	1.0000
K2PIG	58	-1497.6439	3341.1754	0.0000	1.0000
GTRI	64	-1480.7428	3343.1546	0.0000	1.0000
GTRIG	65	-1479.1273	3345.8873	0.0000	1.0000
SYMG	61	-1496.2990	3356.3763	0.0000	1.0000
SYMI	61	-1497.7303	3359.2390	0.0000	1.0000
SYMIG	62	-1496.2986	3362.3392	0.0000	1.0000
HKY	59	-1511.4069	3374.6649	0.0000	1.0000
K2P	56	-1529.8839	3393.7283	0.0000	1.0000
GTR	63	-1510.3880	3396.4815	0.0000	1.0000
F81G	59	-1525.2525	3402.3563	0.0000	1.0000
F81I	59	-1527.4027	3406.6567	0.0000	1.0000
F81IG	60	-1525.3237	3408.4621	0.0000	1.0000
JC69G	56	-1539.8097	3413.5799	0.0000	1.0000
SYM	60	-1528.4873	3414.7894	0.0000	1.0000
JC69I	56	-1541.7924	3417.5452	0.0000	1.0000
JC69IG	57	-1539.8094	3419.5429	0.0000	1.0000
F81	58	-1556.5967	3459.0810	0.0000	1.0000
JC69	55	-1571.5260	3471.0489	0.0000	1.0000



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

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

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


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

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


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

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

