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

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

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


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	-1000.2151	2158.9063	0.5143	0.5143
GTRIG	65	-999.2924	2160.1289	0.2791	0.7934
GTRI	64	-1001.1282	2160.7327	0.2064	0.9998
HKYG	60	-1014.8690	2176.1640	0.0001	0.9999
HKYIG	61	-1013.9323	2177.2703	0.0001	0.9999
HKYI	60	-1015.4953	2177.4167	0.0000	1.0000
SYMG	61	-1016.6450	2182.6958	0.0000	1.0000
SYMIG	62	-1015.6878	2183.7829	0.0000	1.0000
SYMI	61	-1017.4372	2184.2802	0.0000	1.0000
K2PG	57	-1025.9390	2189.4923	0.0000	1.0000
K2PIG	58	-1024.9583	2190.4470	0.0000	1.0000
K2PI	57	-1026.5693	2190.7529	0.0000	1.0000
F81G	59	-1039.3792	2222.2260	0.0000	1.0000
F81IG	60	-1038.4950	2223.4161	0.0000	1.0000
F81I	59	-1040.0073	2223.4822	0.0000	1.0000
GTR	63	-1036.3309	2228.0925	0.0000	1.0000
JC69G	56	-1050.4386	2235.5961	0.0000	1.0000
JC69IG	57	-1049.5878	2236.7899	0.0000	1.0000
JC69I	56	-1051.0760	2236.8708	0.0000	1.0000
HKY	59	-1048.5473	2240.5622	0.0000	1.0000
SYM	60	-1050.1845	2246.7949	0.0000	1.0000
K2P	56	-1059.4120	2253.5428	0.0000	1.0000
F81	58	-1071.8227	2284.1759	0.0000	1.0000
JC69	55	-1082.2868	2296.4175	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1000.2151	2128.4301	0.4298	0.4298
GTRIG	65	-999.2924	2128.5848	0.3978	0.8275
GTRI	64	-1001.1282	2130.2565	0.1724	1.0000
HKYG	60	-1014.8690	2149.7380	0.0000	1.0000
HKYIG	61	-1013.9323	2149.8645	0.0000	1.0000
HKYI	60	-1015.4953	2150.9907	0.0000	1.0000
SYMG	61	-1016.6450	2155.2900	0.0000	1.0000
SYMIG	62	-1015.6878	2155.3756	0.0000	1.0000
SYMI	61	-1017.4372	2156.8744	0.0000	1.0000
K2PG	57	-1025.9390	2165.8780	0.0000	1.0000
K2PIG	58	-1024.9583	2165.9165	0.0000	1.0000
K2PI	57	-1026.5693	2167.1387	0.0000	1.0000
F81G	59	-1039.3792	2196.7584	0.0000	1.0000
F81IG	60	-1038.4950	2196.9901	0.0000	1.0000
F81I	59	-1040.0073	2198.0145	0.0000	1.0000
GTR	63	-1036.3309	2198.6618	0.0000	1.0000
JC69G	56	-1050.4386	2212.8772	0.0000	1.0000
JC69IG	57	-1049.5878	2213.1756	0.0000	1.0000
JC69I	56	-1051.0760	2214.1520	0.0000	1.0000
HKY	59	-1048.5473	2215.0946	0.0000	1.0000
SYM	60	-1050.1845	2220.3689	0.0000	1.0000
K2P	56	-1059.4120	2230.8239	0.0000	1.0000
F81	58	-1071.8227	2259.6455	0.0000	1.0000
JC69	55	-1082.2868	2274.5735	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-1000.2151	2373.1051	0.6137	0.6137
GTRI	64	-1001.1282	2374.9314	0.2462	0.8599
GTRIG	65	-999.2924	2377.0827	0.0840	0.9439
HKYG	60	-1014.8690	2379.1207	0.0303	0.9742
HKYI	60	-1015.4953	2380.3734	0.0162	0.9904
HKYIG	61	-1013.9323	2383.0703	0.0042	0.9946
K2PG	57	-1025.9390	2383.7917	0.0029	0.9976
K2PI	57	-1026.5693	2385.0523	0.0016	0.9991
K2PIG	58	-1024.9583	2387.6532	0.0004	0.9996
SYMG	61	-1016.6450	2388.4958	0.0003	0.9998
SYMI	61	-1017.4372	2390.0802	0.0001	1.0000
SYMIG	62	-1015.6878	2392.4045	0.0000	1.0000
F81G	59	-1039.3792	2422.3181	0.0000	1.0000
F81I	59	-1040.0073	2423.5742	0.0000	1.0000
F81IG	60	-1038.4950	2426.3728	0.0000	1.0000
JC69G	56	-1050.4386	2426.9678	0.0000	1.0000
JC69I	56	-1051.0760	2428.2426	0.0000	1.0000
JC69IG	57	-1049.5878	2431.0892	0.0000	1.0000
GTR	63	-1036.3309	2439.5137	0.0000	1.0000
HKY	59	-1048.5473	2440.6543	0.0000	1.0000
K2P	56	-1059.4120	2444.9145	0.0000	1.0000
SYM	60	-1050.1845	2449.7517	0.0000	1.0000
F81	58	-1071.8227	2481.3821	0.0000	1.0000
JC69	55	-1082.2868	2484.8411	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/chr17_3088.phylip = 1 - 338;
 Partition Dummy = 1:183-phylip/chr17_3088.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/chr17_3088.phylip = 1 - 338;
 Partition Dummy = 1:183-phylip/chr17_3088.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 (GTRG)]

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


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

