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

Input data from file "261-aggressive-phylip/11246_hg18.phylip" (29 taxa, 378 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	-1376.2991	2895.6896	0.5174	0.5174
GTRG	64	-1371.4979	2897.5772	0.2013	0.7187
K2PG	57	-1381.8199	2898.3023	0.1401	0.8588
GTRI	64	-1372.6708	2899.9231	0.0623	0.9211
GTRIG	65	-1371.4645	2900.4289	0.0484	0.9695
HKYIG	61	-1378.6139	2903.1645	0.0123	0.9819
HKYI	60	-1380.3960	2903.8834	0.0086	0.9905
K2PIG	58	-1383.6819	2904.8183	0.0054	0.9958
K2PI	57	-1385.6370	2905.9364	0.0031	0.9989
SYMIG	62	-1379.8245	2908.4489	0.0009	0.9998
SYMG	61	-1382.8173	2911.5713	0.0002	1.0000
SYMI	61	-1385.4164	2916.7696	0.0000	1.0000
HKY	59	-1399.2904	2938.8450	0.0000	1.0000
K2P	56	-1405.1788	2942.2455	0.0000	1.0000
GTR	63	-1396.0084	2943.6983	0.0000	1.0000
SYM	60	-1406.9810	2957.0535	0.0000	1.0000
F81I	59	-1413.0636	2966.3913	0.0000	1.0000
F81IG	60	-1411.8528	2966.7971	0.0000	1.0000
F81G	59	-1415.0749	2970.4139	0.0000	1.0000
JC69G	56	-1419.3546	2970.5971	0.0000	1.0000
JC69I	56	-1422.9499	2977.7878	0.0000	1.0000
JC69IG	57	-1422.4931	2979.6488	0.0000	1.0000
F81	58	-1438.9089	3015.2724	0.0000	1.0000
JC69	55	-1443.1378	3015.4061	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1371.4979	2870.9957	0.4539	0.4539
HKYG	60	-1376.2991	2872.5981	0.2037	0.6575
GTRIG	65	-1371.4645	2872.9289	0.1726	0.8302
GTRI	64	-1372.6708	2873.3416	0.1404	0.9706
K2PG	57	-1381.8199	2877.6398	0.0164	0.9870
HKYIG	61	-1378.6139	2879.2278	0.0074	0.9944
HKYI	60	-1380.3960	2880.7919	0.0034	0.9978
K2PIG	58	-1383.6819	2883.3638	0.0009	0.9987
SYMIG	62	-1379.8245	2883.6489	0.0008	0.9995
K2PI	57	-1385.6370	2885.2739	0.0004	0.9999
SYMG	61	-1382.8173	2887.6346	0.0001	1.0000
SYMI	61	-1385.4164	2892.8329	0.0000	1.0000
HKY	59	-1399.2904	2916.5808	0.0000	1.0000
GTR	63	-1396.0084	2918.0167	0.0000	1.0000
K2P	56	-1405.1788	2922.3577	0.0000	1.0000
SYM	60	-1406.9810	2933.9620	0.0000	1.0000
F81IG	60	-1411.8528	2943.7056	0.0000	1.0000
F81I	59	-1413.0636	2944.1271	0.0000	1.0000
F81G	59	-1415.0749	2948.1497	0.0000	1.0000
JC69G	56	-1419.3546	2950.7093	0.0000	1.0000
JC69I	56	-1422.9499	2957.8999	0.0000	1.0000
JC69IG	57	-1422.4931	2958.9863	0.0000	1.0000
F81	58	-1438.9089	2993.8178	0.0000	1.0000
JC69	55	-1443.1378	2996.2757	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1381.8199	3101.9287	0.9392	0.9392
HKYG	60	-1376.2991	3108.6918	0.0319	0.9711
K2PI	57	-1385.6370	3109.5629	0.0207	0.9918
K2PIG	58	-1383.6819	3111.5876	0.0075	0.9993
HKYI	60	-1380.3960	3116.8856	0.0005	0.9998
HKYIG	61	-1378.6139	3119.2563	0.0002	1.0000
GTRG	64	-1371.4979	3122.8289	0.0000	1.0000
GTRI	64	-1372.6708	3125.1748	0.0000	1.0000
SYMIG	62	-1379.8245	3127.6123	0.0000	1.0000
SYMG	61	-1382.8173	3127.6632	0.0000	1.0000
GTRIG	65	-1371.4645	3128.6971	0.0000	1.0000
SYMI	61	-1385.4164	3132.8614	0.0000	1.0000
K2P	56	-1405.1788	3142.7117	0.0000	1.0000
HKY	59	-1399.2904	3148.7396	0.0000	1.0000
GTR	63	-1396.0084	3165.9151	0.0000	1.0000
SYM	60	-1406.9810	3170.0557	0.0000	1.0000
JC69G	56	-1419.3546	3171.0633	0.0000	1.0000
F81I	59	-1413.0636	3176.2859	0.0000	1.0000
JC69I	56	-1422.9499	3178.2540	0.0000	1.0000
F81IG	60	-1411.8528	3179.7993	0.0000	1.0000
F81G	59	-1415.0749	3180.3085	0.0000	1.0000
JC69IG	57	-1422.4931	3183.2752	0.0000	1.0000
JC69	55	-1443.1378	3212.6948	0.0000	1.0000
F81	58	-1438.9089	3222.0417	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/11246_hg18.phylip = 1 - 378;
 Partition Dummy = 1:261-aggressive-phylip/11246_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/11246_hg18.phylip = 1 - 378;
 Partition Dummy = 1:261-aggressive-phylip/11246_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/11246_hg18.phylip = 1 - 378;
 Partition Dummy = 1:261-aggressive-phylip/11246_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

