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

Input data from file "261-aggressive-phylip/10524_hg18.phylip" (29 taxa, 432 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	-1381.6497	2903.0298	0.8877	0.8877
HKYIG	61	-1382.8212	2908.0857	0.0709	0.9586
HKYI	60	-1385.2070	2910.1444	0.0253	0.9839
GTRG	64	-1380.6966	2912.0636	0.0097	0.9936
GTRIG	65	-1380.0391	2913.5208	0.0047	0.9983
GTRI	64	-1382.4318	2915.5338	0.0017	1.0000
K2PG	57	-1400.8129	2933.3050	0.0000	1.0000
K2PIG	58	-1401.5265	2937.4016	0.0000	1.0000
SYMG	61	-1397.8938	2938.2308	0.0000	1.0000
K2PI	57	-1404.3722	2940.4236	0.0000	1.0000
SYMIG	62	-1398.3824	2941.9355	0.0000	1.0000
SYMI	61	-1400.6661	2943.7755	0.0000	1.0000
F81G	59	-1406.6761	2950.3844	0.0000	1.0000
F81IG	60	-1408.5659	2956.8622	0.0000	1.0000
HKY	59	-1410.5184	2958.0691	0.0000	1.0000
F81I	59	-1415.2860	2967.6043	0.0000	1.0000
GTR	63	-1410.1701	2968.2533	0.0000	1.0000
JC69G	56	-1424.0044	2977.0328	0.0000	1.0000
JC69IG	57	-1425.8750	2983.4292	0.0000	1.0000
K2P	56	-1427.4507	2983.9255	0.0000	1.0000
JC69I	56	-1428.7863	2986.5967	0.0000	1.0000
SYM	60	-1424.0931	2987.9168	0.0000	1.0000
F81	58	-1433.5366	3001.4217	0.0000	1.0000
JC69	55	-1449.7731	3025.9291	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1381.6497	2883.2994	0.8116	0.8116
HKYIG	61	-1382.8212	2887.6424	0.0925	0.9041
GTRG	64	-1380.6966	2889.3933	0.0386	0.9427
GTRIG	65	-1380.0391	2890.0782	0.0274	0.9701
HKYI	60	-1385.2070	2890.4140	0.0231	0.9932
GTRI	64	-1382.4318	2892.8635	0.0068	1.0000
K2PG	57	-1400.8129	2915.6259	0.0000	1.0000
SYMG	61	-1397.8938	2917.7875	0.0000	1.0000
K2PIG	58	-1401.5265	2919.0531	0.0000	1.0000
SYMIG	62	-1398.3824	2920.7647	0.0000	1.0000
K2PI	57	-1404.3722	2922.7444	0.0000	1.0000
SYMI	61	-1400.6661	2923.3323	0.0000	1.0000
F81G	59	-1406.6761	2931.3522	0.0000	1.0000
F81IG	60	-1408.5659	2937.1317	0.0000	1.0000
HKY	59	-1410.5184	2939.0368	0.0000	1.0000
GTR	63	-1410.1701	2946.3403	0.0000	1.0000
F81I	59	-1415.2860	2948.5720	0.0000	1.0000
JC69G	56	-1424.0044	2960.0088	0.0000	1.0000
JC69IG	57	-1425.8750	2965.7501	0.0000	1.0000
K2P	56	-1427.4507	2966.9015	0.0000	1.0000
SYM	60	-1424.0931	2968.1863	0.0000	1.0000
JC69I	56	-1428.7863	2969.5727	0.0000	1.0000
F81	58	-1433.5366	2983.0732	0.0000	1.0000
JC69	55	-1449.7731	3009.5461	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1381.6497	3127.4049	0.9583	0.9583
HKYI	60	-1385.2070	3134.5195	0.0273	0.9857
HKYIG	61	-1382.8212	3135.8164	0.0143	0.9999
K2PG	57	-1400.8129	3147.5261	0.0000	1.0000
GTRG	64	-1380.6966	3149.7725	0.0000	1.0000
GTRI	64	-1382.4318	3153.2428	0.0000	1.0000
GTRIG	65	-1380.0391	3154.5258	0.0000	1.0000
K2PI	57	-1404.3722	3154.6447	0.0000	1.0000
K2PIG	58	-1401.5265	3155.0217	0.0000	1.0000
SYMG	61	-1397.8938	3165.9615	0.0000	1.0000
F81G	59	-1406.6761	3171.3893	0.0000	1.0000
SYMI	61	-1400.6661	3171.5062	0.0000	1.0000
SYMIG	62	-1398.3824	3173.0071	0.0000	1.0000
HKY	59	-1410.5184	3179.0739	0.0000	1.0000
F81IG	60	-1408.5659	3181.2373	0.0000	1.0000
JC69G	56	-1424.0044	3187.8406	0.0000	1.0000
F81I	59	-1415.2860	3188.6091	0.0000	1.0000
K2P	56	-1427.4507	3194.7333	0.0000	1.0000
JC69I	56	-1428.7863	3197.4045	0.0000	1.0000
JC69IG	57	-1425.8750	3197.6503	0.0000	1.0000
GTR	63	-1410.1701	3202.6511	0.0000	1.0000
SYM	60	-1424.0931	3212.2918	0.0000	1.0000
F81	58	-1433.5366	3219.0419	0.0000	1.0000
JC69	55	-1449.7731	3233.3095	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/10524_hg18.phylip = 1 - 432;
 Partition Dummy = 1:261-aggressive-phylip/10524_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/10524_hg18.phylip = 1 - 432;
 Partition Dummy = 1:261-aggressive-phylip/10524_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/10524_hg18.phylip = 1 - 432;
 Partition Dummy = 1:261-aggressive-phylip/10524_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

