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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	-1535.3793	3222.1953	0.7056	0.7056
GTRIG	65	-1534.9050	3224.0473	0.2795	0.9851
GTRI	64	-1539.2405	3229.9177	0.0148	1.0000
HKYG	60	-1551.9445	3244.2789	0.0000	1.0000
HKYIG	61	-1551.6834	3246.4953	0.0000	1.0000
SYMG	61	-1555.7302	3254.5889	0.0000	1.0000
SYMIG	62	-1555.3365	3256.5553	0.0000	1.0000
HKYI	60	-1558.8016	3257.9932	0.0000	1.0000
K2PG	57	-1564.2649	3260.7950	0.0000	1.0000
K2PIG	58	-1563.1707	3261.2998	0.0000	1.0000
SYMI	61	-1560.8207	3264.7700	0.0000	1.0000
K2PI	57	-1568.7839	3269.8330	0.0000	1.0000
F81G	59	-1582.4660	3302.5987	0.0000	1.0000
F81IG	60	-1581.7027	3303.7955	0.0000	1.0000
F81I	59	-1585.7000	3309.0666	0.0000	1.0000
GTR	63	-1582.2837	3313.2190	0.0000	1.0000
JC69G	56	-1594.4628	3318.5124	0.0000	1.0000
JC69IG	57	-1593.6162	3319.4976	0.0000	1.0000
JC69I	56	-1597.9574	3325.5016	0.0000	1.0000
HKY	59	-1595.2650	3328.1967	0.0000	1.0000
SYM	60	-1600.2643	3340.9187	0.0000	1.0000
K2P	56	-1609.1290	3347.8447	0.0000	1.0000
F81	58	-1624.2514	3383.4613	0.0000	1.0000
JC69	55	-1635.8224	3398.5679	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1535.3793	3198.7587	0.6203	0.6203
GTRIG	65	-1534.9050	3199.8101	0.3667	0.9869
GTRI	64	-1539.2405	3206.4811	0.0131	1.0000
HKYG	60	-1551.9445	3223.8889	0.0000	1.0000
HKYIG	61	-1551.6834	3225.3668	0.0000	1.0000
SYMG	61	-1555.7302	3233.4604	0.0000	1.0000
SYMIG	62	-1555.3365	3234.6729	0.0000	1.0000
HKYI	60	-1558.8016	3237.6032	0.0000	1.0000
K2PIG	58	-1563.1707	3242.3414	0.0000	1.0000
K2PG	57	-1564.2649	3242.5298	0.0000	1.0000
SYMI	61	-1560.8207	3243.6415	0.0000	1.0000
K2PI	57	-1568.7839	3251.5678	0.0000	1.0000
F81G	59	-1582.4660	3282.9320	0.0000	1.0000
F81IG	60	-1581.7027	3283.4055	0.0000	1.0000
F81I	59	-1585.7000	3289.3999	0.0000	1.0000
GTR	63	-1582.2837	3290.5673	0.0000	1.0000
JC69G	56	-1594.4628	3300.9256	0.0000	1.0000
JC69IG	57	-1593.6162	3301.2324	0.0000	1.0000
JC69I	56	-1597.9574	3307.9149	0.0000	1.0000
HKY	59	-1595.2650	3308.5301	0.0000	1.0000
SYM	60	-1600.2643	3320.5287	0.0000	1.0000
K2P	56	-1609.1290	3330.2579	0.0000	1.0000
F81	58	-1624.2514	3364.5028	0.0000	1.0000
JC69	55	-1635.8224	3381.6448	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-1535.3793	3457.3350	0.8993	0.8993
GTRIG	65	-1534.9050	3462.4266	0.0705	0.9698
GTRI	64	-1539.2405	3465.0574	0.0189	0.9887
HKYG	60	-1551.9445	3466.3042	0.0101	0.9989
HKYIG	61	-1551.6834	3471.8223	0.0006	0.9995
K2PG	57	-1564.2649	3472.8244	0.0004	0.9999
K2PIG	58	-1563.1707	3476.6762	0.0001	1.0000
SYMG	61	-1555.7302	3479.9160	0.0000	1.0000
HKYI	60	-1558.8016	3480.0185	0.0000	1.0000
K2PI	57	-1568.7839	3481.8623	0.0000	1.0000
SYMIG	62	-1555.3365	3485.1687	0.0000	1.0000
SYMI	61	-1560.8207	3490.0970	0.0000	1.0000
F81G	59	-1582.4660	3521.3070	0.0000	1.0000
F81IG	60	-1581.7027	3525.8208	0.0000	1.0000
JC69G	56	-1594.4628	3527.1799	0.0000	1.0000
F81I	59	-1585.7000	3527.7750	0.0000	1.0000
JC69IG	57	-1593.6162	3531.5269	0.0000	1.0000
JC69I	56	-1597.9574	3534.1691	0.0000	1.0000
GTR	63	-1582.2837	3545.1033	0.0000	1.0000
HKY	59	-1595.2650	3546.9051	0.0000	1.0000
K2P	56	-1609.1290	3556.5122	0.0000	1.0000
SYM	60	-1600.2643	3562.9440	0.0000	1.0000
F81	58	-1624.2514	3598.8376	0.0000	1.0000
JC69	55	-1635.8224	3603.8588	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/7850_hg18.phylip = 1 - 420;
 Partition Dummy = 1:261-aggressive-phylip/7850_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/7850_hg18.phylip = 1 - 420;
 Partition Dummy = 1:261-aggressive-phylip/7850_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 (GTRG)]

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


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

