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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: SYMG


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	64	-960.3576	2078.5360	0.8567	0.8567
SYMG	61	-966.6666	2082.1558	0.1402	0.9969
GTRI	64	-966.2530	2090.3268	0.0024	0.9993
GTRIG	65	-966.3782	2093.6198	0.0005	0.9997
SYMI	61	-973.0836	2094.9899	0.0002	0.9999
SYMIG	62	-973.0620	2097.9247	0.0001	1.0000
HKYG	60	-977.7453	2101.3564	0.0000	1.0000
K2PG	57	-984.1824	2105.4836	0.0000	1.0000
K2PI	57	-984.4523	2106.0234	0.0000	1.0000
F81G	59	-982.5608	2108.0512	0.0000	1.0000
K2PIG	58	-984.0670	2108.1481	0.0000	1.0000
GTR	63	-976.7934	2108.3869	0.0000	1.0000
SYM	60	-981.8437	2109.5531	0.0000	1.0000
JC69G	56	-988.4556	2111.1552	0.0000	1.0000
HKYI	60	-984.8664	2115.5985	0.0000	1.0000
HKYIG	61	-984.7174	2118.2575	0.0000	1.0000
F81I	59	-989.8635	2122.6566	0.0000	1.0000
F81IG	60	-989.7796	2125.4249	0.0000	1.0000
JC69I	56	-996.0528	2126.3495	0.0000	1.0000
JC69IG	57	-995.9083	2128.9355	0.0000	1.0000
HKY	59	-996.0335	2134.9967	0.0000	1.0000
K2P	56	-1001.3710	2136.9859	0.0000	1.0000
F81	58	-1000.0971	2140.2082	0.0000	1.0000
JC69	55	-1005.4461	2142.2810	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-960.3576	2048.7152	0.9613	0.9613
SYMG	61	-966.6666	2055.3331	0.0351	0.9964
GTRI	64	-966.2530	2060.5061	0.0026	0.9991
GTRIG	65	-966.3782	2062.7565	0.0009	0.9999
SYMI	61	-973.0836	2068.1672	0.0001	1.0000
SYMIG	62	-973.0620	2070.1240	0.0000	1.0000
HKYG	60	-977.7453	2075.4906	0.0000	1.0000
GTR	63	-976.7934	2079.5869	0.0000	1.0000
K2PG	57	-984.1824	2082.3647	0.0000	1.0000
K2PI	57	-984.4523	2082.9046	0.0000	1.0000
F81G	59	-982.5608	2083.1216	0.0000	1.0000
SYM	60	-981.8437	2083.6874	0.0000	1.0000
K2PIG	58	-984.0670	2084.1340	0.0000	1.0000
JC69G	56	-988.4556	2088.9113	0.0000	1.0000
HKYI	60	-984.8664	2089.7328	0.0000	1.0000
HKYIG	61	-984.7174	2091.4348	0.0000	1.0000
F81I	59	-989.8635	2097.7270	0.0000	1.0000
F81IG	60	-989.7796	2099.5592	0.0000	1.0000
JC69I	56	-996.0528	2104.1056	0.0000	1.0000
JC69IG	57	-995.9083	2105.8166	0.0000	1.0000
HKY	59	-996.0335	2110.0671	0.0000	1.0000
K2P	56	-1001.3710	2114.7420	0.0000	1.0000
F81	58	-1000.0971	2116.1941	0.0000	1.0000
JC69	55	-1005.4461	2120.8921	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMG	61	-966.6666	2289.6122	0.9136	0.9136
GTRG	64	-960.3576	2294.5163	0.0787	0.9923
K2PG	57	-984.1824	2301.2813	0.0027	0.9950
K2PI	57	-984.4523	2301.8211	0.0020	0.9970
SYMI	61	-973.0836	2302.4463	0.0015	0.9985
JC69G	56	-988.4556	2303.9872	0.0007	0.9992
HKYG	60	-977.7453	2305.9291	0.0003	0.9995
GTRI	64	-966.2530	2306.3071	0.0002	0.9997
K2PIG	58	-984.0670	2306.8912	0.0002	0.9999
SYMIG	62	-973.0620	2308.2438	0.0001	0.9999
F81G	59	-982.5608	2309.7195	0.0000	1.0000
GTRIG	65	-966.3782	2312.3982	0.0000	1.0000
SYM	60	-981.8437	2314.1259	0.0000	1.0000
JC69I	56	-996.0528	2319.1816	0.0000	1.0000
HKYI	60	-984.8664	2320.1713	0.0000	1.0000
GTR	63	-976.7934	2321.5473	0.0000	1.0000
F81I	59	-989.8635	2324.3248	0.0000	1.0000
JC69IG	57	-995.9083	2324.7332	0.0000	1.0000
HKYIG	61	-984.7174	2325.7140	0.0000	1.0000
K2P	56	-1001.3710	2329.8179	0.0000	1.0000
F81IG	60	-989.7796	2329.9977	0.0000	1.0000
JC69	55	-1005.4461	2332.1274	0.0000	1.0000
HKY	59	-996.0335	2336.6649	0.0000	1.0000
F81	58	-1000.0971	2338.9514	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 183-phylip/chr3_23229.phylip = 1 - 344;
 Partition Dummy = 1:183-phylip/chr3_23229.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 183-phylip/chr3_23229.phylip = 1 - 344;
 Partition Dummy = 1:183-phylip/chr3_23229.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 (SYMG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr3_23229.phylip = 1 - 344;
 Partition Dummy = 1:183-phylip/chr3_23229.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=Fixed(Equal);
END;


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

