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

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

Minimum AIC  model: HKYG
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	-2610.1401	5352.1057	0.8234	0.8234
GTRG	64	-2607.0221	5355.5728	0.1455	0.9689
GTRIG	65	-2607.4295	5358.8329	0.0285	0.9974
HKYIG	61	-2614.7034	5363.6463	0.0026	1.0000
HKYI	60	-2620.2928	5372.4111	0.0000	1.0000
GTRI	64	-2617.7128	5376.9541	0.0000	1.0000
K2PG	57	-2628.0428	5380.7158	0.0000	1.0000
K2PIG	58	-2627.9863	5382.9935	0.0000	1.0000
SYMG	61	-2627.3709	5388.9813	0.0000	1.0000
SYMIG	62	-2627.1576	5390.9764	0.0000	1.0000
K2PI	57	-2637.6771	5399.9844	0.0000	1.0000
SYMI	61	-2636.9935	5408.2265	0.0000	1.0000
HKY	59	-2670.4134	5470.2461	0.0000	1.0000
GTR	63	-2669.3879	5477.8668	0.0000	1.0000
K2P	56	-2687.8996	5498.0464	0.0000	1.0000
SYM	60	-2687.0734	5505.9723	0.0000	1.0000
F81G	59	-2699.7805	5528.9804	0.0000	1.0000
F81IG	60	-2700.1979	5532.2213	0.0000	1.0000
F81I	59	-2708.9760	5547.3714	0.0000	1.0000
JC69G	56	-2713.0842	5548.4156	0.0000	1.0000
JC69IG	57	-2713.2496	5551.1294	0.0000	1.0000
JC69I	56	-2722.0853	5566.4178	0.0000	1.0000
F81	58	-2756.6471	5640.3151	0.0000	1.0000
JC69	55	-2769.6024	5659.0767	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-2610.1401	5340.2801	0.6583	0.6583
GTRG	64	-2607.0221	5342.0443	0.2725	0.9307
GTRIG	65	-2607.4295	5344.8589	0.0667	0.9974
HKYIG	61	-2614.7034	5351.4069	0.0025	1.0000
HKYI	60	-2620.2928	5360.5855	0.0000	1.0000
GTRI	64	-2617.7128	5363.4256	0.0000	1.0000
K2PG	57	-2628.0428	5370.0856	0.0000	1.0000
K2PIG	58	-2627.9863	5371.9725	0.0000	1.0000
SYMG	61	-2627.3709	5376.7418	0.0000	1.0000
SYMIG	62	-2627.1576	5378.3152	0.0000	1.0000
K2PI	57	-2637.6771	5389.3541	0.0000	1.0000
SYMI	61	-2636.9935	5395.9870	0.0000	1.0000
HKY	59	-2670.4134	5458.8268	0.0000	1.0000
GTR	63	-2669.3879	5464.7759	0.0000	1.0000
K2P	56	-2687.8996	5487.7992	0.0000	1.0000
SYM	60	-2687.0734	5494.1468	0.0000	1.0000
F81G	59	-2699.7805	5517.5610	0.0000	1.0000
F81IG	60	-2700.1979	5520.3958	0.0000	1.0000
F81I	59	-2708.9760	5535.9521	0.0000	1.0000
JC69G	56	-2713.0842	5538.1684	0.0000	1.0000
JC69IG	57	-2713.2496	5540.4992	0.0000	1.0000
JC69I	56	-2722.0853	5556.1706	0.0000	1.0000
F81	58	-2756.6471	5629.2942	0.0000	1.0000
JC69	55	-2769.6024	5649.2049	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2610.1401	5611.6057	0.9992	0.9992
HKYIG	61	-2614.7034	5627.2545	0.0004	0.9996
K2PG	57	-2628.0428	5627.8449	0.0003	0.9999
GTRG	64	-2607.0221	5631.4582	0.0000	0.9999
HKYI	60	-2620.2928	5631.9111	0.0000	1.0000
K2PIG	58	-2627.9863	5634.2539	0.0000	1.0000
GTRIG	65	-2607.4295	5638.7950	0.0000	1.0000
K2PI	57	-2637.6771	5647.1134	0.0000	1.0000
SYMG	61	-2627.3709	5652.5895	0.0000	1.0000
GTRI	64	-2617.7128	5652.8396	0.0000	1.0000
SYMIG	62	-2627.1576	5658.6849	0.0000	1.0000
SYMI	61	-2636.9935	5671.8347	0.0000	1.0000
HKY	59	-2670.4134	5725.6302	0.0000	1.0000
K2P	56	-2687.8996	5741.0364	0.0000	1.0000
GTR	63	-2669.3879	5749.6677	0.0000	1.0000
SYM	60	-2687.0734	5765.4723	0.0000	1.0000
F81G	59	-2699.7805	5784.3645	0.0000	1.0000
JC69G	56	-2713.0842	5791.4056	0.0000	1.0000
F81IG	60	-2700.1979	5791.7214	0.0000	1.0000
JC69IG	57	-2713.2496	5798.2585	0.0000	1.0000
F81I	59	-2708.9760	5802.7556	0.0000	1.0000
JC69I	56	-2722.0853	5809.4078	0.0000	1.0000
F81	58	-2756.6471	5891.5756	0.0000	1.0000
JC69	55	-2769.6024	5897.9200	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7771_hg18.phylip = 1 - 680;
 Partition Dummy = 1:261-aggressive-phylip/7771_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 AICc model (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7771_hg18.phylip = 1 - 680;
 Partition Dummy = 1:261-aggressive-phylip/7771_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/7771_hg18.phylip = 1 - 680;
 Partition Dummy = 1:261-aggressive-phylip/7771_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

