
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
01:01 02/17/2011

Input data from file "183-phylip/chr8_4221.phylip" (29 taxa, 445 characters)

Minimum AIC  model: HKYIG
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	-1403.5578	2946.1781	0.5013	0.5013
HKYIG	61	-1402.3671	2946.4835	0.4303	0.9316
HKYI	60	-1406.4693	2952.0011	0.0273	0.9588
GTRG	64	-1401.2666	2952.4279	0.0220	0.9808
GTRIG	65	-1400.1006	2952.8397	0.0179	0.9988
GTRI	64	-1404.1538	2958.2024	0.0012	1.0000
K2PG	57	-1424.0986	2979.2825	0.0000	1.0000
K2PIG	58	-1423.1372	2980.0050	0.0000	1.0000
K2PI	57	-1426.2533	2983.5919	0.0000	1.0000
SYMG	61	-1421.7676	2985.2844	0.0000	1.0000
SYMIG	62	-1420.8061	2986.0624	0.0000	1.0000
SYMI	61	-1423.9038	2989.5570	0.0000	1.0000
HKY	59	-1434.5721	3005.5338	0.0000	1.0000
F81G	59	-1437.1471	3010.6837	0.0000	1.0000
F81IG	60	-1436.2837	3011.6298	0.0000	1.0000
GTR	63	-1432.6898	3012.5450	0.0000	1.0000
F81I	59	-1439.3101	3015.0098	0.0000	1.0000
K2P	56	-1454.2909	3037.0353	0.0000	1.0000
JC69G	56	-1456.3891	3041.2317	0.0000	1.0000
JC69IG	57	-1455.7206	3042.5266	0.0000	1.0000
SYM	60	-1452.1763	3043.4152	0.0000	1.0000
JC69I	56	-1458.6158	3045.6851	0.0000	1.0000
F81	58	-1466.7500	3067.2305	0.0000	1.0000
JC69	55	-1485.7018	3097.2391	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-1402.3671	2926.7341	0.4533	0.4533
HKYG	60	-1403.5578	2927.1156	0.3746	0.8279
GTRIG	65	-1400.1006	2930.2012	0.0801	0.9080
GTRG	64	-1401.2666	2930.5331	0.0678	0.9758
HKYI	60	-1406.4693	2932.9386	0.0204	0.9962
GTRI	64	-1404.1538	2936.3076	0.0038	1.0000
K2PG	57	-1424.0986	2962.1972	0.0000	1.0000
K2PIG	58	-1423.1372	2962.2745	0.0000	1.0000
SYMG	61	-1421.7676	2965.5351	0.0000	1.0000
SYMIG	62	-1420.8061	2965.6121	0.0000	1.0000
K2PI	57	-1426.2533	2966.5066	0.0000	1.0000
SYMI	61	-1423.9038	2969.8077	0.0000	1.0000
HKY	59	-1434.5721	2987.1442	0.0000	1.0000
GTR	63	-1432.6898	2991.3796	0.0000	1.0000
F81G	59	-1437.1471	2992.2941	0.0000	1.0000
F81IG	60	-1436.2837	2992.5673	0.0000	1.0000
F81I	59	-1439.3101	2996.6202	0.0000	1.0000
K2P	56	-1454.2909	3020.5817	0.0000	1.0000
SYM	60	-1452.1763	3024.3527	0.0000	1.0000
JC69G	56	-1456.3891	3024.7781	0.0000	1.0000
JC69IG	57	-1455.7206	3025.4413	0.0000	1.0000
JC69I	56	-1458.6158	3029.2315	0.0000	1.0000
F81	58	-1466.7500	3049.5000	0.0000	1.0000
JC69	55	-1485.7018	3081.4036	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1403.5578	3173.0000	0.8262	0.8262
HKYIG	61	-1402.3671	3176.7167	0.1288	0.9550
HKYI	60	-1406.4693	3178.8231	0.0449	0.9999
GTRG	64	-1401.2666	3192.8099	0.0000	1.0000
K2PG	57	-1424.0986	3195.7875	0.0000	1.0000
GTRIG	65	-1400.1006	3196.5760	0.0000	1.0000
GTRI	64	-1404.1538	3198.5844	0.0000	1.0000
K2PIG	58	-1423.1372	3199.9628	0.0000	1.0000
K2PI	57	-1426.2533	3200.0968	0.0000	1.0000
SYMG	61	-1421.7676	3215.5176	0.0000	1.0000
SYMIG	62	-1420.8061	3219.6927	0.0000	1.0000
SYMI	61	-1423.9038	3219.7902	0.0000	1.0000
HKY	59	-1434.5721	3228.9306	0.0000	1.0000
F81G	59	-1437.1471	3234.0805	0.0000	1.0000
F81I	59	-1439.3101	3238.4066	0.0000	1.0000
F81IG	60	-1436.2837	3238.4518	0.0000	1.0000
GTR	63	-1432.6898	3249.5583	0.0000	1.0000
K2P	56	-1454.2909	3250.0739	0.0000	1.0000
JC69G	56	-1456.3891	3254.2703	0.0000	1.0000
JC69I	56	-1458.6158	3258.7237	0.0000	1.0000
JC69IG	57	-1455.7206	3259.0315	0.0000	1.0000
SYM	60	-1452.1763	3270.2371	0.0000	1.0000
F81	58	-1466.7500	3287.1883	0.0000	1.0000
JC69	55	-1485.7018	3306.7977	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYIG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr8_4221.phylip = 1 - 445;
 Partition Dummy = 1:183-phylip/chr8_4221.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr8_4221.phylip = 1 - 445;
 Partition Dummy = 1:183-phylip/chr8_4221.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 183-phylip/chr8_4221.phylip = 1 - 445;
 Partition Dummy = 1:183-phylip/chr8_4221.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

