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

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

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


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	-1269.8750	2682.5280	0.3410	0.3410
GTRI	64	-1269.9922	2682.7623	0.3033	0.6444
GTR	63	-1271.8436	2683.9852	0.1646	0.8089
GTRIG	65	-1269.8677	2685.0023	0.0990	0.9079
HKYG	60	-1277.1098	2687.1296	0.0342	0.9421
HKYI	60	-1277.2452	2687.4005	0.0298	0.9719
HKY	59	-1278.9759	2688.4165	0.0180	0.9898
HKYIG	61	-1277.0990	2689.5620	0.0101	1.0000
F81G	59	-1286.4134	2703.2915	0.0000	1.0000
F81I	59	-1286.5660	2703.5967	0.0000	1.0000
F81	58	-1288.1657	2704.3596	0.0000	1.0000
F81IG	60	-1286.4033	2705.7167	0.0000	1.0000
K2PG	57	-1310.4559	2746.5117	0.0000	1.0000
K2PI	57	-1310.5979	2746.7959	0.0000	1.0000
SYMG	61	-1305.9016	2747.1672	0.0000	1.0000
SYMI	61	-1306.0227	2747.4094	0.0000	1.0000
K2P	56	-1312.4393	2748.0591	0.0000	1.0000
SYM	60	-1307.7746	2748.4592	0.0000	1.0000
K2PIG	58	-1310.4392	2748.9066	0.0000	1.0000
SYMIG	62	-1305.8913	2749.6091	0.0000	1.0000
JC69G	56	-1318.2137	2759.6077	0.0000	1.0000
JC69I	56	-1318.3683	2759.9169	0.0000	1.0000
JC69	55	-1320.0207	2760.8106	0.0000	1.0000
JC69IG	57	-1318.2006	2762.0012	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1269.8750	2667.7500	0.3644	0.3644
GTRI	64	-1269.9922	2667.9843	0.3241	0.6884
GTR	63	-1271.8436	2669.6873	0.1383	0.8267
GTRIG	65	-1269.8677	2669.7354	0.1350	0.9618
HKYG	60	-1277.1098	2674.2195	0.0143	0.9761
HKYI	60	-1277.2452	2674.4904	0.0125	0.9886
HKY	59	-1278.9759	2675.9517	0.0060	0.9947
HKYIG	61	-1277.0990	2676.1980	0.0053	1.0000
F81G	59	-1286.4134	2690.8268	0.0000	1.0000
F81I	59	-1286.5660	2691.1319	0.0000	1.0000
F81	58	-1288.1657	2692.3314	0.0000	1.0000
F81IG	60	-1286.4033	2692.8066	0.0000	1.0000
SYMG	61	-1305.9016	2733.8032	0.0000	1.0000
SYMI	61	-1306.0227	2734.0454	0.0000	1.0000
K2PG	57	-1310.4559	2734.9117	0.0000	1.0000
K2PI	57	-1310.5979	2735.1959	0.0000	1.0000
SYM	60	-1307.7746	2735.5492	0.0000	1.0000
SYMIG	62	-1305.8913	2735.7825	0.0000	1.0000
K2PIG	58	-1310.4392	2736.8785	0.0000	1.0000
K2P	56	-1312.4393	2736.8787	0.0000	1.0000
JC69G	56	-1318.2137	2748.4273	0.0000	1.0000
JC69I	56	-1318.3683	2748.7365	0.0000	1.0000
JC69	55	-1320.0207	2750.0413	0.0000	1.0000
JC69IG	57	-1318.2006	2750.4012	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	59	-1278.9759	2938.0616	0.6657	0.6657
HKYG	60	-1277.1098	2940.7719	0.1717	0.8374
HKYI	60	-1277.2452	2941.0428	0.1499	0.9874
HKYIG	61	-1277.0990	2947.1930	0.0069	0.9943
GTR	63	-1271.8436	2949.5673	0.0021	0.9964
F81	58	-1288.1657	2949.9988	0.0017	0.9981
GTRG	64	-1269.8750	2952.0726	0.0006	0.9987
GTRI	64	-1269.9922	2952.3069	0.0005	0.9992
F81G	59	-1286.4134	2952.9366	0.0004	0.9996
F81I	59	-1286.5660	2953.2418	0.0003	1.0000
GTRIG	65	-1269.8677	2958.5006	0.0000	1.0000
F81IG	60	-1286.4033	2959.3590	0.0000	1.0000
K2P	56	-1312.4393	2985.6609	0.0000	1.0000
K2PG	57	-1310.4559	2988.1365	0.0000	1.0000
K2PI	57	-1310.5979	2988.4207	0.0000	1.0000
JC69	55	-1320.0207	2994.3810	0.0000	1.0000
K2PIG	58	-1310.4392	2994.5458	0.0000	1.0000
JC69G	56	-1318.2137	2997.2096	0.0000	1.0000
JC69I	56	-1318.3683	2997.5188	0.0000	1.0000
SYM	60	-1307.7746	3002.1016	0.0000	1.0000
JC69IG	57	-1318.2006	3003.6260	0.0000	1.0000
SYMG	61	-1305.9016	3004.7982	0.0000	1.0000
SYMI	61	-1306.0227	3005.0404	0.0000	1.0000
SYMIG	62	-1305.8913	3011.2200	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/7524_hg18.phylip = 1 - 628;
 Partition Dummy = 1:261-aggressive-phylip/7524_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/7524_hg18.phylip = 1 - 628;
 Partition Dummy = 1:261-aggressive-phylip/7524_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 (HKY)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7524_hg18.phylip = 1 - 628;
 Partition Dummy = 1:261-aggressive-phylip/7524_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

