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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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	-471.6786	1110.5830	0.4209	0.4209
HKYI	60	-471.7134	1110.6526	0.4065	0.8275
HKYIG	61	-471.6118	1114.3405	0.0643	0.8918
HKY	59	-475.6662	1114.7171	0.0533	0.9450
K2PG	57	-480.6090	1117.0662	0.0165	0.9615
K2PI	57	-480.6125	1117.0730	0.0164	0.9779
GTRG	64	-467.7401	1118.5795	0.0077	0.9856
GTRI	64	-467.8769	1118.8532	0.0067	0.9924
K2PIG	58	-480.5817	1120.7557	0.0026	0.9950
GTR	63	-470.9314	1120.9155	0.0024	0.9974
K2P	56	-485.1118	1122.3746	0.0012	0.9985
GTRIG	65	-467.6587	1122.5174	0.0011	0.9996
SYMG	61	-478.1598	1127.4364	0.0001	0.9997
SYMI	61	-478.2000	1127.5169	0.0001	0.9998
F81G	59	-482.1644	1127.7135	0.0001	0.9999
F81I	59	-482.1908	1127.7662	0.0001	0.9999
SYM	60	-481.9286	1131.0830	0.0000	0.9999
SYMIG	62	-478.0959	1131.2506	0.0000	1.0000
F81IG	60	-482.1213	1131.4685	0.0000	1.0000
F81	58	-486.0769	1131.7462	0.0000	1.0000
JC69G	56	-490.2315	1132.6139	0.0000	1.0000
JC69I	56	-490.2436	1132.6381	0.0000	1.0000
JC69IG	57	-490.2054	1136.2588	0.0000	1.0000
JC69	55	-494.4323	1137.3646	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-471.6786	1063.3572	0.2147	0.2147
HKYI	60	-471.7134	1063.4268	0.2073	0.4220
GTRG	64	-467.7401	1063.4802	0.2019	0.6239
GTRI	64	-467.8769	1063.7538	0.1761	0.7999
HKYIG	61	-471.6118	1065.2236	0.0844	0.8844
GTRIG	65	-467.6587	1065.3174	0.0806	0.9649
GTR	63	-470.9314	1067.8628	0.0226	0.9875
HKY	59	-475.6662	1069.3325	0.0108	0.9983
K2PG	57	-480.6090	1075.2181	0.0006	0.9989
K2PI	57	-480.6125	1075.2249	0.0006	0.9994
K2PIG	58	-480.5817	1077.1633	0.0002	0.9997
SYMG	61	-478.1598	1078.3195	0.0001	0.9998
SYMI	61	-478.2000	1078.4000	0.0001	0.9999
SYMIG	62	-478.0959	1080.1918	0.0000	0.9999
K2P	56	-485.1118	1082.2236	0.0000	1.0000
F81G	59	-482.1644	1082.3289	0.0000	1.0000
F81I	59	-482.1908	1082.3816	0.0000	1.0000
SYM	60	-481.9286	1083.8572	0.0000	1.0000
F81IG	60	-482.1213	1084.2427	0.0000	1.0000
F81	58	-486.0769	1088.1539	0.0000	1.0000
JC69G	56	-490.2315	1092.4630	0.0000	1.0000
JC69I	56	-490.2436	1092.4872	0.0000	1.0000
JC69IG	57	-490.2054	1094.4107	0.0000	1.0000
JC69	55	-494.4323	1098.8646	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-471.6786	1265.8739	0.3075	0.3075
HKYI	60	-471.7134	1265.9435	0.2970	0.6044
K2PG	57	-480.6090	1267.6089	0.1291	0.7336
K2PI	57	-480.6125	1267.6158	0.1287	0.8623
HKY	59	-475.6662	1268.4739	0.0838	0.9461
HKYIG	61	-471.6118	1271.1156	0.0224	0.9684
K2P	56	-485.1118	1271.2392	0.0210	0.9894
K2PIG	58	-480.5817	1272.9295	0.0090	0.9985
GTRG	64	-467.7401	1279.4980	0.0003	0.9988
GTRI	64	-467.8769	1279.7716	0.0003	0.9991
GTR	63	-470.9314	1280.5054	0.0002	0.9993
F81G	59	-482.1644	1281.4703	0.0001	0.9994
JC69G	56	-490.2315	1281.4786	0.0001	0.9996
JC69I	56	-490.2436	1281.5028	0.0001	0.9997
F81I	59	-482.1908	1281.5230	0.0001	0.9998
F81	58	-486.0769	1283.9200	0.0000	0.9999
SYMG	61	-478.1598	1284.2115	0.0000	0.9999
SYMI	61	-478.2000	1284.2920	0.0000	0.9999
JC69	55	-494.4323	1284.5049	0.0000	0.9999
GTRIG	65	-467.6587	1284.7105	0.0000	1.0000
SYM	60	-481.9286	1286.3739	0.0000	1.0000
F81IG	60	-482.1213	1286.7594	0.0000	1.0000
JC69IG	57	-490.2054	1286.8016	0.0000	1.0000
SYMIG	62	-478.0959	1289.4591	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/chr6_4410.phylip = 1 - 216;
 Partition Dummy = 1:183-phylip/chr6_4410.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/chr6_4410.phylip = 1 - 216;
 Partition Dummy = 1:183-phylip/chr6_4410.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 (HKYG)]

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


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

