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

Input data from file "683-tomodel-phylip/chr5_10164.phylip" (7 taxa, 1161 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	-4165.8290	8371.6580	0.7537	0.7537
GTRIG	21	-4166.1035	8374.2071	0.2107	0.9644
GTRI	20	-4168.8911	8377.7822	0.0353	0.9996
HKYG	16	-4177.8661	8387.7322	0.0002	0.9999
HKYIG	17	-4177.8664	8389.7328	0.0001	1.0000
GTR	19	-4177.4608	8392.9216	0.0000	1.0000
HKYI	16	-4180.6663	8393.3327	0.0000	1.0000
HKY	15	-4188.5255	8407.0510	0.0000	1.0000
SYMG	17	-4190.6651	8415.3302	0.0000	1.0000
SYMIG	18	-4190.8106	8417.6212	0.0000	1.0000
K2PG	13	-4197.2123	8420.4246	0.0000	1.0000
SYMI	17	-4193.6047	8421.2095	0.0000	1.0000
K2PIG	14	-4197.1847	8422.3695	0.0000	1.0000
K2PI	13	-4199.8251	8425.6501	0.0000	1.0000
SYM	16	-4200.8669	8433.7338	0.0000	1.0000
K2P	12	-4207.4639	8438.9277	0.0000	1.0000
F81G	15	-4285.1023	8600.2047	0.0000	1.0000
F81I	15	-4286.1307	8602.2613	0.0000	1.0000
F81IG	16	-4287.2173	8606.4347	0.0000	1.0000
F81	14	-4291.4558	8610.9116	0.0000	1.0000
JC69G	12	-4298.4990	8620.9981	0.0000	1.0000
JC69IG	13	-4301.0075	8628.0151	0.0000	1.0000
JC69I	12	-4303.4067	8630.8134	0.0000	1.0000
JC69	11	-4304.7250	8631.4500	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-4177.8661	8468.6448	0.8162	0.8162
GTRG	20	-4165.8290	8472.7988	0.1023	0.9185
HKYI	16	-4180.6663	8474.2453	0.0496	0.9682
HKYIG	17	-4177.8664	8475.7024	0.0239	0.9921
GTRI	20	-4168.8911	8478.9230	0.0048	0.9969
GTRIG	21	-4166.1035	8480.4049	0.0023	0.9992
HKY	15	-4188.5255	8482.9066	0.0007	0.9998
K2PG	13	-4197.2123	8486.1661	0.0001	1.0000
GTR	19	-4177.4608	8489.0053	0.0000	1.0000
K2PI	13	-4199.8251	8491.3916	0.0000	1.0000
K2PIG	14	-4197.1847	8493.1680	0.0000	1.0000
K2P	12	-4207.4639	8499.6122	0.0000	1.0000
SYMG	17	-4190.6651	8501.2998	0.0000	1.0000
SYMI	17	-4193.6047	8507.1791	0.0000	1.0000
SYMIG	18	-4190.8106	8508.6478	0.0000	1.0000
SYM	16	-4200.8669	8514.6464	0.0000	1.0000
F81G	15	-4285.1023	8676.0602	0.0000	1.0000
F81I	15	-4286.1307	8678.1169	0.0000	1.0000
JC69G	12	-4298.4990	8681.6825	0.0000	1.0000
F81	14	-4291.4558	8681.7101	0.0000	1.0000
JC69	11	-4304.7250	8687.0774	0.0000	1.0000
F81IG	16	-4287.2173	8687.3473	0.0000	1.0000
JC69I	12	-4303.4067	8691.4978	0.0000	1.0000
JC69IG	13	-4301.0075	8693.7566	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/chr5_10164.phylip = 1 - 1161;
 Partition Dummy = 1:683-tomodel-phylip/chr5_10164.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/chr5_10164.phylip = 1 - 1161;
 Partition Dummy = 1:683-tomodel-phylip/chr5_10164.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/chr5_10164.phylip = 1 - 1161;
 Partition Dummy = 1:683-tomodel-phylip/chr5_10164.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

