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

Input data from file "261-aggressive-phylip/3412_hg18.phylip" (29 taxa, 577 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	-1438.4366	3011.0592	0.7080	0.7080
HKYIG	61	-1438.4453	3013.5780	0.2009	0.9089
GTRG	64	-1435.7121	3015.6743	0.0704	0.9793
GTRIG	65	-1435.7208	3018.2323	0.0196	0.9989
K2PG	57	-1448.9391	3024.6181	0.0008	0.9997
K2PIG	58	-1448.9479	3027.1082	0.0002	1.0000
SYMG	61	-1447.4327	3031.5527	0.0000	1.0000
SYMIG	62	-1447.4415	3034.0814	0.0000	1.0000
F81G	59	-1473.1797	3078.0537	0.0000	1.0000
F81IG	60	-1473.1883	3080.5627	0.0000	1.0000
JC69G	56	-1481.8274	3087.9316	0.0000	1.0000
JC69IG	57	-1481.8361	3090.4121	0.0000	1.0000
HKY	59	-1490.4910	3112.6764	0.0000	1.0000
HKYI	60	-1490.4951	3115.1762	0.0000	1.0000
GTR	63	-1488.1400	3117.9992	0.0000	1.0000
GTRI	64	-1488.1440	3120.5380	0.0000	1.0000
K2P	56	-1500.8752	3126.0273	0.0000	1.0000
K2PI	57	-1500.8792	3128.4983	0.0000	1.0000
SYM	60	-1499.5086	3133.2032	0.0000	1.0000
SYMI	61	-1499.5126	3135.7126	0.0000	1.0000
F81	58	-1524.8049	3178.8221	0.0000	1.0000
F81I	59	-1524.8090	3181.3124	0.0000	1.0000
JC69	55	-1533.4132	3188.6499	0.0000	1.0000
JC69I	56	-1533.4173	3191.1115	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1438.4366	2996.8731	0.5725	0.5725
HKYIG	61	-1438.4453	2998.8907	0.2088	0.7813
GTRG	64	-1435.7121	2999.4243	0.1599	0.9412
GTRIG	65	-1435.7208	3001.4417	0.0583	0.9995
K2PG	57	-1448.9391	3011.8783	0.0003	0.9998
K2PIG	58	-1448.9479	3013.8958	0.0001	1.0000
SYMG	61	-1447.4327	3016.8653	0.0000	1.0000
SYMIG	62	-1447.4415	3018.8829	0.0000	1.0000
F81G	59	-1473.1797	3064.3593	0.0000	1.0000
F81IG	60	-1473.1883	3066.3767	0.0000	1.0000
JC69G	56	-1481.8274	3075.6547	0.0000	1.0000
JC69IG	57	-1481.8361	3077.6722	0.0000	1.0000
HKY	59	-1490.4910	3098.9820	0.0000	1.0000
HKYI	60	-1490.4951	3100.9902	0.0000	1.0000
GTR	63	-1488.1400	3102.2799	0.0000	1.0000
GTRI	64	-1488.1440	3104.2880	0.0000	1.0000
K2P	56	-1500.8752	3113.7503	0.0000	1.0000
K2PI	57	-1500.8792	3115.7584	0.0000	1.0000
SYM	60	-1499.5086	3119.0172	0.0000	1.0000
SYMI	61	-1499.5126	3121.0252	0.0000	1.0000
F81	58	-1524.8049	3165.6098	0.0000	1.0000
F81I	59	-1524.8090	3167.6180	0.0000	1.0000
JC69	55	-1533.4132	3176.8264	0.0000	1.0000
JC69I	56	-1533.4173	3178.8346	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1438.4366	3258.3437	0.6956	0.6956
K2PG	57	-1448.9391	3260.2753	0.2648	0.9603
HKYIG	61	-1438.4453	3264.7190	0.0287	0.9890
K2PIG	58	-1448.9479	3266.6507	0.0109	1.0000
GTRG	64	-1435.7121	3278.3262	0.0000	1.0000
SYMG	61	-1447.4327	3282.6937	0.0000	1.0000
GTRIG	65	-1435.7208	3284.7014	0.0000	1.0000
SYMIG	62	-1447.4415	3289.0692	0.0000	1.0000
JC69G	56	-1481.8274	3319.6939	0.0000	1.0000
F81G	59	-1473.1797	3321.4720	0.0000	1.0000
JC69IG	57	-1481.8361	3326.0692	0.0000	1.0000
F81IG	60	-1473.1883	3327.8472	0.0000	1.0000
HKY	59	-1490.4910	3356.0947	0.0000	1.0000
K2P	56	-1500.8752	3357.7895	0.0000	1.0000
HKYI	60	-1490.4951	3362.4607	0.0000	1.0000
K2PI	57	-1500.8792	3364.1554	0.0000	1.0000
GTR	63	-1488.1400	3376.8240	0.0000	1.0000
SYM	60	-1499.5086	3380.4877	0.0000	1.0000
GTRI	64	-1488.1440	3383.1899	0.0000	1.0000
SYMI	61	-1499.5126	3386.8536	0.0000	1.0000
JC69	55	-1533.4132	3416.5078	0.0000	1.0000
F81	58	-1524.8049	3418.3646	0.0000	1.0000
JC69I	56	-1533.4173	3422.8738	0.0000	1.0000
F81I	59	-1524.8090	3424.7307	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/3412_hg18.phylip = 1 - 577;
 Partition Dummy = 1:261-aggressive-phylip/3412_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/3412_hg18.phylip = 1 - 577;
 Partition Dummy = 1:261-aggressive-phylip/3412_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/3412_hg18.phylip = 1 - 577;
 Partition Dummy = 1:261-aggressive-phylip/3412_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

