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

Input data from file "591-tomodel-phylip/chr3_19948.phylip" (7 taxa, 490 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	-1478.2106	2989.5713	0.3722	0.3722
HKYI	16	-1478.3673	2989.8846	0.3182	0.6904
HKYIG	17	-1478.2243	2991.7452	0.1255	0.8159
GTRG	20	-1475.3270	2992.4451	0.0885	0.9044
GTRI	20	-1475.6224	2993.0358	0.0658	0.9702
GTRIG	21	-1475.3298	2994.6340	0.0296	0.9998
HKY	15	-1487.1617	3005.3361	0.0001	1.0000
GTR	19	-1484.2547	3008.1264	0.0000	1.0000
K2PG	13	-1493.9138	3014.5923	0.0000	1.0000
K2PI	13	-1494.3724	3015.5096	0.0000	1.0000
K2PIG	14	-1493.9331	3016.7503	0.0000	1.0000
SYMG	17	-1491.6256	3018.5477	0.0000	1.0000
SYMI	17	-1491.9279	3019.1523	0.0000	1.0000
SYMIG	18	-1491.6336	3020.7194	0.0000	1.0000
K2P	12	-1502.1655	3028.9852	0.0000	1.0000
SYM	16	-1499.4444	3032.0388	0.0000	1.0000
F81G	15	-1505.2259	3041.4645	0.0000	1.0000
F81I	15	-1505.4119	3041.8364	0.0000	1.0000
F81IG	16	-1505.2377	3043.6254	0.0000	1.0000
F81	14	-1512.9422	3054.7686	0.0000	1.0000
JC69I	12	-1517.7178	3060.0897	0.0000	1.0000
JC69G	12	-1517.8881	3060.4302	0.0000	1.0000
JC69IG	13	-1517.5355	3061.8357	0.0000	1.0000
JC69	11	-1525.4769	3073.5061	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1478.2106	2988.4212	0.3437	0.3437
HKYI	16	-1478.3673	2988.7345	0.2938	0.6375
HKYIG	17	-1478.2243	2990.4486	0.1247	0.7622
GTRG	20	-1475.3270	2990.6540	0.1125	0.8748
GTRI	20	-1475.6224	2991.2447	0.0838	0.9585
GTRIG	21	-1475.3298	2992.6596	0.0413	0.9998
HKY	15	-1487.1617	3004.3234	0.0001	1.0000
GTR	19	-1484.2547	3006.5094	0.0000	1.0000
K2PG	13	-1493.9138	3013.8276	0.0000	1.0000
K2PI	13	-1494.3724	3014.7449	0.0000	1.0000
K2PIG	14	-1493.9331	3015.8661	0.0000	1.0000
SYMG	17	-1491.6256	3017.2511	0.0000	1.0000
SYMI	17	-1491.9279	3017.8557	0.0000	1.0000
SYMIG	18	-1491.6336	3019.2672	0.0000	1.0000
K2P	12	-1502.1655	3028.3311	0.0000	1.0000
SYM	16	-1499.4444	3030.8887	0.0000	1.0000
F81G	15	-1505.2259	3040.4518	0.0000	1.0000
F81I	15	-1505.4119	3040.8238	0.0000	1.0000
F81IG	16	-1505.2377	3042.4753	0.0000	1.0000
F81	14	-1512.9422	3053.8844	0.0000	1.0000
JC69I	12	-1517.7178	3059.4356	0.0000	1.0000
JC69G	12	-1517.8881	3059.7761	0.0000	1.0000
JC69IG	13	-1517.5355	3061.0710	0.0000	1.0000
JC69	11	-1525.4769	3072.9538	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1478.2106	3055.5317	0.5248	0.5248
HKYI	16	-1478.3673	3055.8450	0.4487	0.9736
HKYIG	17	-1478.2243	3061.7535	0.0234	0.9970
HKY	15	-1487.1617	3067.2395	0.0015	0.9985
K2PG	13	-1493.9138	3068.3549	0.0009	0.9993
K2PI	13	-1494.3724	3069.2721	0.0005	0.9999
GTRG	20	-1475.3270	3074.5421	0.0000	0.9999
K2PIG	14	-1493.9331	3074.5878	0.0000	1.0000
GTRI	20	-1475.6224	3075.1328	0.0000	1.0000
K2P	12	-1502.1655	3078.6640	0.0000	1.0000
GTRIG	21	-1475.3298	3080.7421	0.0000	1.0000
GTR	19	-1484.2547	3086.2031	0.0000	1.0000
SYMG	17	-1491.6256	3088.5560	0.0000	1.0000
SYMI	17	-1491.9279	3089.1606	0.0000	1.0000
SYMIG	18	-1491.6336	3094.7665	0.0000	1.0000
SYM	16	-1499.4444	3097.9992	0.0000	1.0000
F81G	15	-1505.2259	3103.3679	0.0000	1.0000
F81I	15	-1505.4119	3103.7398	0.0000	1.0000
F81IG	16	-1505.2377	3109.5858	0.0000	1.0000
JC69I	12	-1517.7178	3109.7685	0.0000	1.0000
JC69G	12	-1517.8881	3110.1090	0.0000	1.0000
F81	14	-1512.9422	3112.6061	0.0000	1.0000
JC69IG	13	-1517.5355	3115.5983	0.0000	1.0000
JC69	11	-1525.4769	3119.0922	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/chr3_19948.phylip = 1 - 490;
 Partition Dummy = 1:591-tomodel-phylip/chr3_19948.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/chr3_19948.phylip = 1 - 490;
 Partition Dummy = 1:591-tomodel-phylip/chr3_19948.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/chr3_19948.phylip = 1 - 490;
 Partition Dummy = 1:591-tomodel-phylip/chr3_19948.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

