
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:20 04/04/2011

Input data from file "683-tomodel-phylip/chr8_7536.phylip" (7 taxa, 904 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2737.4299	5506.8599	0.4328	0.4328
HKYI	16	-2737.8028	5507.6056	0.2981	0.7309
HKYIG	17	-2737.4233	5508.8465	0.1603	0.8912
GTRG	20	-2735.5558	5511.1116	0.0516	0.9428
GTRI	20	-2735.8660	5511.7319	0.0379	0.9807
GTRIG	21	-2735.5390	5513.0780	0.0193	1.0000
HKY	15	-2754.2996	5538.5992	0.0000	1.0000
GTR	19	-2752.6548	5543.3097	0.0000	1.0000
K2PG	13	-2817.6683	5661.3366	0.0000	1.0000
K2PI	13	-2818.6480	5663.2959	0.0000	1.0000
K2PIG	14	-2817.6831	5663.3662	0.0000	1.0000
SYMG	17	-2815.5249	5665.0498	0.0000	1.0000
SYMIG	18	-2815.5353	5667.0707	0.0000	1.0000
SYMI	17	-2816.8092	5667.6184	0.0000	1.0000
F81G	15	-2832.7180	5695.4360	0.0000	1.0000
F81I	15	-2832.9413	5695.8825	0.0000	1.0000
F81IG	16	-2832.6971	5697.3942	0.0000	1.0000
K2P	12	-2837.6052	5699.2104	0.0000	1.0000
SYM	16	-2836.2712	5704.5424	0.0000	1.0000
F81	14	-2847.9302	5723.8604	0.0000	1.0000
JC69G	12	-2898.7484	5821.4969	0.0000	1.0000
JC69I	12	-2899.4323	5822.8646	0.0000	1.0000
JC69IG	13	-2898.7537	5823.5074	0.0000	1.0000
JC69	11	-2915.8363	5853.6726	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2737.4299	5583.7691	0.5806	0.5806
HKYI	16	-2737.8028	5584.5149	0.3999	0.9806
HKYIG	17	-2737.4233	5590.5626	0.0194	1.0000
GTRG	20	-2735.5558	5607.2482	0.0000	1.0000
GTRI	20	-2735.8660	5607.8685	0.0000	1.0000
HKY	15	-2754.2996	5610.7017	0.0000	1.0000
GTRIG	21	-2735.5390	5614.0214	0.0000	1.0000
GTR	19	-2752.6548	5634.6395	0.0000	1.0000
K2PG	13	-2817.6683	5723.8254	0.0000	1.0000
K2PI	13	-2818.6480	5725.7847	0.0000	1.0000
K2PIG	14	-2817.6831	5730.6618	0.0000	1.0000
SYMG	17	-2815.5249	5746.7659	0.0000	1.0000
SYMI	17	-2816.8092	5749.3345	0.0000	1.0000
SYMIG	18	-2815.5353	5753.5936	0.0000	1.0000
K2P	12	-2837.6052	5756.8924	0.0000	1.0000
F81G	15	-2832.7180	5767.5384	0.0000	1.0000
F81I	15	-2832.9413	5767.9850	0.0000	1.0000
F81IG	16	-2832.6971	5774.3035	0.0000	1.0000
SYM	16	-2836.2712	5781.4516	0.0000	1.0000
F81	14	-2847.9302	5791.1560	0.0000	1.0000
JC69G	12	-2898.7484	5879.1788	0.0000	1.0000
JC69I	12	-2899.4323	5880.5465	0.0000	1.0000
JC69IG	13	-2898.7537	5885.9962	0.0000	1.0000
JC69	11	-2915.8363	5906.5477	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 683-tomodel-phylip/chr8_7536.phylip = 1 - 904;
 Partition Dummy = 1:683-tomodel-phylip/chr8_7536.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 683-tomodel-phylip/chr8_7536.phylip = 1 - 904;
 Partition Dummy = 1:683-tomodel-phylip/chr8_7536.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 683-tomodel-phylip/chr8_7536.phylip = 1 - 904;
 Partition Dummy = 1:683-tomodel-phylip/chr8_7536.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

