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

Input data from file "261-aggressive-phylip/12566_hg18.phylip" (29 taxa, 370 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	-1266.7062	2677.1018	0.4847	0.4847
HKYI	60	-1267.4425	2678.5743	0.2321	0.7169
HKYIG	61	-1266.6955	2679.9495	0.1167	0.8336
GTRG	64	-1262.5887	2680.4561	0.0906	0.9242
GTRI	64	-1263.0874	2681.4534	0.0550	0.9792
GTRIG	65	-1262.5883	2683.4003	0.0208	1.0000
F81G	59	-1282.1495	2705.1378	0.0000	1.0000
F81I	59	-1282.9537	2706.7461	0.0000	1.0000
F81IG	60	-1282.1256	2707.9406	0.0000	1.0000
K2PG	57	-1287.7274	2710.6471	0.0000	1.0000
K2PI	57	-1288.6735	2712.5394	0.0000	1.0000
K2PIG	58	-1287.6583	2713.3230	0.0000	1.0000
SYMG	61	-1283.7809	2714.1202	0.0000	1.0000
SYMI	61	-1284.4814	2715.5212	0.0000	1.0000
SYMIG	62	-1283.7587	2716.9636	0.0000	1.0000
HKY	59	-1294.4478	2729.7344	0.0000	1.0000
JC69G	56	-1300.1084	2732.6129	0.0000	1.0000
JC69I	56	-1301.0818	2734.5597	0.0000	1.0000
JC69IG	57	-1300.0419	2735.2762	0.0000	1.0000
GTR	63	-1291.6504	2735.6537	0.0000	1.0000
F81	58	-1308.0692	2754.1448	0.0000	1.0000
K2P	56	-1314.0369	2760.4700	0.0000	1.0000
SYM	60	-1310.2278	2764.1450	0.0000	1.0000
JC69	55	-1325.6917	2781.0012	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1262.5887	2653.1775	0.2762	0.2762
HKYG	60	-1266.7062	2653.4125	0.2456	0.5217
GTRI	64	-1263.0874	2654.1748	0.1677	0.6895
HKYI	60	-1267.4425	2654.8850	0.1176	0.8071
GTRIG	65	-1262.5883	2655.1766	0.1016	0.9087
HKYIG	61	-1266.6955	2655.3910	0.0913	1.0000
F81G	59	-1282.1495	2682.2991	0.0000	1.0000
F81I	59	-1282.9537	2683.9074	0.0000	1.0000
F81IG	60	-1282.1256	2684.2512	0.0000	1.0000
K2PG	57	-1287.7274	2689.4548	0.0000	1.0000
SYMG	61	-1283.7809	2689.5617	0.0000	1.0000
SYMI	61	-1284.4814	2690.9627	0.0000	1.0000
K2PIG	58	-1287.6583	2691.3166	0.0000	1.0000
K2PI	57	-1288.6735	2691.3471	0.0000	1.0000
SYMIG	62	-1283.7587	2691.5174	0.0000	1.0000
HKY	59	-1294.4478	2706.8957	0.0000	1.0000
GTR	63	-1291.6504	2709.3007	0.0000	1.0000
JC69G	56	-1300.1084	2712.2168	0.0000	1.0000
JC69IG	57	-1300.0419	2714.0839	0.0000	1.0000
JC69I	56	-1301.0818	2714.1635	0.0000	1.0000
F81	58	-1308.0692	2732.1383	0.0000	1.0000
K2P	56	-1314.0369	2740.0739	0.0000	1.0000
SYM	60	-1310.2278	2740.4557	0.0000	1.0000
JC69	55	-1325.6917	2761.3834	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1266.7062	2888.2226	0.6527	0.6527
HKYI	60	-1267.4425	2889.6952	0.3126	0.9652
HKYIG	61	-1266.6955	2894.1147	0.0343	0.9995
GTRG	64	-1262.5887	2903.6417	0.0003	0.9998
GTRI	64	-1263.0874	2904.6390	0.0002	1.0000
GTRIG	65	-1262.5883	2909.5543	0.0000	1.0000
K2PG	57	-1287.7274	2912.5245	0.0000	1.0000
F81G	59	-1282.1495	2913.1958	0.0000	1.0000
K2PI	57	-1288.6735	2914.4167	0.0000	1.0000
F81I	59	-1282.9537	2914.8041	0.0000	1.0000
K2PIG	58	-1287.6583	2918.2997	0.0000	1.0000
F81IG	60	-1282.1256	2919.0614	0.0000	1.0000
SYMG	61	-1283.7809	2928.2854	0.0000	1.0000
SYMI	61	-1284.4814	2929.6864	0.0000	1.0000
JC69G	56	-1300.1084	2931.3729	0.0000	1.0000
JC69I	56	-1301.0818	2933.3197	0.0000	1.0000
SYMIG	62	-1283.7587	2934.1545	0.0000	1.0000
JC69IG	57	-1300.0419	2937.1536	0.0000	1.0000
HKY	59	-1294.4478	2937.7924	0.0000	1.0000
GTR	63	-1291.6504	2955.8514	0.0000	1.0000
F81	58	-1308.0692	2959.1215	0.0000	1.0000
K2P	56	-1314.0369	2959.2300	0.0000	1.0000
SYM	60	-1310.2278	2975.2658	0.0000	1.0000
JC69	55	-1325.6917	2976.6260	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/12566_hg18.phylip = 1 - 370;
 Partition Dummy = 1:261-aggressive-phylip/12566_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/12566_hg18.phylip = 1 - 370;
 Partition Dummy = 1:261-aggressive-phylip/12566_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/12566_hg18.phylip = 1 - 370;
 Partition Dummy = 1:261-aggressive-phylip/12566_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

