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

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

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
Minimum BIC  model: K2PG


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	-663.6973	1361.6426	0.2380	0.2380
HKYI	16	-663.8768	1362.0015	0.1989	0.4369
GTRG	20	-659.6002	1362.7299	0.1382	0.5751
GTRI	20	-659.6004	1362.7303	0.1382	0.7133
HKYIG	17	-663.6733	1363.8860	0.0775	0.7908
GTRIG	21	-659.5316	1364.9619	0.0453	0.8361
K2PG	13	-668.8627	1365.2111	0.0400	0.8760
K2PI	13	-668.9855	1365.4567	0.0353	0.9114
SYMI	17	-664.5554	1365.6502	0.0321	0.9435
SYMG	17	-664.6004	1365.7403	0.0307	0.9741
K2PIG	14	-668.8331	1367.3875	0.0135	0.9876
SYMIG	18	-664.6585	1368.1670	0.0091	0.9967
GTR	19	-665.2354	1371.6507	0.0016	0.9983
HKY	15	-670.1239	1372.2230	0.0012	0.9995
SYM	16	-670.3426	1374.9332	0.0003	0.9998
K2P	12	-675.4297	1376.1276	0.0002	1.0000
F81I	15	-680.1566	1392.2885	0.0000	1.0000
F81G	15	-680.2453	1392.4660	0.0000	1.0000
F81IG	16	-680.4201	1395.0881	0.0000	1.0000
JC69I	12	-686.7777	1398.8236	0.0000	1.0000
JC69G	12	-686.8833	1399.0348	0.0000	1.0000
F81	14	-685.4901	1400.7015	0.0000	1.0000
JC69IG	13	-686.9653	1401.4163	0.0000	1.0000
JC69	11	-691.9871	1407.0430	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-659.6002	1359.2005	0.2009	0.2009
GTRI	20	-659.6004	1359.2009	0.2009	0.4019
HKYG	16	-663.6973	1359.3946	0.1824	0.5842
HKYI	16	-663.8768	1359.7535	0.1524	0.7366
GTRIG	21	-659.5316	1361.0631	0.0792	0.8158
HKYIG	17	-663.6733	1361.3466	0.0687	0.8845
SYMI	17	-664.5554	1363.1108	0.0284	0.9130
SYMG	17	-664.6004	1363.2008	0.0272	0.9401
K2PG	13	-668.8627	1363.7254	0.0209	0.9611
K2PI	13	-668.9855	1363.9710	0.0185	0.9796
SYMIG	18	-664.6585	1365.3170	0.0094	0.9890
K2PIG	14	-668.8331	1365.6662	0.0079	0.9969
GTR	19	-665.2354	1368.4708	0.0020	0.9989
HKY	15	-670.1239	1370.2477	0.0008	0.9997
SYM	16	-670.3426	1372.6852	0.0002	0.9999
K2P	12	-675.4297	1374.8593	0.0001	1.0000
F81I	15	-680.1566	1390.3132	0.0000	1.0000
F81G	15	-680.2453	1390.4907	0.0000	1.0000
F81IG	16	-680.4201	1392.8402	0.0000	1.0000
JC69I	12	-686.7777	1397.5553	0.0000	1.0000
JC69G	12	-686.8833	1397.7665	0.0000	1.0000
F81	14	-685.4901	1398.9802	0.0000	1.0000
JC69IG	13	-686.9653	1399.9306	0.0000	1.0000
JC69	11	-691.9871	1405.9742	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-668.8627	1409.9641	0.4868	0.4868
K2PI	13	-668.9855	1410.2098	0.4305	0.9173
K2PIG	14	-668.8331	1415.4618	0.0312	0.9485
HKYG	16	-663.6973	1416.3039	0.0204	0.9689
HKYI	16	-663.8768	1416.6628	0.0171	0.9860
K2P	12	-675.4297	1417.5413	0.0110	0.9970
HKYIG	17	-663.6733	1421.8126	0.0013	0.9983
SYMI	17	-664.5554	1423.5769	0.0005	0.9989
HKY	15	-670.1239	1423.6002	0.0005	0.9994
SYMG	17	-664.6004	1423.6669	0.0005	0.9999
SYMIG	18	-664.6585	1429.3399	0.0000	0.9999
SYM	16	-670.3426	1429.5945	0.0000	1.0000
GTRG	20	-659.6002	1430.3370	0.0000	1.0000
GTRI	20	-659.6004	1430.3374	0.0000	1.0000
GTRIG	21	-659.5316	1435.7565	0.0000	1.0000
GTR	19	-665.2354	1436.0506	0.0000	1.0000
JC69I	12	-686.7777	1440.2373	0.0000	1.0000
JC69G	12	-686.8833	1440.4485	0.0000	1.0000
F81I	15	-680.1566	1443.6656	0.0000	1.0000
F81G	15	-680.2453	1443.8431	0.0000	1.0000
JC69	11	-691.9871	1445.0993	0.0000	1.0000
JC69IG	13	-686.9653	1446.1693	0.0000	1.0000
F81	14	-685.4901	1448.7758	0.0000	1.0000
F81IG	16	-680.4201	1449.7494	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 591-tomodel-phylip/chr9_6320.phylip = 1 - 259;
 Partition Dummy = 1:591-tomodel-phylip/chr9_6320.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_6320.phylip = 1 - 259;
 Partition Dummy = 1:591-tomodel-phylip/chr9_6320.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 (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_6320.phylip = 1 - 259;
 Partition Dummy = 1:591-tomodel-phylip/chr9_6320.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

