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

Input data from file "261-aggressive-phylip/3762_hg18.phylip" (29 taxa, 364 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	-1521.7660	3187.6905	0.9543	0.9543
GTRG	64	-1519.4376	3194.7013	0.0287	0.9830
HKYIG	61	-1524.9235	3196.8933	0.0096	0.9925
GTRIG	65	-1519.5810	3197.9539	0.0056	0.9982
HKYI	60	-1528.1379	3200.4342	0.0016	0.9998
GTRI	64	-1524.5901	3205.0063	0.0002	1.0000
K2PG	57	-1536.6980	3209.0038	0.0000	1.0000
K2PIG	58	-1536.5152	3211.4698	0.0000	1.0000
K2PI	57	-1540.0443	3215.6965	0.0000	1.0000
SYMG	61	-1535.2716	3217.5896	0.0000	1.0000
SYMIG	62	-1535.6283	3221.2101	0.0000	1.0000
SYMI	61	-1538.9794	3225.0051	0.0000	1.0000
HKY	59	-1568.6371	3278.5637	0.0000	1.0000
GTR	63	-1567.6623	3288.2046	0.0000	1.0000
K2P	56	-1581.0444	3294.8836	0.0000	1.0000
F81G	59	-1580.6801	3302.6496	0.0000	1.0000
SYM	60	-1579.3281	3302.8145	0.0000	1.0000
JC69G	56	-1586.2741	3305.3429	0.0000	1.0000
F81IG	60	-1581.5002	3307.1588	0.0000	1.0000
F81I	59	-1583.9735	3309.2364	0.0000	1.0000
JC69IG	57	-1590.1640	3315.9358	0.0000	1.0000
JC69I	56	-1592.2696	3317.3340	0.0000	1.0000
F81	58	-1626.0244	3390.4882	0.0000	1.0000
JC69	55	-1631.5687	3393.1374	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1521.7660	3163.5321	0.7897	0.7897
GTRG	64	-1519.4376	3166.8752	0.1484	0.9381
GTRIG	65	-1519.5810	3169.1620	0.0473	0.9854
HKYIG	61	-1524.9235	3171.8470	0.0124	0.9978
HKYI	60	-1528.1379	3176.2758	0.0013	0.9991
GTRI	64	-1524.5901	3177.1802	0.0009	1.0000
K2PG	57	-1536.6980	3187.3959	0.0000	1.0000
K2PIG	58	-1536.5152	3189.0305	0.0000	1.0000
SYMG	61	-1535.2716	3192.5433	0.0000	1.0000
K2PI	57	-1540.0443	3194.0886	0.0000	1.0000
SYMIG	62	-1535.6283	3195.2566	0.0000	1.0000
SYMI	61	-1538.9794	3199.9587	0.0000	1.0000
HKY	59	-1568.6371	3255.2743	0.0000	1.0000
GTR	63	-1567.6623	3261.3246	0.0000	1.0000
K2P	56	-1581.0444	3274.0888	0.0000	1.0000
SYM	60	-1579.3281	3278.6561	0.0000	1.0000
F81G	59	-1580.6801	3279.3602	0.0000	1.0000
F81IG	60	-1581.5002	3283.0003	0.0000	1.0000
JC69G	56	-1586.2741	3284.5481	0.0000	1.0000
F81I	59	-1583.9735	3285.9470	0.0000	1.0000
JC69IG	57	-1590.1640	3294.3279	0.0000	1.0000
JC69I	56	-1592.2696	3296.5392	0.0000	1.0000
F81	58	-1626.0244	3368.0489	0.0000	1.0000
JC69	55	-1631.5687	3373.1374	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1521.7660	3397.3613	0.9935	0.9935
K2PG	57	-1536.6980	3409.5337	0.0023	0.9958
HKYIG	61	-1524.9235	3409.5734	0.0022	0.9980
HKYI	60	-1528.1379	3410.1050	0.0017	0.9997
K2PIG	58	-1536.5152	3415.0654	0.0001	0.9998
K2PI	57	-1540.0443	3416.2264	0.0001	0.9999
GTRG	64	-1519.4376	3416.2931	0.0001	1.0000
GTRIG	65	-1519.5810	3422.4770	0.0000	1.0000
GTRI	64	-1524.5901	3426.5981	0.0000	1.0000
SYMG	61	-1535.2716	3430.2697	0.0000	1.0000
SYMIG	62	-1535.6283	3436.8801	0.0000	1.0000
SYMI	61	-1538.9794	3437.6851	0.0000	1.0000
HKY	59	-1568.6371	3485.2063	0.0000	1.0000
K2P	56	-1581.0444	3492.3294	0.0000	1.0000
JC69G	56	-1586.2741	3502.7888	0.0000	1.0000
GTR	63	-1567.6623	3506.8453	0.0000	1.0000
F81G	59	-1580.6801	3509.2922	0.0000	1.0000
SYM	60	-1579.3281	3512.4854	0.0000	1.0000
JC69I	56	-1592.2696	3514.7798	0.0000	1.0000
F81I	59	-1583.9735	3515.8790	0.0000	1.0000
JC69IG	57	-1590.1640	3516.4657	0.0000	1.0000
F81IG	60	-1581.5002	3516.8296	0.0000	1.0000
JC69	55	-1631.5687	3587.4809	0.0000	1.0000
F81	58	-1626.0244	3594.0838	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/3762_hg18.phylip = 1 - 364;
 Partition Dummy = 1:261-aggressive-phylip/3762_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/3762_hg18.phylip = 1 - 364;
 Partition Dummy = 1:261-aggressive-phylip/3762_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/3762_hg18.phylip = 1 - 364;
 Partition Dummy = 1:261-aggressive-phylip/3762_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

