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

Input data from file "183-phylip/chrZ_11420.phylip" (29 taxa, 596 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	60	-1608.7475	3351.1771	0.5722	0.5722
HKYI	60	-1609.9236	3353.5294	0.1765	0.7488
HKYIG	61	-1608.7249	3353.6145	0.1692	0.9179
GTRG	64	-1606.1759	3356.0202	0.0508	0.9687
GTRI	64	-1607.3750	3358.4185	0.0153	0.9840
GTRIG	65	-1606.1624	3358.5136	0.0146	0.9987
HKY	59	-1616.1743	3363.5575	0.0012	0.9998
GTR	63	-1613.0947	3367.3472	0.0002	1.0000
K2PG	57	-1627.2415	3380.7729	0.0000	1.0000
K2PIG	58	-1627.1450	3383.0349	0.0000	1.0000
K2PI	57	-1628.7872	3383.8644	0.0000	1.0000
SYMG	61	-1625.9332	3388.0311	0.0000	1.0000
SYMIG	62	-1625.8493	3390.3552	0.0000	1.0000
SYMI	61	-1627.1887	3390.5421	0.0000	1.0000
K2P	56	-1635.8456	3395.5354	0.0000	1.0000
SYM	60	-1633.9823	3401.6468	0.0000	1.0000
F81G	59	-1637.6826	3406.5742	0.0000	1.0000
F81I	59	-1638.8938	3408.9966	0.0000	1.0000
F81IG	60	-1637.6653	3409.0129	0.0000	1.0000
F81	58	-1644.4642	3417.6734	0.0000	1.0000
JC69G	56	-1653.9624	3431.7689	0.0000	1.0000
JC69IG	57	-1653.9037	3434.0973	0.0000	1.0000
JC69I	56	-1655.5777	3434.9996	0.0000	1.0000
JC69	55	-1661.5453	3444.4980	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1608.7475	3337.4949	0.4789	0.4789
HKYIG	61	-1608.7249	3339.4497	0.1802	0.6590
HKYI	60	-1609.9236	3339.8472	0.1477	0.8067
GTRG	64	-1606.1759	3340.3517	0.1148	0.9215
GTRIG	65	-1606.1624	3342.3249	0.0428	0.9643
GTRI	64	-1607.3750	3342.7499	0.0346	0.9989
HKY	59	-1616.1743	3350.3485	0.0008	0.9997
GTR	63	-1613.0947	3352.1893	0.0003	1.0000
K2PG	57	-1627.2415	3368.4829	0.0000	1.0000
K2PIG	58	-1627.1450	3370.2900	0.0000	1.0000
K2PI	57	-1628.7872	3371.5745	0.0000	1.0000
SYMG	61	-1625.9332	3373.8663	0.0000	1.0000
SYMIG	62	-1625.8493	3375.6985	0.0000	1.0000
SYMI	61	-1627.1887	3376.3773	0.0000	1.0000
K2P	56	-1635.8456	3383.6912	0.0000	1.0000
SYM	60	-1633.9823	3387.9646	0.0000	1.0000
F81G	59	-1637.6826	3393.3652	0.0000	1.0000
F81IG	60	-1637.6653	3395.3306	0.0000	1.0000
F81I	59	-1638.8938	3395.7876	0.0000	1.0000
F81	58	-1644.4642	3404.9285	0.0000	1.0000
JC69G	56	-1653.9624	3419.9248	0.0000	1.0000
JC69IG	57	-1653.9037	3421.8073	0.0000	1.0000
JC69I	56	-1655.5777	3423.1555	0.0000	1.0000
JC69	55	-1661.5453	3433.0906	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1608.7475	3600.9093	0.7325	0.7325
HKYI	60	-1609.9236	3603.2616	0.2260	0.9585
HKYIG	61	-1608.7249	3607.2544	0.0307	0.9892
HKY	59	-1616.1743	3609.3727	0.0106	0.9998
K2PG	57	-1627.2415	3618.7267	0.0001	0.9999
GTRG	64	-1606.1759	3621.3271	0.0000	1.0000
K2PI	57	-1628.7872	3621.8182	0.0000	1.0000
GTRI	64	-1607.3750	3623.7253	0.0000	1.0000
K2PIG	58	-1627.1450	3624.9240	0.0000	1.0000
GTRIG	65	-1606.1624	3627.6905	0.0000	1.0000
GTR	63	-1613.0947	3628.7745	0.0000	1.0000
K2P	56	-1635.8456	3629.5447	0.0000	1.0000
SYMG	61	-1625.9332	3641.6710	0.0000	1.0000
SYMI	61	-1627.1887	3644.1820	0.0000	1.0000
SYMIG	62	-1625.8493	3647.8935	0.0000	1.0000
SYM	60	-1633.9823	3651.3790	0.0000	1.0000
F81G	59	-1637.6826	3652.3894	0.0000	1.0000
F81I	59	-1638.8938	3654.8118	0.0000	1.0000
F81IG	60	-1637.6653	3658.7451	0.0000	1.0000
F81	58	-1644.4642	3659.5624	0.0000	1.0000
JC69G	56	-1653.9624	3665.7783	0.0000	1.0000
JC69I	56	-1655.5777	3669.0090	0.0000	1.0000
JC69IG	57	-1653.9037	3672.0511	0.0000	1.0000
JC69	55	-1661.5453	3674.5539	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 183-phylip/chrZ_11420.phylip = 1 - 596;
 Partition Dummy = 1:183-phylip/chrZ_11420.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 183-phylip/chrZ_11420.phylip = 1 - 596;
 Partition Dummy = 1:183-phylip/chrZ_11420.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 183-phylip/chrZ_11420.phylip = 1 - 596;
 Partition Dummy = 1:183-phylip/chrZ_11420.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

