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

Input data from file "591-tomodel-phylip/chrUn_random_11933.phylip" (7 taxa, 745 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	-2318.4439	4678.0481	0.6163	0.6163
GTRIG	21	-2318.4442	4680.1664	0.2137	0.8299
GTRI	20	-2319.7332	4680.6266	0.1698	0.9997
HKYG	16	-2330.8505	4694.4483	0.0002	0.9999
HKYIG	17	-2330.8470	4696.5357	0.0001	0.9999
HKYI	16	-2332.0630	4696.8732	0.0001	1.0000
GTR	19	-2329.7344	4698.5172	0.0000	1.0000
SYMG	17	-2336.4305	4707.7028	0.0000	1.0000
SYMIG	18	-2336.4308	4709.8038	0.0000	1.0000
SYMI	17	-2337.6686	4710.1791	0.0000	1.0000
K2PG	13	-2343.8545	4714.2069	0.0000	1.0000
HKY	15	-2342.0286	4714.7157	0.0000	1.0000
K2PIG	14	-2343.8523	4716.2799	0.0000	1.0000
K2PI	13	-2345.1498	4716.7975	0.0000	1.0000
SYM	16	-2347.2448	4727.2369	0.0000	1.0000
K2P	12	-2355.4103	4735.2468	0.0000	1.0000
F81G	15	-2364.8745	4760.4075	0.0000	1.0000
F81I	15	-2365.8039	4762.2662	0.0000	1.0000
F81IG	16	-2364.8745	4762.4962	0.0000	1.0000
JC69G	12	-2375.6858	4775.7979	0.0000	1.0000
JC69I	12	-2376.7053	4777.8368	0.0000	1.0000
JC69IG	13	-2375.6872	4777.8724	0.0000	1.0000
F81	14	-2374.7268	4778.0290	0.0000	1.0000
JC69	11	-2385.9563	4794.2728	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2318.4439	4676.8879	0.6084	0.6084
GTRIG	21	-2318.4442	4678.8884	0.2238	0.8322
GTRI	20	-2319.7332	4679.4664	0.1676	0.9998
HKYG	16	-2330.8505	4693.7010	0.0001	0.9999
HKYIG	17	-2330.8470	4695.6939	0.0001	0.9999
HKYI	16	-2332.0630	4696.1259	0.0000	1.0000
GTR	19	-2329.7344	4697.4689	0.0000	1.0000
SYMG	17	-2336.4305	4706.8610	0.0000	1.0000
SYMIG	18	-2336.4308	4708.8616	0.0000	1.0000
SYMI	17	-2337.6686	4709.3373	0.0000	1.0000
K2PG	13	-2343.8545	4713.7090	0.0000	1.0000
HKY	15	-2342.0286	4714.0573	0.0000	1.0000
K2PIG	14	-2343.8523	4715.7046	0.0000	1.0000
K2PI	13	-2345.1498	4716.2996	0.0000	1.0000
SYM	16	-2347.2448	4726.4897	0.0000	1.0000
K2P	12	-2355.4103	4734.8206	0.0000	1.0000
F81G	15	-2364.8745	4759.7491	0.0000	1.0000
F81I	15	-2365.8039	4761.6078	0.0000	1.0000
F81IG	16	-2364.8745	4761.7490	0.0000	1.0000
JC69G	12	-2375.6858	4775.3716	0.0000	1.0000
JC69IG	13	-2375.6872	4777.3744	0.0000	1.0000
JC69I	12	-2376.7053	4777.4106	0.0000	1.0000
F81	14	-2374.7268	4777.4537	0.0000	1.0000
JC69	11	-2385.9563	4793.9127	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2330.8505	4767.5152	0.5069	0.5069
GTRG	20	-2318.4439	4769.1556	0.2232	0.7301
HKYI	16	-2332.0630	4769.9401	0.1508	0.8809
GTRI	20	-2319.7332	4771.7341	0.0615	0.9424
K2PG	13	-2343.8545	4773.6830	0.0232	0.9656
HKYIG	17	-2330.8470	4774.1215	0.0186	0.9842
GTRIG	21	-2318.4442	4775.7695	0.0082	0.9924
K2PI	13	-2345.1498	4776.2736	0.0064	0.9988
K2PIG	14	-2343.8523	4780.2919	0.0009	0.9996
HKY	15	-2342.0286	4783.2580	0.0002	0.9998
GTR	19	-2329.7344	4785.1232	0.0001	0.9999
SYMG	17	-2336.4305	4785.2886	0.0001	1.0000
SYMI	17	-2337.6686	4787.7648	0.0000	1.0000
K2P	12	-2355.4103	4790.1812	0.0000	1.0000
SYMIG	18	-2336.4308	4791.9025	0.0000	1.0000
SYM	16	-2347.2448	4800.3038	0.0000	1.0000
F81G	15	-2364.8745	4828.9498	0.0000	1.0000
JC69G	12	-2375.6858	4830.7323	0.0000	1.0000
F81I	15	-2365.8039	4830.8086	0.0000	1.0000
JC69I	12	-2376.7053	4832.7712	0.0000	1.0000
F81IG	16	-2364.8745	4835.5631	0.0000	1.0000
JC69IG	13	-2375.6872	4837.3484	0.0000	1.0000
F81	14	-2374.7268	4842.0411	0.0000	1.0000
JC69	11	-2385.9563	4844.6599	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/chrUn_random_11933.phylip = 1 - 745;
 Partition Dummy = 1:591-tomodel-phylip/chrUn_random_11933.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 591-tomodel-phylip/chrUn_random_11933.phylip = 1 - 745;
 Partition Dummy = 1:591-tomodel-phylip/chrUn_random_11933.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 591-tomodel-phylip/chrUn_random_11933.phylip = 1 - 745;
 Partition Dummy = 1:591-tomodel-phylip/chrUn_random_11933.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

