
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:19 04/04/2011

Input data from file "683-tomodel-phylip/chr9_3537.phylip" (7 taxa, 1119 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	-3441.0969	6922.1938	0.6583	0.6583
GTRIG	21	-3441.0772	6924.1545	0.2470	0.9053
GTRI	20	-3443.0395	6926.0790	0.0944	0.9996
HKYG	16	-3453.0713	6938.1425	0.0002	0.9999
HKYIG	17	-3453.0558	6940.1116	0.0001	1.0000
HKYI	16	-3454.5850	6941.1700	0.0000	1.0000
GTR	19	-3471.0320	6980.0639	0.0000	1.0000
HKY	15	-3483.1566	6996.3133	0.0000	1.0000
SYMG	17	-3487.8851	7009.7702	0.0000	1.0000
SYMIG	18	-3487.8699	7011.7398	0.0000	1.0000
SYMI	17	-3488.9062	7011.8124	0.0000	1.0000
K2PG	13	-3495.9203	7017.8407	0.0000	1.0000
K2PIG	14	-3495.8976	7019.7952	0.0000	1.0000
K2PI	13	-3496.9635	7019.9270	0.0000	1.0000
SYM	16	-3517.6819	7067.3639	0.0000	1.0000
K2P	12	-3525.8162	7075.6323	0.0000	1.0000
F81G	15	-3538.4963	7106.9927	0.0000	1.0000
F81IG	16	-3538.4780	7108.9560	0.0000	1.0000
F81I	15	-3539.6692	7109.3385	0.0000	1.0000
F81	14	-3565.2914	7158.5829	0.0000	1.0000
JC69G	12	-3573.5801	7171.1603	0.0000	1.0000
JC69IG	13	-3573.5559	7173.1119	0.0000	1.0000
JC69I	12	-3574.5895	7173.1790	0.0000	1.0000
JC69	11	-3600.6146	7223.2292	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3453.0713	7018.4656	0.7147	0.7147
HKYI	16	-3454.5850	7021.4931	0.1573	0.8720
GTRG	20	-3441.0969	7022.5976	0.0905	0.9626
HKYIG	17	-3453.0558	7025.4549	0.0217	0.9843
GTRI	20	-3443.0395	7026.4828	0.0130	0.9972
GTRIG	21	-3441.0772	7029.5785	0.0028	1.0000
HKY	15	-3483.1566	7071.6162	0.0000	1.0000
GTR	19	-3471.0320	7075.4475	0.0000	1.0000
K2PG	13	-3495.9203	7083.1032	0.0000	1.0000
K2PI	13	-3496.9635	7085.1895	0.0000	1.0000
K2PIG	14	-3495.8976	7090.0778	0.0000	1.0000
SYMG	17	-3487.8851	7095.1135	0.0000	1.0000
SYMI	17	-3488.9062	7097.1557	0.0000	1.0000
SYMIG	18	-3487.8699	7102.1032	0.0000	1.0000
K2P	12	-3525.8162	7135.8746	0.0000	1.0000
SYM	16	-3517.6819	7147.6869	0.0000	1.0000
F81G	15	-3538.4963	7182.2955	0.0000	1.0000
F81I	15	-3539.6692	7184.6413	0.0000	1.0000
F81IG	16	-3538.4780	7189.2791	0.0000	1.0000
F81	14	-3565.2914	7228.8655	0.0000	1.0000
JC69G	12	-3573.5801	7231.4025	0.0000	1.0000
JC69I	12	-3574.5895	7233.4213	0.0000	1.0000
JC69IG	13	-3573.5559	7238.3743	0.0000	1.0000
JC69	11	-3600.6146	7278.4513	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 683-tomodel-phylip/chr9_3537.phylip = 1 - 1119;
 Partition Dummy = 1:683-tomodel-phylip/chr9_3537.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 683-tomodel-phylip/chr9_3537.phylip = 1 - 1119;
 Partition Dummy = 1:683-tomodel-phylip/chr9_3537.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 683-tomodel-phylip/chr9_3537.phylip = 1 - 1119;
 Partition Dummy = 1:683-tomodel-phylip/chr9_3537.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

