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

Input data from file "261-aggressive-phylip/10688_hg18.phylip" (29 taxa, 418 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	-1481.1692	3102.8426	0.6486	0.6486
GTRG	64	-1476.6577	3104.8848	0.2336	0.8822
GTRIG	65	-1476.0427	3106.4604	0.1063	0.9885
HKYIG	61	-1484.0139	3111.2749	0.0096	0.9981
GTRI	64	-1482.1653	3115.8999	0.0009	0.9990
HKYI	60	-1488.1960	3116.8962	0.0006	0.9996
K2PG	57	-1493.0403	3118.4473	0.0003	0.9998
K2PIG	58	-1492.5002	3120.0644	0.0001	1.0000
K2PI	57	-1495.2541	3122.8749	0.0000	1.0000
SYMG	61	-1490.9988	3125.2447	0.0000	1.0000
SYMIG	62	-1490.4768	3126.9592	0.0000	1.0000
SYMI	61	-1493.3739	3129.9951	0.0000	1.0000
HKY	59	-1504.8582	3147.4929	0.0000	1.0000
GTR	63	-1502.9588	3154.6973	0.0000	1.0000
K2P	56	-1515.7297	3161.1436	0.0000	1.0000
SYM	60	-1515.2461	3170.9964	0.0000	1.0000
F81G	59	-1520.5988	3178.9741	0.0000	1.0000
F81IG	60	-1520.0905	3180.6852	0.0000	1.0000
F81I	59	-1523.1285	3184.0335	0.0000	1.0000
JC69G	56	-1529.7704	3189.2251	0.0000	1.0000
JC69IG	57	-1529.2667	3190.9001	0.0000	1.0000
JC69I	56	-1532.1525	3193.9892	0.0000	1.0000
F81	58	-1545.9070	3226.8780	0.0000	1.0000
JC69	55	-1554.3478	3235.7122	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1476.6577	3081.3154	0.4352	0.4352
GTRIG	65	-1476.0427	3082.0854	0.2962	0.7314
HKYG	60	-1481.1692	3082.3384	0.2610	0.9923
HKYIG	61	-1484.0139	3090.0277	0.0056	0.9979
GTRI	64	-1482.1653	3092.3305	0.0018	0.9997
HKYI	60	-1488.1960	3096.3920	0.0002	0.9999
K2PG	57	-1493.0403	3100.0807	0.0000	1.0000
K2PIG	58	-1492.5002	3101.0003	0.0000	1.0000
SYMG	61	-1490.9988	3103.9975	0.0000	1.0000
K2PI	57	-1495.2541	3104.5082	0.0000	1.0000
SYMIG	62	-1490.4768	3104.9536	0.0000	1.0000
SYMI	61	-1493.3739	3108.7479	0.0000	1.0000
HKY	59	-1504.8582	3127.7163	0.0000	1.0000
GTR	63	-1502.9588	3131.9176	0.0000	1.0000
K2P	56	-1515.7297	3143.4593	0.0000	1.0000
SYM	60	-1515.2461	3150.4922	0.0000	1.0000
F81G	59	-1520.5988	3159.1975	0.0000	1.0000
F81IG	60	-1520.0905	3160.1810	0.0000	1.0000
F81I	59	-1523.1285	3164.2570	0.0000	1.0000
JC69G	56	-1529.7704	3171.5409	0.0000	1.0000
JC69IG	57	-1529.2667	3172.5334	0.0000	1.0000
JC69I	56	-1532.1525	3176.3050	0.0000	1.0000
F81	58	-1545.9070	3207.8140	0.0000	1.0000
JC69	55	-1554.3478	3218.6957	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1481.1692	3324.4672	0.9297	0.9297
K2PG	57	-1493.0403	3330.1031	0.0555	0.9853
K2PI	57	-1495.2541	3334.5307	0.0061	0.9913
K2PIG	58	-1492.5002	3335.0582	0.0047	0.9960
HKYIG	61	-1484.0139	3336.1921	0.0026	0.9986
HKYI	60	-1488.1960	3338.5209	0.0008	0.9995
GTRG	64	-1476.6577	3339.5862	0.0005	1.0000
GTRIG	65	-1476.0427	3344.3917	0.0000	1.0000
SYMG	61	-1490.9988	3350.1619	0.0000	1.0000
GTRI	64	-1482.1653	3350.6014	0.0000	1.0000
SYMI	61	-1493.3739	3354.9122	0.0000	1.0000
SYMIG	62	-1490.4768	3355.1534	0.0000	1.0000
HKY	59	-1504.8582	3365.8097	0.0000	1.0000
K2P	56	-1515.7297	3369.4463	0.0000	1.0000
GTR	63	-1502.9588	3386.1530	0.0000	1.0000
SYM	60	-1515.2461	3392.6211	0.0000	1.0000
F81G	59	-1520.5988	3397.2909	0.0000	1.0000
JC69G	56	-1529.7704	3397.5278	0.0000	1.0000
JC69I	56	-1532.1525	3402.2919	0.0000	1.0000
F81IG	60	-1520.0905	3402.3099	0.0000	1.0000
F81I	59	-1523.1285	3402.3504	0.0000	1.0000
JC69IG	57	-1529.2667	3402.5559	0.0000	1.0000
JC69	55	-1554.3478	3440.6471	0.0000	1.0000
F81	58	-1545.9070	3441.8719	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/10688_hg18.phylip = 1 - 418;
 Partition Dummy = 1:261-aggressive-phylip/10688_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/10688_hg18.phylip = 1 - 418;
 Partition Dummy = 1:261-aggressive-phylip/10688_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/10688_hg18.phylip = 1 - 418;
 Partition Dummy = 1:261-aggressive-phylip/10688_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

