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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1649.9084	3341.6468	0.4591	0.4591
GTRI	20	-1650.7118	3343.2536	0.2056	0.6647
GTRIG	21	-1649.9084	3343.8343	0.1538	0.8185
HKYG	16	-1656.2217	3345.6184	0.0630	0.8815
GTR	19	-1653.3465	3346.3451	0.0438	0.9254
HKYI	16	-1657.1187	3347.4123	0.0257	0.9511
HKYIG	17	-1656.2218	3347.7684	0.0215	0.9726
SYMG	17	-1657.3385	3350.0016	0.0070	0.9796
HKY	15	-1659.5334	3350.1013	0.0067	0.9863
SYMI	17	-1658.0171	3351.3589	0.0036	0.9899
K2PG	13	-1662.3340	3351.4491	0.0034	0.9933
SYMIG	18	-1657.3385	3352.1608	0.0024	0.9957
K2PI	13	-1663.1451	3353.0714	0.0015	0.9972
K2PIG	14	-1662.3341	3353.5714	0.0012	0.9984
SYM	16	-1660.2590	3353.6930	0.0011	0.9995
K2P	12	-1665.3335	3355.3351	0.0005	1.0000
F81G	15	-1700.6399	3432.3143	0.0000	1.0000
F81I	15	-1701.3108	3433.6561	0.0000	1.0000
F81IG	16	-1700.6400	3434.4550	0.0000	1.0000
JC69G	12	-1704.9729	3434.6138	0.0000	1.0000
F81	14	-1703.0631	3435.0295	0.0000	1.0000
JC69I	12	-1705.6118	3435.8916	0.0000	1.0000
JC69IG	13	-1704.9730	3436.7270	0.0000	1.0000
JC69	11	-1707.3597	3437.2835	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1649.9084	3339.8167	0.4715	0.4715
GTRI	20	-1650.7118	3341.4236	0.2111	0.6826
GTRIG	21	-1649.9084	3341.8168	0.1734	0.8560
HKYG	16	-1656.2217	3344.4435	0.0466	0.9026
GTR	19	-1653.3465	3344.6929	0.0412	0.9438
HKYI	16	-1657.1187	3346.2374	0.0190	0.9628
HKYIG	17	-1656.2218	3346.4437	0.0172	0.9800
SYMG	17	-1657.3385	3348.6769	0.0056	0.9856
HKY	15	-1659.5334	3349.0668	0.0046	0.9902
SYMI	17	-1658.0171	3350.0342	0.0028	0.9931
K2PG	13	-1662.3340	3350.6680	0.0021	0.9951
SYMIG	18	-1657.3385	3350.6771	0.0021	0.9972
K2PI	13	-1663.1451	3352.2903	0.0009	0.9981
SYM	16	-1660.2590	3352.5181	0.0008	0.9990
K2PIG	14	-1662.3341	3352.6682	0.0008	0.9997
K2P	12	-1665.3335	3354.6670	0.0003	1.0000
F81G	15	-1700.6399	3431.2798	0.0000	1.0000
F81I	15	-1701.3108	3432.6216	0.0000	1.0000
F81IG	16	-1700.6400	3433.2800	0.0000	1.0000
JC69G	12	-1704.9729	3433.9457	0.0000	1.0000
F81	14	-1703.0631	3434.1263	0.0000	1.0000
JC69I	12	-1705.6118	3435.2235	0.0000	1.0000
JC69IG	13	-1704.9730	3435.9459	0.0000	1.0000
JC69	11	-1707.3597	3436.7194	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	12	-1665.3335	3404.7525	0.4072	0.4072
K2PG	13	-1662.3340	3404.9272	0.3731	0.7804
K2PI	13	-1663.1451	3406.5495	0.1658	0.9462
K2PIG	14	-1662.3341	3411.1012	0.0170	0.9632
HKYG	16	-1656.2217	3411.2241	0.0160	0.9792
HKY	15	-1659.5334	3411.6736	0.0128	0.9920
HKYI	16	-1657.1187	3413.0180	0.0065	0.9985
HKYIG	17	-1656.2218	3417.3980	0.0007	0.9993
SYM	16	-1660.2590	3419.2987	0.0003	0.9995
SYMG	17	-1657.3385	3419.6313	0.0002	0.9998
SYMI	17	-1658.0171	3420.9885	0.0001	0.9999
GTRG	20	-1649.9084	3423.2924	0.0000	0.9999
GTR	19	-1653.3465	3423.9948	0.0000	1.0000
GTRI	20	-1650.7118	3424.8993	0.0000	1.0000
SYMIG	18	-1657.3385	3425.8052	0.0000	1.0000
GTRIG	21	-1649.9084	3429.4663	0.0000	1.0000
JC69	11	-1707.3597	3482.6310	0.0000	1.0000
JC69G	12	-1704.9729	3484.0311	0.0000	1.0000
JC69I	12	-1705.6118	3485.3089	0.0000	1.0000
JC69IG	13	-1704.9730	3490.2051	0.0000	1.0000
F81	14	-1703.0631	3492.5593	0.0000	1.0000
F81G	15	-1700.6399	3493.8866	0.0000	1.0000
F81I	15	-1701.3108	3495.2284	0.0000	1.0000
F81IG	16	-1700.6400	3500.0606	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/chr3_21547.phylip = 1 - 480;
 Partition Dummy = 1:591-tomodel-phylip/chr3_21547.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 591-tomodel-phylip/chr3_21547.phylip = 1 - 480;
 Partition Dummy = 1:591-tomodel-phylip/chr3_21547.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 (K2P)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_21547.phylip = 1 - 480;
 Partition Dummy = 1:591-tomodel-phylip/chr3_21547.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

