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

Input data from file "591-tomodel-phylip/chr3_496.phylip" (7 taxa, 995 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	-3371.6095	6783.2189	0.5829	0.5829
GTRIG	21	-3371.6124	6785.2248	0.2138	0.7968
HKYG	16	-3377.1955	6786.3911	0.1193	0.9161
HKYIG	17	-3377.1958	6788.3915	0.0439	0.9600
GTRI	20	-3374.4328	6788.8656	0.0346	0.9946
HKYI	16	-3380.2996	6792.5993	0.0054	1.0000
K2PG	13	-3401.2219	6828.4438	0.0000	1.0000
K2PIG	14	-3401.2312	6830.4624	0.0000	1.0000
SYMG	17	-3400.2377	6834.4753	0.0000	1.0000
K2PI	13	-3404.7302	6835.4605	0.0000	1.0000
GTR	19	-3398.9792	6835.9585	0.0000	1.0000
SYMIG	18	-3400.2377	6836.4753	0.0000	1.0000
HKY	15	-3403.8831	6837.7662	0.0000	1.0000
SYMI	17	-3403.6659	6841.3317	0.0000	1.0000
K2P	12	-3428.2702	6880.5404	0.0000	1.0000
SYM	16	-3427.0384	6886.0769	0.0000	1.0000
F81G	15	-3468.6806	6967.3612	0.0000	1.0000
F81IG	16	-3468.6806	6969.3613	0.0000	1.0000
F81I	15	-3471.6388	6973.2775	0.0000	1.0000
JC69G	12	-3486.5968	6997.1937	0.0000	1.0000
JC69IG	13	-3486.5969	6999.1938	0.0000	1.0000
JC69I	12	-3489.9783	7003.9565	0.0000	1.0000
F81	14	-3491.4923	7010.9847	0.0000	1.0000
JC69	11	-3509.8114	7041.6228	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3377.1955	6864.8349	0.9286	0.9286
HKYI	16	-3380.2996	6871.0432	0.0417	0.9703
HKYIG	17	-3377.1958	6871.7382	0.0294	0.9997
GTRG	20	-3371.6095	6881.2738	0.0003	1.0000
GTRI	20	-3374.4328	6886.9204	0.0000	1.0000
GTRIG	21	-3371.6124	6888.1824	0.0000	1.0000
K2PG	13	-3401.2219	6892.1794	0.0000	1.0000
K2PIG	14	-3401.2312	6899.1008	0.0000	1.0000
K2PI	13	-3404.7302	6899.1961	0.0000	1.0000
HKY	15	-3403.8831	6911.3073	0.0000	1.0000
SYMG	17	-3400.2377	6917.8219	0.0000	1.0000
SYMI	17	-3403.6659	6924.6783	0.0000	1.0000
SYMIG	18	-3400.2377	6924.7247	0.0000	1.0000
GTR	19	-3398.9792	6929.1106	0.0000	1.0000
K2P	12	-3428.2702	6939.3733	0.0000	1.0000
SYM	16	-3427.0384	6964.5207	0.0000	1.0000
F81G	15	-3468.6806	7040.9023	0.0000	1.0000
F81I	15	-3471.6388	7046.8187	0.0000	1.0000
F81IG	16	-3468.6806	7047.8052	0.0000	1.0000
JC69G	12	-3486.5968	7056.0266	0.0000	1.0000
JC69I	12	-3489.9783	7062.7894	0.0000	1.0000
JC69IG	13	-3486.5969	7062.9294	0.0000	1.0000
F81	14	-3491.4923	7079.6231	0.0000	1.0000
JC69	11	-3509.8114	7095.5529	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 591-tomodel-phylip/chr3_496.phylip = 1 - 995;
 Partition Dummy = 1:591-tomodel-phylip/chr3_496.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 591-tomodel-phylip/chr3_496.phylip = 1 - 995;
 Partition Dummy = 1:591-tomodel-phylip/chr3_496.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 591-tomodel-phylip/chr3_496.phylip = 1 - 995;
 Partition Dummy = 1:591-tomodel-phylip/chr3_496.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

