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

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

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


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
GTRI	20	-2020.5485	4082.3377	0.4186	0.4186
GTRG	20	-2020.6212	4082.4832	0.3893	0.8079
GTRIG	21	-2020.2646	4083.8960	0.1921	1.0000
HKYI	16	-2038.8354	4110.4695	0.0000	1.0000
HKYG	16	-2038.8827	4110.5642	0.0000	1.0000
HKYIG	17	-2038.6501	4112.2003	0.0000	1.0000
GTR	19	-2042.3619	4123.8448	0.0000	1.0000
HKY	15	-2059.5354	4149.7746	0.0000	1.0000
SYMI	17	-2087.1797	4209.2593	0.0000	1.0000
SYMG	17	-2087.6909	4210.2818	0.0000	1.0000
SYMIG	18	-2087.1836	4211.3746	0.0000	1.0000
K2PI	13	-2094.9780	4216.4881	0.0000	1.0000
K2PG	13	-2095.3037	4217.1395	0.0000	1.0000
K2PIG	14	-2094.9803	4218.5756	0.0000	1.0000
F81I	15	-2100.0904	4230.8847	0.0000	1.0000
F81G	15	-2100.6138	4231.9315	0.0000	1.0000
F81IG	16	-2100.0942	4232.9872	0.0000	1.0000
SYM	16	-2108.3635	4249.5258	0.0000	1.0000
K2P	12	-2116.4345	4257.3244	0.0000	1.0000
F81	14	-2119.1768	4266.9685	0.0000	1.0000
JC69I	12	-2143.7348	4311.9251	0.0000	1.0000
JC69G	12	-2144.3007	4313.0568	0.0000	1.0000
JC69IG	13	-2144.3574	4315.2470	0.0000	1.0000
JC69	11	-2163.4374	4349.2596	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2020.5485	4081.0970	0.4135	0.4135
GTRG	20	-2020.6212	4081.2424	0.3845	0.7979
GTRIG	21	-2020.2646	4082.5291	0.2021	1.0000
HKYI	16	-2038.8354	4109.6707	0.0000	1.0000
HKYG	16	-2038.8827	4109.7654	0.0000	1.0000
HKYIG	17	-2038.6501	4111.3003	0.0000	1.0000
GTR	19	-2042.3619	4122.7239	0.0000	1.0000
HKY	15	-2059.5354	4149.0708	0.0000	1.0000
SYMI	17	-2087.1797	4208.3593	0.0000	1.0000
SYMG	17	-2087.6909	4209.3818	0.0000	1.0000
SYMIG	18	-2087.1836	4210.3672	0.0000	1.0000
K2PI	13	-2094.9780	4215.9560	0.0000	1.0000
K2PG	13	-2095.3037	4216.6073	0.0000	1.0000
K2PIG	14	-2094.9803	4217.9607	0.0000	1.0000
F81I	15	-2100.0904	4230.1809	0.0000	1.0000
F81G	15	-2100.6138	4231.2276	0.0000	1.0000
F81IG	16	-2100.0942	4232.1884	0.0000	1.0000
SYM	16	-2108.3635	4248.7270	0.0000	1.0000
K2P	12	-2116.4345	4256.8689	0.0000	1.0000
F81	14	-2119.1768	4266.3535	0.0000	1.0000
JC69I	12	-2143.7348	4311.4697	0.0000	1.0000
JC69G	12	-2144.3007	4312.6013	0.0000	1.0000
JC69IG	13	-2144.3574	4314.7148	0.0000	1.0000
JC69	11	-2163.4374	4348.8748	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	20	-2020.5485	4172.0613	0.5022	0.5022
GTRG	20	-2020.6212	4172.2068	0.4670	0.9692
GTRIG	21	-2020.2646	4178.0417	0.0252	0.9944
HKYI	16	-2038.8354	4182.4422	0.0028	0.9972
HKYG	16	-2038.8827	4182.5369	0.0027	0.9999
HKYIG	17	-2038.6501	4188.6200	0.0001	1.0000
GTR	19	-2042.3619	4209.1400	0.0000	1.0000
HKY	15	-2059.5354	4217.2941	0.0000	1.0000
K2PI	13	-2094.9780	4275.0828	0.0000	1.0000
K2PG	13	-2095.3037	4275.7342	0.0000	1.0000
K2PIG	14	-2094.9803	4281.6358	0.0000	1.0000
SYMI	17	-2087.1797	4285.6791	0.0000	1.0000
SYMG	17	-2087.6909	4286.7015	0.0000	1.0000
SYMIG	18	-2087.1836	4292.2352	0.0000	1.0000
F81I	15	-2100.0904	4298.4042	0.0000	1.0000
F81G	15	-2100.6138	4299.4509	0.0000	1.0000
F81IG	16	-2100.0942	4304.9599	0.0000	1.0000
K2P	12	-2116.4345	4311.4475	0.0000	1.0000
SYM	16	-2108.3635	4321.4985	0.0000	1.0000
F81	14	-2119.1768	4330.0286	0.0000	1.0000
JC69I	12	-2143.7348	4366.0483	0.0000	1.0000
JC69G	12	-2144.3007	4367.1800	0.0000	1.0000
JC69IG	13	-2144.3574	4373.8417	0.0000	1.0000
JC69	11	-2163.4374	4398.9052	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/chr6_8033.phylip = 1 - 698;
 Partition Dummy = 1:683-tomodel-phylip/chr6_8033.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/chr6_8033.phylip = 1 - 698;
 Partition Dummy = 1:683-tomodel-phylip/chr6_8033.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 (GTRI)]

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


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

