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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
K2PI	13	-1934.3385	3894.6770	0.2513	0.2513
K2PG	13	-1934.3884	3894.7768	0.2391	0.4904
HKYI	16	-1932.3133	3896.6266	0.0948	0.5852
K2PIG	14	-1934.3386	3896.6771	0.0924	0.6777
HKYG	16	-1932.3442	3896.6883	0.0919	0.7696
SYMI	17	-1932.0405	3898.0810	0.0458	0.8154
SYMG	17	-1932.0441	3898.0882	0.0457	0.8611
GTRI	20	-1929.2546	3898.5092	0.0370	0.8980
GTRG	20	-1929.2622	3898.5244	0.0367	0.9348
HKYIG	17	-1932.3130	3898.6260	0.0349	0.9696
SYMIG	18	-1932.0405	3900.0811	0.0169	0.9865
GTRIG	21	-1929.2622	3900.5245	0.0135	1.0000
GTR	19	-1951.3279	3940.6559	0.0000	1.0000
K2P	12	-1958.8111	3941.6222	0.0000	1.0000
SYM	16	-1955.1841	3942.3683	0.0000	1.0000
HKY	15	-1956.2706	3942.5412	0.0000	1.0000
JC69I	12	-1968.1980	3960.3961	0.0000	1.0000
JC69G	12	-1968.2090	3960.4181	0.0000	1.0000
F81I	15	-1966.0562	3962.1124	0.0000	1.0000
F81G	15	-1966.0789	3962.1577	0.0000	1.0000
JC69IG	13	-1968.2091	3962.4182	0.0000	1.0000
F81IG	16	-1966.0789	3964.1578	0.0000	1.0000
JC69	11	-1989.9202	4001.8404	0.0000	1.0000
F81	14	-1987.3940	4002.7881	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-1934.3385	3956.4261	0.5035	0.5035
K2PG	13	-1934.3884	3956.5259	0.4790	0.9825
K2PIG	14	-1934.3386	3963.1762	0.0172	0.9997
HKYI	16	-1932.3133	3972.6255	0.0002	0.9998
HKYG	16	-1932.3442	3972.6872	0.0001	1.0000
SYMI	17	-1932.0405	3978.8298	0.0000	1.0000
SYMG	17	-1932.0441	3978.8370	0.0000	1.0000
HKYIG	17	-1932.3130	3979.3748	0.0000	1.0000
SYMIG	18	-1932.0405	3985.5798	0.0000	1.0000
GTRI	20	-1929.2546	3993.5078	0.0000	1.0000
GTRG	20	-1929.2622	3993.5230	0.0000	1.0000
K2P	12	-1958.8111	3998.6213	0.0000	1.0000
GTRIG	21	-1929.2622	4000.2731	0.0000	1.0000
HKY	15	-1956.2706	4013.7902	0.0000	1.0000
JC69I	12	-1968.1980	4017.3952	0.0000	1.0000
JC69G	12	-1968.2090	4017.4173	0.0000	1.0000
SYM	16	-1955.1841	4018.3672	0.0000	1.0000
JC69IG	13	-1968.2091	4024.1673	0.0000	1.0000
GTR	19	-1951.3279	4030.9046	0.0000	1.0000
F81I	15	-1966.0562	4033.3614	0.0000	1.0000
F81G	15	-1966.0789	4033.4067	0.0000	1.0000
F81IG	16	-1966.0789	4040.1567	0.0000	1.0000
JC69	11	-1989.9202	4054.0896	0.0000	1.0000
F81	14	-1987.3940	4069.2871	0.0000	1.0000



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

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

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


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

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


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

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

