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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-2401.5626	4835.1251	0.4540	0.4540
HKYG	16	-2402.3178	4836.6357	0.2133	0.6673
HKYIG	17	-2401.5786	4837.1572	0.1644	0.8317
GTRG	20	-2399.0914	4838.1827	0.0984	0.9301
GTRI	20	-2399.7468	4839.4936	0.0511	0.9812
GTRIG	21	-2399.7495	4841.4990	0.0187	1.0000
HKY	15	-2411.9145	4853.8290	0.0000	1.0000
GTR	19	-2410.5770	4859.1539	0.0000	1.0000
K2PI	13	-2428.7606	4883.5213	0.0000	1.0000
K2PG	13	-2429.2702	4884.5403	0.0000	1.0000
K2PIG	14	-2429.8165	4887.6329	0.0000	1.0000
SYMI	17	-2427.4889	4888.9778	0.0000	1.0000
SYMG	17	-2428.0795	4890.1589	0.0000	1.0000
SYMIG	18	-2428.1540	4892.3080	0.0000	1.0000
K2P	12	-2440.7142	4905.4283	0.0000	1.0000
SYM	16	-2439.7003	4911.4007	0.0000	1.0000
F81I	15	-2444.0397	4918.0795	0.0000	1.0000
F81G	15	-2444.1724	4918.3448	0.0000	1.0000
F81IG	16	-2444.0646	4920.1293	0.0000	1.0000
F81	14	-2452.6207	4933.2414	0.0000	1.0000
JC69I	12	-2468.7802	4961.5604	0.0000	1.0000
JC69G	12	-2468.9044	4961.8089	0.0000	1.0000
JC69IG	13	-2468.5022	4963.0044	0.0000	1.0000
JC69	11	-2477.7194	4977.4388	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2401.5626	4911.3288	0.6648	0.6648
HKYG	16	-2402.3178	4912.8394	0.3124	0.9771
HKYIG	17	-2401.5786	4918.1236	0.0222	0.9994
HKY	15	-2411.9145	4925.2699	0.0006	1.0000
GTRG	20	-2399.0914	4933.4373	0.0000	1.0000
GTRI	20	-2399.7468	4934.7482	0.0000	1.0000
GTRIG	21	-2399.7495	4941.5163	0.0000	1.0000
K2PI	13	-2428.7606	4945.4367	0.0000	1.0000
K2PG	13	-2429.2702	4946.4558	0.0000	1.0000
GTR	19	-2410.5770	4949.6458	0.0000	1.0000
K2PIG	14	-2429.8165	4954.3111	0.0000	1.0000
K2P	12	-2440.7142	4962.5811	0.0000	1.0000
SYMI	17	-2427.4889	4969.9442	0.0000	1.0000
SYMG	17	-2428.0795	4971.1253	0.0000	1.0000
SYMIG	18	-2428.1540	4978.0372	0.0000	1.0000
SYM	16	-2439.7003	4987.6043	0.0000	1.0000
F81I	15	-2444.0397	4989.5204	0.0000	1.0000
F81G	15	-2444.1724	4989.7857	0.0000	1.0000
F81IG	16	-2444.0646	4996.3330	0.0000	1.0000
F81	14	-2452.6207	4999.9196	0.0000	1.0000
JC69I	12	-2468.7802	5018.7132	0.0000	1.0000
JC69G	12	-2468.9044	5018.9616	0.0000	1.0000
JC69IG	13	-2468.5022	5024.9198	0.0000	1.0000
JC69	11	-2477.7194	5029.8288	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_18148.phylip = 1 - 865;
 Partition Dummy = 1:591-tomodel-phylip/chr3_18148.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_18148.phylip = 1 - 865;
 Partition Dummy = 1:591-tomodel-phylip/chr3_18148.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_18148.phylip = 1 - 865;
 Partition Dummy = 1:591-tomodel-phylip/chr3_18148.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

