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

Input data from file "591-tomodel-phylip/chr2_21265.phylip" (7 taxa, 1035 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	-4089.9368	8219.8737	0.7263	0.7263
GTRIG	21	-4089.9390	8221.8779	0.2666	0.9929
GTRI	20	-4094.5803	8229.1605	0.0070	0.9999
HKYG	16	-4103.6406	8239.2812	0.0000	1.0000
HKYIG	17	-4103.5715	8241.1431	0.0000	1.0000
HKYI	16	-4106.2216	8244.4433	0.0000	1.0000
SYMG	17	-4115.8793	8265.7586	0.0000	1.0000
SYMIG	18	-4115.9139	8267.8278	0.0000	1.0000
SYMI	17	-4118.3000	8270.6001	0.0000	1.0000
K2PG	13	-4123.5894	8273.1787	0.0000	1.0000
K2PIG	14	-4123.4220	8274.8440	0.0000	1.0000
K2PI	13	-4125.2491	8276.4981	0.0000	1.0000
GTR	19	-4135.5650	8309.1300	0.0000	1.0000
HKY	15	-4147.2613	8324.5226	0.0000	1.0000
SYM	16	-4160.6990	8353.3981	0.0000	1.0000
K2P	12	-4167.6311	8359.2622	0.0000	1.0000
F81G	15	-4190.7583	8411.5166	0.0000	1.0000
F81IG	16	-4190.6763	8413.3525	0.0000	1.0000
F81I	15	-4192.4622	8414.9243	0.0000	1.0000
JC69G	12	-4205.5546	8435.1092	0.0000	1.0000
JC69IG	13	-4205.4135	8436.8271	0.0000	1.0000
JC69I	12	-4206.7012	8437.4023	0.0000	1.0000
F81	14	-4229.8208	8487.6417	0.0000	1.0000
JC69	11	-4245.0505	8512.1011	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-4103.6406	8318.3557	0.5056	0.5056
GTRG	20	-4089.9368	8318.7168	0.4221	0.9277
HKYI	16	-4106.2216	8323.5178	0.0383	0.9660
HKYIG	17	-4103.5715	8325.1597	0.0168	0.9828
GTRIG	21	-4089.9390	8325.6632	0.0131	0.9959
GTRI	20	-4094.5803	8328.0037	0.0041	1.0000
K2PG	13	-4123.5894	8337.4268	0.0000	1.0000
K2PI	13	-4125.2491	8340.7462	0.0000	1.0000
K2PIG	14	-4123.4220	8344.0342	0.0000	1.0000
SYMG	17	-4115.8793	8349.7753	0.0000	1.0000
SYMI	17	-4118.3000	8354.6167	0.0000	1.0000
SYMIG	18	-4115.9139	8356.7866	0.0000	1.0000
HKY	15	-4147.2613	8398.6550	0.0000	1.0000
GTR	19	-4135.5650	8403.0310	0.0000	1.0000
K2P	12	-4167.6311	8418.5681	0.0000	1.0000
SYM	16	-4160.6990	8432.4726	0.0000	1.0000
F81G	15	-4190.7583	8485.6489	0.0000	1.0000
F81I	15	-4192.4622	8489.0567	0.0000	1.0000
F81IG	16	-4190.6763	8492.4270	0.0000	1.0000
JC69G	12	-4205.5546	8494.4151	0.0000	1.0000
JC69I	12	-4206.7012	8496.7082	0.0000	1.0000
JC69IG	13	-4205.4135	8501.0751	0.0000	1.0000
F81	14	-4229.8208	8556.8319	0.0000	1.0000
JC69	11	-4245.0505	8566.4648	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/chr2_21265.phylip = 1 - 1035;
 Partition Dummy = 1:591-tomodel-phylip/chr2_21265.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/chr2_21265.phylip = 1 - 1035;
 Partition Dummy = 1:591-tomodel-phylip/chr2_21265.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/chr2_21265.phylip = 1 - 1035;
 Partition Dummy = 1:591-tomodel-phylip/chr2_21265.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

