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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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
GTRG	20	-2322.8144	4686.8253	0.5670	0.5670
GTRI	20	-2323.6285	4688.4535	0.2512	0.8182
GTRIG	21	-2322.8915	4689.1012	0.1817	1.0000
HKYG	16	-2337.4884	4707.7473	0.0000	1.0000
HKYI	16	-2337.6788	4708.1282	0.0000	1.0000
HKYIG	17	-2337.5987	4710.0654	0.0000	1.0000
SYMG	17	-2338.7675	4712.4032	0.0000	1.0000
SYMI	17	-2338.8479	4712.5638	0.0000	1.0000
SYMIG	18	-2338.8783	4714.7282	0.0000	1.0000
K2PI	13	-2347.4688	4721.4509	0.0000	1.0000
K2PG	13	-2347.4969	4721.5072	0.0000	1.0000
K2PIG	14	-2347.6170	4723.8272	0.0000	1.0000
GTR	19	-2358.4221	4755.9253	0.0000	1.0000
HKY	15	-2371.4310	4773.5408	0.0000	1.0000
SYM	16	-2372.0903	4776.9512	0.0000	1.0000
K2P	12	-2380.9166	4786.2726	0.0000	1.0000
F81I	15	-2380.2940	4791.2670	0.0000	1.0000
F81G	15	-2380.6998	4792.0786	0.0000	1.0000
F81IG	16	-2380.8025	4794.3756	0.0000	1.0000
JC69I	12	-2387.1779	4798.7952	0.0000	1.0000
JC69G	12	-2387.6689	4799.7773	0.0000	1.0000
JC69IG	13	-2387.7773	4802.0680	0.0000	1.0000
F81	14	-2409.0532	4846.6996	0.0000	1.0000
JC69	11	-2416.2066	4854.7846	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2322.8144	4685.6287	0.5607	0.5607
GTRI	20	-2323.6285	4687.2570	0.2484	0.8090
GTRIG	21	-2322.8915	4687.7831	0.1909	1.0000
HKYG	16	-2337.4884	4706.9768	0.0000	1.0000
HKYI	16	-2337.6788	4707.3577	0.0000	1.0000
HKYIG	17	-2337.5987	4709.1973	0.0000	1.0000
SYMG	17	-2338.7675	4711.5351	0.0000	1.0000
SYMI	17	-2338.8479	4711.6957	0.0000	1.0000
SYMIG	18	-2338.8783	4713.7566	0.0000	1.0000
K2PI	13	-2347.4688	4720.9375	0.0000	1.0000
K2PG	13	-2347.4969	4720.9938	0.0000	1.0000
K2PIG	14	-2347.6170	4723.2340	0.0000	1.0000
GTR	19	-2358.4221	4754.8442	0.0000	1.0000
HKY	15	-2371.4310	4772.8619	0.0000	1.0000
SYM	16	-2372.0903	4776.1807	0.0000	1.0000
K2P	12	-2380.9166	4785.8332	0.0000	1.0000
F81I	15	-2380.2940	4790.5881	0.0000	1.0000
F81G	15	-2380.6998	4791.3996	0.0000	1.0000
F81IG	16	-2380.8025	4793.6051	0.0000	1.0000
JC69I	12	-2387.1779	4798.3557	0.0000	1.0000
JC69G	12	-2387.6689	4799.3379	0.0000	1.0000
JC69IG	13	-2387.7773	4801.5546	0.0000	1.0000
F81	14	-2409.0532	4846.1064	0.0000	1.0000
JC69	11	-2416.2066	4854.4133	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-2322.8144	4777.2969	0.4360	0.4360
GTRI	20	-2323.6285	4778.9251	0.1932	0.6292
HKYG	16	-2337.4884	4780.3113	0.0966	0.7257
K2PI	13	-2347.4688	4780.5219	0.0869	0.8127
K2PG	13	-2347.4969	4780.5781	0.0845	0.8972
HKYI	16	-2337.6788	4780.6922	0.0798	0.9770
GTRIG	21	-2322.8915	4784.0347	0.0150	0.9920
HKYIG	17	-2337.5987	4787.1153	0.0032	0.9953
K2PIG	14	-2347.6170	4787.4017	0.0028	0.9980
SYMG	17	-2338.7675	4789.4530	0.0010	0.9990
SYMI	17	-2338.8479	4789.6137	0.0009	1.0000
SYMIG	18	-2338.8783	4796.2579	0.0000	1.0000
K2P	12	-2380.9166	4840.8341	0.0000	1.0000
HKY	15	-2371.4310	4841.6131	0.0000	1.0000
GTR	19	-2358.4221	4841.9290	0.0000	1.0000
SYM	16	-2372.0903	4849.5152	0.0000	1.0000
JC69I	12	-2387.1779	4853.3566	0.0000	1.0000
JC69G	12	-2387.6689	4854.3388	0.0000	1.0000
F81I	15	-2380.2940	4859.3392	0.0000	1.0000
F81G	15	-2380.6998	4860.1508	0.0000	1.0000
JC69IG	13	-2387.7773	4861.1389	0.0000	1.0000
F81IG	16	-2380.8025	4866.9396	0.0000	1.0000
JC69	11	-2416.2066	4904.8308	0.0000	1.0000
F81	14	-2409.0532	4910.2741	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr6_8088.phylip = 1 - 723;
 Partition Dummy = 1:683-tomodel-phylip/chr6_8088.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr6_8088.phylip = 1 - 723;
 Partition Dummy = 1:683-tomodel-phylip/chr6_8088.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr6_8088.phylip = 1 - 723;
 Partition Dummy = 1:683-tomodel-phylip/chr6_8088.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

