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

Input data from file "261-aggressive-phylip/9230_hg18.phylip" (29 taxa, 631 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	64	-2469.4248	5081.5492	0.9912	0.9912
GTRIG	65	-2472.9107	5091.0072	0.0088	0.9999
HKYIG	61	-2482.8012	5100.8959	0.0001	1.0000
HKYG	60	-2484.9663	5102.7747	0.0000	1.0000
SYMG	61	-2485.7199	5106.7334	0.0000	1.0000
GTRI	64	-2485.9946	5114.6888	0.0000	1.0000
SYMIG	62	-2489.3185	5116.3905	0.0000	1.0000
K2PIG	58	-2498.5673	5125.0995	0.0000	1.0000
K2PG	57	-2500.5931	5126.7254	0.0000	1.0000
SYMI	61	-2496.7649	5128.8234	0.0000	1.0000
HKYI	60	-2499.8517	5132.5455	0.0000	1.0000
K2PI	57	-2513.3313	5152.2018	0.0000	1.0000
GTR	63	-2556.1892	5252.6007	0.0000	1.0000
SYM	60	-2573.4548	5279.7517	0.0000	1.0000
HKY	59	-2581.6137	5293.6268	0.0000	1.0000
F81G	59	-2589.8865	5310.1723	0.0000	1.0000
K2P	56	-2599.5623	5322.2466	0.0000	1.0000
F81IG	60	-2595.1930	5323.2281	0.0000	1.0000
F81I	59	-2606.0439	5342.4872	0.0000	1.0000
JC69G	56	-2611.0399	5345.2018	0.0000	1.0000
JC69IG	57	-2613.4187	5352.3767	0.0000	1.0000
JC69I	56	-2623.0345	5369.1909	0.0000	1.0000
F81	58	-2682.3248	5492.6146	0.0000	1.0000
JC69	55	-2698.3350	5517.3830	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2469.4248	5066.8496	0.9888	0.9888
GTRIG	65	-2472.9107	5075.8213	0.0111	1.0000
HKYIG	61	-2482.8012	5087.6024	0.0000	1.0000
HKYG	60	-2484.9663	5089.9326	0.0000	1.0000
SYMG	61	-2485.7199	5093.4399	0.0000	1.0000
GTRI	64	-2485.9946	5099.9891	0.0000	1.0000
SYMIG	62	-2489.3185	5102.6370	0.0000	1.0000
K2PIG	58	-2498.5673	5113.1345	0.0000	1.0000
K2PG	57	-2500.5931	5115.1861	0.0000	1.0000
SYMI	61	-2496.7649	5115.5299	0.0000	1.0000
HKYI	60	-2499.8517	5119.7034	0.0000	1.0000
K2PI	57	-2513.3313	5140.6625	0.0000	1.0000
GTR	63	-2556.1892	5238.3784	0.0000	1.0000
SYM	60	-2573.4548	5266.9096	0.0000	1.0000
HKY	59	-2581.6137	5281.2275	0.0000	1.0000
F81G	59	-2589.8865	5297.7730	0.0000	1.0000
F81IG	60	-2595.1930	5310.3860	0.0000	1.0000
K2P	56	-2599.5623	5311.1247	0.0000	1.0000
F81I	59	-2606.0439	5330.0879	0.0000	1.0000
JC69G	56	-2611.0399	5334.0798	0.0000	1.0000
JC69IG	57	-2613.4187	5340.8375	0.0000	1.0000
JC69I	56	-2623.0345	5358.0689	0.0000	1.0000
F81	58	-2682.3248	5480.6496	0.0000	1.0000
JC69	55	-2698.3350	5506.6699	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-2469.4248	5351.4772	0.9105	0.9105
HKYG	60	-2484.9663	5356.7710	0.0645	0.9751
HKYIG	61	-2482.8012	5358.8881	0.0224	0.9975
SYMG	61	-2485.7199	5364.7255	0.0012	0.9987
GTRIG	65	-2472.9107	5364.8962	0.0011	0.9998
K2PG	57	-2500.5931	5368.6826	0.0002	0.9999
K2PIG	58	-2498.5673	5371.0782	0.0001	1.0000
SYMIG	62	-2489.3185	5378.3699	0.0000	1.0000
GTRI	64	-2485.9946	5384.6167	0.0000	1.0000
HKYI	60	-2499.8517	5386.5418	0.0000	1.0000
SYMI	61	-2496.7649	5386.8155	0.0000	1.0000
K2PI	57	-2513.3313	5394.1590	0.0000	1.0000
GTR	63	-2556.1892	5518.5587	0.0000	1.0000
SYM	60	-2573.4548	5533.7480	0.0000	1.0000
HKY	59	-2581.6137	5543.6185	0.0000	1.0000
F81G	59	-2589.8865	5560.1641	0.0000	1.0000
K2P	56	-2599.5623	5560.1738	0.0000	1.0000
F81IG	60	-2595.1930	5577.2243	0.0000	1.0000
JC69G	56	-2611.0399	5583.1290	0.0000	1.0000
F81I	59	-2606.0439	5592.4789	0.0000	1.0000
JC69IG	57	-2613.4187	5594.3339	0.0000	1.0000
JC69I	56	-2623.0345	5607.1181	0.0000	1.0000
F81	58	-2682.3248	5738.5933	0.0000	1.0000
JC69	55	-2698.3350	5751.2717	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/9230_hg18.phylip = 1 - 631;
 Partition Dummy = 1:261-aggressive-phylip/9230_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/9230_hg18.phylip = 1 - 631;
 Partition Dummy = 1:261-aggressive-phylip/9230_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 (GTRG)]

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


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

