
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:30 08/01/2011

Input data from file "261-aggressive-phylip/3166_hg18.phylip" (29 taxa, 735 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-2628.6006	5397.6191	0.5331	0.5331
GTRIG	65	-2627.6289	5398.0830	0.4227	0.9558
HKYG	60	-2636.4885	5403.8375	0.0238	0.9796
HKYIG	61	-2635.5354	5404.3101	0.0188	0.9984
GTRI	64	-2634.4321	5409.2821	0.0016	0.9999
HKYI	60	-2642.4118	5415.6842	0.0001	1.0000
K2PG	57	-2649.4153	5422.5972	0.0000	1.0000
K2PIG	58	-2648.5122	5423.1487	0.0000	1.0000
SYMG	61	-2645.6772	5424.5936	0.0000	1.0000
SYMIG	62	-2644.7604	5425.1459	0.0000	1.0000
K2PI	57	-2654.1165	5431.9996	0.0000	1.0000
SYMI	61	-2651.7863	5436.8119	0.0000	1.0000
GTR	63	-2684.0113	5506.0405	0.0000	1.0000
HKY	59	-2690.9263	5510.3415	0.0000	1.0000
K2P	56	-2701.1032	5523.6222	0.0000	1.0000
SYM	60	-2698.6589	5528.1783	0.0000	1.0000
F81IG	60	-2699.0717	5529.0040	0.0000	1.0000
F81G	59	-2700.9442	5530.3773	0.0000	1.0000
F81I	59	-2704.4990	5537.4869	0.0000	1.0000
JC69G	56	-2710.5578	5542.5314	0.0000	1.0000
JC69IG	57	-2709.7706	5543.3078	0.0000	1.0000
JC69I	56	-2713.5823	5548.5805	0.0000	1.0000
F81	58	-2750.7715	5627.6673	0.0000	1.0000
JC69	55	-2757.6284	5634.3289	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2628.6006	5385.2012	0.4963	0.4963
GTRIG	65	-2627.6289	5385.2579	0.4824	0.9786
HKYG	60	-2636.4885	5392.9770	0.0102	0.9888
HKYIG	61	-2635.5354	5393.0709	0.0097	0.9985
GTRI	64	-2634.4321	5396.8642	0.0015	1.0000
HKYI	60	-2642.4118	5404.8236	0.0000	1.0000
K2PG	57	-2649.4153	5412.8306	0.0000	1.0000
K2PIG	58	-2648.5122	5413.0245	0.0000	1.0000
SYMG	61	-2645.6772	5413.3544	0.0000	1.0000
SYMIG	62	-2644.7604	5413.5209	0.0000	1.0000
K2PI	57	-2654.1165	5422.2329	0.0000	1.0000
SYMI	61	-2651.7863	5425.5726	0.0000	1.0000
GTR	63	-2684.0113	5494.0227	0.0000	1.0000
HKY	59	-2690.9263	5499.8526	0.0000	1.0000
K2P	56	-2701.1032	5514.2063	0.0000	1.0000
SYM	60	-2698.6589	5517.3178	0.0000	1.0000
F81IG	60	-2699.0717	5518.1434	0.0000	1.0000
F81G	59	-2700.9442	5519.8884	0.0000	1.0000
F81I	59	-2704.4990	5526.9980	0.0000	1.0000
JC69G	56	-2710.5578	5533.1155	0.0000	1.0000
JC69IG	57	-2709.7706	5533.5412	0.0000	1.0000
JC69I	56	-2713.5823	5539.1645	0.0000	1.0000
F81	58	-2750.7715	5617.5430	0.0000	1.0000
JC69	55	-2757.6284	5625.2568	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2636.4885	5668.9692	0.8642	0.8642
HKYIG	61	-2635.5354	5673.6630	0.0827	0.9469
K2PG	57	-2649.4153	5675.0232	0.0419	0.9888
GTRG	64	-2628.6006	5679.5929	0.0043	0.9931
K2PIG	58	-2648.5122	5679.8169	0.0038	0.9969
HKYI	60	-2642.4118	5680.8158	0.0023	0.9992
GTRIG	65	-2627.6289	5684.2494	0.0004	0.9996
K2PI	57	-2654.1165	5684.4256	0.0004	1.0000
GTRI	64	-2634.4321	5691.2559	0.0000	1.0000
SYMG	61	-2645.6772	5693.9465	0.0000	1.0000
SYMIG	62	-2644.7604	5698.7129	0.0000	1.0000
SYMI	61	-2651.7863	5706.1647	0.0000	1.0000
HKY	59	-2690.9263	5771.2450	0.0000	1.0000
K2P	56	-2701.1032	5771.7991	0.0000	1.0000
GTR	63	-2684.0113	5783.8145	0.0000	1.0000
JC69G	56	-2710.5578	5790.7083	0.0000	1.0000
F81G	59	-2700.9442	5791.2807	0.0000	1.0000
SYM	60	-2698.6589	5793.3100	0.0000	1.0000
F81IG	60	-2699.0717	5794.1357	0.0000	1.0000
JC69IG	57	-2709.7706	5795.7338	0.0000	1.0000
JC69I	56	-2713.5823	5796.7573	0.0000	1.0000
F81I	59	-2704.4990	5798.3904	0.0000	1.0000
JC69	55	-2757.6284	5878.2496	0.0000	1.0000
F81	58	-2750.7715	5884.3355	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 261-aggressive-phylip/3166_hg18.phylip = 1 - 735;
 Partition Dummy = 1:261-aggressive-phylip/3166_hg18.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 261-aggressive-phylip/3166_hg18.phylip = 1 - 735;
 Partition Dummy = 1:261-aggressive-phylip/3166_hg18.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 (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/3166_hg18.phylip = 1 - 735;
 Partition Dummy = 1:261-aggressive-phylip/3166_hg18.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

