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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3450.4160	6940.8319	0.5997	0.5997
GTRIG	21	-3450.1980	6942.3961	0.2743	0.8740
GTRI	20	-3451.9758	6943.9515	0.1260	1.0000
HKYG	16	-3472.3112	6976.6223	0.0000	1.0000
HKYIG	17	-3472.2817	6978.5635	0.0000	1.0000
HKYI	16	-3473.4378	6978.8756	0.0000	1.0000
SYMG	17	-3475.3941	6984.7882	0.0000	1.0000
SYMIG	18	-3475.6156	6987.2312	0.0000	1.0000
SYMI	17	-3476.7810	6987.5619	0.0000	1.0000
K2PG	13	-3493.0608	7012.1217	0.0000	1.0000
K2PI	13	-3494.0414	7014.0827	0.0000	1.0000
K2PIG	14	-3493.0583	7014.1166	0.0000	1.0000
GTR	19	-3501.0168	7040.0336	0.0000	1.0000
F81G	15	-3521.6411	7073.2822	0.0000	1.0000
F81I	15	-3522.2200	7074.4401	0.0000	1.0000
F81IG	16	-3521.6923	7075.3847	0.0000	1.0000
HKY	15	-3523.1086	7076.2173	0.0000	1.0000
SYM	16	-3525.4885	7082.9771	0.0000	1.0000
JC69IG	13	-3539.2956	7104.5913	0.0000	1.0000
JC69G	12	-3541.8602	7107.7203	0.0000	1.0000
JC69I	12	-3542.9183	7109.8367	0.0000	1.0000
K2P	12	-3545.8562	7115.7123	0.0000	1.0000
F81	14	-3572.5017	7173.0034	0.0000	1.0000
JC69	11	-3590.5322	7203.0644	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-3450.4160	7041.7474	0.8016	0.8016
GTRI	20	-3451.9758	7044.8671	0.1685	0.9701
GTRIG	21	-3450.1980	7048.3574	0.0294	0.9996
HKYG	16	-3472.3112	7057.3547	0.0003	0.9999
HKYI	16	-3473.4378	7059.6081	0.0001	1.0000
HKYIG	17	-3472.2817	7064.3417	0.0000	1.0000
SYMG	17	-3475.3941	7070.5664	0.0000	1.0000
SYMI	17	-3476.7810	7073.3401	0.0000	1.0000
K2PG	13	-3493.0608	7077.7168	0.0000	1.0000
SYMIG	18	-3475.6156	7078.0551	0.0000	1.0000
K2PI	13	-3494.0414	7079.6778	0.0000	1.0000
K2PIG	14	-3493.0583	7084.7575	0.0000	1.0000
GTR	19	-3501.0168	7135.9033	0.0000	1.0000
F81G	15	-3521.6411	7148.9689	0.0000	1.0000
F81I	15	-3522.2200	7150.1267	0.0000	1.0000
HKY	15	-3523.1086	7151.9039	0.0000	1.0000
F81IG	16	-3521.6923	7156.1171	0.0000	1.0000
SYM	16	-3525.4885	7163.7095	0.0000	1.0000
JC69G	12	-3541.8602	7168.2697	0.0000	1.0000
JC69IG	13	-3539.2956	7170.1864	0.0000	1.0000
JC69I	12	-3542.9183	7170.3860	0.0000	1.0000
K2P	12	-3545.8562	7176.2616	0.0000	1.0000
F81	14	-3572.5017	7243.6443	0.0000	1.0000
JC69	11	-3590.5322	7258.5679	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/chr1_12757.phylip = 1 - 1148;
 Partition Dummy = 1:591-tomodel-phylip/chr1_12757.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/chr1_12757.phylip = 1 - 1148;
 Partition Dummy = 1:591-tomodel-phylip/chr1_12757.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 (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_12757.phylip = 1 - 1148;
 Partition Dummy = 1:591-tomodel-phylip/chr1_12757.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;


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

