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

Input data from file "591-tomodel-phylip/chr2_24792.phylip" (7 taxa, 1043 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2611.7495	5263.4991	0.4296	0.4296
GTRI	20	-2611.8066	5263.6132	0.4058	0.8355
GTRIG	21	-2611.7497	5265.4993	0.1580	0.9935
HKYI	16	-2620.7094	5273.4187	0.0030	0.9965
HKYG	16	-2620.8771	5273.7541	0.0025	0.9991
HKYIG	17	-2620.8771	5275.7541	0.0009	1.0000
GTR	19	-2631.5040	5301.0080	0.0000	1.0000
HKY	15	-2639.6983	5309.3966	0.0000	1.0000
F81I	15	-2650.8853	5331.7707	0.0000	1.0000
F81G	15	-2651.1315	5332.2631	0.0000	1.0000
F81IG	16	-2651.1315	5334.2631	0.0000	1.0000
SYMG	17	-2657.1593	5348.3186	0.0000	1.0000
SYMI	17	-2657.2235	5348.4469	0.0000	1.0000
SYMIG	18	-2657.1758	5350.3516	0.0000	1.0000
K2PI	13	-2662.2180	5350.4360	0.0000	1.0000
K2PG	13	-2662.3529	5350.7058	0.0000	1.0000
K2PIG	14	-2662.4208	5352.8417	0.0000	1.0000
F81	14	-2668.6950	5365.3899	0.0000	1.0000
SYM	16	-2676.4983	5384.9967	0.0000	1.0000
K2P	12	-2681.5476	5387.0952	0.0000	1.0000
JC69I	12	-2689.2356	5402.4712	0.0000	1.0000
JC69G	12	-2689.4715	5402.9430	0.0000	1.0000
JC69IG	13	-2689.4715	5404.9431	0.0000	1.0000
JC69	11	-2707.1812	5436.3625	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2620.7094	5352.6164	0.5302	0.5302
HKYG	16	-2620.8771	5352.9518	0.4484	0.9786
HKYIG	17	-2620.8771	5359.9017	0.0139	0.9925
GTRG	20	-2611.7495	5362.4962	0.0038	0.9963
GTRI	20	-2611.8066	5362.6103	0.0036	0.9999
GTRIG	21	-2611.7497	5369.4463	0.0001	1.0000
HKY	15	-2639.6983	5383.6445	0.0000	1.0000
GTR	19	-2631.5040	5395.0553	0.0000	1.0000
F81I	15	-2650.8853	5406.0185	0.0000	1.0000
F81G	15	-2651.1315	5406.5109	0.0000	1.0000
F81IG	16	-2651.1315	5413.4608	0.0000	1.0000
K2PI	13	-2662.2180	5414.7841	0.0000	1.0000
K2PG	13	-2662.3529	5415.0539	0.0000	1.0000
K2PIG	14	-2662.4208	5422.1397	0.0000	1.0000
SYMG	17	-2657.1593	5432.4661	0.0000	1.0000
SYMI	17	-2657.2235	5432.5945	0.0000	1.0000
F81	14	-2668.6950	5434.6879	0.0000	1.0000
SYMIG	18	-2657.1758	5439.4491	0.0000	1.0000
K2P	12	-2681.5476	5446.4935	0.0000	1.0000
JC69I	12	-2689.2356	5461.8695	0.0000	1.0000
JC69G	12	-2689.4715	5462.3413	0.0000	1.0000
SYM	16	-2676.4983	5464.1944	0.0000	1.0000
JC69IG	13	-2689.4715	5469.2912	0.0000	1.0000
JC69	11	-2707.1812	5490.8109	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_24792.phylip = 1 - 1043;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24792.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_24792.phylip = 1 - 1043;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24792.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 (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_24792.phylip = 1 - 1043;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24792.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

