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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-1528.1904	3198.9734	0.6432	0.6432
HKYIG	61	-1528.1902	3201.7984	0.1566	0.7999
HKYI	60	-1530.0787	3202.7500	0.0973	0.8972
GTRG	64	-1524.6221	3203.2442	0.0760	0.9732
GTRIG	65	-1524.6207	3206.1379	0.0179	0.9911
GTRI	64	-1526.7697	3207.5395	0.0089	1.0000
HKY	59	-1542.3013	3224.3872	0.0000	1.0000
K2PG	57	-1546.1933	3226.6067	0.0000	1.0000
GTR	63	-1538.9495	3229.0204	0.0000	1.0000
K2PIG	58	-1546.1944	3229.3827	0.0000	1.0000
K2PI	57	-1547.8819	3229.9839	0.0000	1.0000
SYMG	61	-1542.8212	3231.0604	0.0000	1.0000
SYMIG	62	-1542.8223	3233.9054	0.0000	1.0000
SYMI	61	-1544.4328	3234.2837	0.0000	1.0000
K2P	56	-1560.8981	3253.2596	0.0000	1.0000
SYM	60	-1556.7837	3256.1599	0.0000	1.0000
F81G	59	-1561.9718	3263.7283	0.0000	1.0000
F81IG	60	-1561.7878	3266.1681	0.0000	1.0000
F81I	59	-1563.2158	3266.2162	0.0000	1.0000
F81	58	-1575.8865	3288.7669	0.0000	1.0000
JC69G	56	-1580.5645	3292.5925	0.0000	1.0000
JC69I	56	-1581.8220	3295.1074	0.0000	1.0000
JC69IG	57	-1580.5680	3295.3561	0.0000	1.0000
JC69	55	-1594.4455	3317.6145	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1528.1904	3176.3808	0.4026	0.4026
GTRG	64	-1524.6221	3177.2442	0.2615	0.6641
HKYIG	61	-1528.1902	3178.3804	0.1481	0.8122
GTRIG	65	-1524.6207	3179.2414	0.0963	0.9085
HKYI	60	-1530.0787	3180.1574	0.0609	0.9695
GTRI	64	-1526.7697	3181.5395	0.0305	1.0000
HKY	59	-1542.3013	3202.6025	0.0000	1.0000
GTR	63	-1538.9495	3203.8989	0.0000	1.0000
K2PG	57	-1546.1933	3206.3865	0.0000	1.0000
SYMG	61	-1542.8212	3207.6425	0.0000	1.0000
K2PIG	58	-1546.1944	3208.3888	0.0000	1.0000
SYMIG	62	-1542.8223	3209.6445	0.0000	1.0000
K2PI	57	-1547.8819	3209.7637	0.0000	1.0000
SYMI	61	-1544.4328	3210.8657	0.0000	1.0000
SYM	60	-1556.7837	3233.5673	0.0000	1.0000
K2P	56	-1560.8981	3233.7962	0.0000	1.0000
F81G	59	-1561.9718	3241.9437	0.0000	1.0000
F81IG	60	-1561.7878	3243.5755	0.0000	1.0000
F81I	59	-1563.2158	3244.4316	0.0000	1.0000
F81	58	-1575.8865	3267.7730	0.0000	1.0000
JC69G	56	-1580.5645	3273.1291	0.0000	1.0000
JC69IG	57	-1580.5680	3275.1359	0.0000	1.0000
JC69I	56	-1581.8220	3275.6440	0.0000	1.0000
JC69	55	-1594.4455	3298.8911	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1528.1904	3413.5754	0.8314	0.8314
HKYI	60	-1530.0787	3417.3520	0.1258	0.9573
HKYIG	61	-1528.1902	3419.5283	0.0424	0.9996
GTRG	64	-1524.6221	3430.2518	0.0002	0.9998
K2PG	57	-1546.1933	3431.7214	0.0001	0.9999
GTRI	64	-1526.7697	3434.5471	0.0000	1.0000
K2PI	57	-1547.8819	3435.0986	0.0000	1.0000
HKY	59	-1542.3013	3435.8439	0.0000	1.0000
GTRIG	65	-1524.6207	3436.2022	0.0000	1.0000
K2PIG	58	-1546.1944	3437.6769	0.0000	1.0000
SYMG	61	-1542.8212	3448.7903	0.0000	1.0000
SYMI	61	-1544.4328	3452.0135	0.0000	1.0000
GTR	63	-1538.9495	3452.9532	0.0000	1.0000
SYMIG	62	-1542.8223	3454.7456	0.0000	1.0000
K2P	56	-1560.8981	3455.1778	0.0000	1.0000
SYM	60	-1556.7837	3470.7619	0.0000	1.0000
F81G	59	-1561.9718	3475.1850	0.0000	1.0000
F81I	59	-1563.2158	3477.6729	0.0000	1.0000
F81IG	60	-1561.7878	3480.7701	0.0000	1.0000
JC69G	56	-1580.5645	3494.5107	0.0000	1.0000
JC69I	56	-1581.8220	3497.0256	0.0000	1.0000
F81	58	-1575.8865	3497.0611	0.0000	1.0000
JC69IG	57	-1580.5680	3500.4708	0.0000	1.0000
JC69	55	-1594.4455	3516.3195	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9746_hg18.phylip = 1 - 385;
 Partition Dummy = 1:261-aggressive-phylip/9746_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 AICc model (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9746_hg18.phylip = 1 - 385;
 Partition Dummy = 1:261-aggressive-phylip/9746_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 (HKYG)]

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

