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

Input data from file "591-tomodel-phylip/chr7_6598.phylip" (7 taxa, 993 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	-3532.4520	7096.9040	0.3826	0.3826
GTRG	20	-3528.9904	7097.9807	0.2233	0.6060
HKYIG	17	-3532.4581	7098.9161	0.1399	0.7459
HKYI	16	-3533.7027	7099.4054	0.1095	0.8554
GTRIG	21	-3528.9967	7099.9934	0.0816	0.9371
GTRI	20	-3530.2568	7100.5135	0.0629	1.0000
HKY	15	-3556.1799	7142.3597	0.0000	1.0000
GTR	19	-3553.3422	7144.6844	0.0000	1.0000
SYMG	17	-3617.0363	7268.0725	0.0000	1.0000
K2PG	13	-3621.9064	7269.8129	0.0000	1.0000
SYMIG	18	-3617.0596	7270.1192	0.0000	1.0000
SYMI	17	-3618.7923	7271.5847	0.0000	1.0000
K2PIG	14	-3621.9271	7271.8543	0.0000	1.0000
K2PI	13	-3623.4271	7272.8542	0.0000	1.0000
F81G	15	-3629.3200	7288.6399	0.0000	1.0000
F81IG	16	-3629.3431	7290.6862	0.0000	1.0000
F81I	15	-3630.4749	7290.9499	0.0000	1.0000
SYM	16	-3639.7886	7311.5773	0.0000	1.0000
K2P	12	-3644.0347	7312.0694	0.0000	1.0000
F81	14	-3649.4492	7326.8984	0.0000	1.0000
JC69G	12	-3697.0787	7418.1575	0.0000	1.0000
JC69IG	13	-3697.0994	7420.1987	0.0000	1.0000
JC69I	12	-3698.2638	7420.5275	0.0000	1.0000
JC69	11	-3716.9193	7455.8386	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3532.4520	7175.3157	0.7588	0.7588
HKYI	16	-3533.7027	7177.8171	0.2172	0.9760
HKYIG	17	-3532.4581	7182.2285	0.0239	1.0000
GTRG	20	-3528.9904	7195.9953	0.0000	1.0000
GTRI	20	-3530.2568	7198.5281	0.0000	1.0000
GTRIG	21	-3528.9967	7202.9087	0.0000	1.0000
HKY	15	-3556.1799	7215.8707	0.0000	1.0000
GTR	19	-3553.3422	7237.7983	0.0000	1.0000
K2PG	13	-3621.9064	7333.5224	0.0000	1.0000
K2PI	13	-3623.4271	7336.5637	0.0000	1.0000
K2PIG	14	-3621.9271	7340.4645	0.0000	1.0000
SYMG	17	-3617.0363	7351.3850	0.0000	1.0000
SYMI	17	-3618.7923	7354.8971	0.0000	1.0000
SYMIG	18	-3617.0596	7358.3324	0.0000	1.0000
F81G	15	-3629.3200	7362.1509	0.0000	1.0000
F81I	15	-3630.4749	7364.4608	0.0000	1.0000
F81IG	16	-3629.3431	7369.0979	0.0000	1.0000
K2P	12	-3644.0347	7370.8782	0.0000	1.0000
SYM	16	-3639.7886	7389.9890	0.0000	1.0000
F81	14	-3649.4492	7395.5086	0.0000	1.0000
JC69G	12	-3697.0787	7476.9662	0.0000	1.0000
JC69I	12	-3698.2638	7479.3363	0.0000	1.0000
JC69IG	13	-3697.0994	7483.9082	0.0000	1.0000
JC69	11	-3716.9193	7509.7466	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 591-tomodel-phylip/chr7_6598.phylip = 1 - 993;
 Partition Dummy = 1:591-tomodel-phylip/chr7_6598.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 591-tomodel-phylip/chr7_6598.phylip = 1 - 993;
 Partition Dummy = 1:591-tomodel-phylip/chr7_6598.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 591-tomodel-phylip/chr7_6598.phylip = 1 - 993;
 Partition Dummy = 1:591-tomodel-phylip/chr7_6598.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

