
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:52 02/17/2011

Input data from file "183-phylip/chr4_15987.phylip" (29 taxa, 555 characters)

Minimum AIC  model: HKYIG
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	60	-1428.5344	2991.8866	0.3863	0.3863
HKYIG	61	-1427.3407	2992.0242	0.3606	0.7468
HKYI	60	-1429.6375	2994.0929	0.1282	0.8750
GTRIG	65	-1424.1183	2995.7826	0.0551	0.9301
GTRG	64	-1425.4586	2995.8969	0.0520	0.9821
GTRI	64	-1426.5352	2998.0501	0.0177	0.9998
K2PG	57	-1440.5847	3008.4732	0.0001	0.9999
K2PIG	58	-1439.4054	3008.6093	0.0001	1.0000
K2PI	57	-1441.7603	3010.8244	0.0000	1.0000
SYMG	61	-1439.6885	3016.7199	0.0000	1.0000
SYMIG	62	-1438.4692	3016.8164	0.0000	1.0000
SYMI	61	-1440.8727	3019.0882	0.0000	1.0000
F81G	59	-1446.5192	3025.3414	0.0000	1.0000
F81IG	60	-1445.4270	3025.6718	0.0000	1.0000
F81I	59	-1447.6051	3027.5133	0.0000	1.0000
HKY	59	-1451.0439	3034.3909	0.0000	1.0000
JC69G	56	-1456.4254	3037.6700	0.0000	1.0000
JC69IG	57	-1455.3126	3037.9290	0.0000	1.0000
GTR	63	-1448.7001	3039.8238	0.0000	1.0000
JC69I	56	-1457.5519	3039.9231	0.0000	1.0000
K2P	56	-1462.6438	3050.1068	0.0000	1.0000
SYM	60	-1461.5034	3057.8247	0.0000	1.0000
F81	58	-1467.7277	3065.2538	0.0000	1.0000
JC69	55	-1477.6659	3077.6765	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-1427.3407	2976.6814	0.3338	0.3338
HKYG	60	-1428.5344	2977.0688	0.2751	0.6089
GTRIG	65	-1424.1183	2978.2366	0.1534	0.7623
GTRG	64	-1425.4586	2978.9173	0.1092	0.8715
HKYI	60	-1429.6375	2979.2751	0.0913	0.9627
GTRI	64	-1426.5352	2981.0705	0.0372	0.9999
K2PIG	58	-1439.4054	2994.8109	0.0000	1.0000
K2PG	57	-1440.5847	2995.1694	0.0000	1.0000
K2PI	57	-1441.7603	2997.5206	0.0000	1.0000
SYMIG	62	-1438.4692	3000.9384	0.0000	1.0000
SYMG	61	-1439.6885	3001.3771	0.0000	1.0000
SYMI	61	-1440.8727	3003.7454	0.0000	1.0000
F81IG	60	-1445.4270	3010.8540	0.0000	1.0000
F81G	59	-1446.5192	3011.0384	0.0000	1.0000
F81I	59	-1447.6051	3013.2103	0.0000	1.0000
HKY	59	-1451.0439	3020.0879	0.0000	1.0000
GTR	63	-1448.7001	3023.4002	0.0000	1.0000
JC69IG	57	-1455.3126	3024.6252	0.0000	1.0000
JC69G	56	-1456.4254	3024.8507	0.0000	1.0000
JC69I	56	-1457.5519	3027.1038	0.0000	1.0000
K2P	56	-1462.6438	3037.2875	0.0000	1.0000
SYM	60	-1461.5034	3043.0069	0.0000	1.0000
F81	58	-1467.7277	3051.4554	0.0000	1.0000
JC69	55	-1477.6659	3065.3318	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1428.5344	3236.2068	0.6319	0.6319
HKYI	60	-1429.6375	3238.4131	0.2097	0.8416
HKYIG	61	-1427.3407	3240.1384	0.0885	0.9301
K2PG	57	-1440.5847	3241.3506	0.0483	0.9784
K2PI	57	-1441.7603	3243.7018	0.0149	0.9933
K2PIG	58	-1439.4054	3245.3110	0.0067	0.9999
GTRG	64	-1425.4586	3255.3312	0.0000	1.0000
GTRI	64	-1426.5352	3257.4845	0.0000	1.0000
GTRIG	65	-1424.1183	3258.9695	0.0000	1.0000
SYMG	61	-1439.6885	3264.8341	0.0000	1.0000
F81G	59	-1446.5192	3265.8575	0.0000	1.0000
JC69G	56	-1456.4254	3266.7129	0.0000	1.0000
SYMI	61	-1440.8727	3267.2025	0.0000	1.0000
F81I	59	-1447.6051	3268.0294	0.0000	1.0000
SYMIG	62	-1438.4692	3268.7144	0.0000	1.0000
JC69I	56	-1457.5519	3268.9660	0.0000	1.0000
F81IG	60	-1445.4270	3269.9920	0.0000	1.0000
JC69IG	57	-1455.3126	3270.8064	0.0000	1.0000
HKY	59	-1451.0439	3274.9070	0.0000	1.0000
K2P	56	-1462.6438	3279.1497	0.0000	1.0000
GTR	63	-1448.7001	3295.4952	0.0000	1.0000
F81	58	-1467.7277	3301.9555	0.0000	1.0000
SYM	60	-1461.5034	3302.1449	0.0000	1.0000
JC69	55	-1477.6659	3302.8750	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr4_15987.phylip = 1 - 555;
 Partition Dummy = 1:183-phylip/chr4_15987.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr4_15987.phylip = 1 - 555;
 Partition Dummy = 1:183-phylip/chr4_15987.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 183-phylip/chr4_15987.phylip = 1 - 555;
 Partition Dummy = 1:183-phylip/chr4_15987.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

