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

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

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


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	-2524.3820	5089.9866	0.4286	0.4286
GTRG	20	-2524.4619	5090.1465	0.3957	0.8243
GTRIG	21	-2524.2458	5091.8384	0.1698	0.9941
HKYI	16	-2533.7452	5100.2776	0.0025	0.9966
HKYG	16	-2533.7958	5100.3789	0.0024	0.9990
HKYIG	17	-2533.6214	5102.1298	0.0010	1.0000
GTR	19	-2538.2638	5115.6323	0.0000	1.0000
HKY	15	-2546.3303	5123.3542	0.0000	1.0000
SYMI	17	-2550.6544	5136.1957	0.0000	1.0000
SYMG	17	-2550.8610	5136.6090	0.0000	1.0000
SYMIG	18	-2550.6060	5138.2048	0.0000	1.0000
K2PI	13	-2556.6128	5139.7501	0.0000	1.0000
K2PG	13	-2556.7995	5140.1236	0.0000	1.0000
K2PIG	14	-2556.5601	5141.7263	0.0000	1.0000
SYM	16	-2563.1954	5159.1780	0.0000	1.0000
K2P	12	-2568.6549	5161.7586	0.0000	1.0000
F81G	15	-2588.0137	5206.7210	0.0000	1.0000
F81I	15	-2588.0808	5206.8552	0.0000	1.0000
F81IG	16	-2587.9422	5208.6717	0.0000	1.0000
F81	14	-2598.5979	5225.8019	0.0000	1.0000
JC69G	12	-2605.8907	5236.2302	0.0000	1.0000
JC69I	12	-2605.9534	5236.3558	0.0000	1.0000
JC69IG	13	-2605.8215	5238.1675	0.0000	1.0000
JC69	11	-2616.4337	5255.2467	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2524.3820	5088.7639	0.4245	0.4245
GTRG	20	-2524.4619	5088.9238	0.3919	0.8163
GTRIG	21	-2524.2458	5090.4915	0.1789	0.9953
HKYI	16	-2533.7452	5099.4903	0.0020	0.9973
HKYG	16	-2533.7958	5099.5916	0.0019	0.9992
HKYIG	17	-2533.6214	5101.2428	0.0008	1.0000
GTR	19	-2538.2638	5114.5277	0.0000	1.0000
HKY	15	-2546.3303	5122.6605	0.0000	1.0000
SYMI	17	-2550.6544	5135.3087	0.0000	1.0000
SYMG	17	-2550.8610	5135.7221	0.0000	1.0000
SYMIG	18	-2550.6060	5137.2121	0.0000	1.0000
K2PI	13	-2556.6128	5139.2256	0.0000	1.0000
K2PG	13	-2556.7995	5139.5991	0.0000	1.0000
K2PIG	14	-2556.5601	5141.1203	0.0000	1.0000
SYM	16	-2563.1954	5158.3908	0.0000	1.0000
K2P	12	-2568.6549	5161.3097	0.0000	1.0000
F81G	15	-2588.0137	5206.0274	0.0000	1.0000
F81I	15	-2588.0808	5206.1616	0.0000	1.0000
F81IG	16	-2587.9422	5207.8845	0.0000	1.0000
F81	14	-2598.5979	5225.1958	0.0000	1.0000
JC69G	12	-2605.8907	5235.7813	0.0000	1.0000
JC69I	12	-2605.9534	5235.9068	0.0000	1.0000
JC69IG	13	-2605.8215	5237.6430	0.0000	1.0000
JC69	11	-2616.4337	5254.8674	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2533.7452	5172.4894	0.4905	0.4905
HKYG	16	-2533.7958	5172.5907	0.4662	0.9567
HKYIG	17	-2533.6214	5178.8044	0.0209	0.9775
GTRI	20	-2524.3820	5180.0128	0.0114	0.9889
GTRG	20	-2524.4619	5180.1726	0.0105	0.9995
GTRIG	21	-2524.2458	5186.3028	0.0005	1.0000
HKY	15	-2546.3303	5191.0972	0.0000	1.0000
K2PI	13	-2556.6128	5198.5374	0.0000	1.0000
K2PG	13	-2556.7995	5198.9109	0.0000	1.0000
GTR	19	-2538.2638	5201.2141	0.0000	1.0000
K2PIG	14	-2556.5601	5204.9945	0.0000	1.0000
SYMI	17	-2550.6544	5212.8703	0.0000	1.0000
SYMG	17	-2550.8610	5213.2836	0.0000	1.0000
K2P	12	-2568.6549	5216.0590	0.0000	1.0000
SYMIG	18	-2550.6060	5219.3361	0.0000	1.0000
SYM	16	-2563.1954	5231.3899	0.0000	1.0000
F81G	15	-2588.0137	5274.4640	0.0000	1.0000
F81I	15	-2588.0808	5274.5982	0.0000	1.0000
F81IG	16	-2587.9422	5280.8836	0.0000	1.0000
F81	14	-2598.5979	5289.0701	0.0000	1.0000
JC69G	12	-2605.8907	5290.5306	0.0000	1.0000
JC69I	12	-2605.9534	5290.6562	0.0000	1.0000
JC69IG	13	-2605.8215	5296.9548	0.0000	1.0000
JC69	11	-2616.4337	5305.0543	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 591-tomodel-phylip/chr7_10302.phylip = 1 - 708;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10302.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 591-tomodel-phylip/chr7_10302.phylip = 1 - 708;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10302.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 591-tomodel-phylip/chr7_10302.phylip = 1 - 708;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10302.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

