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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3004.2842	6048.5684	0.5425	0.5425
GTRIG	21	-3004.2483	6050.4966	0.2069	0.7493
GTRI	20	-3005.8993	6051.7986	0.1079	0.8572
HKYG	16	-3010.1894	6052.3787	0.0807	0.9380
HKYIG	17	-3009.9516	6053.9031	0.0377	0.9756
HKYI	16	-3011.3868	6054.7736	0.0244	1.0000
GTR	19	-3035.2622	6108.5243	0.0000	1.0000
HKY	15	-3041.1294	6112.2588	0.0000	1.0000
K2PG	13	-3061.5231	6149.0463	0.0000	1.0000
K2PI	13	-3062.4433	6150.8866	0.0000	1.0000
K2PIG	14	-3061.9815	6151.9630	0.0000	1.0000
SYMG	17	-3060.0674	6154.1348	0.0000	1.0000
SYMI	17	-3060.7752	6155.5503	0.0000	1.0000
SYMIG	18	-3060.0522	6156.1045	0.0000	1.0000
K2P	12	-3093.0228	6210.0456	0.0000	1.0000
SYM	16	-3090.5287	6213.0573	0.0000	1.0000
F81G	15	-3092.7399	6215.4798	0.0000	1.0000
F81I	15	-3093.1793	6216.3585	0.0000	1.0000
F81IG	16	-3092.7370	6217.4739	0.0000	1.0000
F81	14	-3119.3551	6266.7103	0.0000	1.0000
JC69G	12	-3134.2053	6292.4105	0.0000	1.0000
JC69I	12	-3134.5448	6293.0895	0.0000	1.0000
JC69IG	13	-3134.2070	6294.4140	0.0000	1.0000
JC69	11	-3161.9553	6345.9107	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3010.1894	6130.0653	0.7442	0.7442
HKYI	16	-3011.3868	6132.4601	0.2247	0.9690
HKYIG	17	-3009.9516	6136.4451	0.0306	0.9996
GTRG	20	-3004.2842	6145.6766	0.0003	0.9999
GTRI	20	-3005.8993	6148.9068	0.0001	1.0000
GTRIG	21	-3004.2483	6152.4602	0.0000	1.0000
HKY	15	-3041.1294	6185.0900	0.0000	1.0000
GTR	19	-3035.2622	6200.7771	0.0000	1.0000
K2PG	13	-3061.5231	6212.1666	0.0000	1.0000
K2PI	13	-3062.4433	6214.0069	0.0000	1.0000
K2PIG	14	-3061.9815	6219.9387	0.0000	1.0000
SYMG	17	-3060.0674	6236.6767	0.0000	1.0000
SYMI	17	-3060.7752	6238.0923	0.0000	1.0000
SYMIG	18	-3060.0522	6243.5018	0.0000	1.0000
K2P	12	-3093.0228	6268.3105	0.0000	1.0000
F81G	15	-3092.7399	6288.3109	0.0000	1.0000
F81I	15	-3093.1793	6289.1896	0.0000	1.0000
SYM	16	-3090.5287	6290.7439	0.0000	1.0000
F81IG	16	-3092.7370	6295.1605	0.0000	1.0000
F81	14	-3119.3551	6334.6860	0.0000	1.0000
JC69G	12	-3134.2053	6350.6754	0.0000	1.0000
JC69I	12	-3134.5448	6351.3544	0.0000	1.0000
JC69IG	13	-3134.2070	6357.5343	0.0000	1.0000
JC69	11	-3161.9553	6399.3202	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 683-tomodel-phylip/chr3_245.phylip = 1 - 949;
 Partition Dummy = 1:683-tomodel-phylip/chr3_245.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 683-tomodel-phylip/chr3_245.phylip = 1 - 949;
 Partition Dummy = 1:683-tomodel-phylip/chr3_245.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 683-tomodel-phylip/chr3_245.phylip = 1 - 949;
 Partition Dummy = 1:683-tomodel-phylip/chr3_245.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

