
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
10:55 08/01/2011

Input data from file "261-aggressive-phylip/12088_hg18.phylip" (29 taxa, 751 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-1050.6061	2241.3404	0.5195	0.5195
HKYG	60	-1056.1833	2242.9754	0.2294	0.7488
GTRIG	65	-1050.6445	2243.8145	0.1508	0.8996
HKYIG	61	-1056.2218	2245.4219	0.0675	0.9671
GTR	63	-1055.7640	2249.2661	0.0099	0.9770
HKY	59	-1060.5128	2249.2717	0.0098	0.9868
K2PG	57	-1063.8395	2251.2201	0.0037	0.9905
HKYI	60	-1060.5502	2251.7091	0.0029	0.9935
GTRI	64	-1055.8011	2251.7304	0.0029	0.9963
SYMG	61	-1059.9044	2252.7869	0.0017	0.9980
K2PIG	58	-1063.8779	2253.6460	0.0011	0.9991
SYMIG	62	-1059.9426	2255.2399	0.0005	0.9996
K2P	56	-1067.9016	2257.0020	0.0002	0.9998
SYM	60	-1064.2503	2259.1092	0.0001	0.9999
K2PI	57	-1067.9390	2259.4191	0.0001	1.0000
SYMI	61	-1064.2875	2261.5532	0.0000	1.0000
F81G	59	-1069.0225	2266.2910	0.0000	1.0000
F81IG	60	-1069.0611	2268.7308	0.0000	1.0000
F81	58	-1072.9541	2271.7983	0.0000	1.0000
JC69G	56	-1075.3865	2271.9718	0.0000	1.0000
F81I	59	-1072.9915	2274.2290	0.0000	1.0000
JC69IG	57	-1075.4249	2274.3909	0.0000	1.0000
JC69	55	-1079.2333	2277.3299	0.0000	1.0000
JC69I	56	-1079.2707	2279.7402	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1050.6061	2229.2121	0.5987	0.5987
GTRIG	65	-1050.6445	2231.2889	0.2120	0.8107
HKYG	60	-1056.1833	2232.3667	0.1237	0.9343
HKYIG	61	-1056.2218	2234.4436	0.0438	0.9781
GTR	63	-1055.7640	2237.5281	0.0094	0.9875
HKY	59	-1060.5128	2239.0257	0.0044	0.9919
GTRI	64	-1055.8011	2239.6021	0.0033	0.9952
HKYI	60	-1060.5502	2241.1004	0.0016	0.9968
K2PG	57	-1063.8395	2241.6790	0.0012	0.9980
SYMG	61	-1059.9044	2241.8087	0.0011	0.9991
K2PIG	58	-1063.8779	2243.7558	0.0004	0.9995
SYMIG	62	-1059.9426	2243.8852	0.0004	0.9999
K2P	56	-1067.9016	2247.8032	0.0001	0.9999
SYM	60	-1064.2503	2248.5005	0.0000	1.0000
K2PI	57	-1067.9390	2249.8780	0.0000	1.0000
SYMI	61	-1064.2875	2250.5750	0.0000	1.0000
F81G	59	-1069.0225	2256.0450	0.0000	1.0000
F81IG	60	-1069.0611	2258.1221	0.0000	1.0000
F81	58	-1072.9541	2261.9081	0.0000	1.0000
JC69G	56	-1075.3865	2262.7730	0.0000	1.0000
F81I	59	-1072.9915	2263.9829	0.0000	1.0000
JC69IG	57	-1075.4249	2264.8497	0.0000	1.0000
JC69	55	-1079.2333	2268.4666	0.0000	1.0000
JC69I	56	-1079.2707	2270.5414	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1063.8395	2505.0991	0.5993	0.5993
K2P	56	-1067.9016	2506.6019	0.2827	0.8819
HKYG	60	-1056.1833	2509.6510	0.0615	0.9435
HKY	59	-1060.5128	2511.6886	0.0222	0.9657
K2PIG	58	-1063.8779	2511.7974	0.0210	0.9867
K2PI	57	-1067.9390	2513.2981	0.0099	0.9967
HKYIG	61	-1056.2218	2516.3494	0.0022	0.9988
HKYI	60	-1060.5502	2518.3847	0.0008	0.9996
JC69G	56	-1075.3865	2521.5717	0.0002	0.9998
JC69	55	-1079.2333	2522.6439	0.0001	0.9999
SYMG	61	-1059.9044	2523.7144	0.0001	0.9999
GTRG	64	-1050.6061	2524.9821	0.0000	1.0000
SYM	60	-1064.2503	2525.7849	0.0000	1.0000
JC69IG	57	-1075.4249	2528.2699	0.0000	1.0000
GTR	63	-1055.7640	2528.6767	0.0000	1.0000
F81G	59	-1069.0225	2528.7080	0.0000	1.0000
JC69I	56	-1079.2707	2529.3401	0.0000	1.0000
F81	58	-1072.9541	2529.9496	0.0000	1.0000
SYMIG	62	-1059.9426	2530.4124	0.0000	1.0000
GTRIG	65	-1050.6445	2531.6803	0.0000	1.0000
SYMI	61	-1064.2875	2532.4807	0.0000	1.0000
GTRI	64	-1055.8011	2535.3721	0.0000	1.0000
F81IG	60	-1069.0611	2535.4064	0.0000	1.0000
F81I	59	-1072.9915	2536.6459	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 261-aggressive-phylip/12088_hg18.phylip = 1 - 751;
 Partition Dummy = 1:261-aggressive-phylip/12088_hg18.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 261-aggressive-phylip/12088_hg18.phylip = 1 - 751;
 Partition Dummy = 1:261-aggressive-phylip/12088_hg18.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 (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/12088_hg18.phylip = 1 - 751;
 Partition Dummy = 1:261-aggressive-phylip/12088_hg18.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

