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

Input data from file "591-tomodel-phylip/chr7_4601.phylip" (7 taxa, 1150 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	-4124.8402	8281.6804	0.6506	0.6506
HKYIG	17	-4124.8080	8283.6159	0.2472	0.8978
GTRG	20	-4123.3557	8286.7114	0.0526	0.9503
HKYI	16	-4127.9917	8287.9834	0.0278	0.9782
GTRIG	21	-4123.3205	8288.6410	0.0200	0.9982
GTRI	20	-4126.7331	8293.4662	0.0018	1.0000
HKY	15	-4163.7080	8357.4160	0.0000	1.0000
GTR	19	-4161.6446	8361.2892	0.0000	1.0000
SYMG	17	-4194.8746	8423.7492	0.0000	1.0000
SYMIG	18	-4194.8833	8425.7665	0.0000	1.0000
SYMI	17	-4198.0443	8430.0885	0.0000	1.0000
K2PG	13	-4203.2255	8432.4510	0.0000	1.0000
K2PIG	14	-4203.2216	8434.4431	0.0000	1.0000
K2PI	13	-4206.0759	8438.1517	0.0000	1.0000
F81G	15	-4220.2829	8470.5658	0.0000	1.0000
F81IG	16	-4220.2853	8472.5705	0.0000	1.0000
F81I	15	-4222.5373	8475.0745	0.0000	1.0000
SYM	16	-4232.2533	8496.5067	0.0000	1.0000
K2P	12	-4240.9509	8505.9018	0.0000	1.0000
F81	14	-4253.7422	8535.4844	0.0000	1.0000
JC69G	12	-4281.7504	8587.5007	0.0000	1.0000
JC69IG	13	-4281.7595	8589.5190	0.0000	1.0000
JC69I	12	-4283.9708	8591.9417	0.0000	1.0000
JC69	11	-4315.1748	8652.3495	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-4124.8402	8362.4407	0.9318	0.9318
HKYI	16	-4127.9917	8368.7437	0.0399	0.9716
HKYIG	17	-4124.8080	8369.4237	0.0284	1.0000
GTRG	20	-4123.3557	8387.6617	0.0000	1.0000
GTRI	20	-4126.7331	8394.4166	0.0000	1.0000
GTRIG	21	-4123.3205	8394.6388	0.0000	1.0000
HKY	15	-4163.7080	8433.1287	0.0000	1.0000
GTR	19	-4161.6446	8457.1920	0.0000	1.0000
K2PG	13	-4203.2255	8498.0688	0.0000	1.0000
K2PI	13	-4206.0759	8503.7694	0.0000	1.0000
K2PIG	14	-4203.2216	8505.1084	0.0000	1.0000
SYMG	17	-4194.8746	8509.5570	0.0000	1.0000
SYMI	17	-4198.0443	8515.8963	0.0000	1.0000
SYMIG	18	-4194.8833	8516.6218	0.0000	1.0000
F81G	15	-4220.2829	8546.2785	0.0000	1.0000
F81I	15	-4222.5373	8550.7873	0.0000	1.0000
F81IG	16	-4220.2853	8553.3308	0.0000	1.0000
K2P	12	-4240.9509	8566.4720	0.0000	1.0000
SYM	16	-4232.2533	8577.2669	0.0000	1.0000
F81	14	-4253.7422	8606.1496	0.0000	1.0000
JC69G	12	-4281.7504	8648.0709	0.0000	1.0000
JC69I	12	-4283.9708	8652.5119	0.0000	1.0000
JC69IG	13	-4281.7595	8655.1367	0.0000	1.0000
JC69	11	-4315.1748	8707.8722	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_4601.phylip = 1 - 1150;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4601.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_4601.phylip = 1 - 1150;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4601.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_4601.phylip = 1 - 1150;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4601.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

