
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:09 04/04/2011

Input data from file "683-tomodel-phylip/chr2_18686.phylip" (7 taxa, 216 characters)

Minimum AIC  model: K2PI
Minimum AICc model: K2P
Minimum BIC  model: K2P


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
K2P	12	-477.3141	980.1652	0.2656	0.2656
K2PI	13	-476.3075	980.4169	0.2342	0.4998
K2PG	13	-476.3106	980.4232	0.2334	0.7332
K2PIG	14	-476.3114	982.7123	0.0743	0.8075
SYMG	17	-473.7530	984.5969	0.0290	0.8365
HKY	15	-476.1123	984.6246	0.0286	0.8650
SYMI	17	-473.7781	984.6472	0.0282	0.8933
SYM	16	-474.9942	984.7220	0.0272	0.9205
HKYI	16	-475.0279	984.7895	0.0263	0.9468
HKYG	16	-475.0331	984.7999	0.0262	0.9730
SYMIG	18	-473.7469	986.9659	0.0089	0.9818
HKYIG	17	-475.0344	987.1597	0.0080	0.9899
GTRG	20	-472.8104	989.9284	0.0020	0.9919
GTRI	20	-472.8272	989.9620	0.0020	0.9939
GTR	19	-474.1170	990.1115	0.0018	0.9957
JC69	11	-483.8060	990.9061	0.0012	0.9969
JC69G	12	-482.9983	991.5335	0.0009	0.9978
JC69I	12	-483.0007	991.5383	0.0009	0.9987
GTRIG	21	-472.8053	992.3735	0.0006	0.9993
JC69IG	13	-482.9982	993.7984	0.0003	0.9996
F81	14	-482.6168	995.3232	0.0001	0.9998
F81I	15	-481.7608	995.9215	0.0001	0.9999
F81G	15	-481.7617	995.9233	0.0001	1.0000
F81IG	16	-481.7618	998.2573	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PI	13	-476.3075	978.6149	0.2088	0.2088
K2PG	13	-476.3106	978.6212	0.2082	0.4170
K2P	12	-477.3141	978.6282	0.2074	0.6245
K2PIG	14	-476.3114	980.6228	0.0765	0.7010
SYMG	17	-473.7530	981.5060	0.0492	0.7502
SYMI	17	-473.7781	981.5563	0.0480	0.7982
SYM	16	-474.9942	981.9883	0.0387	0.8368
HKYI	16	-475.0279	982.0558	0.0374	0.8742
HKYG	16	-475.0331	982.0662	0.0372	0.9114
HKY	15	-476.1123	982.2246	0.0344	0.9458
SYMIG	18	-473.7469	983.4938	0.0182	0.9640
HKYIG	17	-475.0344	984.0688	0.0137	0.9776
GTRG	20	-472.8104	985.6207	0.0063	0.9839
GTRI	20	-472.8272	985.6543	0.0062	0.9901
GTR	19	-474.1170	986.2340	0.0046	0.9947
GTRIG	21	-472.8053	987.6106	0.0023	0.9971
JC69	11	-483.8060	989.6120	0.0009	0.9979
JC69G	12	-482.9983	989.9966	0.0007	0.9986
JC69I	12	-483.0007	990.0014	0.0007	0.9993
JC69IG	13	-482.9982	991.9964	0.0003	0.9996
F81	14	-482.6168	993.2336	0.0001	0.9997
F81I	15	-481.7608	993.5215	0.0001	0.9998
F81G	15	-481.7617	993.5233	0.0001	1.0000
F81IG	16	-481.7618	995.5236	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	12	-477.3141	1019.1316	0.7064	0.7064
K2PI	13	-476.3075	1022.4935	0.1315	0.8380
K2PG	13	-476.3106	1022.4998	0.1311	0.9691
JC69	11	-483.8060	1026.7401	0.0157	0.9848
K2PIG	14	-476.3114	1027.8767	0.0089	0.9937
JC69G	12	-482.9983	1030.4999	0.0024	0.9961
JC69I	12	-483.0007	1030.5047	0.0024	0.9985
HKY	15	-476.1123	1032.8538	0.0007	0.9993
JC69IG	13	-482.9982	1035.8750	0.0002	0.9994
SYM	16	-474.9942	1035.9928	0.0002	0.9996
HKYI	16	-475.0279	1036.0603	0.0001	0.9997
HKYG	16	-475.0331	1036.0707	0.0001	0.9999
SYMG	17	-473.7530	1038.8857	0.0000	0.9999
SYMI	17	-473.7781	1038.9360	0.0000	1.0000
F81	14	-482.6168	1040.4875	0.0000	1.0000
HKYIG	17	-475.0344	1041.4486	0.0000	1.0000
F81I	15	-481.7608	1044.1507	0.0000	1.0000
F81G	15	-481.7617	1044.1525	0.0000	1.0000
SYMIG	18	-473.7469	1044.2488	0.0000	1.0000
F81IG	16	-481.7618	1049.5281	0.0000	1.0000
GTR	19	-474.1170	1050.3643	0.0000	1.0000
GTRG	20	-472.8104	1053.1263	0.0000	1.0000
GTRI	20	-472.8272	1053.1599	0.0000	1.0000
GTRIG	21	-472.8053	1058.4914	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_18686.phylip = 1 - 216;
 Partition Dummy = 1:683-tomodel-phylip/chr2_18686.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_18686.phylip = 1 - 216;
 Partition Dummy = 1:683-tomodel-phylip/chr2_18686.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_18686.phylip = 1 - 216;
 Partition Dummy = 1:683-tomodel-phylip/chr2_18686.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

