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

Input data from file "591-tomodel-phylip/chr7_5228.phylip" (7 taxa, 713 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-1794.3072	3621.3959	0.2817	0.2817
HKYG	16	-1794.5537	3621.8890	0.2202	0.5019
GTRI	20	-1790.3912	3621.9964	0.2087	0.7106
GTRG	20	-1790.6707	3622.5553	0.1578	0.8684
HKYIG	17	-1794.5538	3623.9881	0.0771	0.9454
GTRIG	21	-1790.6708	3624.6789	0.0546	1.0000
HKY	15	-1820.7058	3672.1002	0.0000	1.0000
GTR	19	-1816.5322	3672.1610	0.0000	1.0000
K2PI	13	-1823.1943	3672.9094	0.0000	1.0000
K2PG	13	-1823.3670	3673.2547	0.0000	1.0000
K2PIG	14	-1823.3670	3675.3358	0.0000	1.0000
SYMI	17	-1821.1145	3677.1096	0.0000	1.0000
SYMG	17	-1821.3189	3677.5184	0.0000	1.0000
SYMIG	18	-1821.3190	3679.6236	0.0000	1.0000
F81I	15	-1834.3111	3699.3109	0.0000	1.0000
F81G	15	-1834.7174	3700.1234	0.0000	1.0000
F81IG	16	-1834.7176	3702.2168	0.0000	1.0000
K2P	12	-1850.4049	3725.2556	0.0000	1.0000
SYM	16	-1847.5875	3727.9567	0.0000	1.0000
JC69I	12	-1859.7229	3743.8916	0.0000	1.0000
JC69G	12	-1860.0729	3744.5914	0.0000	1.0000
F81	14	-1858.8344	3746.2704	0.0000	1.0000
JC69IG	13	-1860.0730	3746.6668	0.0000	1.0000
JC69	11	-1884.9898	3792.3562	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1794.3072	3620.6143	0.2539	0.2539
GTRI	20	-1790.3912	3620.7825	0.2334	0.4872
HKYG	16	-1794.5537	3621.1074	0.1984	0.6856
GTRG	20	-1790.6707	3621.3414	0.1765	0.8621
HKYIG	17	-1794.5538	3623.1075	0.0730	0.9351
GTRIG	21	-1790.6708	3623.3417	0.0649	1.0000
GTR	19	-1816.5322	3671.0643	0.0000	1.0000
HKY	15	-1820.7058	3671.4115	0.0000	1.0000
K2PI	13	-1823.1943	3672.3887	0.0000	1.0000
K2PG	13	-1823.3670	3672.7339	0.0000	1.0000
K2PIG	14	-1823.3670	3674.7341	0.0000	1.0000
SYMI	17	-1821.1145	3676.2290	0.0000	1.0000
SYMG	17	-1821.3189	3676.6379	0.0000	1.0000
SYMIG	18	-1821.3190	3678.6380	0.0000	1.0000
F81I	15	-1834.3111	3698.6223	0.0000	1.0000
F81G	15	-1834.7174	3699.4348	0.0000	1.0000
F81IG	16	-1834.7176	3701.4352	0.0000	1.0000
K2P	12	-1850.4049	3724.8099	0.0000	1.0000
SYM	16	-1847.5875	3727.1751	0.0000	1.0000
JC69I	12	-1859.7229	3743.4458	0.0000	1.0000
JC69G	12	-1860.0729	3744.1457	0.0000	1.0000
F81	14	-1858.8344	3745.6687	0.0000	1.0000
JC69IG	13	-1860.0730	3746.1461	0.0000	1.0000
JC69	11	-1884.9898	3791.9796	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1794.3072	3693.7260	0.5522	0.5522
HKYG	16	-1794.5537	3694.2191	0.4315	0.9837
HKYIG	17	-1794.5538	3700.7887	0.0162	0.9999
GTRI	20	-1790.3912	3712.1721	0.0001	1.0000
GTRG	20	-1790.6707	3712.7310	0.0000	1.0000
GTRIG	21	-1790.6708	3719.3008	0.0000	1.0000
K2PI	13	-1823.1943	3731.7919	0.0000	1.0000
K2PG	13	-1823.3670	3732.1372	0.0000	1.0000
K2PIG	14	-1823.3670	3738.7068	0.0000	1.0000
HKY	15	-1820.7058	3739.9537	0.0000	1.0000
SYMI	17	-1821.1145	3753.9102	0.0000	1.0000
SYMG	17	-1821.3189	3754.3190	0.0000	1.0000
GTR	19	-1816.5322	3757.8845	0.0000	1.0000
SYMIG	18	-1821.3190	3760.8887	0.0000	1.0000
F81I	15	-1834.3111	3767.1645	0.0000	1.0000
F81G	15	-1834.7174	3767.9770	0.0000	1.0000
F81IG	16	-1834.7176	3774.5469	0.0000	1.0000
K2P	12	-1850.4049	3779.6436	0.0000	1.0000
JC69I	12	-1859.7229	3798.2796	0.0000	1.0000
JC69G	12	-1860.0729	3798.9795	0.0000	1.0000
SYM	16	-1847.5875	3800.2868	0.0000	1.0000
JC69IG	13	-1860.0730	3805.5493	0.0000	1.0000
F81	14	-1858.8344	3809.6414	0.0000	1.0000
JC69	11	-1884.9898	3842.2439	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_5228.phylip = 1 - 713;
 Partition Dummy = 1:591-tomodel-phylip/chr7_5228.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_5228.phylip = 1 - 713;
 Partition Dummy = 1:591-tomodel-phylip/chr7_5228.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_5228.phylip = 1 - 713;
 Partition Dummy = 1:591-tomodel-phylip/chr7_5228.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

