
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:33 04/04/2011

Input data from file "591-tomodel-phylip/chr2_18714.phylip" (7 taxa, 564 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	16	-2369.3691	4771.7328	0.3172	0.3172
GTRG	20	-2365.3088	4772.1645	0.2556	0.5729
HKYI	16	-2370.4104	4773.8152	0.1120	0.6849
HKYIG	17	-2369.3576	4773.8361	0.1108	0.7957
GTRIG	21	-2365.3009	4774.3067	0.0876	0.8833
GTRI	20	-2366.6371	4774.8211	0.0677	0.9510
K2PG	13	-2375.0953	4776.8523	0.0245	0.9756
K2PIG	14	-2375.0997	4778.9644	0.0085	0.9841
K2PI	13	-2376.2183	4779.0983	0.0080	0.9921
SYMG	17	-2372.4738	4780.0685	0.0049	0.9970
SYMIG	18	-2372.4884	4782.2318	0.0017	0.9987
SYMI	17	-2373.7759	4782.6726	0.0013	1.0000
HKY	15	-2389.5702	4810.0164	0.0000	1.0000
GTR	19	-2386.4221	4812.2413	0.0000	1.0000
K2P	12	-2394.2937	4813.1536	0.0000	1.0000
SYM	16	-2391.6962	4816.3870	0.0000	1.0000
F81G	15	-2437.8365	4906.5489	0.0000	1.0000
JC69G	12	-2441.2182	4907.0026	0.0000	1.0000
F81I	15	-2438.7167	4908.3093	0.0000	1.0000
F81IG	16	-2437.8477	4908.6899	0.0000	1.0000
JC69I	12	-2442.1230	4908.8122	0.0000	1.0000
JC69IG	13	-2441.2309	4909.1237	0.0000	1.0000
JC69	11	-2457.8512	4938.1807	0.0000	1.0000
F81	14	-2454.9265	4938.6180	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2365.3088	4770.6176	0.2950	0.2950
HKYG	16	-2369.3691	4770.7383	0.2778	0.5728
GTRIG	21	-2365.3009	4772.6019	0.1094	0.6822
HKYIG	17	-2369.3576	4772.7153	0.1034	0.7855
HKYI	16	-2370.4104	4772.8207	0.0981	0.8836
GTRI	20	-2366.6371	4773.2741	0.0782	0.9618
K2PG	13	-2375.0953	4776.1905	0.0182	0.9799
K2PIG	14	-2375.0997	4778.1994	0.0067	0.9866
K2PI	13	-2376.2183	4778.4365	0.0059	0.9925
SYMG	17	-2372.4738	4778.9477	0.0046	0.9971
SYMIG	18	-2372.4884	4780.9768	0.0017	0.9988
SYMI	17	-2373.7759	4781.5517	0.0012	1.0000
HKY	15	-2389.5702	4809.1405	0.0000	1.0000
GTR	19	-2386.4221	4810.8443	0.0000	1.0000
K2P	12	-2394.2937	4812.5874	0.0000	1.0000
SYM	16	-2391.6962	4815.3925	0.0000	1.0000
F81G	15	-2437.8365	4905.6730	0.0000	1.0000
JC69G	12	-2441.2182	4906.4363	0.0000	1.0000
F81I	15	-2438.7167	4907.4334	0.0000	1.0000
F81IG	16	-2437.8477	4907.6954	0.0000	1.0000
JC69I	12	-2442.1230	4908.2459	0.0000	1.0000
JC69IG	13	-2441.2309	4908.4619	0.0000	1.0000
JC69	11	-2457.8512	4937.7025	0.0000	1.0000
F81	14	-2454.9265	4937.8530	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2375.0953	4832.5462	0.7147	0.7147
K2PI	13	-2376.2183	4834.7922	0.2325	0.9472
K2PIG	14	-2375.0997	4838.8901	0.0300	0.9771
HKYG	16	-2369.3691	4840.0991	0.0164	0.9935
HKYI	16	-2370.4104	4842.1816	0.0058	0.9993
HKYIG	17	-2369.3576	4846.4112	0.0007	1.0000
SYMG	17	-2372.4738	4852.6436	0.0000	1.0000
SYMI	17	-2373.7759	4855.2477	0.0000	1.0000
GTRG	20	-2365.3088	4857.3187	0.0000	1.0000
SYMIG	18	-2372.4884	4859.0078	0.0000	1.0000
GTRI	20	-2366.6371	4859.9752	0.0000	1.0000
GTRIG	21	-2365.3009	4863.6380	0.0000	1.0000
K2P	12	-2394.2937	4864.6081	0.0000	1.0000
HKY	15	-2389.5702	4874.1663	0.0000	1.0000
SYM	16	-2391.6962	4884.7533	0.0000	1.0000
GTR	19	-2386.4221	4893.2103	0.0000	1.0000
JC69G	12	-2441.2182	4958.4570	0.0000	1.0000
JC69I	12	-2442.1230	4960.2666	0.0000	1.0000
JC69IG	13	-2441.2309	4964.8176	0.0000	1.0000
F81G	15	-2437.8365	4970.6988	0.0000	1.0000
F81I	15	-2438.7167	4972.4592	0.0000	1.0000
F81IG	16	-2437.8477	4977.0562	0.0000	1.0000
JC69	11	-2457.8512	4985.3881	0.0000	1.0000
F81	14	-2454.9265	4998.5437	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 591-tomodel-phylip/chr2_18714.phylip = 1 - 564;
 Partition Dummy = 1:591-tomodel-phylip/chr2_18714.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 591-tomodel-phylip/chr2_18714.phylip = 1 - 564;
 Partition Dummy = 1:591-tomodel-phylip/chr2_18714.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 591-tomodel-phylip/chr2_18714.phylip = 1 - 564;
 Partition Dummy = 1:591-tomodel-phylip/chr2_18714.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

