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

Input data from file "261-aggressive-phylip/13347_hg18.phylip" (29 taxa, 609 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	-1515.6278	3174.5497	0.3630	0.3630
GTRI	64	-1516.2380	3175.7702	0.1972	0.5601
HKYG	60	-1521.3366	3176.0308	0.1731	0.7332
GTRIG	65	-1515.1233	3176.0477	0.1716	0.9048
HKYI	60	-1521.9622	3177.2820	0.0926	0.9974
HKYIG	61	-1524.2932	3184.4146	0.0026	1.0000
GTR	63	-1529.2108	3199.2180	0.0000	1.0000
HKY	59	-1534.8862	3200.6686	0.0000	1.0000
K2PG	57	-1548.1028	3222.2056	0.0000	1.0000
K2PI	57	-1548.5428	3223.0856	0.0000	1.0000
SYMG	61	-1544.8123	3225.4528	0.0000	1.0000
F81G	59	-1547.3474	3225.5910	0.0000	1.0000
F81I	59	-1547.6518	3226.1997	0.0000	1.0000
SYMI	61	-1545.2908	3226.4098	0.0000	1.0000
F81IG	60	-1547.0404	3227.4384	0.0000	1.0000
K2PIG	58	-1551.2358	3230.9152	0.0000	1.0000
SYMIG	62	-1547.9256	3234.1590	0.0000	1.0000
K2P	56	-1561.8909	3247.3470	0.0000	1.0000
F81	58	-1560.4738	3249.3911	0.0000	1.0000
SYM	60	-1558.5111	3250.3799	0.0000	1.0000
JC69G	56	-1571.3240	3266.2132	0.0000	1.0000
JC69I	56	-1571.7524	3267.0700	0.0000	1.0000
JC69IG	57	-1575.1342	3276.2684	0.0000	1.0000
JC69	55	-1584.9080	3290.9553	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1515.6278	3159.2556	0.4109	0.4109
GTRIG	65	-1515.1233	3160.2466	0.2503	0.6612
GTRI	64	-1516.2380	3160.4761	0.2232	0.8844
HKYG	60	-1521.3366	3162.6732	0.0744	0.9588
HKYI	60	-1521.9622	3163.9243	0.0398	0.9986
HKYIG	61	-1524.2932	3170.5864	0.0014	1.0000
GTR	63	-1529.2108	3184.4216	0.0000	1.0000
HKY	59	-1534.8862	3187.7724	0.0000	1.0000
K2PG	57	-1548.1028	3210.2056	0.0000	1.0000
K2PI	57	-1548.5428	3211.0856	0.0000	1.0000
SYMG	61	-1544.8123	3211.6246	0.0000	1.0000
SYMI	61	-1545.2908	3212.5817	0.0000	1.0000
F81G	59	-1547.3474	3212.6948	0.0000	1.0000
F81I	59	-1547.6518	3213.3036	0.0000	1.0000
F81IG	60	-1547.0404	3214.0807	0.0000	1.0000
K2PIG	58	-1551.2358	3218.4715	0.0000	1.0000
SYMIG	62	-1547.9256	3219.8513	0.0000	1.0000
K2P	56	-1561.8909	3235.7818	0.0000	1.0000
F81	58	-1560.4738	3236.9475	0.0000	1.0000
SYM	60	-1558.5111	3237.0223	0.0000	1.0000
JC69G	56	-1571.3240	3254.6480	0.0000	1.0000
JC69I	56	-1571.7524	3255.5048	0.0000	1.0000
JC69IG	57	-1575.1342	3264.2684	0.0000	1.0000
JC69	55	-1584.9080	3279.8161	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1521.3366	3427.3823	0.6500	0.6500
HKYI	60	-1521.9622	3428.6334	0.3477	0.9978
HKYIG	61	-1524.2932	3439.7073	0.0014	0.9991
GTRG	64	-1515.6278	3441.6119	0.0005	0.9997
GTRI	64	-1516.2380	3442.8324	0.0003	0.9999
GTRIG	65	-1515.1233	3447.0148	0.0000	1.0000
HKY	59	-1534.8862	3448.0697	0.0000	1.0000
K2PG	57	-1548.1028	3461.6792	0.0000	1.0000
GTR	63	-1529.2108	3462.3662	0.0000	1.0000
K2PI	57	-1548.5428	3462.5592	0.0000	1.0000
F81G	59	-1547.3474	3472.9921	0.0000	1.0000
F81I	59	-1547.6518	3473.6008	0.0000	1.0000
K2PIG	58	-1551.2358	3474.3570	0.0000	1.0000
F81IG	60	-1547.0404	3478.7898	0.0000	1.0000
SYMG	61	-1544.8123	3480.7455	0.0000	1.0000
SYMI	61	-1545.2908	3481.7026	0.0000	1.0000
K2P	56	-1561.8909	3482.8436	0.0000	1.0000
F81	58	-1560.4738	3492.8330	0.0000	1.0000
SYMIG	62	-1547.9256	3493.3840	0.0000	1.0000
JC69G	56	-1571.3240	3501.7098	0.0000	1.0000
SYM	60	-1558.5111	3501.7314	0.0000	1.0000
JC69I	56	-1571.7524	3502.5666	0.0000	1.0000
JC69IG	57	-1575.1342	3515.7420	0.0000	1.0000
JC69	55	-1584.9080	3522.4661	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/13347_hg18.phylip = 1 - 609;
 Partition Dummy = 1:261-aggressive-phylip/13347_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/13347_hg18.phylip = 1 - 609;
 Partition Dummy = 1:261-aggressive-phylip/13347_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/13347_hg18.phylip = 1 - 609;
 Partition Dummy = 1:261-aggressive-phylip/13347_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

