
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:14 04/04/2011

Input data from file "683-tomodel-phylip/chr2_9136.phylip" (7 taxa, 781 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-2365.2012	4771.5077	0.4631	0.4631
GTRI	20	-2365.6099	4772.3250	0.3077	0.7708
GTRIG	21	-2365.2075	4773.6324	0.1601	0.9309
HKYG	16	-2372.3257	4777.3634	0.0248	0.9557
HKYI	16	-2372.4785	4777.6691	0.0213	0.9769
GTR	19	-2369.8376	4778.6739	0.0129	0.9898
HKYIG	17	-2372.3246	4779.4513	0.0087	0.9985
HKY	15	-2376.1819	4782.9913	0.0015	1.0000
K2PI	13	-2391.5313	4809.5372	0.0000	1.0000
K2PG	13	-2391.5814	4809.6374	0.0000	1.0000
SYMG	17	-2388.0696	4810.9412	0.0000	1.0000
SYMI	17	-2388.1861	4811.1744	0.0000	1.0000
K2PIG	14	-2391.5168	4811.5819	0.0000	1.0000
SYMIG	18	-2388.0549	4813.0074	0.0000	1.0000
K2P	12	-2395.3699	4815.1460	0.0000	1.0000
SYM	16	-2392.1068	4816.9257	0.0000	1.0000
F81I	15	-2442.0367	4914.7009	0.0000	1.0000
F81G	15	-2442.1180	4914.8635	0.0000	1.0000
F81IG	16	-2442.0773	4916.8666	0.0000	1.0000
F81	14	-2445.2253	4918.9990	0.0000	1.0000
JC69I	12	-2457.8256	4940.0574	0.0000	1.0000
JC69G	12	-2457.9942	4940.3946	0.0000	1.0000
JC69IG	13	-2457.8287	4942.1320	0.0000	1.0000
JC69	11	-2461.0563	4944.4559	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2365.2012	4770.4024	0.4636	0.4636
GTRI	20	-2365.6099	4771.2198	0.3081	0.7717
GTRIG	21	-2365.2075	4772.4150	0.1695	0.9412
HKYG	16	-2372.3257	4776.6514	0.0204	0.9616
HKYI	16	-2372.4785	4776.9571	0.0175	0.9791
GTR	19	-2369.8376	4777.6752	0.0122	0.9913
HKYIG	17	-2372.3246	4778.6492	0.0075	0.9988
HKY	15	-2376.1819	4782.3639	0.0012	1.0000
K2PI	13	-2391.5313	4809.0626	0.0000	1.0000
K2PG	13	-2391.5814	4809.1628	0.0000	1.0000
SYMG	17	-2388.0696	4810.1391	0.0000	1.0000
SYMI	17	-2388.1861	4810.3723	0.0000	1.0000
K2PIG	14	-2391.5168	4811.0336	0.0000	1.0000
SYMIG	18	-2388.0549	4812.1098	0.0000	1.0000
K2P	12	-2395.3699	4814.7397	0.0000	1.0000
SYM	16	-2392.1068	4816.2136	0.0000	1.0000
F81I	15	-2442.0367	4914.0734	0.0000	1.0000
F81G	15	-2442.1180	4914.2361	0.0000	1.0000
F81IG	16	-2442.0773	4916.1545	0.0000	1.0000
F81	14	-2445.2253	4918.4507	0.0000	1.0000
JC69I	12	-2457.8256	4939.6512	0.0000	1.0000
JC69G	12	-2457.9942	4939.9884	0.0000	1.0000
JC69IG	13	-2457.8287	4941.6574	0.0000	1.0000
JC69	11	-2461.0563	4944.1126	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2372.3257	4851.2206	0.4017	0.4017
HKYI	16	-2372.4785	4851.5263	0.3448	0.7465
HKY	15	-2376.1819	4852.2725	0.2374	0.9839
HKYIG	17	-2372.3246	4857.8789	0.0144	0.9983
GTRG	20	-2365.2012	4863.6139	0.0008	0.9991
GTRI	20	-2365.6099	4864.4313	0.0005	0.9996
GTR	19	-2369.8376	4866.2261	0.0002	0.9999
K2PI	13	-2391.5313	4869.6501	0.0000	0.9999
K2PG	13	-2391.5814	4869.7503	0.0000	0.9999
GTRIG	21	-2365.2075	4870.2871	0.0000	1.0000
K2P	12	-2395.3699	4870.6666	0.0000	1.0000
K2PIG	14	-2391.5168	4876.2817	0.0000	1.0000
SYMG	17	-2388.0696	4889.3689	0.0000	1.0000
SYMI	17	-2388.1861	4889.6020	0.0000	1.0000
SYM	16	-2392.1068	4890.7828	0.0000	1.0000
SYMIG	18	-2388.0549	4896.0002	0.0000	1.0000
F81	14	-2445.2253	4983.6987	0.0000	1.0000
F81I	15	-2442.0367	4983.9820	0.0000	1.0000
F81G	15	-2442.1180	4984.1447	0.0000	1.0000
F81IG	16	-2442.0773	4990.7237	0.0000	1.0000
JC69	11	-2461.0563	4995.3789	0.0000	1.0000
JC69I	12	-2457.8256	4995.5781	0.0000	1.0000
JC69G	12	-2457.9942	4995.9153	0.0000	1.0000
JC69IG	13	-2457.8287	5002.2449	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 683-tomodel-phylip/chr2_9136.phylip = 1 - 781;
 Partition Dummy = 1:683-tomodel-phylip/chr2_9136.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 (GTRG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_9136.phylip = 1 - 781;
 Partition Dummy = 1:683-tomodel-phylip/chr2_9136.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 BIC model (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_9136.phylip = 1 - 781;
 Partition Dummy = 1:683-tomodel-phylip/chr2_9136.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

