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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1016.9688	2076.2580	0.3245	0.3245
GTRI	20	-1017.4568	2077.2341	0.1992	0.5237
HKYI	16	-1022.1119	2077.7101	0.1570	0.6808
HKYG	16	-1022.1357	2077.7578	0.1533	0.8341
GTRIG	21	-1016.8834	2078.3263	0.1154	0.9495
HKYIG	17	-1022.1505	2079.9777	0.0505	1.0000
GTR	19	-1031.2889	2102.6714	0.0000	1.0000
HKY	15	-1035.7073	2102.7225	0.0000	1.0000
K2PI	13	-1039.0267	2105.0397	0.0000	1.0000
K2PG	13	-1039.1151	2105.2167	0.0000	1.0000
SYMI	17	-1035.1262	2105.9290	0.0000	1.0000
SYMG	17	-1035.3131	2106.3029	0.0000	1.0000
K2PIG	14	-1039.1151	2107.3715	0.0000	1.0000
SYMIG	18	-1035.3131	2108.5053	0.0000	1.0000
F81I	15	-1043.0545	2117.4169	0.0000	1.0000
F81G	15	-1043.1744	2117.6568	0.0000	1.0000
F81IG	16	-1043.1745	2119.8353	0.0000	1.0000
K2P	12	-1051.5425	2127.9282	0.0000	1.0000
SYM	16	-1047.4249	2128.3361	0.0000	1.0000
F81	14	-1054.5479	2138.2371	0.0000	1.0000
JC69I	12	-1059.3980	2143.6391	0.0000	1.0000
JC69G	12	-1059.4724	2143.7880	0.0000	1.0000
JC69IG	13	-1059.4984	2145.9832	0.0000	1.0000
JC69	11	-1070.2981	2163.3077	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1016.9688	2073.9376	0.3626	0.3626
GTRI	20	-1017.4568	2074.9136	0.2226	0.5852
GTRIG	21	-1016.8834	2075.7667	0.1453	0.7306
HKYI	16	-1022.1119	2076.2237	0.1156	0.8462
HKYG	16	-1022.1357	2076.2714	0.1129	0.9591
HKYIG	17	-1022.1505	2078.3010	0.0409	1.0000
GTR	19	-1031.2889	2100.5777	0.0000	1.0000
HKY	15	-1035.7073	2101.4146	0.0000	1.0000
K2PI	13	-1039.0267	2104.0533	0.0000	1.0000
K2PG	13	-1039.1151	2104.2302	0.0000	1.0000
SYMI	17	-1035.1262	2104.2523	0.0000	1.0000
SYMG	17	-1035.3131	2104.6262	0.0000	1.0000
K2PIG	14	-1039.1151	2106.2302	0.0000	1.0000
SYMIG	18	-1035.3131	2106.6262	0.0000	1.0000
F81I	15	-1043.0545	2116.1090	0.0000	1.0000
F81G	15	-1043.1744	2116.3489	0.0000	1.0000
F81IG	16	-1043.1745	2118.3489	0.0000	1.0000
SYM	16	-1047.4249	2126.8498	0.0000	1.0000
K2P	12	-1051.5425	2127.0849	0.0000	1.0000
F81	14	-1054.5479	2137.0958	0.0000	1.0000
JC69I	12	-1059.3980	2142.7959	0.0000	1.0000
JC69G	12	-1059.4724	2142.9447	0.0000	1.0000
JC69IG	13	-1059.4984	2144.9967	0.0000	1.0000
JC69	11	-1070.2981	2162.5961	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1022.1119	2139.3923	0.4930	0.4930
HKYG	16	-1022.1357	2139.4400	0.4814	0.9745
HKYIG	17	-1022.1505	2145.4176	0.0242	0.9987
GTRG	20	-1016.9688	2152.8983	0.0006	0.9993
GTRI	20	-1017.4568	2153.8743	0.0004	0.9996
K2PI	13	-1039.0267	2155.3778	0.0002	0.9998
K2PG	13	-1039.1151	2155.5547	0.0002	0.9999
GTRIG	21	-1016.8834	2158.6755	0.0000	1.0000
HKY	15	-1035.7073	2160.6351	0.0000	1.0000
K2PIG	14	-1039.1151	2161.5027	0.0000	1.0000
SYMI	17	-1035.1262	2171.3689	0.0000	1.0000
SYMG	17	-1035.3131	2171.7428	0.0000	1.0000
K2P	12	-1051.5425	2174.4614	0.0000	1.0000
F81I	15	-1043.0545	2175.3296	0.0000	1.0000
F81G	15	-1043.1744	2175.5694	0.0000	1.0000
GTR	19	-1031.2889	2175.5904	0.0000	1.0000
SYMIG	18	-1035.3131	2177.6908	0.0000	1.0000
F81IG	16	-1043.1745	2181.5175	0.0000	1.0000
SYM	16	-1047.4249	2190.0183	0.0000	1.0000
JC69I	12	-1059.3980	2190.1723	0.0000	1.0000
JC69G	12	-1059.4724	2190.3211	0.0000	1.0000
F81	14	-1054.5479	2192.3683	0.0000	1.0000
JC69IG	13	-1059.4984	2196.3212	0.0000	1.0000
JC69	11	-1070.2981	2206.0245	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_10135.phylip = 1 - 383;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10135.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_10135.phylip = 1 - 383;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10135.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

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

