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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-1859.0633	3759.4269	0.4989	0.4989
GTRI	20	-1859.6008	3760.5018	0.2915	0.7904
GTRIG	21	-1859.0265	3761.4856	0.1782	0.9686
SYMG	17	-1865.7763	3766.4955	0.0146	0.9831
SYMI	17	-1866.0122	3766.9675	0.0115	0.9946
SYMIG	18	-1865.7804	3768.6165	0.0050	0.9997
HKYG	16	-1871.5817	3776.0002	0.0001	0.9998
HKYI	16	-1871.6914	3776.2197	0.0001	0.9999
HKYIG	17	-1871.5816	3778.1062	0.0000	1.0000
K2PG	13	-1876.9341	3780.4256	0.0000	1.0000
K2PI	13	-1876.9486	3780.4545	0.0000	1.0000
K2PIG	14	-1876.9340	3782.5122	0.0000	1.0000
GTR	19	-1888.6739	3816.5224	0.0000	1.0000
SYM	16	-1894.5491	3821.9351	0.0000	1.0000
HKY	15	-1902.1297	3834.9966	0.0000	1.0000
K2P	12	-1907.0584	3838.5939	0.0000	1.0000
F81I	15	-1904.1019	3838.9411	0.0000	1.0000
F81G	15	-1904.1184	3838.9742	0.0000	1.0000
JC69I	12	-1908.2384	3840.9539	0.0000	1.0000
JC69G	12	-1908.2866	3841.0502	0.0000	1.0000
F81IG	16	-1904.1184	3841.0738	0.0000	1.0000
JC69IG	13	-1908.2865	3843.1305	0.0000	1.0000
F81	14	-1933.4966	3895.6373	0.0000	1.0000
JC69	11	-1937.5441	3897.4912	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1859.0633	3758.1266	0.4953	0.4953
GTRI	20	-1859.6008	3759.2015	0.2894	0.7847
GTRIG	21	-1859.0265	3760.0531	0.1890	0.9737
SYMG	17	-1865.7763	3765.5525	0.0121	0.9858
SYMI	17	-1866.0122	3766.0245	0.0095	0.9953
SYMIG	18	-1865.7804	3767.5609	0.0044	0.9998
HKYG	16	-1871.5817	3775.1633	0.0001	0.9999
HKYI	16	-1871.6914	3775.3828	0.0001	0.9999
HKYIG	17	-1871.5816	3777.1633	0.0000	1.0000
K2PG	13	-1876.9341	3779.8681	0.0000	1.0000
K2PI	13	-1876.9486	3779.8971	0.0000	1.0000
K2PIG	14	-1876.9340	3781.8680	0.0000	1.0000
GTR	19	-1888.6739	3815.3477	0.0000	1.0000
SYM	16	-1894.5491	3821.0982	0.0000	1.0000
HKY	15	-1902.1297	3834.2593	0.0000	1.0000
K2P	12	-1907.0584	3838.1169	0.0000	1.0000
F81I	15	-1904.1019	3838.2038	0.0000	1.0000
F81G	15	-1904.1184	3838.2369	0.0000	1.0000
F81IG	16	-1904.1184	3840.2369	0.0000	1.0000
JC69I	12	-1908.2384	3840.4769	0.0000	1.0000
JC69G	12	-1908.2866	3840.5731	0.0000	1.0000
JC69IG	13	-1908.2865	3842.5731	0.0000	1.0000
F81	14	-1933.4966	3894.9931	0.0000	1.0000
JC69	11	-1937.5441	3897.0882	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1876.9341	3838.4044	0.4258	0.4258
K2PI	13	-1876.9486	3838.4334	0.4197	0.8455
SYMG	17	-1865.7763	3842.1000	0.0671	0.9126
SYMI	17	-1866.0122	3842.5719	0.0530	0.9656
K2PIG	14	-1876.9340	3844.9071	0.0165	0.9821
HKYG	16	-1871.5817	3847.2080	0.0052	0.9873
HKYI	16	-1871.6914	3847.4274	0.0047	0.9920
GTRG	20	-1859.0633	3848.1824	0.0032	0.9952
SYMIG	18	-1865.7804	3848.6111	0.0026	0.9978
GTRI	20	-1859.6008	3849.2573	0.0019	0.9997
HKYIG	17	-1871.5816	3853.7107	0.0002	0.9999
GTRIG	21	-1859.0265	3854.6117	0.0001	1.0000
K2P	12	-1907.0584	3892.1503	0.0000	1.0000
SYM	16	-1894.5491	3893.1428	0.0000	1.0000
JC69I	12	-1908.2384	3894.5103	0.0000	1.0000
JC69G	12	-1908.2866	3894.6066	0.0000	1.0000
GTR	19	-1888.6739	3900.9008	0.0000	1.0000
JC69IG	13	-1908.2865	3901.1093	0.0000	1.0000
HKY	15	-1902.1297	3901.8012	0.0000	1.0000
F81I	15	-1904.1019	3905.7457	0.0000	1.0000
F81G	15	-1904.1184	3905.7787	0.0000	1.0000
F81IG	16	-1904.1184	3912.2815	0.0000	1.0000
JC69	11	-1937.5441	3946.6189	0.0000	1.0000
F81	14	-1933.4966	3958.0322	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_17005.phylip = 1 - 667;
 Partition Dummy = 1:591-tomodel-phylip/chr2_17005.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_17005.phylip = 1 - 667;
 Partition Dummy = 1:591-tomodel-phylip/chr2_17005.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 (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_17005.phylip = 1 - 667;
 Partition Dummy = 1:591-tomodel-phylip/chr2_17005.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

