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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-639.8136	1440.5210	0.3253	0.3253
HKY	59	-641.7031	1440.7396	0.2916	0.6169
HKYI	60	-640.4604	1441.8146	0.1704	0.7873
K2P	56	-648.5516	1443.9884	0.0575	0.8447
K2PG	57	-646.8456	1444.0208	0.0565	0.9013
HKYIG	61	-639.8133	1444.1210	0.0538	0.9551
K2PI	57	-647.4885	1445.3067	0.0297	0.9848
K2PIG	58	-646.8503	1447.5128	0.0099	0.9946
GTR	63	-639.5453	1450.9088	0.0018	0.9964
GTRG	64	-637.6848	1450.9124	0.0018	0.9983
GTRI	64	-638.3251	1452.1930	0.0010	0.9992
GTRIG	65	-637.6839	1454.6781	0.0003	0.9995
SYM	60	-647.1976	1455.2890	0.0002	0.9997
SYMG	61	-645.4523	1455.3990	0.0002	0.9999
SYMI	61	-646.0933	1456.6810	0.0001	1.0000
SYMIG	62	-645.4520	1459.0395	0.0000	1.0000
F81	58	-665.3979	1484.6080	0.0000	1.0000
F81G	59	-663.8261	1484.9856	0.0000	1.0000
F81I	59	-664.4365	1486.2063	0.0000	1.0000
F81IG	60	-663.8297	1488.5532	0.0000	1.0000
JC69	55	-672.6697	1488.8176	0.0000	1.0000
JC69G	56	-671.0871	1489.0594	0.0000	1.0000
JC69I	56	-671.7132	1490.3116	0.0000	1.0000
JC69IG	57	-671.0867	1492.5031	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-639.8136	1399.6271	0.3697	0.3697
HKYI	60	-640.4604	1400.9208	0.1936	0.5633
HKY	59	-641.7031	1401.4062	0.1519	0.7152
HKYIG	61	-639.8133	1401.6266	0.1360	0.8513
GTRG	64	-637.6848	1403.3696	0.0569	0.9082
GTRI	64	-638.3251	1404.6502	0.0300	0.9382
GTR	63	-639.5453	1405.0906	0.0241	0.9623
GTRIG	65	-637.6839	1405.3678	0.0210	0.9832
K2PG	57	-646.8456	1407.6911	0.0066	0.9898
K2PI	57	-647.4885	1408.9770	0.0034	0.9932
K2P	56	-648.5516	1409.1032	0.0032	0.9965
K2PIG	58	-646.8503	1409.7006	0.0024	0.9989
SYMG	61	-645.4523	1412.9046	0.0005	0.9993
SYMI	61	-646.0933	1414.1866	0.0003	0.9996
SYM	60	-647.1976	1414.3952	0.0002	0.9998
SYMIG	62	-645.4520	1414.9039	0.0002	1.0000
F81G	59	-663.8261	1445.6522	0.0000	1.0000
F81	58	-665.3979	1446.7959	0.0000	1.0000
F81I	59	-664.4365	1446.8730	0.0000	1.0000
F81IG	60	-663.8297	1447.6593	0.0000	1.0000
JC69G	56	-671.0871	1454.1742	0.0000	1.0000
JC69	55	-672.6697	1455.3393	0.0000	1.0000
JC69I	56	-671.7132	1455.4264	0.0000	1.0000
JC69IG	57	-671.0867	1456.1734	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	56	-648.5516	1604.0190	0.5023	0.5023
K2PG	57	-646.8456	1606.0875	0.1786	0.6809
HKY	59	-641.7031	1606.7639	0.1273	0.8082
K2PI	57	-647.4885	1607.3734	0.0939	0.9021
HKYG	60	-639.8136	1608.4655	0.0544	0.9565
HKYI	60	-640.4604	1609.7591	0.0285	0.9849
K2PIG	58	-646.8503	1611.5777	0.0115	0.9964
HKYIG	61	-639.8133	1613.9456	0.0035	0.9999
SYM	60	-647.1976	1623.2335	0.0000	0.9999
GTR	63	-639.5453	1624.3709	0.0000	1.0000
SYMG	61	-645.4523	1625.2236	0.0000	1.0000
GTRG	64	-637.6848	1626.1305	0.0000	1.0000
SYMI	61	-646.0933	1626.5056	0.0000	1.0000
GTRI	64	-638.3251	1627.4111	0.0000	1.0000
SYMIG	62	-645.4520	1630.7036	0.0000	1.0000
GTRIG	65	-637.6839	1631.6093	0.0000	1.0000
JC69	55	-672.6697	1646.7745	0.0000	1.0000
F81	58	-665.3979	1648.6729	0.0000	1.0000
JC69G	56	-671.0871	1649.0900	0.0000	1.0000
JC69I	56	-671.7132	1650.3422	0.0000	1.0000
F81G	59	-663.8261	1651.0099	0.0000	1.0000
F81I	59	-664.4365	1652.2307	0.0000	1.0000
JC69IG	57	-671.0867	1654.5698	0.0000	1.0000
F81IG	60	-663.8297	1656.4976	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 183-phylip/chr3_9700.phylip = 1 - 240;
 Partition Dummy = 1:183-phylip/chr3_9700.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 183-phylip/chr3_9700.phylip = 1 - 240;
 Partition Dummy = 1:183-phylip/chr3_9700.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 (K2P)]

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

