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

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

Minimum AIC  model: K2P
Minimum AICc model: K2P
Minimum BIC  model: K2P


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
K2P	12	-365.9905	758.0892	0.2446	0.2446
K2PI	13	-365.2819	759.0400	0.1521	0.3967
K2PG	13	-365.3106	759.0974	0.1478	0.5445
SYM	16	-362.2778	760.3333	0.0797	0.6242
HKY	15	-363.8330	760.9763	0.0578	0.6819
K2PIG	14	-365.3116	761.5000	0.0445	0.7264
SYMI	17	-361.7148	761.7093	0.0400	0.7664
SYMG	17	-361.7896	761.8589	0.0371	0.8036
HKYI	16	-363.1234	762.0247	0.0342	0.8377
HKYG	16	-363.1420	762.0617	0.0336	0.8713
GTR	19	-359.6752	762.7405	0.0239	0.8952
JC69	11	-369.6478	763.0674	0.0203	0.9155
JC69I	12	-368.9710	764.0502	0.0124	0.9279
GTRI	20	-359.0342	764.0685	0.0123	0.9402
JC69G	12	-369.0168	764.1416	0.0119	0.9521
GTRG	20	-359.0869	764.1739	0.0117	0.9638
SYMIG	18	-361.7847	764.3864	0.0105	0.9743
HKYIG	17	-363.1428	764.5654	0.0096	0.9839
F81	14	-367.7053	766.2873	0.0041	0.9879
JC69IG	13	-369.0171	766.5104	0.0036	0.9916
GTRIG	21	-359.0863	766.8201	0.0031	0.9947
F81I	15	-367.0258	767.3619	0.0024	0.9971
F81G	15	-367.0652	767.4407	0.0023	0.9993
F81IG	16	-367.0657	769.9091	0.0007	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2P	12	-365.9905	755.9811	0.1362	0.1362
SYM	16	-362.2778	756.5555	0.1022	0.2384
K2PI	13	-365.2819	756.5638	0.1018	0.3402
K2PG	13	-365.3106	756.6212	0.0989	0.4391
GTR	19	-359.6752	757.3504	0.0687	0.5077
SYMI	17	-361.7148	757.4296	0.0660	0.5737
SYMG	17	-361.7896	757.5792	0.0613	0.6350
HKY	15	-363.8330	757.6660	0.0587	0.6937
GTRI	20	-359.0342	758.0685	0.0480	0.7416
GTRG	20	-359.0869	758.1739	0.0455	0.7871
HKYI	16	-363.1234	758.2469	0.0439	0.8310
HKYG	16	-363.1420	758.2839	0.0431	0.8740
K2PIG	14	-365.3116	758.6233	0.0363	0.9104
SYMIG	18	-361.7847	759.5695	0.0226	0.9330
GTRIG	21	-359.0863	760.1726	0.0167	0.9498
HKYIG	17	-363.1428	760.2857	0.0158	0.9656
JC69	11	-369.6478	761.2956	0.0096	0.9752
JC69I	12	-368.9710	761.9421	0.0069	0.9821
JC69G	12	-369.0168	762.0335	0.0066	0.9887
F81	14	-367.7053	763.4105	0.0033	0.9920
JC69IG	13	-369.0171	764.0342	0.0024	0.9944
F81I	15	-367.0258	764.0516	0.0024	0.9968
F81G	15	-367.0652	764.1303	0.0023	0.9991
F81IG	16	-367.0657	766.1313	0.0009	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	12	-365.9905	792.9580	0.5658	0.5658
JC69	11	-369.6478	795.1910	0.1852	0.7510
K2PI	13	-365.2819	796.6220	0.0906	0.8416
K2PG	13	-365.3106	796.6795	0.0880	0.9296
JC69I	12	-368.9710	798.9189	0.0287	0.9583
JC69G	12	-369.0168	799.0104	0.0274	0.9857
K2PIG	14	-365.3116	801.7629	0.0069	0.9927
HKY	15	-363.8330	803.8870	0.0024	0.9951
JC69IG	13	-369.0171	804.0925	0.0022	0.9972
SYM	16	-362.2778	805.8580	0.0009	0.9981
F81	14	-367.7053	806.5502	0.0006	0.9988
HKYI	16	-363.1234	807.5494	0.0004	0.9991
HKYG	16	-363.1420	807.5864	0.0004	0.9995
SYMI	17	-361.7148	809.8135	0.0001	0.9996
SYMG	17	-361.7896	809.9631	0.0001	0.9998
F81I	15	-367.0258	810.2726	0.0001	0.9999
F81G	15	-367.0652	810.3514	0.0001	0.9999
HKYIG	17	-363.1428	812.6695	0.0000	1.0000
SYMIG	18	-361.7847	815.0348	0.0000	1.0000
F81IG	16	-367.0657	815.4338	0.0000	1.0000
GTR	19	-359.6752	815.8971	0.0000	1.0000
GTRI	20	-359.0342	819.6966	0.0000	1.0000
GTRG	20	-359.0869	819.8020	0.0000	1.0000
GTRIG	21	-359.0863	824.8821	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_8651.phylip = 1 - 161;
 Partition Dummy = 1:591-tomodel-phylip/chr2_8651.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_8651.phylip = 1 - 161;
 Partition Dummy = 1:591-tomodel-phylip/chr2_8651.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_8651.phylip = 1 - 161;
 Partition Dummy = 1:591-tomodel-phylip/chr2_8651.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


-------------------------------------------------------------
End of Output

