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

Input data from file "261-aggressive-phylip/9235_hg18.phylip" (29 taxa, 688 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	-2463.5778	5068.5104	0.8772	0.8772
GTRIG	65	-2464.3817	5072.5577	0.1159	0.9931
GTRI	64	-2468.4291	5078.2129	0.0069	1.0000
HKYG	60	-2486.3766	5104.4278	0.0000	1.0000
HKYIG	61	-2485.8219	5105.7268	0.0000	1.0000
HKYI	60	-2491.3834	5114.4414	0.0000	1.0000
SYMG	61	-2513.6948	5161.4727	0.0000	1.0000
SYMIG	62	-2513.1285	5162.7562	0.0000	1.0000
K2PG	57	-2521.9268	5168.3488	0.0000	1.0000
K2PIG	58	-2521.0710	5169.0227	0.0000	1.0000
SYMI	61	-2518.0992	5170.2814	0.0000	1.0000
K2PI	57	-2524.4210	5173.3372	0.0000	1.0000
GTR	63	-2527.3400	5193.6030	0.0000	1.0000
HKY	59	-2549.1336	5227.5410	0.0000	1.0000
SYM	60	-2568.3324	5268.3395	0.0000	1.0000
K2P	56	-2578.4401	5278.9974	0.0000	1.0000
F81G	59	-2579.3975	5288.0688	0.0000	1.0000
F81IG	60	-2578.8981	5289.4709	0.0000	1.0000
F81I	59	-2583.7201	5296.7141	0.0000	1.0000
JC69G	56	-2603.6818	5329.4808	0.0000	1.0000
JC69IG	57	-2603.0148	5330.5248	0.0000	1.0000
JC69I	56	-2604.7726	5331.6625	0.0000	1.0000
F81	58	-2636.1201	5399.1209	0.0000	1.0000
JC69	55	-2656.1994	5432.1456	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2463.5778	5055.1557	0.8529	0.8529
GTRIG	65	-2464.3817	5058.7634	0.1404	0.9933
GTRI	64	-2468.4291	5064.8582	0.0067	1.0000
HKYG	60	-2486.3766	5092.7532	0.0000	1.0000
HKYIG	61	-2485.8219	5093.6437	0.0000	1.0000
HKYI	60	-2491.3834	5102.7667	0.0000	1.0000
SYMG	61	-2513.6948	5149.3896	0.0000	1.0000
SYMIG	62	-2513.1285	5150.2570	0.0000	1.0000
K2PG	57	-2521.9268	5157.8535	0.0000	1.0000
K2PIG	58	-2521.0710	5158.1419	0.0000	1.0000
SYMI	61	-2518.0992	5158.1983	0.0000	1.0000
K2PI	57	-2524.4210	5162.8419	0.0000	1.0000
GTR	63	-2527.3400	5180.6799	0.0000	1.0000
HKY	59	-2549.1336	5216.2671	0.0000	1.0000
SYM	60	-2568.3324	5256.6649	0.0000	1.0000
K2P	56	-2578.4401	5268.8801	0.0000	1.0000
F81G	59	-2579.3975	5276.7950	0.0000	1.0000
F81IG	60	-2578.8981	5277.7963	0.0000	1.0000
F81I	59	-2583.7201	5285.4402	0.0000	1.0000
JC69G	56	-2603.6818	5319.3635	0.0000	1.0000
JC69IG	57	-2603.0148	5320.0295	0.0000	1.0000
JC69I	56	-2604.7726	5321.5452	0.0000	1.0000
F81	58	-2636.1201	5388.2402	0.0000	1.0000
JC69	55	-2656.1994	5422.3988	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-2463.5778	5345.3182	0.9757	0.9757
GTRIG	65	-2464.3817	5353.4597	0.0166	0.9923
GTRI	64	-2468.4291	5355.0206	0.0076	0.9999
HKYG	60	-2486.3766	5364.7805	0.0001	1.0000
HKYIG	61	-2485.8219	5370.2049	0.0000	1.0000
HKYI	60	-2491.3834	5374.7941	0.0000	1.0000
K2PG	57	-2521.9268	5416.2795	0.0000	1.0000
K2PIG	58	-2521.0710	5421.1017	0.0000	1.0000
K2PI	57	-2524.4210	5421.2679	0.0000	1.0000
SYMG	61	-2513.6948	5425.9508	0.0000	1.0000
SYMIG	62	-2513.1285	5431.3519	0.0000	1.0000
SYMI	61	-2518.0992	5434.7595	0.0000	1.0000
GTR	63	-2527.3400	5466.3086	0.0000	1.0000
HKY	59	-2549.1336	5483.7607	0.0000	1.0000
K2P	56	-2578.4401	5522.7723	0.0000	1.0000
SYM	60	-2568.3324	5528.6922	0.0000	1.0000
F81G	59	-2579.3975	5544.2885	0.0000	1.0000
F81IG	60	-2578.8981	5549.8236	0.0000	1.0000
F81I	59	-2583.7201	5552.9338	0.0000	1.0000
JC69G	56	-2603.6818	5573.2557	0.0000	1.0000
JC69I	56	-2604.7726	5575.4374	0.0000	1.0000
JC69IG	57	-2603.0148	5578.4555	0.0000	1.0000
F81	58	-2636.1201	5651.1999	0.0000	1.0000
JC69	55	-2656.1994	5671.7571	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/9235_hg18.phylip = 1 - 688;
 Partition Dummy = 1:261-aggressive-phylip/9235_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/9235_hg18.phylip = 1 - 688;
 Partition Dummy = 1:261-aggressive-phylip/9235_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/9235_hg18.phylip = 1 - 688;
 Partition Dummy = 1:261-aggressive-phylip/9235_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

