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

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

Minimum AIC  model: K2PG
Minimum AICc model: K2PG
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
K2PG	57	-583.9207	1334.3176	0.5465	0.5465
K2PI	57	-584.4685	1335.4132	0.3160	0.8625
K2PIG	58	-583.7420	1338.2361	0.0770	0.9395
JC69G	56	-589.4843	1341.2363	0.0172	0.9567
HKYG	60	-581.0520	1341.6161	0.0142	0.9709
JC69I	56	-590.0173	1342.3024	0.0101	0.9810
HKYI	60	-581.5068	1342.5257	0.0090	0.9900
K2P	56	-591.0078	1344.2833	0.0037	0.9938
JC69IG	57	-589.3061	1345.0884	0.0025	0.9963
HKYIG	61	-580.8909	1345.7818	0.0018	0.9980
SYMG	61	-581.9588	1347.9176	0.0006	0.9986
SYMI	61	-582.3664	1348.7327	0.0004	0.9990
F81G	59	-586.9658	1349.0285	0.0003	0.9994
F81I	59	-587.4798	1350.0563	0.0002	0.9996
HKY	59	-588.1076	1351.3119	0.0001	0.9997
JC69	55	-596.6141	1351.3532	0.0001	0.9998
SYMIG	62	-581.8089	1352.1798	0.0001	0.9999
F81IG	60	-586.7972	1353.1067	0.0000	0.9999
GTRG	64	-578.1084	1354.1327	0.0000	1.0000
GTRI	64	-578.4734	1354.8627	0.0000	1.0000
SYM	60	-588.6070	1356.7261	0.0000	1.0000
F81	58	-593.8823	1358.5166	0.0000	1.0000
GTRIG	65	-577.9556	1358.6230	0.0000	1.0000
GTR	63	-584.6885	1362.5770	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	57	-583.9207	1281.8414	0.2087	0.2087
HKYG	60	-581.0520	1282.1039	0.1831	0.3918
K2PI	57	-584.4685	1282.9370	0.1207	0.5125
HKYI	60	-581.5068	1283.0135	0.1162	0.6287
K2PIG	58	-583.7420	1283.4841	0.0918	0.7205
HKYIG	61	-580.8909	1283.7818	0.0791	0.7996
GTRG	64	-578.1084	1284.2167	0.0637	0.8632
GTRI	64	-578.4734	1284.9468	0.0442	0.9074
GTRIG	65	-577.9556	1285.9111	0.0273	0.9347
SYMG	61	-581.9588	1285.9176	0.0272	0.9619
SYMI	61	-582.3664	1286.7327	0.0181	0.9800
SYMIG	62	-581.8089	1287.6178	0.0116	0.9916
JC69G	56	-589.4843	1290.9686	0.0022	0.9938
F81G	59	-586.9658	1291.9317	0.0013	0.9951
JC69I	56	-590.0173	1292.0347	0.0013	0.9964
JC69IG	57	-589.3061	1292.6122	0.0010	0.9974
F81I	59	-587.4798	1292.9596	0.0008	0.9982
F81IG	60	-586.7972	1293.5945	0.0006	0.9988
K2P	56	-591.0078	1294.0155	0.0005	0.9992
HKY	59	-588.1076	1294.2152	0.0004	0.9997
GTR	63	-584.6885	1295.3770	0.0002	0.9999
SYM	60	-588.6070	1297.2139	0.0001	1.0000
JC69	55	-596.6141	1303.2282	0.0000	1.0000
F81	58	-593.8823	1303.7646	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-583.9207	1465.0927	0.5623	0.5623
K2PI	57	-584.4685	1466.1883	0.3251	0.8874
K2PIG	58	-583.7420	1469.9503	0.0496	0.9370
JC69G	56	-589.4843	1471.0050	0.0292	0.9662
JC69I	56	-590.0173	1472.0711	0.0172	0.9834
K2P	56	-591.0078	1474.0519	0.0064	0.9898
HKYG	60	-581.0520	1475.0001	0.0040	0.9937
JC69IG	57	-589.3061	1475.8635	0.0026	0.9963
HKYI	60	-581.5068	1475.9097	0.0025	0.9988
HKYIG	61	-580.8909	1479.8929	0.0003	0.9992
JC69	55	-596.6141	1480.0496	0.0003	0.9995
F81G	59	-586.9658	1481.6129	0.0001	0.9996
SYMG	61	-581.9588	1482.0286	0.0001	0.9998
F81I	59	-587.4798	1482.6408	0.0001	0.9998
SYMI	61	-582.3664	1482.8438	0.0001	0.9999
HKY	59	-588.1076	1483.8964	0.0000	1.0000
F81IG	60	-586.7972	1486.4906	0.0000	1.0000
SYMIG	62	-581.8089	1486.9439	0.0000	1.0000
GTRG	64	-578.1084	1489.9726	0.0000	1.0000
SYM	60	-588.6070	1490.1101	0.0000	1.0000
F81	58	-593.8823	1490.2308	0.0000	1.0000
GTRI	64	-578.4734	1490.7027	0.0000	1.0000
GTRIG	65	-577.9556	1494.8820	0.0000	1.0000
GTR	63	-584.6885	1497.9179	0.0000	1.0000



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

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

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


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

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


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

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

