
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:39 04/04/2011

Input data from file "591-tomodel-phylip/chr3_6135.phylip" (7 taxa, 893 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1768.6482	3577.2964	0.4939	0.4939
GTRIG	21	-1768.5882	3579.1764	0.1929	0.6868
HKYI	16	-1773.6541	3579.3083	0.1806	0.8674
HKYIG	17	-1773.4539	3580.9078	0.0812	0.9486
GTRG	20	-1771.1640	3582.3281	0.0399	0.9885
HKYG	16	-1776.4120	3584.8239	0.0115	0.9999
SYMI	17	-1781.2365	3596.4730	0.0000	0.9999
K2PI	13	-1785.6122	3597.2245	0.0000	1.0000
SYMIG	18	-1781.1586	3598.3173	0.0000	1.0000
K2PIG	14	-1785.4474	3598.8949	0.0000	1.0000
SYMG	17	-1784.0232	3602.0464	0.0000	1.0000
K2PG	13	-1788.2674	3602.5348	0.0000	1.0000
GTR	19	-1794.6023	3627.2047	0.0000	1.0000
HKY	15	-1800.1527	3630.3054	0.0000	1.0000
SYM	16	-1807.7708	3647.5417	0.0000	1.0000
K2P	12	-1811.8334	3647.6667	0.0000	1.0000
F81I	15	-1809.1501	3648.3002	0.0000	1.0000
F81IG	16	-1809.1500	3650.3001	0.0000	1.0000
F81G	15	-1811.4792	3652.9584	0.0000	1.0000
JC69I	12	-1821.6440	3667.2881	0.0000	1.0000
JC69IG	13	-1821.6440	3669.2880	0.0000	1.0000
JC69G	12	-1824.0759	3672.1518	0.0000	1.0000
F81	14	-1834.6301	3697.2603	0.0000	1.0000
JC69	11	-1847.3637	3716.7273	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1773.6541	3656.0217	0.7726	0.7726
K2PI	13	-1785.6122	3659.5541	0.1321	0.9047
HKYG	16	-1776.4120	3661.5373	0.0490	0.9537
HKYIG	17	-1773.4539	3662.4158	0.0316	0.9853
K2PG	13	-1788.2674	3664.8644	0.0093	0.9946
K2PIG	14	-1785.4474	3666.0191	0.0052	0.9998
GTRI	20	-1768.6482	3673.1881	0.0001	1.0000
SYMI	17	-1781.2365	3677.9810	0.0000	1.0000
GTRG	20	-1771.1640	3678.2198	0.0000	1.0000
GTRIG	21	-1768.5882	3679.8627	0.0000	1.0000
SYMG	17	-1784.0232	3683.5543	0.0000	1.0000
SYMIG	18	-1781.1586	3684.6198	0.0000	1.0000
HKY	15	-1800.1527	3702.2242	0.0000	1.0000
K2P	12	-1811.8334	3705.2018	0.0000	1.0000
GTR	19	-1794.6023	3718.3018	0.0000	1.0000
F81I	15	-1809.1501	3720.2190	0.0000	1.0000
SYM	16	-1807.7708	3724.2551	0.0000	1.0000
JC69I	12	-1821.6440	3724.8231	0.0000	1.0000
F81G	15	-1811.4792	3724.8772	0.0000	1.0000
F81IG	16	-1809.1500	3727.0134	0.0000	1.0000
JC69G	12	-1824.0759	3729.6868	0.0000	1.0000
JC69IG	13	-1821.6440	3731.6176	0.0000	1.0000
F81	14	-1834.6301	3764.3845	0.0000	1.0000
JC69	11	-1847.3637	3769.4678	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_6135.phylip = 1 - 893;
 Partition Dummy = 1:591-tomodel-phylip/chr3_6135.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_6135.phylip = 1 - 893;
 Partition Dummy = 1:591-tomodel-phylip/chr3_6135.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_6135.phylip = 1 - 893;
 Partition Dummy = 1:591-tomodel-phylip/chr3_6135.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

