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

Input data from file "183-phylip/chr1_8696.phylip" (29 taxa, 222 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	-621.1742	1404.0520	0.6712	0.6712
HKYG	60	-620.7919	1407.0496	0.1499	0.8211
HKYI	60	-620.8005	1407.0669	0.1486	0.9697
HKYIG	61	-620.7920	1410.8589	0.0223	0.9921
GTR	63	-618.4733	1413.9845	0.0047	0.9967
GTRI	64	-618.1365	1417.2666	0.0009	0.9976
GTRG	64	-618.1786	1417.3508	0.0009	0.9985
K2P	56	-633.3505	1417.3919	0.0009	0.9994
K2PG	57	-632.7973	1419.9117	0.0002	0.9996
K2PI	57	-632.8141	1419.9452	0.0002	0.9998
GTRIG	65	-618.1386	1421.2773	0.0001	1.0000
K2PIG	58	-632.7988	1423.5852	0.0000	1.0000
SYM	60	-631.9609	1429.3877	0.0000	1.0000
SYMI	61	-631.4132	1432.1013	0.0000	1.0000
SYMG	61	-631.4236	1432.1222	0.0000	1.0000
SYMIG	62	-631.4199	1435.9719	0.0000	1.0000
F81	58	-640.1139	1438.2156	0.0000	1.0000
F81G	59	-639.7032	1441.1101	0.0000	1.0000
F81I	59	-639.7165	1441.1367	0.0000	1.0000
F81IG	60	-639.7046	1444.8750	0.0000	1.0000
JC69	55	-651.2129	1449.5342	0.0000	1.0000
JC69G	56	-650.6778	1452.0466	0.0000	1.0000
JC69I	56	-650.7040	1452.0989	0.0000	1.0000
JC69IG	57	-650.6813	1455.6797	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKY	59	-621.1742	1360.3483	0.3482	0.3482
HKYG	60	-620.7919	1361.5837	0.1877	0.5359
HKYI	60	-620.8005	1361.6011	0.1861	0.7221
GTR	63	-618.4733	1362.9465	0.0950	0.8170
HKYIG	61	-620.7920	1363.5839	0.0691	0.8861
GTRI	64	-618.1365	1364.2730	0.0489	0.9350
GTRG	64	-618.1786	1364.3571	0.0469	0.9819
GTRIG	65	-618.1386	1366.2773	0.0180	0.9999
K2P	56	-633.3505	1378.7010	0.0000	0.9999
K2PG	57	-632.7973	1379.5947	0.0000	1.0000
K2PI	57	-632.8141	1379.6282	0.0000	1.0000
K2PIG	58	-632.7988	1381.5975	0.0000	1.0000
SYM	60	-631.9609	1383.9218	0.0000	1.0000
SYMI	61	-631.4132	1384.8263	0.0000	1.0000
SYMG	61	-631.4236	1384.8472	0.0000	1.0000
SYMIG	62	-631.4199	1386.8398	0.0000	1.0000
F81	58	-640.1139	1396.2279	0.0000	1.0000
F81G	59	-639.7032	1397.4064	0.0000	1.0000
F81I	59	-639.7165	1397.4330	0.0000	1.0000
F81IG	60	-639.7046	1399.4092	0.0000	1.0000
JC69	55	-651.2129	1412.4257	0.0000	1.0000
JC69G	56	-650.6778	1413.3557	0.0000	1.0000
JC69I	56	-650.7040	1413.4080	0.0000	1.0000
JC69IG	57	-650.6813	1415.3626	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	59	-621.1742	1561.1063	0.8170	0.8170
HKYG	60	-620.7919	1565.7444	0.0804	0.8974
HKYI	60	-620.8005	1565.7617	0.0797	0.9771
K2P	56	-633.3505	1569.2509	0.0139	0.9910
HKYIG	61	-620.7920	1571.1472	0.0054	0.9964
K2PG	57	-632.7973	1573.5473	0.0016	0.9980
K2PI	57	-632.8141	1573.5808	0.0016	0.9996
GTR	63	-618.4733	1577.3152	0.0002	0.9998
K2PIG	58	-632.7988	1578.9528	0.0001	1.0000
GTRI	64	-618.1365	1582.0444	0.0000	1.0000
GTRG	64	-618.1786	1582.1285	0.0000	1.0000
GTRIG	65	-618.1386	1587.4513	0.0000	1.0000
SYM	60	-631.9609	1588.0825	0.0000	1.0000
SYMI	61	-631.4132	1592.3896	0.0000	1.0000
SYMG	61	-631.4236	1592.4105	0.0000	1.0000
F81	58	-640.1139	1593.5832	0.0000	1.0000
SYMIG	62	-631.4199	1597.8058	0.0000	1.0000
F81G	59	-639.7032	1598.1644	0.0000	1.0000
F81I	59	-639.7165	1598.1910	0.0000	1.0000
JC69	55	-651.2129	1599.5730	0.0000	1.0000
F81IG	60	-639.7046	1603.5698	0.0000	1.0000
JC69G	56	-650.6778	1603.9056	0.0000	1.0000
JC69I	56	-650.7040	1603.9579	0.0000	1.0000
JC69IG	57	-650.6813	1609.3152	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/chr1_8696.phylip = 1 - 222;
 Partition Dummy = 1:183-phylip/chr1_8696.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/chr1_8696.phylip = 1 - 222;
 Partition Dummy = 1:183-phylip/chr1_8696.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/chr1_8696.phylip = 1 - 222;
 Partition Dummy = 1:183-phylip/chr1_8696.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

