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

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

Minimum AIC  model: GTRG
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	-1324.0364	2786.7941	0.4156	0.4156
HKYIG	61	-1323.3387	2788.0722	0.2193	0.6349
GTRG	64	-1319.6820	2788.8626	0.1477	0.7826
HKYI	60	-1325.3723	2789.4658	0.1093	0.8919
GTRIG	65	-1318.9226	2790.0731	0.0807	0.9725
GTRI	64	-1321.3644	2792.2275	0.0275	1.0000
K2PG	57	-1346.5306	2823.8428	0.0000	1.0000
K2PIG	58	-1345.7459	2824.9066	0.0000	1.0000
K2PI	57	-1347.4219	2825.6255	0.0000	1.0000
HKY	59	-1347.2008	2830.4629	0.0000	1.0000
GTR	63	-1342.8058	2832.3952	0.0000	1.0000
SYMG	61	-1345.5336	2832.4620	0.0000	1.0000
SYMIG	62	-1344.7298	2833.5419	0.0000	1.0000
SYMI	61	-1346.5086	2834.4121	0.0000	1.0000
F81G	59	-1367.5854	2871.2320	0.0000	1.0000
F81IG	60	-1366.8924	2872.5061	0.0000	1.0000
K2P	56	-1372.2294	2872.6207	0.0000	1.0000
F81I	59	-1368.7255	2873.5123	0.0000	1.0000
SYM	60	-1370.8547	2880.4306	0.0000	1.0000
JC69G	56	-1386.8377	2901.8373	0.0000	1.0000
JC69IG	57	-1386.0927	2902.9672	0.0000	1.0000
JC69I	56	-1387.7977	2903.7574	0.0000	1.0000
F81	58	-1391.4710	2916.3568	0.0000	1.0000
JC69	55	-1411.9597	2949.4749	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1319.6820	2767.3639	0.2962	0.2962
GTRIG	65	-1318.9226	2767.8451	0.2328	0.5290
HKYG	60	-1324.0364	2768.0729	0.2078	0.7367
HKYIG	61	-1323.3387	2768.6773	0.1536	0.8903
GTRI	64	-1321.3644	2770.7288	0.0551	0.9454
HKYI	60	-1325.3723	2770.7446	0.0546	1.0000
K2PG	57	-1346.5306	2807.0611	0.0000	1.0000
K2PIG	58	-1345.7459	2807.4919	0.0000	1.0000
K2PI	57	-1347.4219	2808.8438	0.0000	1.0000
GTR	63	-1342.8058	2811.6117	0.0000	1.0000
HKY	59	-1347.2008	2812.4017	0.0000	1.0000
SYMG	61	-1345.5336	2813.0671	0.0000	1.0000
SYMIG	62	-1344.7298	2813.4596	0.0000	1.0000
SYMI	61	-1346.5086	2815.0172	0.0000	1.0000
F81G	59	-1367.5854	2853.1707	0.0000	1.0000
F81IG	60	-1366.8924	2853.7849	0.0000	1.0000
F81I	59	-1368.7255	2855.4511	0.0000	1.0000
K2P	56	-1372.2294	2856.4587	0.0000	1.0000
SYM	60	-1370.8547	2861.7094	0.0000	1.0000
JC69G	56	-1386.8377	2885.6753	0.0000	1.0000
JC69IG	57	-1386.0927	2886.1855	0.0000	1.0000
JC69I	56	-1387.7977	2887.5954	0.0000	1.0000
F81	58	-1391.4710	2898.9421	0.0000	1.0000
JC69	55	-1411.9597	2933.9193	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1324.0364	3014.8938	0.7364	0.7364
HKYI	60	-1325.3723	3017.5655	0.1936	0.9300
HKYIG	61	-1323.3387	3019.6119	0.0696	0.9996
GTRG	64	-1319.6820	3030.6396	0.0003	0.9999
GTRI	64	-1321.3644	3034.0045	0.0001	1.0000
GTRIG	65	-1318.9226	3035.2345	0.0000	1.0000
K2PG	57	-1346.5306	3041.5410	0.0000	1.0000
K2PI	57	-1347.4219	3043.3237	0.0000	1.0000
K2PIG	58	-1345.7459	3046.0854	0.0000	1.0000
HKY	59	-1347.2008	3055.1089	0.0000	1.0000
SYMG	61	-1345.5336	3064.0017	0.0000	1.0000
SYMI	61	-1346.5086	3065.9518	0.0000	1.0000
SYMIG	62	-1344.7298	3068.5079	0.0000	1.0000
GTR	63	-1342.8058	3070.7737	0.0000	1.0000
K2P	56	-1372.2294	3086.8249	0.0000	1.0000
F81G	59	-1367.5854	3095.8780	0.0000	1.0000
F81I	59	-1368.7255	3098.1583	0.0000	1.0000
F81IG	60	-1366.8924	3100.6058	0.0000	1.0000
SYM	60	-1370.8547	3108.5303	0.0000	1.0000
JC69G	56	-1386.8377	3116.0415	0.0000	1.0000
JC69I	56	-1387.7977	3117.9616	0.0000	1.0000
JC69IG	57	-1386.0927	3120.6653	0.0000	1.0000
F81	58	-1391.4710	3137.5356	0.0000	1.0000
JC69	55	-1411.9597	3160.1718	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/12378_hg18.phylip = 1 - 452;
 Partition Dummy = 1:261-aggressive-phylip/12378_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/12378_hg18.phylip = 1 - 452;
 Partition Dummy = 1:261-aggressive-phylip/12378_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/12378_hg18.phylip = 1 - 452;
 Partition Dummy = 1:261-aggressive-phylip/12378_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

