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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1705.4832	3452.1698	0.4323	0.4323
GTRI	20	-1705.5777	3452.3589	0.3933	0.8256
GTRIG	21	-1705.4606	3454.2469	0.1530	0.9786
HKYI	16	-1713.5438	3459.8625	0.0092	0.9878
HKYG	16	-1713.5798	3459.9346	0.0089	0.9967
HKYIG	17	-1713.5709	3462.0149	0.0031	0.9999
GTR	19	-1714.5681	3468.2235	0.0001	1.0000
HKY	15	-1722.6276	3475.9379	0.0000	1.0000
F81I	15	-1749.0555	3528.7937	0.0000	1.0000
F81G	15	-1749.1380	3528.9587	0.0000	1.0000
F81IG	16	-1749.1363	3531.0474	0.0000	1.0000
F81	14	-1756.2655	3541.1276	0.0000	1.0000
SYMI	17	-1800.6992	3636.2715	0.0000	1.0000
SYMG	17	-1800.8874	3636.6478	0.0000	1.0000
SYMIG	18	-1800.7015	3638.3802	0.0000	1.0000
K2PI	13	-1810.6755	3647.8672	0.0000	1.0000
K2PG	13	-1810.8992	3648.3147	0.0000	1.0000
K2PIG	14	-1810.9086	3650.4137	0.0000	1.0000
SYM	16	-1809.4587	3651.6923	0.0000	1.0000
K2P	12	-1819.3375	3663.1169	0.0000	1.0000
JC69I	12	-1853.2535	3730.9489	0.0000	1.0000
JC69G	12	-1853.5361	3731.5142	0.0000	1.0000
JC69IG	13	-1853.2606	3733.0376	0.0000	1.0000
JC69	11	-1860.5179	3743.4092	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1705.4832	3450.9664	0.4298	0.4298
GTRI	20	-1705.5777	3451.1554	0.3911	0.8209
GTRIG	21	-1705.4606	3452.9212	0.1617	0.9827
HKYI	16	-1713.5438	3459.0875	0.0074	0.9901
HKYG	16	-1713.5798	3459.1597	0.0071	0.9972
HKYIG	17	-1713.5709	3461.1419	0.0027	0.9999
GTR	19	-1714.5681	3467.1362	0.0001	1.0000
HKY	15	-1722.6276	3475.2551	0.0000	1.0000
F81I	15	-1749.0555	3528.1110	0.0000	1.0000
F81G	15	-1749.1380	3528.2759	0.0000	1.0000
F81IG	16	-1749.1363	3530.2725	0.0000	1.0000
F81	14	-1756.2655	3540.5310	0.0000	1.0000
SYMI	17	-1800.6992	3635.3985	0.0000	1.0000
SYMG	17	-1800.8874	3635.7748	0.0000	1.0000
SYMIG	18	-1800.7015	3637.4030	0.0000	1.0000
K2PI	13	-1810.6755	3647.3509	0.0000	1.0000
K2PG	13	-1810.8992	3647.7983	0.0000	1.0000
K2PIG	14	-1810.9086	3649.8171	0.0000	1.0000
SYM	16	-1809.4587	3650.9174	0.0000	1.0000
K2P	12	-1819.3375	3662.6750	0.0000	1.0000
JC69I	12	-1853.2535	3730.5070	0.0000	1.0000
JC69G	12	-1853.5361	3731.0723	0.0000	1.0000
JC69IG	13	-1853.2606	3732.5213	0.0000	1.0000
JC69	11	-1860.5179	3743.0358	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1713.5438	3532.3333	0.4961	0.4961
HKYG	16	-1713.5798	3532.4055	0.4785	0.9746
HKYIG	17	-1713.5709	3538.9655	0.0180	0.9926
GTRG	20	-1705.4832	3542.5236	0.0030	0.9956
GTRI	20	-1705.5777	3542.7127	0.0028	0.9984
HKY	15	-1722.6276	3543.9231	0.0015	0.9999
GTRIG	21	-1705.4606	3549.0563	0.0001	1.0000
GTR	19	-1714.5681	3554.1156	0.0000	1.0000
F81I	15	-1749.0555	3596.7789	0.0000	1.0000
F81G	15	-1749.1380	3596.9438	0.0000	1.0000
F81IG	16	-1749.1363	3603.5183	0.0000	1.0000
F81	14	-1756.2655	3604.6210	0.0000	1.0000
K2PI	13	-1810.6755	3706.8631	0.0000	1.0000
K2PG	13	-1810.8992	3707.3105	0.0000	1.0000
SYMI	17	-1800.6992	3713.2221	0.0000	1.0000
SYMG	17	-1800.8874	3713.5984	0.0000	1.0000
K2PIG	14	-1810.9086	3713.9072	0.0000	1.0000
K2P	12	-1819.3375	3717.6093	0.0000	1.0000
SYMIG	18	-1800.7015	3719.8045	0.0000	1.0000
SYM	16	-1809.4587	3724.1631	0.0000	1.0000
JC69I	12	-1853.2535	3785.4414	0.0000	1.0000
JC69G	12	-1853.5361	3786.0066	0.0000	1.0000
JC69IG	13	-1853.2606	3792.0335	0.0000	1.0000
JC69	11	-1860.5179	3793.3923	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/chrZ_6110.phylip = 1 - 719;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_6110.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/chrZ_6110.phylip = 1 - 719;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_6110.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 (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_6110.phylip = 1 - 719;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_6110.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

