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

Input data from file "591-tomodel-phylip/chr7_4777.phylip" (7 taxa, 755 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	16	-2131.7750	4296.2870	0.5488	0.5488
HKYI	16	-2132.6469	4298.0310	0.2295	0.7783
HKYIG	17	-2131.7779	4298.3862	0.1921	0.9705
GTRG	20	-2131.0547	4303.2538	0.0169	0.9873
GTRI	20	-2131.9621	4305.0687	0.0068	0.9941
GTRIG	21	-2131.0559	4305.3723	0.0058	1.0000
HKY	15	-2142.1799	4315.0094	0.0000	1.0000
GTR	19	-2141.7252	4322.4845	0.0000	1.0000
K2PG	13	-2157.7594	4342.0100	0.0000	1.0000
K2PI	13	-2158.5996	4343.6904	0.0000	1.0000
K2PIG	14	-2157.7615	4344.0906	0.0000	1.0000
SYMG	17	-2157.1818	4349.1941	0.0000	1.0000
SYMI	17	-2157.9883	4350.8070	0.0000	1.0000
SYMIG	18	-2157.1833	4351.2959	0.0000	1.0000
K2P	12	-2167.5462	4359.5128	0.0000	1.0000
SYM	16	-2167.0990	4366.9351	0.0000	1.0000
F81G	15	-2175.1746	4380.9987	0.0000	1.0000
F81I	15	-2176.0594	4382.7684	0.0000	1.0000
F81IG	16	-2175.1752	4383.0875	0.0000	1.0000
F81	14	-2184.6649	4397.8974	0.0000	1.0000
JC69G	12	-2197.2083	4418.8372	0.0000	1.0000
JC69I	12	-2198.0445	4420.5095	0.0000	1.0000
JC69IG	13	-2197.2083	4420.9079	0.0000	1.0000
JC69	11	-2206.4794	4435.3142	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2131.7750	4295.5499	0.5401	0.5401
HKYI	16	-2132.6469	4297.2939	0.2258	0.7659
HKYIG	17	-2131.7779	4297.5558	0.1981	0.9640
GTRG	20	-2131.0547	4302.1094	0.0203	0.9843
GTRI	20	-2131.9621	4303.9243	0.0082	0.9925
GTRIG	21	-2131.0559	4304.1118	0.0075	1.0000
HKY	15	-2142.1799	4314.3598	0.0000	1.0000
GTR	19	-2141.7252	4321.4505	0.0000	1.0000
K2PG	13	-2157.7594	4341.5188	0.0000	1.0000
K2PI	13	-2158.5996	4343.1992	0.0000	1.0000
K2PIG	14	-2157.7615	4343.5231	0.0000	1.0000
SYMG	17	-2157.1818	4348.3637	0.0000	1.0000
SYMI	17	-2157.9883	4349.9766	0.0000	1.0000
SYMIG	18	-2157.1833	4350.3665	0.0000	1.0000
K2P	12	-2167.5462	4359.0923	0.0000	1.0000
SYM	16	-2167.0990	4366.1980	0.0000	1.0000
F81G	15	-2175.1746	4380.3492	0.0000	1.0000
F81I	15	-2176.0594	4382.1189	0.0000	1.0000
F81IG	16	-2175.1752	4382.3504	0.0000	1.0000
F81	14	-2184.6649	4397.3298	0.0000	1.0000
JC69G	12	-2197.2083	4418.4167	0.0000	1.0000
JC69I	12	-2198.0445	4420.0890	0.0000	1.0000
JC69IG	13	-2197.2083	4420.4167	0.0000	1.0000
JC69	11	-2206.4794	4434.9589	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2131.7750	4369.5774	0.6872	0.6872
HKYI	16	-2132.6469	4371.3214	0.2873	0.9745
HKYIG	17	-2131.7779	4376.2100	0.0249	0.9994
HKY	15	-2142.1799	4383.7606	0.0006	1.0000
GTRG	20	-2131.0547	4394.6437	0.0000	1.0000
GTRI	20	-2131.9621	4396.4586	0.0000	1.0000
GTRIG	21	-2131.0559	4401.2728	0.0000	1.0000
K2PG	13	-2157.7594	4401.6661	0.0000	1.0000
K2PI	13	-2158.5996	4403.3465	0.0000	1.0000
K2PIG	14	-2157.7615	4408.2971	0.0000	1.0000
GTR	19	-2141.7252	4409.3581	0.0000	1.0000
K2P	12	-2167.5462	4414.6130	0.0000	1.0000
SYMG	17	-2157.1818	4427.0179	0.0000	1.0000
SYMI	17	-2157.9883	4428.6308	0.0000	1.0000
SYMIG	18	-2157.1833	4433.6475	0.0000	1.0000
SYM	16	-2167.0990	4440.2255	0.0000	1.0000
F81G	15	-2175.1746	4449.7499	0.0000	1.0000
F81I	15	-2176.0594	4451.5196	0.0000	1.0000
F81IG	16	-2175.1752	4456.3779	0.0000	1.0000
F81	14	-2184.6649	4462.1039	0.0000	1.0000
JC69G	12	-2197.2083	4473.9373	0.0000	1.0000
JC69I	12	-2198.0445	4475.6096	0.0000	1.0000
JC69IG	13	-2197.2083	4480.5640	0.0000	1.0000
JC69	11	-2206.4794	4485.8528	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 591-tomodel-phylip/chr7_4777.phylip = 1 - 755;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4777.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 591-tomodel-phylip/chr7_4777.phylip = 1 - 755;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4777.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 591-tomodel-phylip/chr7_4777.phylip = 1 - 755;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4777.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

