
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:51 04/04/2011

Input data from file "591-tomodel-phylip/chr8_9158.phylip" (7 taxa, 389 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-986.1839	2014.6504	0.3307	0.3307
GTRG	20	-986.6456	2015.5738	0.2084	0.5391
GTRIG	21	-986.1975	2016.9127	0.1067	0.6458
HKYI	16	-991.9013	2017.2650	0.0895	0.7353
HKYG	16	-992.2029	2017.8682	0.0662	0.8014
GTR	19	-989.1711	2018.4018	0.0507	0.8521
SYMI	17	-991.5646	2018.7788	0.0420	0.8941
HKYIG	17	-991.9105	2019.4705	0.0297	0.9238
SYMG	17	-991.9794	2019.6084	0.0277	0.9515
SYMIG	18	-991.5768	2021.0022	0.0138	0.9653
HKY	15	-995.0677	2021.4223	0.0112	0.9765
SYM	16	-994.1166	2021.6956	0.0098	0.9863
K2PI	13	-997.9819	2022.9344	0.0053	0.9915
K2PG	13	-998.2675	2023.5057	0.0039	0.9955
K2PIG	14	-997.9903	2025.1036	0.0018	0.9972
K2P	12	-1000.8216	2026.4730	0.0009	0.9981
F81I	15	-997.8650	2027.0168	0.0007	0.9988
F81G	15	-998.0725	2027.4318	0.0006	0.9994
F81IG	16	-997.8712	2029.2048	0.0002	0.9996
F81	14	-1000.0804	2029.2839	0.0002	0.9998
JC69I	12	-1003.3452	2031.5202	0.0001	0.9999
JC69G	12	-1003.5474	2031.9246	0.0001	0.9999
JC69	11	-1005.4179	2033.5360	0.0000	1.0000
JC69IG	13	-1003.3515	2033.6736	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-986.1839	2012.3678	0.3627	0.3627
GTRG	20	-986.6456	2013.2912	0.2285	0.5912
GTRIG	21	-986.1975	2014.3950	0.1316	0.7228
HKYI	16	-991.9013	2015.8026	0.0651	0.7879
GTR	19	-989.1711	2016.3422	0.0497	0.8376
HKYG	16	-992.2029	2016.4059	0.0482	0.8858
SYMI	17	-991.5646	2017.1292	0.0335	0.9193
HKYIG	17	-991.9105	2017.8209	0.0237	0.9431
SYMG	17	-991.9794	2017.9588	0.0222	0.9652
SYMIG	18	-991.5768	2019.1536	0.0122	0.9774
HKY	15	-995.0677	2020.1355	0.0075	0.9849
SYM	16	-994.1166	2020.2332	0.0071	0.9920
K2PI	13	-997.9819	2021.9638	0.0030	0.9950
K2PG	13	-998.2675	2022.5350	0.0022	0.9972
K2PIG	14	-997.9903	2023.9806	0.0011	0.9983
K2P	12	-1000.8216	2025.6432	0.0005	0.9988
F81I	15	-997.8650	2025.7300	0.0005	0.9992
F81G	15	-998.0725	2026.1450	0.0004	0.9996
F81IG	16	-997.8712	2027.7424	0.0002	0.9998
F81	14	-1000.0804	2028.1609	0.0001	0.9999
JC69I	12	-1003.3452	2030.6904	0.0000	0.9999
JC69G	12	-1003.5474	2031.0948	0.0000	1.0000
JC69IG	13	-1003.3515	2032.7030	0.0000	1.0000
JC69	11	-1005.4179	2032.8358	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	12	-1000.8216	2073.2062	0.3265	0.3265
K2PI	13	-997.9819	2073.4903	0.2832	0.6097
K2PG	13	-998.2675	2074.0616	0.2129	0.8225
JC69	11	-1005.4179	2076.4352	0.0650	0.8875
JC69I	12	-1003.3452	2078.2534	0.0262	0.9137
JC69G	12	-1003.5474	2078.6578	0.0214	0.9350
HKYI	16	-991.9013	2079.2199	0.0161	0.9512
K2PIG	14	-997.9903	2079.4707	0.0142	0.9654
HKY	15	-995.0677	2079.5892	0.0134	0.9788
HKYG	16	-992.2029	2079.8231	0.0119	0.9908
SYM	16	-994.1166	2083.6505	0.0018	0.9925
F81	14	-1000.0804	2083.6510	0.0018	0.9943
JC69IG	13	-1003.3515	2084.2295	0.0013	0.9956
SYMI	17	-991.5646	2084.5100	0.0011	0.9968
F81I	15	-997.8650	2085.1837	0.0008	0.9976
HKYIG	17	-991.9105	2085.2018	0.0008	0.9984
SYMG	17	-991.9794	2085.3396	0.0008	0.9991
F81G	15	-998.0725	2085.5987	0.0007	0.9998
SYMIG	18	-991.5768	2090.4980	0.0001	0.9999
F81IG	16	-997.8712	2091.1597	0.0000	0.9999
GTRI	20	-986.1839	2091.6393	0.0000	0.9999
GTR	19	-989.1711	2091.6502	0.0000	1.0000
GTRG	20	-986.6456	2092.5628	0.0000	1.0000
GTRIG	21	-986.1975	2097.6301	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 591-tomodel-phylip/chr8_9158.phylip = 1 - 389;
 Partition Dummy = 1:591-tomodel-phylip/chr8_9158.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 (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr8_9158.phylip = 1 - 389;
 Partition Dummy = 1:591-tomodel-phylip/chr8_9158.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 BIC model (K2P)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr8_9158.phylip = 1 - 389;
 Partition Dummy = 1:591-tomodel-phylip/chr8_9158.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

