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

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

Minimum AIC  model: GTRG
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	-1485.3354	3108.7895	0.3750	0.3750
GTRG	64	-1480.6002	3110.0004	0.2047	0.5797
HKYI	60	-1486.0337	3110.1863	0.1865	0.7662
HKYIG	61	-1485.2868	3111.3429	0.1046	0.8708
GTRI	64	-1481.6517	3112.1034	0.0715	0.9423
GTRIG	65	-1480.5149	3112.5336	0.0577	1.0000
K2PG	57	-1502.2954	3134.8365	0.0000	1.0000
K2PI	57	-1502.9702	3136.1861	0.0000	1.0000
K2PIG	58	-1502.2510	3137.3591	0.0000	1.0000
HKY	59	-1501.3895	3138.2605	0.0000	1.0000
GTR	63	-1497.1128	3140.3353	0.0000	1.0000
SYMG	61	-1499.9282	3140.6257	0.0000	1.0000
SYMI	61	-1500.5908	3141.9508	0.0000	1.0000
SYMIG	62	-1499.8850	3143.2028	0.0000	1.0000
F81G	59	-1513.9087	3163.2988	0.0000	1.0000
K2P	56	-1518.2436	3164.1344	0.0000	1.0000
F81IG	60	-1513.8539	3165.8267	0.0000	1.0000
F81I	59	-1515.9890	3167.4595	0.0000	1.0000
SYM	60	-1515.9799	3170.0786	0.0000	1.0000
JC69G	56	-1530.4412	3188.5295	0.0000	1.0000
JC69I	56	-1531.1854	3190.0180	0.0000	1.0000
JC69IG	57	-1530.4005	3191.0466	0.0000	1.0000
F81	58	-1530.5566	3193.9704	0.0000	1.0000
JC69	55	-1547.3431	3219.7473	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1480.6002	3089.2004	0.3769	0.3769
HKYG	60	-1485.3354	3090.6707	0.1807	0.5576
GTRIG	65	-1480.5149	3091.0299	0.1510	0.7086
GTRI	64	-1481.6517	3091.3034	0.1317	0.8403
HKYI	60	-1486.0337	3092.0675	0.0899	0.9302
HKYIG	61	-1485.2868	3092.5737	0.0698	1.0000
K2PG	57	-1502.2954	3118.5908	0.0000	1.0000
K2PI	57	-1502.9702	3119.9404	0.0000	1.0000
GTR	63	-1497.1128	3120.2256	0.0000	1.0000
K2PIG	58	-1502.2510	3120.5020	0.0000	1.0000
HKY	59	-1501.3895	3120.7790	0.0000	1.0000
SYMG	61	-1499.9282	3121.8565	0.0000	1.0000
SYMI	61	-1500.5908	3123.1816	0.0000	1.0000
SYMIG	62	-1499.8850	3123.7699	0.0000	1.0000
F81G	59	-1513.9087	3145.8173	0.0000	1.0000
F81IG	60	-1513.8539	3147.7079	0.0000	1.0000
K2P	56	-1518.2436	3148.4873	0.0000	1.0000
F81I	59	-1515.9890	3149.9781	0.0000	1.0000
SYM	60	-1515.9799	3151.9598	0.0000	1.0000
JC69G	56	-1530.4412	3172.8825	0.0000	1.0000
JC69I	56	-1531.1854	3174.3709	0.0000	1.0000
JC69IG	57	-1530.4005	3174.8009	0.0000	1.0000
F81	58	-1530.5566	3177.1132	0.0000	1.0000
JC69	55	-1547.3431	3204.6862	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1485.3354	3339.1930	0.6462	0.6462
HKYI	60	-1486.0337	3340.5897	0.3214	0.9676
HKYIG	61	-1485.2868	3345.2379	0.0315	0.9991
GTRG	64	-1480.6002	3354.2908	0.0003	0.9994
K2PG	57	-1502.2954	3354.6870	0.0003	0.9997
K2PI	57	-1502.9702	3356.0365	0.0001	0.9998
GTRI	64	-1481.6517	3356.3938	0.0001	1.0000
GTRIG	65	-1480.5149	3360.2623	0.0000	1.0000
K2PIG	58	-1502.2510	3360.7402	0.0000	1.0000
HKY	59	-1501.3895	3365.1592	0.0000	1.0000
SYMG	61	-1499.9282	3374.5207	0.0000	1.0000
SYMI	61	-1500.5908	3375.8459	0.0000	1.0000
K2P	56	-1518.2436	3380.4414	0.0000	1.0000
SYMIG	62	-1499.8850	3380.5762	0.0000	1.0000
GTR	63	-1497.1128	3381.1739	0.0000	1.0000
F81G	59	-1513.9087	3390.1975	0.0000	1.0000
F81I	59	-1515.9890	3394.3583	0.0000	1.0000
F81IG	60	-1513.8539	3396.2301	0.0000	1.0000
SYM	60	-1515.9799	3400.4820	0.0000	1.0000
JC69G	56	-1530.4412	3404.8366	0.0000	1.0000
JC69I	56	-1531.1854	3406.3250	0.0000	1.0000
JC69IG	57	-1530.4005	3410.8971	0.0000	1.0000
F81	58	-1530.5566	3417.3514	0.0000	1.0000
JC69	55	-1547.3431	3432.4982	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/9088_hg18.phylip = 1 - 465;
 Partition Dummy = 1:261-aggressive-phylip/9088_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 (HKYG)]

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

