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

Input data from file "591-tomodel-phylip/chr2_13502.phylip" (7 taxa, 524 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	-1840.0262	3721.7225	0.4041	0.4041
GTRI	20	-1840.1487	3721.9674	0.3575	0.7615
GTRIG	21	-1839.7894	3723.4194	0.1730	0.9345
HKYG	16	-1846.8154	3726.7038	0.0335	0.9680
HKYI	16	-1847.3477	3727.7685	0.0197	0.9876
HKYIG	17	-1846.7444	3728.6982	0.0124	1.0000
GTR	19	-1861.0997	3761.7074	0.0000	1.0000
HKY	15	-1868.2595	3767.4639	0.0000	1.0000
K2PG	13	-1872.8501	3772.4138	0.0000	1.0000
K2PI	13	-1873.5526	3773.8189	0.0000	1.0000
K2PIG	14	-1872.8482	3774.5216	0.0000	1.0000
SYMG	17	-1871.8082	3778.8258	0.0000	1.0000
SYMI	17	-1872.4746	3780.1586	0.0000	1.0000
SYMIG	18	-1871.8071	3780.9687	0.0000	1.0000
F81G	15	-1887.9492	3806.8432	0.0000	1.0000
F81I	15	-1888.0035	3806.9518	0.0000	1.0000
F81IG	16	-1887.8400	3808.7529	0.0000	1.0000
K2P	12	-1896.3520	3817.3145	0.0000	1.0000
SYM	16	-1894.4145	3821.9021	0.0000	1.0000
JC69G	12	-1909.1261	3842.8627	0.0000	1.0000
F81	14	-1907.1943	3843.2137	0.0000	1.0000
JC69I	12	-1909.4574	3843.5253	0.0000	1.0000
JC69IG	13	-1909.1499	3845.0134	0.0000	1.0000
JC69	11	-1929.6950	3881.9055	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1840.0262	3720.0525	0.4044	0.4044
GTRI	20	-1840.1487	3720.2974	0.3578	0.7622
GTRIG	21	-1839.7894	3721.5788	0.1885	0.9507
HKYG	16	-1846.8154	3725.6308	0.0249	0.9756
HKYI	16	-1847.3477	3726.6955	0.0146	0.9902
HKYIG	17	-1846.7444	3727.4888	0.0098	1.0000
GTR	19	-1861.0997	3760.1994	0.0000	1.0000
HKY	15	-1868.2595	3766.5190	0.0000	1.0000
K2PG	13	-1872.8501	3771.7001	0.0000	1.0000
K2PI	13	-1873.5526	3773.1052	0.0000	1.0000
K2PIG	14	-1872.8482	3773.6964	0.0000	1.0000
SYMG	17	-1871.8082	3777.6163	0.0000	1.0000
SYMI	17	-1872.4746	3778.9491	0.0000	1.0000
SYMIG	18	-1871.8071	3779.6142	0.0000	1.0000
F81G	15	-1887.9492	3805.8983	0.0000	1.0000
F81I	15	-1888.0035	3806.0069	0.0000	1.0000
F81IG	16	-1887.8400	3807.6799	0.0000	1.0000
K2P	12	-1896.3520	3816.7039	0.0000	1.0000
SYM	16	-1894.4145	3820.8291	0.0000	1.0000
JC69G	12	-1909.1261	3842.2522	0.0000	1.0000
F81	14	-1907.1943	3842.3886	0.0000	1.0000
JC69I	12	-1909.4574	3842.9147	0.0000	1.0000
JC69IG	13	-1909.1499	3844.2997	0.0000	1.0000
JC69	11	-1929.6950	3881.3899	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1846.8154	3793.8147	0.6096	0.6096
HKYI	16	-1847.3477	3794.8793	0.3580	0.9676
HKYIG	17	-1846.7444	3799.9341	0.0286	0.9962
GTRG	20	-1840.0262	3805.2823	0.0020	0.9981
GTRI	20	-1840.1487	3805.5273	0.0017	0.9999
GTRIG	21	-1839.7894	3811.0701	0.0001	1.0000
K2PG	13	-1872.8501	3827.0995	0.0000	1.0000
K2PI	13	-1873.5526	3828.5046	0.0000	1.0000
HKY	15	-1868.2595	3830.4414	0.0000	1.0000
K2PIG	14	-1872.8482	3833.3573	0.0000	1.0000
GTR	19	-1861.0997	3841.1678	0.0000	1.0000
SYMG	17	-1871.8082	3850.0617	0.0000	1.0000
SYMI	17	-1872.4746	3851.3945	0.0000	1.0000
SYMIG	18	-1871.8071	3856.3211	0.0000	1.0000
K2P	12	-1896.3520	3867.8418	0.0000	1.0000
F81G	15	-1887.9492	3869.8207	0.0000	1.0000
F81I	15	-1888.0035	3869.9293	0.0000	1.0000
F81IG	16	-1887.8400	3875.8638	0.0000	1.0000
SYM	16	-1894.4145	3889.0129	0.0000	1.0000
JC69G	12	-1909.1261	3893.3901	0.0000	1.0000
JC69I	12	-1909.4574	3894.0526	0.0000	1.0000
JC69IG	13	-1909.1499	3899.6991	0.0000	1.0000
F81	14	-1907.1943	3902.0495	0.0000	1.0000
JC69	11	-1929.6950	3928.2663	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_13502.phylip = 1 - 524;
 Partition Dummy = 1:591-tomodel-phylip/chr2_13502.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/chr2_13502.phylip = 1 - 524;
 Partition Dummy = 1:591-tomodel-phylip/chr2_13502.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/chr2_13502.phylip = 1 - 524;
 Partition Dummy = 1:591-tomodel-phylip/chr2_13502.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

