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

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

Minimum AIC  model: HKYI
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	-872.2656	1894.2873	0.4628	0.4628
HKYG	60	-872.4523	1894.6607	0.3840	0.8468
HKYIG	61	-872.4340	1897.7415	0.0823	0.9291
GTRG	64	-868.6350	1899.6502	0.0317	0.9608
GTRI	64	-868.7307	1899.8416	0.0288	0.9896
GTRIG	65	-868.6220	1902.8458	0.0064	0.9960
HKY	59	-878.6641	1903.9921	0.0036	0.9996
GTR	63	-874.5587	1908.3026	0.0004	1.0000
F81I	59	-887.1857	1921.0353	0.0000	1.0000
F81G	59	-887.2845	1921.2329	0.0000	1.0000
F81IG	60	-887.2622	1924.2804	0.0000	1.0000
K2PG	57	-894.1273	1928.8088	0.0000	1.0000
K2PI	57	-894.2921	1929.1383	0.0000	1.0000
F81	58	-893.6043	1930.8054	0.0000	1.0000
K2PIG	58	-894.1026	1931.8019	0.0000	1.0000
SYMG	61	-892.6848	1938.2430	0.0000	1.0000
SYMI	61	-892.7408	1938.3551	0.0000	1.0000
K2P	56	-901.7393	1941.0146	0.0000	1.0000
SYMIG	62	-892.6760	1941.3684	0.0000	1.0000
SYM	60	-899.9879	1949.7320	0.0000	1.0000
JC69G	56	-906.4315	1950.3990	0.0000	1.0000
JC69I	56	-906.6796	1950.8951	0.0000	1.0000
JC69IG	57	-906.3995	1953.3532	0.0000	1.0000
JC69	55	-913.6791	1961.9000	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-872.2656	1864.5312	0.2683	0.2683
HKYG	60	-872.4523	1864.9046	0.2226	0.4909
GTRG	64	-868.6350	1865.2701	0.1854	0.6764
GTRI	64	-868.7307	1865.4614	0.1685	0.8449
HKYIG	61	-872.4340	1866.8680	0.0834	0.9283
GTRIG	65	-868.6220	1867.2441	0.0691	0.9974
GTR	63	-874.5587	1875.1174	0.0013	0.9988
HKY	59	-878.6641	1875.3282	0.0012	1.0000
F81I	59	-887.1857	1892.3713	0.0000	1.0000
F81G	59	-887.2845	1892.5689	0.0000	1.0000
F81IG	60	-887.2622	1894.5243	0.0000	1.0000
K2PG	57	-894.1273	1902.2546	0.0000	1.0000
K2PI	57	-894.2921	1902.5841	0.0000	1.0000
F81	58	-893.6043	1903.2087	0.0000	1.0000
K2PIG	58	-894.1026	1904.2052	0.0000	1.0000
SYMG	61	-892.6848	1907.3696	0.0000	1.0000
SYMI	61	-892.7408	1907.4816	0.0000	1.0000
SYMIG	62	-892.6760	1909.3520	0.0000	1.0000
K2P	56	-901.7393	1915.4786	0.0000	1.0000
SYM	60	-899.9879	1919.9759	0.0000	1.0000
JC69G	56	-906.4315	1924.8630	0.0000	1.0000
JC69I	56	-906.6796	1925.3591	0.0000	1.0000
JC69IG	57	-906.3995	1926.7990	0.0000	1.0000
JC69	55	-913.6791	1937.3581	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-872.2656	2088.1421	0.5241	0.5241
HKYG	60	-872.4523	2088.5154	0.4349	0.9590
HKYIG	61	-872.4340	2094.2057	0.0253	0.9843
HKY	59	-878.6641	2095.2122	0.0153	0.9996
GTRG	64	-868.6350	2103.7883	0.0002	0.9998
GTRI	64	-868.7307	2103.9797	0.0002	1.0000
GTRIG	65	-868.6220	2109.4892	0.0000	1.0000
GTR	63	-874.5587	2109.9088	0.0000	1.0000
F81I	59	-887.1857	2112.2553	0.0000	1.0000
F81G	59	-887.2845	2112.4529	0.0000	1.0000
K2PG	57	-894.1273	2114.6849	0.0000	1.0000
K2PI	57	-894.2921	2115.0144	0.0000	1.0000
F81IG	60	-887.2622	2118.1352	0.0000	1.0000
F81	58	-893.6043	2119.3658	0.0000	1.0000
K2PIG	58	-894.1026	2120.3623	0.0000	1.0000
K2P	56	-901.7393	2124.1821	0.0000	1.0000
JC69G	56	-906.4315	2133.5664	0.0000	1.0000
JC69I	56	-906.6796	2134.0626	0.0000	1.0000
SYMG	61	-892.6848	2134.7073	0.0000	1.0000
SYMI	61	-892.7408	2134.8193	0.0000	1.0000
JC69IG	57	-906.3995	2139.2293	0.0000	1.0000
SYMIG	62	-892.6760	2140.4166	0.0000	1.0000
JC69	55	-913.6791	2142.3348	0.0000	1.0000
SYM	60	-899.9879	2143.5868	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr6_4403.phylip = 1 - 307;
 Partition Dummy = 1:183-phylip/chr6_4403.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 AICc model (HKYI)]

BEGIN MRBAYES;
 Charset 183-phylip/chr6_4403.phylip = 1 - 307;
 Partition Dummy = 1:183-phylip/chr6_4403.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/chr6_4403.phylip = 1 - 307;
 Partition Dummy = 1:183-phylip/chr6_4403.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

