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

Input data from file "683-tomodel-phylip/chr7_11588.phylip" (7 taxa, 879 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: K2PI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2370.6516	4781.3031	0.3182	0.3182
GTRG	20	-2370.7791	4781.5582	0.2801	0.5983
HKYI	16	-2375.5876	4783.1753	0.1248	0.7231
HKYG	16	-2375.6549	4783.3097	0.1167	0.8398
GTRIG	21	-2370.8350	4783.6701	0.0974	0.9372
HKYIG	17	-2375.6748	4785.3496	0.0421	0.9793
K2PI	13	-2381.8381	4789.6762	0.0048	0.9841
K2PG	13	-2381.8885	4789.7770	0.0046	0.9887
SYMI	17	-2377.9238	4789.8476	0.0044	0.9932
SYMG	17	-2378.0784	4790.1569	0.0038	0.9970
K2PIG	14	-2381.9121	4791.8241	0.0017	0.9986
SYMIG	18	-2378.0949	4792.1898	0.0014	1.0000
GTR	19	-2392.0839	4822.1678	0.0000	1.0000
HKY	15	-2396.5810	4823.1620	0.0000	1.0000
K2P	12	-2402.3951	4828.7901	0.0000	1.0000
SYM	16	-2398.5358	4829.0717	0.0000	1.0000
F81G	15	-2429.1950	4888.3899	0.0000	1.0000
F81I	15	-2429.2681	4888.5361	0.0000	1.0000
F81IG	16	-2429.2346	4890.4693	0.0000	1.0000
JC69G	12	-2434.0283	4892.0567	0.0000	1.0000
JC69I	12	-2434.1085	4892.2171	0.0000	1.0000
JC69IG	13	-2434.0615	4894.1230	0.0000	1.0000
F81	14	-2447.9270	4923.8539	0.0000	1.0000
JC69	11	-2452.6398	4927.2796	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-2381.8381	4851.8004	0.4947	0.4947
K2PG	13	-2381.8885	4851.9012	0.4704	0.9651
K2PIG	14	-2381.9121	4858.7271	0.0155	0.9806
HKYI	16	-2375.5876	4859.6358	0.0098	0.9904
HKYG	16	-2375.6549	4859.7703	0.0092	0.9996
HKYIG	17	-2375.6748	4866.5889	0.0003	0.9999
SYMI	17	-2377.9238	4871.0869	0.0000	1.0000
SYMG	17	-2378.0784	4871.3962	0.0000	1.0000
GTRI	20	-2370.6516	4876.8788	0.0000	1.0000
GTRG	20	-2370.7791	4877.1339	0.0000	1.0000
SYMIG	18	-2378.0949	4878.2079	0.0000	1.0000
GTRIG	21	-2370.8350	4884.0245	0.0000	1.0000
K2P	12	-2402.3951	4886.1355	0.0000	1.0000
HKY	15	-2396.5810	4894.8438	0.0000	1.0000
SYM	16	-2398.5358	4905.5322	0.0000	1.0000
GTR	19	-2392.0839	4912.9647	0.0000	1.0000
JC69G	12	-2434.0283	4949.4021	0.0000	1.0000
JC69I	12	-2434.1085	4949.5625	0.0000	1.0000
JC69IG	13	-2434.0615	4956.2472	0.0000	1.0000
F81G	15	-2429.1950	4960.0717	0.0000	1.0000
F81I	15	-2429.2681	4960.2179	0.0000	1.0000
F81IG	16	-2429.2346	4966.9298	0.0000	1.0000
JC69	11	-2452.6398	4979.8462	0.0000	1.0000
F81	14	-2447.9270	4990.7569	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_11588.phylip = 1 - 879;
 Partition Dummy = 1:683-tomodel-phylip/chr7_11588.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_11588.phylip = 1 - 879;
 Partition Dummy = 1:683-tomodel-phylip/chr7_11588.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (K2PI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_11588.phylip = 1 - 879;
 Partition Dummy = 1:683-tomodel-phylip/chr7_11588.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

