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

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

Minimum AIC  model: HKY
Minimum AICc model: HKY
Minimum BIC  model: HKY


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
HKY	59	-653.8837	1457.9493	0.5152	0.5152
HKYI	60	-653.6024	1460.6294	0.1349	0.6501
HKYG	60	-653.6699	1460.7644	0.1261	0.7762
F81	58	-656.9460	1460.8603	0.1202	0.8964
F81I	59	-656.6638	1463.5094	0.0320	0.9284
F81G	59	-656.7250	1463.6318	0.0301	0.9584
HKYIG	61	-653.6042	1463.9056	0.0262	0.9846
F81IG	60	-656.6658	1466.7564	0.0063	0.9909
GTR	63	-652.2925	1467.9184	0.0035	0.9945
K2P	56	-664.4382	1469.5042	0.0016	0.9961
GTRI	64	-651.8446	1470.3868	0.0010	0.9971
GTRG	64	-651.9282	1470.5540	0.0009	0.9980
JC69	55	-667.2021	1471.9041	0.0005	0.9985
K2PI	57	-664.1363	1472.0564	0.0004	0.9990
K2PG	57	-664.2018	1472.1874	0.0004	0.9994
GTRIG	65	-651.8476	1473.7886	0.0002	0.9996
JC69I	56	-666.8953	1474.4185	0.0001	0.9997
JC69G	56	-666.9609	1474.5495	0.0001	0.9998
K2PIG	58	-664.1384	1475.2451	0.0001	0.9999
JC69IG	57	-666.8977	1477.5792	0.0000	1.0000
SYM	60	-662.1157	1477.6561	0.0000	1.0000
SYMI	61	-661.5929	1479.8830	0.0000	1.0000
SYMG	61	-661.6681	1480.0335	0.0000	1.0000
SYMIG	62	-661.5946	1483.1893	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKY	59	-653.8837	1425.7674	0.3843	0.3843
HKYI	60	-653.6024	1427.2048	0.1873	0.5716
HKYG	60	-653.6699	1427.3397	0.1751	0.7467
HKYIG	61	-653.6042	1429.2084	0.0688	0.8155
F81	58	-656.9460	1429.8920	0.0489	0.8643
GTR	63	-652.2925	1430.5851	0.0346	0.8989
F81I	59	-656.6638	1431.3275	0.0238	0.9227
F81G	59	-656.7250	1431.4500	0.0224	0.9452
GTRI	64	-651.8446	1431.6891	0.0199	0.9651
GTRG	64	-651.9282	1431.8563	0.0183	0.9834
F81IG	60	-656.6658	1433.3317	0.0088	0.9921
GTRIG	65	-651.8476	1433.6952	0.0073	0.9994
K2P	56	-664.4382	1440.8764	0.0002	0.9996
K2PI	57	-664.1363	1442.2726	0.0001	0.9997
K2PG	57	-664.2018	1442.4036	0.0001	0.9998
SYM	60	-662.1157	1444.2314	0.0000	0.9998
K2PIG	58	-664.1384	1444.2768	0.0000	0.9999
JC69	55	-667.2021	1444.4041	0.0000	0.9999
SYMI	61	-661.5929	1445.1858	0.0000	0.9999
SYMG	61	-661.6681	1445.3362	0.0000	1.0000
JC69I	56	-666.8953	1445.7907	0.0000	1.0000
JC69G	56	-666.9609	1445.9217	0.0000	1.0000
SYMIG	62	-661.5946	1447.1893	0.0000	1.0000
JC69IG	57	-666.8977	1447.7954	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	59	-653.8837	1640.2200	0.4243	0.4243
F81	58	-656.9460	1640.7098	0.3322	0.7565
JC69	55	-667.2021	1644.3176	0.0547	0.8112
K2P	56	-664.4382	1644.4247	0.0518	0.8631
HKYI	60	-653.6024	1645.2922	0.0336	0.8967
HKYG	60	-653.6699	1645.4271	0.0314	0.9281
F81I	59	-656.6638	1645.7801	0.0263	0.9544
F81G	59	-656.7250	1645.9025	0.0248	0.9792
JC69I	56	-666.8953	1649.3389	0.0044	0.9836
K2PI	57	-664.1363	1649.4556	0.0042	0.9878
JC69G	56	-666.9609	1649.4699	0.0042	0.9919
K2PG	57	-664.2018	1649.5866	0.0039	0.9959
HKYIG	61	-653.6042	1650.9306	0.0020	0.9979
F81IG	60	-656.6658	1651.4191	0.0016	0.9994
JC69IG	57	-666.8977	1654.9784	0.0003	0.9997
K2PIG	58	-664.1384	1655.0946	0.0002	1.0000
GTR	63	-652.2925	1659.5768	0.0000	1.0000
SYM	60	-662.1157	1662.3188	0.0000	1.0000
GTRI	64	-651.8446	1664.3157	0.0000	1.0000
GTRG	64	-651.9282	1664.4829	0.0000	1.0000
SYMI	61	-661.5929	1666.9080	0.0000	1.0000
SYMG	61	-661.6681	1667.0584	0.0000	1.0000
GTRIG	65	-651.8476	1669.9565	0.0000	1.0000
SYMIG	62	-661.5946	1672.5462	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_22628.phylip = 1 - 280;
 Partition Dummy = 1:183-phylip/chr2_22628.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_22628.phylip = 1 - 280;
 Partition Dummy = 1:183-phylip/chr2_22628.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_22628.phylip = 1 - 280;
 Partition Dummy = 1:183-phylip/chr2_22628.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

