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

Input data from file "261-aggressive-phylip/3357_hg18.phylip" (29 taxa, 425 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	-1568.8322	3288.7756	0.6848	0.6848
GTRIG	65	-1568.7404	3291.3805	0.1862	0.8710
HKYG	60	-1576.3637	3292.8373	0.0899	0.9608
HKYIG	61	-1576.2891	3295.4156	0.0248	0.9856
GTRI	64	-1572.9248	3296.9607	0.0114	0.9970
HKYI	60	-1579.7722	3299.6544	0.0030	1.0000
GTR	63	-1586.7051	3321.7482	0.0000	1.0000
HKY	59	-1594.0520	3325.5013	0.0000	1.0000
K2PG	57	-1604.5023	3341.0208	0.0000	1.0000
K2PIG	58	-1604.4253	3343.5500	0.0000	1.0000
K2PI	57	-1607.0368	3346.0900	0.0000	1.0000
SYMG	61	-1602.7313	3348.3000	0.0000	1.0000
F81G	59	-1606.5556	3350.5084	0.0000	1.0000
SYMIG	62	-1602.6559	3350.8919	0.0000	1.0000
F81IG	60	-1606.4463	3353.0025	0.0000	1.0000
SYMI	61	-1605.1574	3353.1522	0.0000	1.0000
F81I	59	-1609.5282	3356.4536	0.0000	1.0000
K2P	56	-1624.1464	3377.6405	0.0000	1.0000
SYM	60	-1621.5978	3383.3054	0.0000	1.0000
F81	58	-1624.5685	3383.8364	0.0000	1.0000
JC69G	56	-1630.4226	3390.1931	0.0000	1.0000
JC69IG	57	-1630.3483	3392.7129	0.0000	1.0000
JC69I	56	-1632.9588	3395.2655	0.0000	1.0000
JC69	55	-1649.7072	3426.1082	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1568.8322	3265.6645	0.6841	0.6841
GTRIG	65	-1568.7404	3267.4808	0.2759	0.9600
HKYG	60	-1576.3637	3272.7274	0.0200	0.9800
GTRI	64	-1572.9248	3273.8496	0.0114	0.9914
HKYIG	61	-1576.2891	3274.5782	0.0079	0.9993
HKYI	60	-1579.7722	3279.5445	0.0007	1.0000
GTR	63	-1586.7051	3299.4103	0.0000	1.0000
HKY	59	-1594.0520	3306.1041	0.0000	1.0000
K2PG	57	-1604.5023	3323.0045	0.0000	1.0000
K2PIG	58	-1604.4253	3324.8505	0.0000	1.0000
SYMG	61	-1602.7313	3327.4625	0.0000	1.0000
K2PI	57	-1607.0368	3328.0737	0.0000	1.0000
SYMIG	62	-1602.6559	3329.3118	0.0000	1.0000
F81G	59	-1606.5556	3331.1111	0.0000	1.0000
SYMI	61	-1605.1574	3332.3147	0.0000	1.0000
F81IG	60	-1606.4463	3332.8926	0.0000	1.0000
F81I	59	-1609.5282	3337.0563	0.0000	1.0000
K2P	56	-1624.1464	3360.2927	0.0000	1.0000
SYM	60	-1621.5978	3363.1955	0.0000	1.0000
F81	58	-1624.5685	3365.1370	0.0000	1.0000
JC69G	56	-1630.4226	3372.8453	0.0000	1.0000
JC69IG	57	-1630.3483	3374.6965	0.0000	1.0000
JC69I	56	-1632.9588	3377.9176	0.0000	1.0000
JC69	55	-1649.7072	3409.4144	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1576.3637	3515.8527	0.9121	0.9121
HKYIG	61	-1576.2891	3521.7556	0.0477	0.9597
HKYI	60	-1579.7722	3522.6698	0.0302	0.9899
GTRG	64	-1568.8322	3524.9982	0.0094	0.9993
GTRIG	65	-1568.7404	3530.8666	0.0005	0.9998
GTRI	64	-1572.9248	3533.1833	0.0002	1.0000
HKY	59	-1594.0520	3545.1773	0.0000	1.0000
K2PG	57	-1604.5023	3553.9736	0.0000	1.0000
GTR	63	-1586.7051	3554.6919	0.0000	1.0000
K2PI	57	-1607.0368	3559.0428	0.0000	1.0000
K2PIG	58	-1604.4253	3559.8717	0.0000	1.0000
F81G	59	-1606.5556	3570.1844	0.0000	1.0000
SYMG	61	-1602.7313	3574.6400	0.0000	1.0000
F81IG	60	-1606.4463	3576.0180	0.0000	1.0000
F81I	59	-1609.5282	3576.1296	0.0000	1.0000
SYMI	61	-1605.1574	3579.4922	0.0000	1.0000
SYMIG	62	-1602.6559	3580.5413	0.0000	1.0000
K2P	56	-1624.1464	3587.2097	0.0000	1.0000
JC69G	56	-1630.4226	3599.7623	0.0000	1.0000
F81	58	-1624.5685	3600.1582	0.0000	1.0000
JC69I	56	-1632.9588	3604.8346	0.0000	1.0000
JC69IG	57	-1630.3483	3605.6656	0.0000	1.0000
SYM	60	-1621.5978	3606.3209	0.0000	1.0000
JC69	55	-1649.7072	3632.2793	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/3357_hg18.phylip = 1 - 425;
 Partition Dummy = 1:261-aggressive-phylip/3357_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/3357_hg18.phylip = 1 - 425;
 Partition Dummy = 1:261-aggressive-phylip/3357_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/3357_hg18.phylip = 1 - 425;
 Partition Dummy = 1:261-aggressive-phylip/3357_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

