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

Input data from file "591-tomodel-phylip/chr13_2812.phylip" (7 taxa, 397 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	-1573.5535	3180.5386	0.6225	0.6225
HKYIG	17	-1573.5552	3182.7252	0.2086	0.8311
GTRG	20	-1571.1767	3184.5875	0.0822	0.9133
HKYI	16	-1576.0057	3185.4430	0.0536	0.9669
GTRIG	21	-1571.1767	3186.8175	0.0270	0.9939
GTRI	20	-1573.7785	3189.7911	0.0061	1.0000
HKY	15	-1591.5776	3214.4151	0.0000	1.0000
GTR	19	-1589.5649	3219.1458	0.0000	1.0000
K2PG	13	-1598.0685	3223.0875	0.0000	1.0000
K2PIG	14	-1598.0634	3225.2262	0.0000	1.0000
K2PI	13	-1599.6112	3226.1728	0.0000	1.0000
SYMG	17	-1597.2138	3230.0424	0.0000	1.0000
SYMIG	18	-1597.2013	3232.2120	0.0000	1.0000
SYMI	17	-1598.6136	3232.8420	0.0000	1.0000
K2P	12	-1614.9565	3254.7255	0.0000	1.0000
SYM	16	-1614.1634	3261.7583	0.0000	1.0000
F81G	15	-1617.5544	3266.3686	0.0000	1.0000
F81IG	16	-1617.5577	3268.5471	0.0000	1.0000
F81I	15	-1619.4651	3270.1901	0.0000	1.0000
F81	14	-1633.5369	3296.1734	0.0000	1.0000
JC69G	12	-1643.3190	3311.4504	0.0000	1.0000
JC69IG	13	-1643.3181	3313.5867	0.0000	1.0000
JC69I	12	-1644.7239	3314.2604	0.0000	1.0000
JC69	11	-1658.2379	3339.1614	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1573.5535	3179.1070	0.5754	0.5754
HKYIG	17	-1573.5552	3181.1105	0.2113	0.7868
GTRG	20	-1571.1767	3182.3534	0.1135	0.9003
HKYI	16	-1576.0057	3184.0114	0.0495	0.9498
GTRIG	21	-1571.1767	3184.3535	0.0418	0.9916
GTRI	20	-1573.7785	3187.5571	0.0084	1.0000
HKY	15	-1591.5776	3213.1552	0.0000	1.0000
GTR	19	-1589.5649	3217.1299	0.0000	1.0000
K2PG	13	-1598.0685	3222.1371	0.0000	1.0000
K2PIG	14	-1598.0634	3224.1268	0.0000	1.0000
K2PI	13	-1599.6112	3225.2224	0.0000	1.0000
SYMG	17	-1597.2138	3228.4276	0.0000	1.0000
SYMIG	18	-1597.2013	3230.4025	0.0000	1.0000
SYMI	17	-1598.6136	3231.2272	0.0000	1.0000
K2P	12	-1614.9565	3253.9130	0.0000	1.0000
SYM	16	-1614.1634	3260.3267	0.0000	1.0000
F81G	15	-1617.5544	3265.1088	0.0000	1.0000
F81IG	16	-1617.5577	3267.1155	0.0000	1.0000
F81I	15	-1619.4651	3268.9303	0.0000	1.0000
F81	14	-1633.5369	3295.0739	0.0000	1.0000
JC69G	12	-1643.3190	3310.6379	0.0000	1.0000
JC69IG	13	-1643.3181	3312.6363	0.0000	1.0000
JC69I	12	-1644.7239	3313.4479	0.0000	1.0000
JC69	11	-1658.2379	3338.4757	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1573.5535	3242.8500	0.8801	0.8801
HKYI	16	-1576.0057	3247.7544	0.0758	0.9558
HKYIG	17	-1573.5552	3248.8374	0.0441	0.9999
GTRG	20	-1571.1767	3262.0321	0.0001	1.0000
GTRI	20	-1573.7785	3267.2358	0.0000	1.0000
GTRIG	21	-1571.1767	3268.0161	0.0000	1.0000
HKY	15	-1591.5776	3272.9143	0.0000	1.0000
K2PG	13	-1598.0685	3273.9283	0.0000	1.0000
K2PI	13	-1599.6112	3277.0136	0.0000	1.0000
K2PIG	14	-1598.0634	3279.9019	0.0000	1.0000
GTR	19	-1589.5649	3292.8247	0.0000	1.0000
SYMG	17	-1597.2138	3296.1545	0.0000	1.0000
SYMI	17	-1598.6136	3298.9541	0.0000	1.0000
K2P	12	-1614.9565	3301.7202	0.0000	1.0000
SYMIG	18	-1597.2013	3302.1134	0.0000	1.0000
SYM	16	-1614.1634	3324.0697	0.0000	1.0000
F81G	15	-1617.5544	3324.8678	0.0000	1.0000
F81I	15	-1619.4651	3328.6893	0.0000	1.0000
F81IG	16	-1617.5577	3330.8585	0.0000	1.0000
F81	14	-1633.5369	3350.8490	0.0000	1.0000
JC69G	12	-1643.3190	3358.4452	0.0000	1.0000
JC69I	12	-1644.7239	3361.2551	0.0000	1.0000
JC69IG	13	-1643.3181	3364.4274	0.0000	1.0000
JC69	11	-1658.2379	3382.2990	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/chr13_2812.phylip = 1 - 397;
 Partition Dummy = 1:591-tomodel-phylip/chr13_2812.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/chr13_2812.phylip = 1 - 397;
 Partition Dummy = 1:591-tomodel-phylip/chr13_2812.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/chr13_2812.phylip = 1 - 397;
 Partition Dummy = 1:591-tomodel-phylip/chr13_2812.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

