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

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

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


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
K2PG	57	-1054.8630	2248.6770	0.3074	0.3074
HKYG	60	-1050.7312	2249.4013	0.2140	0.5213
K2PI	57	-1055.3663	2249.6834	0.1858	0.7072
HKYI	60	-1051.2024	2250.3437	0.1336	0.8407
K2PIG	58	-1054.9086	2251.7415	0.0664	0.9072
HKYIG	61	-1050.7730	2252.5268	0.0448	0.9520
SYMG	61	-1051.9799	2254.9406	0.0134	0.9654
GTRG	64	-1047.4420	2255.1320	0.0122	0.9776
SYMI	61	-1052.4379	2255.8567	0.0085	0.9861
GTRI	64	-1047.9507	2256.1496	0.0073	0.9934
SYMIG	62	-1051.8462	2257.7386	0.0033	0.9967
GTRIG	65	-1047.1892	2257.7636	0.0033	1.0000
JC69G	56	-1068.1287	2272.2573	0.0000	1.0000
F81G	59	-1064.0570	2273.0342	0.0000	1.0000
JC69I	56	-1068.6160	2273.2319	0.0000	1.0000
F81I	59	-1064.5726	2274.0653	0.0000	1.0000
JC69IG	57	-1067.7032	2274.3573	0.0000	1.0000
F81IG	60	-1063.5558	2275.0506	0.0000	1.0000
K2P	56	-1089.6544	2315.3088	0.0000	1.0000
HKY	59	-1085.7924	2316.5050	0.0000	1.0000
SYM	60	-1087.4765	2322.8919	0.0000	1.0000
GTR	63	-1083.2420	2323.6192	0.0000	1.0000
JC69	55	-1102.7505	2338.5723	0.0000	1.0000
F81	58	-1098.7355	2339.3953	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1050.7312	2221.4624	0.2602	0.2602
HKYI	60	-1051.2024	2222.4048	0.1624	0.4226
GTRG	64	-1047.4420	2222.8840	0.1278	0.5504
HKYIG	61	-1050.7730	2223.5460	0.0918	0.6422
K2PG	57	-1054.8630	2223.7261	0.0839	0.7260
GTRI	64	-1047.9507	2223.9015	0.0768	0.8029
GTRIG	65	-1047.1892	2224.3783	0.0605	0.8634
K2PI	57	-1055.3663	2224.7325	0.0507	0.9141
K2PIG	58	-1054.9086	2225.8172	0.0295	0.9436
SYMG	61	-1051.9799	2225.9597	0.0275	0.9711
SYMI	61	-1052.4379	2226.8758	0.0174	0.9885
SYMIG	62	-1051.8462	2227.6924	0.0115	1.0000
F81G	59	-1064.0570	2246.1141	0.0000	1.0000
F81IG	60	-1063.5558	2247.1116	0.0000	1.0000
F81I	59	-1064.5726	2247.1452	0.0000	1.0000
JC69G	56	-1068.1287	2248.2573	0.0000	1.0000
JC69I	56	-1068.6160	2249.2319	0.0000	1.0000
JC69IG	57	-1067.7032	2249.4063	0.0000	1.0000
HKY	59	-1085.7924	2289.5848	0.0000	1.0000
K2P	56	-1089.6544	2291.3088	0.0000	1.0000
GTR	63	-1083.2420	2292.4841	0.0000	1.0000
SYM	60	-1087.4765	2294.9530	0.0000	1.0000
F81	58	-1098.7355	2313.4711	0.0000	1.0000
JC69	55	-1102.7505	2315.5011	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1054.8630	2439.0523	0.5966	0.5966
K2PI	57	-1055.3663	2440.0587	0.3607	0.9573
K2PIG	58	-1054.9086	2444.9211	0.0317	0.9891
HKYG	60	-1050.7312	2448.1215	0.0064	0.9955
HKYI	60	-1051.2024	2449.0639	0.0040	0.9995
HKYIG	61	-1050.7730	2453.9828	0.0003	0.9998
SYMG	61	-1051.9799	2456.3965	0.0001	0.9999
SYMI	61	-1052.4379	2457.3126	0.0001	1.0000
JC69G	56	-1068.1287	2459.8059	0.0000	1.0000
JC69I	56	-1068.6160	2460.7804	0.0000	1.0000
SYMIG	62	-1051.8462	2461.9069	0.0000	1.0000
GTRG	64	-1047.4420	2464.6537	0.0000	1.0000
JC69IG	57	-1067.7032	2464.7325	0.0000	1.0000
GTRI	64	-1047.9507	2465.6712	0.0000	1.0000
F81G	59	-1064.0570	2468.9955	0.0000	1.0000
GTRIG	65	-1047.1892	2469.9257	0.0000	1.0000
F81I	59	-1064.5726	2470.0266	0.0000	1.0000
F81IG	60	-1063.5558	2473.7708	0.0000	1.0000
K2P	56	-1089.6544	2502.8574	0.0000	1.0000
HKY	59	-1085.7924	2512.4663	0.0000	1.0000
SYM	60	-1087.4765	2521.6121	0.0000	1.0000
JC69	55	-1102.7505	2523.2720	0.0000	1.0000
GTR	63	-1083.2420	2530.4762	0.0000	1.0000
F81	58	-1098.7355	2532.5749	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/chr7_4732.phylip = 1 - 323;
 Partition Dummy = 1:183-phylip/chr7_4732.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 (K2PG)]

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


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

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


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

