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

Input data from file "591-tomodel-phylip/chr7_10687.phylip" (7 taxa, 610 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	-1568.5419	3170.0012	0.4694	0.4694
HKYI	16	-1569.1243	3171.1659	0.2622	0.7315
HKYIG	17	-1568.3102	3171.6542	0.2054	0.9369
GTRG	20	-1566.9799	3175.3860	0.0318	0.9687
GTRI	20	-1567.6172	3176.6605	0.0168	0.9855
GTRIG	21	-1566.7305	3177.0324	0.0140	0.9995
K2PG	13	-1579.2203	3185.0512	0.0003	0.9997
K2PI	13	-1579.7852	3186.1812	0.0001	0.9999
K2PIG	14	-1578.9923	3186.6904	0.0001	1.0000
SYMG	17	-1578.9240	3192.8818	0.0000	1.0000
SYMI	17	-1579.5004	3194.0345	0.0000	1.0000
SYMIG	18	-1578.6914	3194.5401	0.0000	1.0000
HKY	15	-1584.6608	3200.1297	0.0000	1.0000
GTR	19	-1583.3824	3206.0530	0.0000	1.0000
F81G	15	-1590.2599	3211.3280	0.0000	1.0000
F81I	15	-1590.7325	3212.2731	0.0000	1.0000
F81IG	16	-1590.1053	3213.1280	0.0000	1.0000
K2P	12	-1595.1145	3214.7516	0.0000	1.0000
SYM	16	-1594.7883	3222.4939	0.0000	1.0000
JC69G	12	-1599.5888	3223.7002	0.0000	1.0000
JC69I	12	-1600.0533	3224.6292	0.0000	1.0000
JC69IG	13	-1599.4350	3225.4806	0.0000	1.0000
F81	14	-1605.7542	3240.2144	0.0000	1.0000
JC69	11	-1614.9949	3252.4312	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1568.5419	3169.0838	0.4549	0.4549
HKYI	16	-1569.1243	3170.2485	0.2541	0.7091
HKYIG	17	-1568.3102	3170.6204	0.2110	0.9201
GTRG	20	-1566.9799	3173.9599	0.0397	0.9598
GTRI	20	-1567.6172	3175.2344	0.0210	0.9808
GTRIG	21	-1566.7305	3175.4610	0.0188	0.9996
K2PG	13	-1579.2203	3184.4405	0.0002	0.9998
K2PI	13	-1579.7852	3185.5705	0.0001	0.9999
K2PIG	14	-1578.9923	3185.9845	0.0001	1.0000
SYMG	17	-1578.9240	3191.8480	0.0000	1.0000
SYMI	17	-1579.5004	3193.0007	0.0000	1.0000
SYMIG	18	-1578.6914	3193.3827	0.0000	1.0000
HKY	15	-1584.6608	3199.3216	0.0000	1.0000
GTR	19	-1583.3824	3204.7648	0.0000	1.0000
F81G	15	-1590.2599	3210.5199	0.0000	1.0000
F81I	15	-1590.7325	3211.4650	0.0000	1.0000
F81IG	16	-1590.1053	3212.2106	0.0000	1.0000
K2P	12	-1595.1145	3214.2289	0.0000	1.0000
SYM	16	-1594.7883	3221.5766	0.0000	1.0000
JC69G	12	-1599.5888	3223.1776	0.0000	1.0000
JC69I	12	-1600.0533	3224.1066	0.0000	1.0000
JC69IG	13	-1599.4350	3224.8699	0.0000	1.0000
F81	14	-1605.7542	3239.5085	0.0000	1.0000
JC69	11	-1614.9949	3251.9898	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1568.5419	3239.6992	0.4605	0.4605
HKYI	16	-1569.1243	3240.8638	0.2572	0.7177
K2PG	13	-1579.2203	3241.8155	0.1598	0.8775
K2PI	13	-1579.7852	3242.9455	0.0908	0.9684
HKYIG	17	-1568.3102	3245.6492	0.0235	0.9919
K2PIG	14	-1578.9923	3247.7729	0.0081	1.0000
GTRG	20	-1566.9799	3262.2290	0.0000	1.0000
GTRI	20	-1567.6172	3263.5036	0.0000	1.0000
HKY	15	-1584.6608	3265.5235	0.0000	1.0000
SYMG	17	-1578.9240	3266.8768	0.0000	1.0000
K2P	12	-1595.1145	3267.1904	0.0000	1.0000
SYMI	17	-1579.5004	3268.0295	0.0000	1.0000
GTRIG	21	-1566.7305	3268.1436	0.0000	1.0000
SYMIG	18	-1578.6914	3272.8250	0.0000	1.0000
JC69G	12	-1599.5888	3276.1391	0.0000	1.0000
F81G	15	-1590.2599	3276.7218	0.0000	1.0000
JC69I	12	-1600.0533	3277.0681	0.0000	1.0000
F81I	15	-1590.7325	3277.6669	0.0000	1.0000
JC69IG	13	-1599.4350	3282.2449	0.0000	1.0000
F81IG	16	-1590.1053	3282.8260	0.0000	1.0000
GTR	19	-1583.3824	3288.6205	0.0000	1.0000
SYM	16	-1594.7883	3292.1919	0.0000	1.0000
JC69	11	-1614.9949	3300.5378	0.0000	1.0000
F81	14	-1605.7542	3301.2969	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/chr7_10687.phylip = 1 - 610;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10687.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/chr7_10687.phylip = 1 - 610;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10687.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/chr7_10687.phylip = 1 - 610;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10687.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

