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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-1419.5273	2992.8930	0.7436	0.7436
GTRIG	65	-1419.2486	2995.2263	0.2316	0.9751
HKYG	60	-1429.1101	3000.6742	0.0152	0.9903
HKYIG	61	-1428.8511	3002.9761	0.0048	0.9951
GTRI	64	-1424.6301	3003.0986	0.0045	0.9996
HKYI	60	-1432.8271	3008.1082	0.0004	1.0000
GTR	63	-1444.1676	3039.3011	0.0000	1.0000
K2PG	57	-1454.0668	3042.2309	0.0000	1.0000
K2PIG	58	-1453.7901	3044.4460	0.0000	1.0000
HKY	59	-1453.4188	3046.4890	0.0000	1.0000
K2PI	57	-1457.0121	3048.1215	0.0000	1.0000
SYMG	61	-1452.1197	3049.5131	0.0000	1.0000
SYMIG	62	-1451.8480	3051.8070	0.0000	1.0000
SYMI	61	-1455.2736	3055.8211	0.0000	1.0000
F81G	59	-1459.7171	3059.0857	0.0000	1.0000
F81IG	60	-1459.4949	3061.4437	0.0000	1.0000
F81I	59	-1462.8206	3065.2925	0.0000	1.0000
JC69G	56	-1479.4940	3090.3335	0.0000	1.0000
K2P	56	-1480.4341	3092.2136	0.0000	1.0000
JC69IG	57	-1479.2537	3092.6046	0.0000	1.0000
JC69I	56	-1482.3453	3096.0361	0.0000	1.0000
SYM	60	-1477.4690	3097.3920	0.0000	1.0000
F81	58	-1484.2068	3105.2795	0.0000	1.0000
JC69	55	-1505.2391	3139.0886	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1419.5273	2967.0545	0.6676	0.6676
GTRIG	65	-1419.2486	2968.4973	0.3245	0.9922
GTRI	64	-1424.6301	2977.2601	0.0041	0.9962
HKYG	60	-1429.1101	2978.2202	0.0025	0.9987
HKYIG	61	-1428.8511	2979.7022	0.0012	0.9999
HKYI	60	-1432.8271	2985.6542	0.0001	1.0000
GTR	63	-1444.1676	3014.3352	0.0000	1.0000
K2PG	57	-1454.0668	3022.1337	0.0000	1.0000
K2PIG	58	-1453.7901	3023.5801	0.0000	1.0000
HKY	59	-1453.4188	3024.8376	0.0000	1.0000
SYMG	61	-1452.1197	3026.2393	0.0000	1.0000
SYMIG	62	-1451.8480	3027.6959	0.0000	1.0000
K2PI	57	-1457.0121	3028.0243	0.0000	1.0000
SYMI	61	-1455.2736	3032.5473	0.0000	1.0000
F81G	59	-1459.7171	3037.4343	0.0000	1.0000
F81IG	60	-1459.4949	3038.9898	0.0000	1.0000
F81I	59	-1462.8206	3043.6412	0.0000	1.0000
JC69G	56	-1479.4940	3070.9880	0.0000	1.0000
JC69IG	57	-1479.2537	3072.5074	0.0000	1.0000
K2P	56	-1480.4341	3072.8681	0.0000	1.0000
SYM	60	-1477.4690	3074.9381	0.0000	1.0000
JC69I	56	-1482.3453	3076.6906	0.0000	1.0000
F81	58	-1484.2068	3084.4137	0.0000	1.0000
JC69	55	-1505.2391	3120.4783	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1429.1101	3215.7257	0.8374	0.8374
GTRG	64	-1419.5273	3220.3937	0.0812	0.9186
HKYIG	61	-1428.8511	3221.1661	0.0552	0.9737
HKYI	60	-1432.8271	3223.1597	0.0204	0.9941
GTRIG	65	-1419.2486	3225.7949	0.0055	0.9995
GTRI	64	-1424.6301	3230.5993	0.0005	1.0000
K2PG	57	-1454.0668	3247.7639	0.0000	1.0000
K2PIG	58	-1453.7901	3253.1687	0.0000	1.0000
K2PI	57	-1457.0121	3253.6545	0.0000	1.0000
HKY	59	-1453.4188	3258.3847	0.0000	1.0000
GTR	63	-1444.1676	3263.7159	0.0000	1.0000
SYMG	61	-1452.1197	3267.7032	0.0000	1.0000
F81G	59	-1459.7171	3270.9813	0.0000	1.0000
SYMIG	62	-1451.8480	3273.1183	0.0000	1.0000
SYMI	61	-1455.2736	3274.0112	0.0000	1.0000
F81IG	60	-1459.4949	3276.4952	0.0000	1.0000
F81I	59	-1462.8206	3277.1882	0.0000	1.0000
JC69G	56	-1479.4940	3292.6598	0.0000	1.0000
K2P	56	-1480.4341	3294.5399	0.0000	1.0000
JC69IG	57	-1479.2537	3298.1376	0.0000	1.0000
JC69I	56	-1482.3453	3298.3624	0.0000	1.0000
SYM	60	-1477.4690	3312.4435	0.0000	1.0000
F81	58	-1484.2068	3314.0023	0.0000	1.0000
JC69	55	-1505.2391	3338.1916	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr1_32424.phylip = 1 - 387;
 Partition Dummy = 1:183-phylip/chr1_32424.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr1_32424.phylip = 1 - 387;
 Partition Dummy = 1:183-phylip/chr1_32424.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

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

