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

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

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
Minimum BIC  model: K2PG


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
HKYG	60	-1705.1635	3542.9696	0.4612	0.4612
GTRG	64	-1700.7195	3543.9085	0.2884	0.7495
HKYIG	61	-1705.1679	3545.4223	0.1353	0.8848
GTRIG	65	-1700.7240	3546.3958	0.0831	0.9680
K2PG	57	-1711.8091	3548.9790	0.0229	0.9908
K2PIG	58	-1711.8134	3551.4065	0.0068	0.9976
SYMG	61	-1709.4659	3554.0182	0.0018	0.9995
SYMIG	62	-1709.4703	3556.4795	0.0005	1.0000
GTR	63	-1722.0473	3584.0945	0.0000	1.0000
HKY	59	-1727.4444	3585.0957	0.0000	1.0000
GTRI	64	-1722.0495	3586.5686	0.0000	1.0000
HKYI	60	-1727.4463	3587.5351	0.0000	1.0000
F81G	59	-1729.2358	3588.6784	0.0000	1.0000
K2P	56	-1734.0721	3591.0945	0.0000	1.0000
F81IG	60	-1729.2401	3591.1228	0.0000	1.0000
JC69G	56	-1734.8156	3592.5814	0.0000	1.0000
K2PI	57	-1734.0738	3593.5085	0.0000	1.0000
SYM	60	-1730.4695	3593.5814	0.0000	1.0000
JC69IG	57	-1734.8199	3595.0006	0.0000	1.0000
SYMI	61	-1730.4715	3596.0296	0.0000	1.0000
F81	58	-1750.9478	3629.6754	0.0000	1.0000
F81I	59	-1750.9496	3632.1062	0.0000	1.0000
JC69	55	-1756.6522	3633.8523	0.0000	1.0000
JC69I	56	-1756.6540	3636.2582	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1700.7195	3529.4390	0.4406	0.4406
HKYG	60	-1705.1635	3530.3271	0.2826	0.7232
GTRIG	65	-1700.7240	3531.4481	0.1613	0.8845
HKYIG	61	-1705.1679	3532.3358	0.1035	0.9880
K2PG	57	-1711.8091	3537.6182	0.0074	0.9954
K2PIG	58	-1711.8134	3539.6268	0.0027	0.9981
SYMG	61	-1709.4659	3540.9317	0.0014	0.9995
SYMIG	62	-1709.4703	3542.9405	0.0005	1.0000
GTR	63	-1722.0473	3570.0945	0.0000	1.0000
GTRI	64	-1722.0495	3572.0990	0.0000	1.0000
HKY	59	-1727.4444	3572.8888	0.0000	1.0000
HKYI	60	-1727.4463	3574.8927	0.0000	1.0000
F81G	59	-1729.2358	3576.4715	0.0000	1.0000
F81IG	60	-1729.2401	3578.4803	0.0000	1.0000
K2P	56	-1734.0721	3580.1443	0.0000	1.0000
SYM	60	-1730.4695	3580.9389	0.0000	1.0000
JC69G	56	-1734.8156	3581.6311	0.0000	1.0000
K2PI	57	-1734.0738	3582.1477	0.0000	1.0000
SYMI	61	-1730.4715	3582.9431	0.0000	1.0000
JC69IG	57	-1734.8199	3583.6398	0.0000	1.0000
F81	58	-1750.9478	3617.8957	0.0000	1.0000
F81I	59	-1750.9496	3619.8993	0.0000	1.0000
JC69	55	-1756.6522	3623.3044	0.0000	1.0000
JC69I	56	-1756.6540	3625.3079	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1711.8091	3791.9219	0.9185	0.9185
HKYG	60	-1705.1635	3798.0152	0.0436	0.9621
K2PIG	58	-1711.8134	3798.3920	0.0361	0.9982
HKYIG	61	-1705.1679	3804.4853	0.0017	1.0000
SYMG	61	-1709.4659	3813.0813	0.0000	1.0000
GTRG	64	-1700.7195	3814.9729	0.0000	1.0000
SYMIG	62	-1709.4703	3819.5515	0.0000	1.0000
GTRIG	65	-1700.7240	3821.4435	0.0000	1.0000
K2P	56	-1734.0721	3829.9865	0.0000	1.0000
JC69G	56	-1734.8156	3831.4733	0.0000	1.0000
HKY	59	-1727.4444	3836.1154	0.0000	1.0000
K2PI	57	-1734.0738	3836.4514	0.0000	1.0000
JC69IG	57	-1734.8199	3837.9435	0.0000	1.0000
F81G	59	-1729.2358	3839.6981	0.0000	1.0000
HKYI	60	-1727.4463	3842.5808	0.0000	1.0000
F81IG	60	-1729.2401	3846.1684	0.0000	1.0000
SYM	60	-1730.4695	3848.6270	0.0000	1.0000
GTR	63	-1722.0473	3851.1670	0.0000	1.0000
SYMI	61	-1730.4715	3855.0926	0.0000	1.0000
GTRI	64	-1722.0495	3857.6330	0.0000	1.0000
JC69	55	-1756.6522	3868.6851	0.0000	1.0000
JC69I	56	-1756.6540	3875.1502	0.0000	1.0000
F81	58	-1750.9478	3876.6608	0.0000	1.0000
F81I	59	-1750.9496	3883.1259	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/7258_hg18.phylip = 1 - 640;
 Partition Dummy = 1:261-aggressive-phylip/7258_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 (HKYG)]

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


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7258_hg18.phylip = 1 - 640;
 Partition Dummy = 1:261-aggressive-phylip/7258_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

