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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
HKYG	16	-2101.0352	4234.9353	0.3940	0.3940
HKYI	16	-2101.4363	4235.7375	0.2638	0.6577
HKYIG	17	-2100.4872	4235.9489	0.2373	0.8951
GTRG	20	-2099.0445	4239.4329	0.0416	0.9367
GTRI	20	-2099.2944	4239.9329	0.0324	0.9690
GTRIG	21	-2098.2710	4240.0228	0.0310	1.0000
HKY	15	-2113.0443	4256.8505	0.0000	1.0000
GTR	19	-2111.2303	4261.6747	0.0000	1.0000
F81G	15	-2149.1803	4329.1225	0.0000	1.0000
F81I	15	-2149.5260	4329.8138	0.0000	1.0000
F81IG	16	-2149.1377	4331.1403	0.0000	1.0000
F81	14	-2160.1501	4348.9659	0.0000	1.0000
SYMG	17	-2187.0621	4409.0987	0.0000	1.0000
SYMI	17	-2187.3269	4409.6283	0.0000	1.0000
SYMIG	18	-2186.5759	4410.2427	0.0000	1.0000
K2PG	13	-2192.5586	4411.6932	0.0000	1.0000
K2PI	13	-2192.7983	4412.1726	0.0000	1.0000
K2PIG	14	-2191.9545	4412.5746	0.0000	1.0000
SYM	16	-2199.8321	4432.5291	0.0000	1.0000
K2P	12	-2204.6722	4433.8373	0.0000	1.0000
JC69G	12	-2224.0628	4472.6186	0.0000	1.0000
JC69I	12	-2224.2733	4473.0396	0.0000	1.0000
JC69IG	13	-2223.9590	4474.4940	0.0000	1.0000
JC69	11	-2235.8298	4494.0760	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2101.0352	4234.0705	0.3772	0.3772
HKYI	16	-2101.4363	4234.8727	0.2525	0.6297
HKYIG	17	-2100.4872	4234.9743	0.2400	0.8697
GTRG	20	-2099.0445	4238.0889	0.0506	0.9203
GTRIG	21	-2098.2710	4238.5421	0.0403	0.9606
GTRI	20	-2099.2944	4238.5889	0.0394	1.0000
HKY	15	-2113.0443	4256.0886	0.0000	1.0000
GTR	19	-2111.2303	4260.4607	0.0000	1.0000
F81G	15	-2149.1803	4328.3606	0.0000	1.0000
F81I	15	-2149.5260	4329.0519	0.0000	1.0000
F81IG	16	-2149.1377	4330.2754	0.0000	1.0000
F81	14	-2160.1501	4348.3002	0.0000	1.0000
SYMG	17	-2187.0621	4408.1242	0.0000	1.0000
SYMI	17	-2187.3269	4408.6538	0.0000	1.0000
SYMIG	18	-2186.5759	4409.1518	0.0000	1.0000
K2PG	13	-2192.5586	4411.1172	0.0000	1.0000
K2PI	13	-2192.7983	4411.5966	0.0000	1.0000
K2PIG	14	-2191.9545	4411.9090	0.0000	1.0000
SYM	16	-2199.8321	4431.6643	0.0000	1.0000
K2P	12	-2204.6722	4433.3444	0.0000	1.0000
JC69G	12	-2224.0628	4472.1257	0.0000	1.0000
JC69I	12	-2224.2733	4472.5467	0.0000	1.0000
JC69IG	13	-2223.9590	4473.9181	0.0000	1.0000
JC69	11	-2235.8298	4493.6596	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2101.0352	4305.6033	0.5754	0.5754
HKYI	16	-2101.4363	4306.4055	0.3853	0.9607
HKYIG	17	-2100.4872	4310.9779	0.0392	0.9999
HKY	15	-2113.0443	4323.1506	0.0001	1.0000
GTRG	20	-2099.0445	4327.5049	0.0000	1.0000
GTRI	20	-2099.2944	4328.0049	0.0000	1.0000
GTRIG	21	-2098.2710	4332.4289	0.0000	1.0000
GTR	19	-2111.2303	4345.4059	0.0000	1.0000
F81G	15	-2149.1803	4395.4226	0.0000	1.0000
F81I	15	-2149.5260	4396.1139	0.0000	1.0000
F81IG	16	-2149.1377	4401.8082	0.0000	1.0000
F81	14	-2160.1501	4410.8914	0.0000	1.0000
K2PG	13	-2192.5586	4469.2376	0.0000	1.0000
K2PI	13	-2192.7983	4469.7170	0.0000	1.0000
K2PIG	14	-2191.9545	4474.5002	0.0000	1.0000
SYMG	17	-2187.0621	4484.1278	0.0000	1.0000
SYMI	17	-2187.3269	4484.6574	0.0000	1.0000
K2P	12	-2204.6722	4486.9940	0.0000	1.0000
SYMIG	18	-2186.5759	4489.6262	0.0000	1.0000
SYM	16	-2199.8321	4503.1971	0.0000	1.0000
JC69G	12	-2224.0628	4525.7753	0.0000	1.0000
JC69I	12	-2224.2733	4526.1963	0.0000	1.0000
JC69IG	13	-2223.9590	4532.0385	0.0000	1.0000
JC69	11	-2235.8298	4542.8384	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_7148.phylip = 1 - 646;
 Partition Dummy = 1:591-tomodel-phylip/chr9_7148.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_7148.phylip = 1 - 646;
 Partition Dummy = 1:591-tomodel-phylip/chr9_7148.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_7148.phylip = 1 - 646;
 Partition Dummy = 1:591-tomodel-phylip/chr9_7148.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

