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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2387.8031	4815.6062	0.4082	0.4082
GTRG	20	-2388.1397	4816.2795	0.2915	0.6997
GTRIG	21	-2387.8030	4817.6060	0.1502	0.8499
HKYI	16	-2393.5329	4819.0658	0.0724	0.9222
HKYG	16	-2393.8807	4819.7614	0.0511	0.9734
HKYIG	17	-2393.5327	4821.0653	0.0266	1.0000
GTR	19	-2416.2465	4870.4930	0.0000	1.0000
HKY	15	-2421.0156	4872.0313	0.0000	1.0000
F81I	15	-2429.9611	4889.9221	0.0000	1.0000
F81G	15	-2430.1712	4890.3424	0.0000	1.0000
F81IG	16	-2429.9606	4891.9212	0.0000	1.0000
F81	14	-2455.8146	4939.6291	0.0000	1.0000
K2PI	13	-2480.9289	4987.8578	0.0000	1.0000
K2PG	13	-2481.1376	4988.2752	0.0000	1.0000
K2PIG	14	-2480.8202	4989.6404	0.0000	1.0000
SYMI	17	-2479.1446	4992.2892	0.0000	1.0000
SYMG	17	-2479.3229	4992.6457	0.0000	1.0000
SYMIG	18	-2479.0744	4994.1488	0.0000	1.0000
K2P	12	-2505.9621	5035.9242	0.0000	1.0000
SYM	16	-2504.1257	5040.2514	0.0000	1.0000
JC69I	12	-2509.1270	5042.2541	0.0000	1.0000
JC69G	12	-2509.2545	5042.5089	0.0000	1.0000
JC69IG	13	-2509.1270	5044.2541	0.0000	1.0000
JC69	11	-2533.9234	5089.8468	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2393.5329	4899.7703	0.5760	0.5760
HKYG	16	-2393.8807	4900.4659	0.4068	0.9827
HKYIG	17	-2393.5327	4906.8139	0.0170	0.9998
GTRI	20	-2387.8031	4916.4869	0.0001	0.9999
GTRG	20	-2388.1397	4917.1601	0.0001	1.0000
GTRIG	21	-2387.8030	4923.5307	0.0000	1.0000
HKY	15	-2421.0156	4947.6918	0.0000	1.0000
F81I	15	-2429.9611	4965.5826	0.0000	1.0000
F81G	15	-2430.1712	4966.0029	0.0000	1.0000
GTR	19	-2416.2465	4966.3296	0.0000	1.0000
F81IG	16	-2429.9606	4972.6257	0.0000	1.0000
F81	14	-2455.8146	5010.2456	0.0000	1.0000
K2PI	13	-2480.9289	5053.4302	0.0000	1.0000
K2PG	13	-2481.1376	5053.8476	0.0000	1.0000
K2PIG	14	-2480.8202	5060.2569	0.0000	1.0000
SYMI	17	-2479.1446	5078.0377	0.0000	1.0000
SYMG	17	-2479.3229	5078.3943	0.0000	1.0000
SYMIG	18	-2479.0744	5084.9414	0.0000	1.0000
K2P	12	-2505.9621	5096.4526	0.0000	1.0000
JC69I	12	-2509.1270	5102.7825	0.0000	1.0000
JC69G	12	-2509.2545	5103.0373	0.0000	1.0000
JC69IG	13	-2509.1270	5109.8265	0.0000	1.0000
SYM	16	-2504.1257	5120.9559	0.0000	1.0000
JC69	11	-2533.9234	5145.3311	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_4911.phylip = 1 - 1146;
 Partition Dummy = 1:683-tomodel-phylip/chr7_4911.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_4911.phylip = 1 - 1146;
 Partition Dummy = 1:683-tomodel-phylip/chr7_4911.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 (HKYI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_4911.phylip = 1 - 1146;
 Partition Dummy = 1:683-tomodel-phylip/chr7_4911.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

