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

Input data from file "591-tomodel-phylip/chr9_5218.phylip" (7 taxa, 899 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	-3115.3421	6270.6842	0.4300	0.4300
GTRI	20	-3115.4411	6270.8821	0.3895	0.8194
GTRIG	21	-3115.2098	6272.4196	0.1806	1.0000
GTR	19	-3127.0904	6292.1808	0.0000	1.0000
HKYG	16	-3134.3786	6300.7571	0.0000	1.0000
HKYI	16	-3134.9948	6301.9896	0.0000	1.0000
HKYIG	17	-3134.3904	6302.7807	0.0000	1.0000
SYMG	17	-3135.6656	6305.3313	0.0000	1.0000
SYMI	17	-3136.1992	6306.3983	0.0000	1.0000
SYMIG	18	-3135.6532	6307.3063	0.0000	1.0000
HKY	15	-3145.5303	6321.0607	0.0000	1.0000
K2PG	13	-3148.6438	6323.2877	0.0000	1.0000
K2PI	13	-3149.4776	6324.9552	0.0000	1.0000
SYM	16	-3146.5615	6325.1230	0.0000	1.0000
K2PIG	14	-3148.6372	6325.2744	0.0000	1.0000
K2P	12	-3159.6990	6343.3980	0.0000	1.0000
F81G	15	-3182.3375	6394.6749	0.0000	1.0000
F81I	15	-3182.4058	6394.8116	0.0000	1.0000
F81IG	16	-3182.2657	6396.5313	0.0000	1.0000
F81	14	-3191.1366	6410.2731	0.0000	1.0000
JC69G	12	-3193.9646	6411.9292	0.0000	1.0000
JC69I	12	-3194.2273	6412.4547	0.0000	1.0000
JC69IG	13	-3193.9482	6413.8963	0.0000	1.0000
JC69	11	-3203.0354	6428.0708	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-3115.3421	6366.7098	0.5125	0.5125
GTRI	20	-3115.4411	6366.9078	0.4642	0.9768
GTRIG	21	-3115.2098	6373.2465	0.0195	0.9963
HKYG	16	-3134.3786	6377.5777	0.0022	0.9985
HKYI	16	-3134.9948	6378.8102	0.0012	0.9997
GTR	19	-3127.0904	6383.4052	0.0001	0.9998
HKYIG	17	-3134.3904	6384.4025	0.0001	0.9999
K2PG	13	-3148.6438	6385.7044	0.0000	0.9999
SYMG	17	-3135.6656	6386.9531	0.0000	1.0000
K2PI	13	-3149.4776	6387.3718	0.0000	1.0000
SYMI	17	-3136.1992	6388.0202	0.0000	1.0000
K2PIG	14	-3148.6372	6392.4923	0.0000	1.0000
HKY	15	-3145.5303	6393.0799	0.0000	1.0000
SYMIG	18	-3135.6532	6393.7294	0.0000	1.0000
K2P	12	-3159.6990	6401.0134	0.0000	1.0000
SYM	16	-3146.5615	6401.9436	0.0000	1.0000
F81G	15	-3182.3375	6466.6942	0.0000	1.0000
F81I	15	-3182.4058	6466.8309	0.0000	1.0000
JC69G	12	-3193.9646	6469.5446	0.0000	1.0000
JC69I	12	-3194.2273	6470.0701	0.0000	1.0000
F81IG	16	-3182.2657	6473.3519	0.0000	1.0000
JC69IG	13	-3193.9482	6476.3130	0.0000	1.0000
F81	14	-3191.1366	6477.4911	0.0000	1.0000
JC69	11	-3203.0354	6480.8849	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/chr9_5218.phylip = 1 - 899;
 Partition Dummy = 1:591-tomodel-phylip/chr9_5218.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/chr9_5218.phylip = 1 - 899;
 Partition Dummy = 1:591-tomodel-phylip/chr9_5218.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/chr9_5218.phylip = 1 - 899;
 Partition Dummy = 1:591-tomodel-phylip/chr9_5218.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

