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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	60	-976.0749	2099.3617	0.7474	0.7474
HKYIG	61	-976.2167	2102.6573	0.1438	0.8912
HKYI	60	-978.2193	2103.6504	0.0875	0.9788
GTRG	64	-973.5758	2106.5478	0.0206	0.9993
K2PG	57	-987.9973	2114.3035	0.0004	0.9997
GTRIG	65	-976.8863	2116.2726	0.0002	0.9999
SYMG	61	-984.4331	2119.0900	0.0000	0.9999
GTRI	64	-980.1048	2119.6058	0.0000	1.0000
K2PIG	58	-989.6534	2120.5613	0.0000	1.0000
K2PI	57	-992.0795	2122.4678	0.0000	1.0000
SYMI	61	-987.0545	2124.3330	0.0000	1.0000
SYMIG	62	-988.3508	2129.9599	0.0000	1.0000
HKY	59	-997.6430	2139.5083	0.0000	1.0000
F81G	59	-998.1886	2140.5994	0.0000	1.0000
JC69G	56	-1002.9155	2141.2156	0.0000	1.0000
GTR	63	-995.5675	2147.4507	0.0000	1.0000
F81IG	60	-1000.3547	2147.9213	0.0000	1.0000
F81I	59	-1003.3330	2150.8883	0.0000	1.0000
K2P	56	-1008.1392	2151.6630	0.0000	1.0000
JC69IG	57	-1007.5328	2153.3744	0.0000	1.0000
JC69I	56	-1010.2707	2155.9261	0.0000	1.0000
SYM	60	-1004.6692	2156.5503	0.0000	1.0000
F81	58	-1018.0982	2177.4511	0.0000	1.0000
JC69	55	-1023.1820	2178.8457	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-976.0749	2072.1498	0.6014	0.6014
HKYIG	61	-976.2167	2074.4334	0.1920	0.7934
GTRG	64	-973.5758	2075.1516	0.1341	0.9274
HKYI	60	-978.2193	2076.4385	0.0704	0.9979
GTRIG	65	-976.8863	2083.7726	0.0018	0.9997
GTRI	64	-980.1048	2088.2095	0.0002	0.9999
K2PG	57	-987.9973	2089.9946	0.0001	0.9999
SYMG	61	-984.4331	2090.8661	0.0001	1.0000
K2PIG	58	-989.6534	2095.3067	0.0000	1.0000
SYMI	61	-987.0545	2096.1091	0.0000	1.0000
K2PI	57	-992.0795	2098.1590	0.0000	1.0000
SYMIG	62	-988.3508	2100.7015	0.0000	1.0000
HKY	59	-997.6430	2113.2861	0.0000	1.0000
F81G	59	-998.1886	2114.3771	0.0000	1.0000
GTR	63	-995.5675	2117.1349	0.0000	1.0000
JC69G	56	-1002.9155	2117.8310	0.0000	1.0000
F81IG	60	-1000.3547	2120.7094	0.0000	1.0000
F81I	59	-1003.3330	2124.6661	0.0000	1.0000
K2P	56	-1008.1392	2128.2784	0.0000	1.0000
JC69IG	57	-1007.5328	2129.0656	0.0000	1.0000
SYM	60	-1004.6692	2129.3384	0.0000	1.0000
JC69I	56	-1010.2707	2132.5415	0.0000	1.0000
F81	58	-1018.0982	2152.1965	0.0000	1.0000
JC69	55	-1023.1820	2156.3640	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-976.0749	2300.0954	0.8292	0.8292
HKYI	60	-978.2193	2304.3841	0.0971	0.9264
HKYIG	61	-976.2167	2306.1781	0.0396	0.9660
K2PG	57	-987.9973	2306.5429	0.0330	0.9990
K2PI	57	-992.0795	2314.7073	0.0006	0.9995
K2PIG	58	-989.6534	2315.6541	0.0003	0.9999
GTRG	64	-973.5758	2318.2935	0.0001	1.0000
SYMG	61	-984.4331	2322.6108	0.0000	1.0000
SYMI	61	-987.0545	2327.8537	0.0000	1.0000
JC69G	56	-1002.9155	2330.5802	0.0000	1.0000
GTRIG	65	-976.8863	2330.7137	0.0000	1.0000
GTRI	64	-980.1048	2331.3515	0.0000	1.0000
SYMIG	62	-988.3508	2336.2453	0.0000	1.0000
HKY	59	-997.6430	2337.4325	0.0000	1.0000
F81G	59	-998.1886	2338.5236	0.0000	1.0000
K2P	56	-1008.1392	2341.0276	0.0000	1.0000
JC69I	56	-1010.2707	2345.2906	0.0000	1.0000
JC69IG	57	-1007.5328	2345.6139	0.0000	1.0000
F81IG	60	-1000.3547	2348.6550	0.0000	1.0000
F81I	59	-1003.3330	2348.8126	0.0000	1.0000
GTR	63	-995.5675	2356.4777	0.0000	1.0000
SYM	60	-1004.6692	2357.2840	0.0000	1.0000
JC69	55	-1023.1820	2365.3141	0.0000	1.0000
F81	58	-1018.0982	2372.5438	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr6_3291.phylip = 1 - 330;
 Partition Dummy = 1:183-phylip/chr6_3291.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr6_3291.phylip = 1 - 330;
 Partition Dummy = 1:183-phylip/chr6_3291.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr6_3291.phylip = 1 - 330;
 Partition Dummy = 1:183-phylip/chr6_3291.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

