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

Input data from file "261-aggressive-phylip/7863_hg18.phylip" (29 taxa, 811 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	-2453.6548	5046.4624	0.7646	0.7646
GTRIG	65	-2453.6619	5048.8406	0.2328	0.9974
HKYG	60	-2464.3133	5058.3866	0.0020	0.9994
HKYIG	61	-2464.3206	5060.7401	0.0006	1.0000
K2PG	57	-2499.1291	5121.0390	0.0000	1.0000
SYMG	61	-2495.8539	5123.8066	0.0000	1.0000
K2PIG	58	-2501.4809	5128.0629	0.0000	1.0000
SYMIG	62	-2498.0624	5130.5686	0.0000	1.0000
F81G	59	-2513.7984	5155.0243	0.0000	1.0000
F81IG	60	-2515.2683	5160.2966	0.0000	1.0000
JC69G	56	-2544.0027	5208.4722	0.0000	1.0000
JC69IG	57	-2544.0101	5210.8011	0.0000	1.0000
GTRI	64	-2539.3499	5217.8527	0.0000	1.0000
GTR	63	-2543.7455	5224.2862	0.0000	1.0000
HKYI	60	-2550.7040	5231.1680	0.0000	1.0000
HKY	59	-2554.0044	5235.4363	0.0000	1.0000
SYM	60	-2584.7693	5299.2987	0.0000	1.0000
K2P	56	-2589.7245	5299.9159	0.0000	1.0000
SYMI	61	-2584.8184	5301.7355	0.0000	1.0000
K2PI	57	-2589.6635	5302.1078	0.0000	1.0000
F81	58	-2602.7417	5330.5844	0.0000	1.0000
F81I	59	-2602.9024	5333.2322	0.0000	1.0000
JC69	55	-2631.5625	5381.2840	0.0000	1.0000
JC69I	56	-2631.2823	5383.0315	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2453.6548	5035.3096	0.7315	0.7315
GTRIG	65	-2453.6619	5037.3239	0.2672	0.9987
HKYG	60	-2464.3133	5048.6266	0.0009	0.9997
HKYIG	61	-2464.3206	5050.6413	0.0003	1.0000
K2PG	57	-2499.1291	5112.2581	0.0000	1.0000
SYMG	61	-2495.8539	5113.7078	0.0000	1.0000
K2PIG	58	-2501.4809	5118.9618	0.0000	1.0000
SYMIG	62	-2498.0624	5120.1247	0.0000	1.0000
F81G	59	-2513.7984	5145.5969	0.0000	1.0000
F81IG	60	-2515.2683	5150.5366	0.0000	1.0000
JC69G	56	-2544.0027	5200.0054	0.0000	1.0000
JC69IG	57	-2544.0101	5202.0202	0.0000	1.0000
GTRI	64	-2539.3499	5206.6999	0.0000	1.0000
GTR	63	-2543.7455	5213.4910	0.0000	1.0000
HKYI	60	-2550.7040	5221.4080	0.0000	1.0000
HKY	59	-2554.0044	5226.0089	0.0000	1.0000
SYM	60	-2584.7693	5289.5387	0.0000	1.0000
K2P	56	-2589.7245	5291.4491	0.0000	1.0000
SYMI	61	-2584.8184	5291.6367	0.0000	1.0000
K2PI	57	-2589.6635	5293.3269	0.0000	1.0000
F81	58	-2602.7417	5321.4833	0.0000	1.0000
F81I	59	-2602.9024	5323.8047	0.0000	1.0000
JC69	55	-2631.5625	5373.1250	0.0000	1.0000
JC69I	56	-2631.2823	5374.5647	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2464.3133	5330.5227	0.9076	0.9076
GTRG	64	-2453.6548	5335.9987	0.0587	0.9663
HKYIG	61	-2464.3206	5337.2356	0.0316	0.9980
GTRIG	65	-2453.6619	5342.7113	0.0020	1.0000
K2PG	57	-2499.1291	5380.0594	0.0000	1.0000
K2PIG	58	-2501.4809	5391.4614	0.0000	1.0000
SYMG	61	-2495.8539	5400.3022	0.0000	1.0000
SYMIG	62	-2498.0624	5411.4173	0.0000	1.0000
F81G	59	-2513.7984	5422.7947	0.0000	1.0000
F81IG	60	-2515.2683	5432.4327	0.0000	1.0000
JC69G	56	-2544.0027	5463.1084	0.0000	1.0000
JC69IG	57	-2544.0101	5469.8215	0.0000	1.0000
HKY	59	-2554.0044	5503.2067	0.0000	1.0000
HKYI	60	-2550.7040	5503.3041	0.0000	1.0000
GTRI	64	-2539.3499	5507.3890	0.0000	1.0000
GTR	63	-2543.7455	5509.4819	0.0000	1.0000
K2P	56	-2589.7245	5554.5521	0.0000	1.0000
K2PI	57	-2589.6635	5561.1282	0.0000	1.0000
SYM	60	-2584.7693	5571.4347	0.0000	1.0000
SYMI	61	-2584.8184	5578.2311	0.0000	1.0000
F81	58	-2602.7417	5593.9829	0.0000	1.0000
F81I	59	-2602.9024	5601.0025	0.0000	1.0000
JC69	55	-2631.5625	5631.5298	0.0000	1.0000
JC69I	56	-2631.2823	5637.6677	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/7863_hg18.phylip = 1 - 811;
 Partition Dummy = 1:261-aggressive-phylip/7863_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/7863_hg18.phylip = 1 - 811;
 Partition Dummy = 1:261-aggressive-phylip/7863_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/7863_hg18.phylip = 1 - 811;
 Partition Dummy = 1:261-aggressive-phylip/7863_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

