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

Input data from file "591-tomodel-phylip/chr6_7323.phylip" (7 taxa, 905 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	-3541.9508	7123.9016	0.4477	0.4477
GTRI	20	-3542.5090	7125.0180	0.2562	0.7039
GTRIG	21	-3541.8358	7125.6716	0.1848	0.8887
HKYG	16	-3548.0325	7128.0650	0.0558	0.9446
HKYI	16	-3548.5739	7129.1479	0.0325	0.9771
HKYIG	17	-3547.9283	7129.8566	0.0228	0.9999
SYMG	17	-3553.8474	7141.6948	0.0001	0.9999
SYMI	17	-3554.3534	7142.7069	0.0000	1.0000
SYMIG	18	-3553.7215	7143.4430	0.0000	1.0000
GTR	19	-3555.0501	7148.1001	0.0000	1.0000
K2PG	13	-3561.7791	7149.5582	0.0000	1.0000
K2PI	13	-3562.2047	7150.4094	0.0000	1.0000
K2PIG	14	-3561.6515	7151.3031	0.0000	1.0000
HKY	15	-3561.9680	7153.9359	0.0000	1.0000
SYM	16	-3566.9788	7165.9575	0.0000	1.0000
K2P	12	-3575.3542	7174.7085	0.0000	1.0000
F81I	15	-3633.8636	7297.7273	0.0000	1.0000
F81G	15	-3634.1065	7298.2131	0.0000	1.0000
F81IG	16	-3633.7164	7299.4328	0.0000	1.0000
JC69I	12	-3645.8537	7315.7074	0.0000	1.0000
JC69G	12	-3646.1167	7316.2334	0.0000	1.0000
JC69IG	13	-3645.7139	7317.4277	0.0000	1.0000
F81	14	-3645.3897	7318.7794	0.0000	1.0000
JC69	11	-3657.0992	7336.1984	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3548.0325	7204.9919	0.5991	0.5991
HKYI	16	-3548.5739	7206.0748	0.3486	0.9477
HKYIG	17	-3547.9283	7211.5915	0.0221	0.9698
K2PG	13	-3561.7791	7212.0613	0.0175	0.9873
K2PI	13	-3562.2047	7212.9126	0.0114	0.9987
K2PIG	14	-3561.6515	7218.6142	0.0007	0.9994
GTRG	20	-3541.9508	7220.0603	0.0003	0.9997
GTRI	20	-3542.5090	7221.1767	0.0002	0.9999
SYMG	17	-3553.8474	7223.4297	0.0001	0.9999
SYMI	17	-3554.3534	7224.4418	0.0000	1.0000
HKY	15	-3561.9680	7226.0550	0.0000	1.0000
GTRIG	21	-3541.8358	7226.6382	0.0000	1.0000
SYMIG	18	-3553.7215	7229.9858	0.0000	1.0000
K2P	12	-3575.3542	7232.4037	0.0000	1.0000
GTR	19	-3555.0501	7239.4509	0.0000	1.0000
SYM	16	-3566.9788	7242.8845	0.0000	1.0000
F81I	15	-3633.8636	7369.8463	0.0000	1.0000
F81G	15	-3634.1065	7370.3321	0.0000	1.0000
JC69I	12	-3645.8537	7373.4026	0.0000	1.0000
JC69G	12	-3646.1167	7373.9286	0.0000	1.0000
F81IG	16	-3633.7164	7376.3597	0.0000	1.0000
JC69IG	13	-3645.7139	7379.9309	0.0000	1.0000
F81	14	-3645.3897	7386.0905	0.0000	1.0000
JC69	11	-3657.0992	7389.0857	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/chr6_7323.phylip = 1 - 905;
 Partition Dummy = 1:591-tomodel-phylip/chr6_7323.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/chr6_7323.phylip = 1 - 905;
 Partition Dummy = 1:591-tomodel-phylip/chr6_7323.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/chr6_7323.phylip = 1 - 905;
 Partition Dummy = 1:591-tomodel-phylip/chr6_7323.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

