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

Input data from file "591-tomodel-phylip/chrZ_7889.phylip" (7 taxa, 780 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	16	-1716.6567	3466.0264	0.3693	0.3693
HKYI	16	-1716.9207	3466.5544	0.2836	0.6530
HKYIG	17	-1716.5698	3467.9427	0.1417	0.7946
GTRG	20	-1713.8099	3468.7265	0.0957	0.8904
GTRI	20	-1714.0743	3469.2554	0.0735	0.9639
GTRIG	21	-1713.7377	3470.6944	0.0358	0.9997
HKY	15	-1725.0196	3480.6676	0.0002	0.9999
GTR	19	-1721.9522	3482.9044	0.0001	1.0000
K2PG	13	-1733.3585	3493.1922	0.0000	1.0000
K2PI	13	-1733.6236	3493.7225	0.0000	1.0000
K2PIG	14	-1733.2238	3494.9966	0.0000	1.0000
SYMG	17	-1731.8791	3498.5613	0.0000	1.0000
SYMI	17	-1732.1491	3499.1014	0.0000	1.0000
SYMIG	18	-1731.7567	3500.4123	0.0000	1.0000
F81G	15	-1737.1267	3504.8817	0.0000	1.0000
F81I	15	-1737.3705	3505.3692	0.0000	1.0000
F81IG	16	-1737.0498	3506.8126	0.0000	1.0000
K2P	12	-1742.2480	3508.9028	0.0000	1.0000
SYM	16	-1740.5718	3513.8565	0.0000	1.0000
F81	14	-1745.4008	3519.3506	0.0000	1.0000
JC69G	12	-1752.9490	3530.3048	0.0000	1.0000
JC69I	12	-1753.1948	3530.7964	0.0000	1.0000
JC69IG	13	-1752.8435	3532.1623	0.0000	1.0000
JC69	11	-1761.4908	3545.3253	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1716.6567	3465.3134	0.3505	0.3505
HKYI	16	-1716.9207	3465.8415	0.2692	0.6197
HKYIG	17	-1716.5698	3467.1396	0.1407	0.7604
GTRG	20	-1713.8099	3467.6198	0.1106	0.8710
GTRI	20	-1714.0743	3468.1487	0.0849	0.9559
GTRIG	21	-1713.7377	3469.4754	0.0437	0.9997
HKY	15	-1725.0196	3480.0393	0.0002	0.9999
GTR	19	-1721.9522	3481.9044	0.0001	1.0000
K2PG	13	-1733.3585	3492.7170	0.0000	1.0000
K2PI	13	-1733.6236	3493.2473	0.0000	1.0000
K2PIG	14	-1733.2238	3494.4476	0.0000	1.0000
SYMG	17	-1731.8791	3497.7582	0.0000	1.0000
SYMI	17	-1732.1491	3498.2982	0.0000	1.0000
SYMIG	18	-1731.7567	3499.5135	0.0000	1.0000
F81G	15	-1737.1267	3504.2534	0.0000	1.0000
F81I	15	-1737.3705	3504.7410	0.0000	1.0000
F81IG	16	-1737.0498	3506.0996	0.0000	1.0000
K2P	12	-1742.2480	3508.4961	0.0000	1.0000
SYM	16	-1740.5718	3513.1435	0.0000	1.0000
F81	14	-1745.4008	3518.8016	0.0000	1.0000
JC69G	12	-1752.9490	3529.8980	0.0000	1.0000
JC69I	12	-1753.1948	3530.3896	0.0000	1.0000
JC69IG	13	-1752.8435	3531.6871	0.0000	1.0000
JC69	11	-1761.4908	3544.9815	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1716.6567	3539.8621	0.5507	0.5507
HKYI	16	-1716.9207	3540.3902	0.4229	0.9737
HKYIG	17	-1716.5698	3546.3476	0.0215	0.9952
HKY	15	-1725.0196	3549.9287	0.0036	0.9988
K2PG	13	-1733.3585	3553.2878	0.0007	0.9994
K2PI	13	-1733.6236	3553.8181	0.0005	0.9999
K2PIG	14	-1733.2238	3559.6777	0.0000	1.0000
GTRG	20	-1713.8099	3560.8057	0.0000	1.0000
GTRI	20	-1714.0743	3561.3345	0.0000	1.0000
K2P	12	-1742.2480	3564.4076	0.0000	1.0000
GTRIG	21	-1713.7377	3567.3206	0.0000	1.0000
GTR	19	-1721.9522	3570.4310	0.0000	1.0000
F81G	15	-1737.1267	3574.1428	0.0000	1.0000
F81I	15	-1737.3705	3574.6304	0.0000	1.0000
SYMG	17	-1731.8791	3576.9662	0.0000	1.0000
SYMI	17	-1732.1491	3577.5062	0.0000	1.0000
F81IG	16	-1737.0498	3580.6483	0.0000	1.0000
SYMIG	18	-1731.7567	3583.3808	0.0000	1.0000
F81	14	-1745.4008	3584.0317	0.0000	1.0000
JC69G	12	-1752.9490	3585.8096	0.0000	1.0000
JC69I	12	-1753.1948	3586.3011	0.0000	1.0000
SYM	16	-1740.5718	3587.6922	0.0000	1.0000
JC69IG	13	-1752.8435	3592.2579	0.0000	1.0000
JC69	11	-1761.4908	3596.2338	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 591-tomodel-phylip/chrZ_7889.phylip = 1 - 780;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_7889.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 591-tomodel-phylip/chrZ_7889.phylip = 1 - 780;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_7889.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 591-tomodel-phylip/chrZ_7889.phylip = 1 - 780;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_7889.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

