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

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

Minimum AIC  model: GTRI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	60	-519.7716	1229.2574	0.3649	0.3649
HKYG	60	-520.0584	1229.8310	0.2740	0.6389
K2PI	57	-527.8732	1230.9687	0.1551	0.7940
K2PG	57	-528.0201	1231.2624	0.1339	0.9279
HKYIG	61	-519.8158	1234.3624	0.0284	0.9564
HKY	59	-525.2924	1235.3772	0.0171	0.9735
K2PIG	58	-528.0109	1235.9845	0.0126	0.9861
K2P	56	-534.0109	1238.5906	0.0034	0.9895
GTRI	64	-514.2345	1238.8453	0.0030	0.9926
SYMI	61	-522.1942	1239.1191	0.0026	0.9952
SYMG	61	-522.6091	1239.9489	0.0017	0.9969
GTRG	64	-514.7918	1239.9599	0.0017	0.9987
F81I	59	-529.2926	1243.3776	0.0003	0.9990
F81G	59	-529.6659	1244.1243	0.0002	0.9992
JC69I	56	-537.0167	1244.6022	0.0002	0.9994
JC69G	56	-537.1213	1244.8114	0.0002	0.9995
SYMIG	62	-522.6084	1245.0615	0.0001	0.9997
GTR	63	-520.0313	1245.1215	0.0001	0.9998
GTRIG	65	-514.7910	1245.3819	0.0001	0.9999
SYM	60	-528.9511	1247.6164	0.0000	0.9999
F81IG	60	-529.4056	1248.5255	0.0000	1.0000
F81	58	-534.5709	1249.1045	0.0000	1.0000
JC69IG	57	-537.1200	1249.4622	0.0000	1.0000
JC69	55	-542.4800	1250.9600	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	64	-514.2345	1156.4690	0.4420	0.4420
GTRG	64	-514.7918	1157.5836	0.2531	0.6951
HKYI	60	-519.7716	1159.5432	0.0950	0.7901
GTRIG	65	-514.7910	1159.5819	0.0932	0.8833
HKYG	60	-520.0584	1160.1167	0.0713	0.9547
HKYIG	61	-519.8158	1161.6316	0.0334	0.9881
GTR	63	-520.0313	1166.0627	0.0036	0.9917
SYMI	61	-522.1942	1166.3883	0.0031	0.9948
SYMG	61	-522.6091	1167.2181	0.0020	0.9969
HKY	59	-525.2924	1168.5847	0.0010	0.9979
SYMIG	62	-522.6084	1169.2169	0.0008	0.9987
K2PI	57	-527.8732	1169.7465	0.0006	0.9993
K2PG	57	-528.0201	1170.0402	0.0005	0.9998
K2PIG	58	-528.0109	1172.0218	0.0002	0.9999
F81I	59	-529.2926	1176.5852	0.0000	1.0000
F81G	59	-529.6659	1177.3319	0.0000	1.0000
SYM	60	-528.9511	1177.9022	0.0000	1.0000
F81IG	60	-529.4056	1178.8112	0.0000	1.0000
K2P	56	-534.0109	1180.0218	0.0000	1.0000
F81	58	-534.5709	1185.1419	0.0000	1.0000
JC69I	56	-537.0167	1186.0334	0.0000	1.0000
JC69G	56	-537.1213	1186.2426	0.0000	1.0000
JC69IG	57	-537.1200	1188.2400	0.0000	1.0000
JC69	55	-542.4800	1194.9600	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-519.7716	1346.2624	0.3148	0.3148
HKYG	60	-520.0584	1346.8360	0.2363	0.5511
K2PI	57	-527.8732	1347.1298	0.2040	0.7551
K2PG	57	-528.0201	1347.4235	0.1762	0.9313
HKYIG	61	-519.8158	1351.4629	0.0234	0.9546
HKY	59	-525.2924	1352.1920	0.0162	0.9709
K2PIG	58	-528.0109	1352.5171	0.0138	0.9847
K2P	56	-534.0109	1354.2931	0.0057	0.9903
GTRI	64	-514.2345	1355.6362	0.0029	0.9932
SYMI	61	-522.1942	1356.2196	0.0022	0.9954
GTRG	64	-514.7918	1356.7508	0.0017	0.9971
SYMG	61	-522.6091	1357.0494	0.0014	0.9985
F81I	59	-529.2926	1360.1924	0.0003	0.9988
JC69I	56	-537.0167	1360.3048	0.0003	0.9991
JC69G	56	-537.1213	1360.5139	0.0003	0.9993
F81G	59	-529.6659	1360.9392	0.0002	0.9995
GTRIG	65	-514.7910	1361.8611	0.0001	0.9997
GTR	63	-520.0313	1362.1179	0.0001	0.9998
SYMIG	62	-522.6084	1362.1601	0.0001	0.9999
SYM	60	-528.9511	1364.6214	0.0000	0.9999
F81IG	60	-529.4056	1365.5304	0.0000	0.9999
JC69IG	57	-537.1200	1365.6233	0.0000	1.0000
F81	58	-534.5709	1365.6372	0.0000	1.0000
JC69	55	-542.4800	1366.1194	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_17712.phylip = 1 - 166;
 Partition Dummy = 1:183-phylip/chr3_17712.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_17712.phylip = 1 - 166;
 Partition Dummy = 1:183-phylip/chr3_17712.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_17712.phylip = 1 - 166;
 Partition Dummy = 1:183-phylip/chr3_17712.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

