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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-1595.4924	3232.3288	0.4797	0.4797
GTRG	20	-1595.8953	3233.1347	0.3206	0.8003
GTRIG	21	-1595.5035	3234.4878	0.1630	0.9632
GTR	19	-1599.3803	3237.9747	0.0285	0.9917
HKYI	16	-1604.5573	3241.9795	0.0038	0.9956
HKYG	16	-1604.8937	3242.6522	0.0027	0.9983
HKYIG	17	-1604.5672	3244.1088	0.0013	0.9997
HKY	15	-1608.0767	3246.9153	0.0003	1.0000
SYMI	17	-1622.4813	3279.9371	0.0000	1.0000
K2PI	13	-1626.7869	3280.1498	0.0000	1.0000
SYMG	17	-1622.8080	3280.5906	0.0000	1.0000
K2PG	13	-1627.0701	3280.7161	0.0000	1.0000
SYMIG	18	-1622.4912	3282.0734	0.0000	1.0000
K2PIG	14	-1626.7956	3282.2568	0.0000	1.0000
SYM	16	-1625.9820	3284.8288	0.0000	1.0000
K2P	12	-1630.4524	3285.3978	0.0000	1.0000
F81I	15	-1637.4423	3305.6465	0.0000	1.0000
F81G	15	-1637.8205	3306.4028	0.0000	1.0000
F81IG	16	-1637.4543	3307.7735	0.0000	1.0000
F81	14	-1640.4607	3309.5870	0.0000	1.0000
JC69I	12	-1656.5468	3337.5865	0.0000	1.0000
JC69G	12	-1656.8884	3338.2697	0.0000	1.0000
JC69IG	13	-1656.5575	3339.6910	0.0000	1.0000
JC69	11	-1659.6892	3341.7948	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1595.4924	3230.9848	0.4759	0.4759
GTRG	20	-1595.8953	3231.7907	0.3180	0.7939
GTRIG	21	-1595.5035	3233.0070	0.1731	0.9670
GTR	19	-1599.3803	3236.7607	0.0265	0.9935
HKYI	16	-1604.5573	3241.1146	0.0030	0.9965
HKYG	16	-1604.8937	3241.7874	0.0021	0.9987
HKYIG	17	-1604.5672	3243.1343	0.0011	0.9998
HKY	15	-1608.0767	3246.1534	0.0002	1.0000
SYMI	17	-1622.4813	3278.9626	0.0000	1.0000
K2PI	13	-1626.7869	3279.5738	0.0000	1.0000
SYMG	17	-1622.8080	3279.6161	0.0000	1.0000
K2PG	13	-1627.0701	3280.1402	0.0000	1.0000
SYMIG	18	-1622.4912	3280.9825	0.0000	1.0000
K2PIG	14	-1626.7956	3281.5912	0.0000	1.0000
SYM	16	-1625.9820	3283.9639	0.0000	1.0000
K2P	12	-1630.4524	3284.9049	0.0000	1.0000
F81I	15	-1637.4423	3304.8846	0.0000	1.0000
F81G	15	-1637.8205	3305.6409	0.0000	1.0000
F81IG	16	-1637.4543	3306.9087	0.0000	1.0000
F81	14	-1640.4607	3308.9214	0.0000	1.0000
JC69I	12	-1656.5468	3337.0936	0.0000	1.0000
JC69G	12	-1656.8884	3337.7768	0.0000	1.0000
JC69IG	13	-1656.5575	3339.1151	0.0000	1.0000
JC69	11	-1659.6892	3341.3784	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1604.5573	3312.6474	0.3918	0.3918
HKY	15	-1608.0767	3313.2154	0.2949	0.6868
HKYG	16	-1604.8937	3313.3202	0.2799	0.9666
HKYIG	17	-1604.5672	3319.1379	0.0153	0.9819
GTRI	20	-1595.4924	3320.4008	0.0081	0.9900
GTRG	20	-1595.8953	3321.2067	0.0054	0.9955
GTR	19	-1599.3803	3321.7059	0.0042	0.9997
GTRIG	21	-1595.5035	3326.8938	0.0003	1.0000
K2PI	13	-1626.7869	3337.6942	0.0000	1.0000
K2PG	13	-1627.0701	3338.2606	0.0000	1.0000
K2P	12	-1630.4524	3338.5545	0.0000	1.0000
K2PIG	14	-1626.7956	3344.1824	0.0000	1.0000
SYMI	17	-1622.4813	3354.9662	0.0000	1.0000
SYM	16	-1625.9820	3355.4967	0.0000	1.0000
SYMG	17	-1622.8080	3355.6197	0.0000	1.0000
SYMIG	18	-1622.4912	3361.4569	0.0000	1.0000
F81	14	-1640.4607	3371.5126	0.0000	1.0000
F81I	15	-1637.4423	3371.9466	0.0000	1.0000
F81G	15	-1637.8205	3372.7029	0.0000	1.0000
F81IG	16	-1637.4543	3378.4415	0.0000	1.0000
JC69	11	-1659.6892	3390.5572	0.0000	1.0000
JC69I	12	-1656.5468	3390.7432	0.0000	1.0000
JC69G	12	-1656.8884	3391.4264	0.0000	1.0000
JC69IG	13	-1656.5575	3397.2355	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_8801.phylip = 1 - 646;
 Partition Dummy = 1:591-tomodel-phylip/chr6_8801.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_8801.phylip = 1 - 646;
 Partition Dummy = 1:591-tomodel-phylip/chr6_8801.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

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

