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

Input data from file "183-phylip/chr5_11304.phylip" (29 taxa, 161 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	-520.2221	1218.6383	0.6980	0.6980
K2PI	57	-521.4456	1221.0854	0.2053	0.9034
K2PIG	58	-519.9675	1223.0331	0.0775	0.9809
K2P	56	-526.5260	1226.4366	0.0141	0.9950
HKYG	60	-518.6476	1230.4953	0.0019	0.9969
SYMG	61	-516.6184	1231.6408	0.0010	0.9979
HKYI	60	-519.8839	1232.9678	0.0005	0.9985
JC69G	56	-529.8469	1233.0785	0.0005	0.9990
SYMI	61	-517.5223	1233.4487	0.0004	0.9994
HKYIG	61	-518.4036	1235.2113	0.0002	0.9996
JC69I	56	-531.0388	1235.4622	0.0002	0.9998
SYMIG	62	-516.4274	1236.5692	0.0001	0.9998
JC69IG	57	-529.6167	1237.4275	0.0001	0.9999
HKY	59	-524.8201	1237.7391	0.0000	0.9999
SYM	60	-522.6033	1238.4066	0.0000	1.0000
JC69	55	-536.0086	1240.6839	0.0000	1.0000
F81G	59	-527.4548	1243.0086	0.0000	1.0000
F81I	59	-528.7089	1245.5169	0.0000	1.0000
GTRG	64	-515.8911	1246.4488	0.0000	1.0000
F81IG	60	-527.1991	1247.5981	0.0000	1.0000
GTRI	64	-516.8335	1248.3336	0.0000	1.0000
F81	58	-533.8068	1250.7117	0.0000	1.0000
GTRIG	65	-515.6917	1251.6991	0.0000	1.0000
GTR	63	-521.8777	1252.8894	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	57	-520.2221	1154.4441	0.2796	0.2796
SYMG	61	-516.6184	1155.2368	0.1881	0.4678
K2PIG	58	-519.9675	1155.9350	0.1327	0.6005
SYMIG	62	-516.4274	1156.8549	0.0838	0.6843
K2PI	57	-521.4456	1156.8912	0.0823	0.7665
SYMI	61	-517.5223	1157.0447	0.0762	0.8427
HKYG	60	-518.6476	1157.2953	0.0672	0.9099
HKYIG	61	-518.4036	1158.8073	0.0316	0.9415
HKYI	60	-519.8839	1159.7678	0.0195	0.9610
GTRG	64	-515.8911	1159.7822	0.0194	0.9804
GTRIG	65	-515.6917	1161.3834	0.0087	0.9891
GTRI	64	-516.8335	1161.6669	0.0076	0.9967
K2P	56	-526.5260	1165.0520	0.0014	0.9981
SYM	60	-522.6033	1165.2066	0.0013	0.9993
HKY	59	-524.8201	1167.6401	0.0004	0.9997
GTR	63	-521.8777	1169.7554	0.0001	0.9999
JC69G	56	-529.8469	1171.6939	0.0001	0.9999
F81G	59	-527.4548	1172.9096	0.0000	0.9999
JC69IG	57	-529.6167	1173.2333	0.0000	1.0000
JC69I	56	-531.0388	1174.0776	0.0000	1.0000
F81IG	60	-527.1991	1174.3981	0.0000	1.0000
F81I	59	-528.7089	1175.4179	0.0000	1.0000
JC69	55	-536.0086	1182.0173	0.0000	1.0000
F81	58	-533.8068	1183.6136	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-520.2221	1330.0842	0.7014	0.7014
K2PI	57	-521.4456	1332.5312	0.2063	0.9077
K2PIG	58	-519.9675	1334.6565	0.0713	0.9790
K2P	56	-526.5260	1337.6106	0.0163	0.9953
HKYG	60	-518.6476	1342.1796	0.0017	0.9969
SYMG	61	-516.6184	1343.2024	0.0010	0.9979
JC69G	56	-529.8469	1344.2525	0.0006	0.9985
HKYI	60	-519.8839	1344.6520	0.0005	0.9990
SYMI	61	-517.5223	1345.0103	0.0004	0.9994
JC69I	56	-531.0388	1346.6362	0.0002	0.9996
HKYIG	61	-518.4036	1346.7730	0.0002	0.9997
SYMIG	62	-516.4274	1347.9020	0.0001	0.9998
JC69IG	57	-529.6167	1348.8734	0.0001	0.9999
HKY	59	-524.8201	1349.4430	0.0000	0.9999
SYM	60	-522.6033	1350.0909	0.0000	1.0000
JC69	55	-536.0086	1351.4945	0.0000	1.0000
F81G	59	-527.4548	1354.7124	0.0000	1.0000
GTRG	64	-515.8911	1356.9921	0.0000	1.0000
F81I	59	-528.7089	1357.2207	0.0000	1.0000
GTRI	64	-516.8335	1358.8768	0.0000	1.0000
F81IG	60	-527.1991	1359.2824	0.0000	1.0000
GTRIG	65	-515.6917	1361.6746	0.0000	1.0000
F81	58	-533.8068	1362.3351	0.0000	1.0000
GTR	63	-521.8777	1363.8838	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/chr5_11304.phylip = 1 - 161;
 Partition Dummy = 1:183-phylip/chr5_11304.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/chr5_11304.phylip = 1 - 161;
 Partition Dummy = 1:183-phylip/chr5_11304.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/chr5_11304.phylip = 1 - 161;
 Partition Dummy = 1:183-phylip/chr5_11304.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

