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

Input data from file "591-tomodel-phylip/chr3_9734.phylip" (7 taxa, 543 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	-1575.0877	3183.2096	0.4094	0.4094
HKYI	16	-1575.0994	3183.2330	0.4046	0.8140
HKYIG	17	-1575.0878	3185.3413	0.1410	0.9550
GTRG	20	-1573.8364	3189.2820	0.0197	0.9746
GTRI	20	-1573.8847	3189.3787	0.0187	0.9933
GTRIG	21	-1573.8365	3191.4465	0.0067	1.0000
HKY	15	-1592.0892	3215.0892	0.0000	1.0000
GTR	19	-1591.1729	3221.7990	0.0000	1.0000
F81I	15	-1595.6979	3222.3065	0.0000	1.0000
F81G	15	-1595.7856	3222.4820	0.0000	1.0000
F81IG	16	-1595.7856	3224.6053	0.0000	1.0000
F81	14	-1609.8530	3248.5014	0.0000	1.0000
K2PG	13	-1614.1528	3254.9937	0.0000	1.0000
K2PI	13	-1614.4934	3255.6749	0.0000	1.0000
K2PIG	14	-1614.1597	3257.1149	0.0000	1.0000
SYMG	17	-1611.6062	3258.3781	0.0000	1.0000
SYMI	17	-1612.2272	3259.6201	0.0000	1.0000
SYMIG	18	-1611.6057	3260.5167	0.0000	1.0000
JC69G	12	-1628.5436	3281.6759	0.0000	1.0000
K2P	12	-1628.8095	3282.2077	0.0000	1.0000
JC69I	12	-1629.2854	3283.1595	0.0000	1.0000
JC69IG	13	-1629.2319	3285.1519	0.0000	1.0000
SYM	16	-1626.3330	3285.7003	0.0000	1.0000
JC69	11	-1641.1973	3304.8918	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1575.0877	3182.1754	0.3992	0.3992
HKYI	16	-1575.0994	3182.1988	0.3946	0.7938
HKYIG	17	-1575.0878	3184.1756	0.1469	0.9407
GTRG	20	-1573.8364	3187.6728	0.0256	0.9662
GTRI	20	-1573.8847	3187.7695	0.0243	0.9906
GTRIG	21	-1573.8365	3189.6730	0.0094	1.0000
HKY	15	-1592.0892	3214.1784	0.0000	1.0000
GTR	19	-1591.1729	3220.3458	0.0000	1.0000
F81I	15	-1595.6979	3221.3957	0.0000	1.0000
F81G	15	-1595.7856	3221.5711	0.0000	1.0000
F81IG	16	-1595.7856	3223.5711	0.0000	1.0000
F81	14	-1609.8530	3247.7060	0.0000	1.0000
K2PG	13	-1614.1528	3254.3056	0.0000	1.0000
K2PI	13	-1614.4934	3254.9868	0.0000	1.0000
K2PIG	14	-1614.1597	3256.3195	0.0000	1.0000
SYMG	17	-1611.6062	3257.2123	0.0000	1.0000
SYMI	17	-1612.2272	3258.4544	0.0000	1.0000
SYMIG	18	-1611.6057	3259.2113	0.0000	1.0000
JC69G	12	-1628.5436	3281.0873	0.0000	1.0000
K2P	12	-1628.8095	3281.6190	0.0000	1.0000
JC69I	12	-1629.2854	3282.5708	0.0000	1.0000
JC69IG	13	-1629.2319	3284.4638	0.0000	1.0000
SYM	16	-1626.3330	3284.6661	0.0000	1.0000
JC69	11	-1641.1973	3304.3946	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1575.0877	3250.9292	0.4923	0.4923
HKYI	16	-1575.0994	3250.9526	0.4866	0.9789
HKYIG	17	-1575.0878	3257.2264	0.0211	1.0000
GTRG	20	-1573.8364	3273.6150	0.0000	1.0000
GTRI	20	-1573.8847	3273.7117	0.0000	1.0000
HKY	15	-1592.0892	3278.6350	0.0000	1.0000
GTRIG	21	-1573.8365	3279.9123	0.0000	1.0000
F81I	15	-1595.6979	3285.8523	0.0000	1.0000
F81G	15	-1595.7856	3286.0278	0.0000	1.0000
F81IG	16	-1595.7856	3292.3248	0.0000	1.0000
GTR	19	-1591.1729	3301.9909	0.0000	1.0000
F81	14	-1609.8530	3307.8655	0.0000	1.0000
K2PG	13	-1614.1528	3310.1680	0.0000	1.0000
K2PI	13	-1614.4934	3310.8492	0.0000	1.0000
K2PIG	14	-1614.1597	3316.4790	0.0000	1.0000
SYMG	17	-1611.6062	3330.2632	0.0000	1.0000
SYMI	17	-1612.2272	3331.5053	0.0000	1.0000
JC69G	12	-1628.5436	3332.6526	0.0000	1.0000
K2P	12	-1628.8095	3333.1843	0.0000	1.0000
JC69I	12	-1629.2854	3334.1361	0.0000	1.0000
SYMIG	18	-1611.6057	3336.5593	0.0000	1.0000
JC69IG	13	-1629.2319	3340.3263	0.0000	1.0000
JC69	11	-1641.1973	3351.6628	0.0000	1.0000
SYM	16	-1626.3330	3353.4198	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/chr3_9734.phylip = 1 - 543;
 Partition Dummy = 1:591-tomodel-phylip/chr3_9734.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/chr3_9734.phylip = 1 - 543;
 Partition Dummy = 1:591-tomodel-phylip/chr3_9734.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/chr3_9734.phylip = 1 - 543;
 Partition Dummy = 1:591-tomodel-phylip/chr3_9734.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

