
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:54 02/17/2011

Input data from file "183-phylip/chr6_8064.phylip" (29 taxa, 629 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKY


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	64	-1821.0878	3784.9274	0.3720	0.3720
GTRI	64	-1821.7830	3786.3178	0.1856	0.5577
HKYG	60	-1826.8963	3786.6799	0.1549	0.7126
GTRIG	65	-1821.0872	3787.4142	0.1073	0.8199
HKYI	60	-1827.7552	3788.3978	0.0656	0.8855
HKYIG	61	-1826.8958	3789.1320	0.0455	0.9309
GTR	63	-1824.4867	3789.2459	0.0429	0.9739
HKY	59	-1829.8983	3790.2395	0.0261	1.0000
K2PG	57	-1853.0621	3831.7039	0.0000	1.0000
K2PI	57	-1854.0704	3833.7206	0.0000	1.0000
K2PIG	58	-1853.0615	3834.1300	0.0000	1.0000
K2P	56	-1856.2841	3835.7290	0.0000	1.0000
SYMG	61	-1851.5701	3838.4806	0.0000	1.0000
SYMI	61	-1852.4710	3840.2824	0.0000	1.0000
SYMIG	62	-1851.5696	3840.9413	0.0000	1.0000
SYM	60	-1854.7319	3842.3510	0.0000	1.0000
F81G	59	-1858.6852	3847.8133	0.0000	1.0000
F81I	59	-1859.4648	3849.3724	0.0000	1.0000
F81IG	60	-1858.6848	3850.2569	0.0000	1.0000
F81	58	-1861.2522	3850.5114	0.0000	1.0000
JC69G	56	-1881.0276	3885.2161	0.0000	1.0000
JC69I	56	-1881.8653	3886.8914	0.0000	1.0000
JC69IG	57	-1881.0271	3887.6339	0.0000	1.0000
JC69	55	-1883.7917	3888.3338	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1821.0878	3770.1756	0.4398	0.4398
GTRI	64	-1821.7830	3771.5661	0.2195	0.6593
GTRIG	65	-1821.0872	3772.1744	0.1619	0.8212
HKYG	60	-1826.8963	3773.7926	0.0721	0.8933
GTR	63	-1824.4867	3774.9733	0.0399	0.9332
HKYI	60	-1827.7552	3775.5105	0.0305	0.9637
HKYIG	61	-1826.8958	3775.7916	0.0265	0.9903
HKY	59	-1829.8983	3777.7966	0.0097	1.0000
K2PG	57	-1853.0621	3820.1242	0.0000	1.0000
K2PIG	58	-1853.0615	3822.1230	0.0000	1.0000
K2PI	57	-1854.0704	3822.1409	0.0000	1.0000
K2P	56	-1856.2841	3824.5682	0.0000	1.0000
SYMG	61	-1851.5701	3825.1402	0.0000	1.0000
SYMI	61	-1852.4710	3826.9420	0.0000	1.0000
SYMIG	62	-1851.5696	3827.1391	0.0000	1.0000
SYM	60	-1854.7319	3829.4637	0.0000	1.0000
F81G	59	-1858.6852	3835.3704	0.0000	1.0000
F81I	59	-1859.4648	3836.9295	0.0000	1.0000
F81IG	60	-1858.6848	3837.3696	0.0000	1.0000
F81	58	-1861.2522	3838.5044	0.0000	1.0000
JC69G	56	-1881.0276	3874.0552	0.0000	1.0000
JC69I	56	-1881.8653	3875.7305	0.0000	1.0000
JC69IG	57	-1881.0271	3876.0543	0.0000	1.0000
JC69	55	-1883.7917	3877.5834	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	59	-1829.8983	4040.0004	0.4596	0.4596
HKYG	60	-1826.8963	4040.4405	0.3688	0.8283
HKYI	60	-1827.7552	4042.1584	0.1562	0.9846
HKYIG	61	-1826.8958	4046.8836	0.0147	0.9993
GTRG	64	-1821.0878	4054.6000	0.0003	0.9996
GTR	63	-1824.4867	4054.9536	0.0003	0.9998
GTRI	64	-1821.7830	4055.9905	0.0002	1.0000
GTRIG	65	-1821.0872	4061.0429	0.0000	1.0000
K2P	56	-1856.2841	4073.4395	0.0000	1.0000
K2PG	57	-1853.0621	4073.4397	0.0000	1.0000
K2PI	57	-1854.0704	4075.4564	0.0000	1.0000
K2PIG	58	-1853.0615	4079.8826	0.0000	1.0000
SYM	60	-1854.7319	4096.1116	0.0000	1.0000
SYMG	61	-1851.5701	4096.2322	0.0000	1.0000
F81	58	-1861.2522	4096.2640	0.0000	1.0000
F81G	59	-1858.6852	4097.5741	0.0000	1.0000
SYMI	61	-1852.4710	4098.0340	0.0000	1.0000
F81I	59	-1859.4648	4099.1333	0.0000	1.0000
SYMIG	62	-1851.5696	4102.6753	0.0000	1.0000
F81IG	60	-1858.6848	4104.0174	0.0000	1.0000
JC69	55	-1883.7917	4122.0106	0.0000	1.0000
JC69G	56	-1881.0276	4122.9266	0.0000	1.0000
JC69I	56	-1881.8653	4124.6019	0.0000	1.0000
JC69IG	57	-1881.0271	4129.3697	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 183-phylip/chr6_8064.phylip = 1 - 629;
 Partition Dummy = 1:183-phylip/chr6_8064.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 183-phylip/chr6_8064.phylip = 1 - 629;
 Partition Dummy = 1:183-phylip/chr6_8064.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 (HKY)]

BEGIN MRBAYES;
 Charset 183-phylip/chr6_8064.phylip = 1 - 629;
 Partition Dummy = 1:183-phylip/chr6_8064.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

