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

Input data from file "183-phylip/chr7_4994.phylip" (29 taxa, 297 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	-1136.4862	2436.8345	0.4143	0.4143
HKYG	60	-1143.7738	2438.5645	0.1744	0.5887
K2PG	57	-1148.5689	2438.8031	0.1548	0.7435
GTRIG	65	-1135.9609	2439.0646	0.1358	0.8794
GTRI	64	-1138.6471	2441.1563	0.0477	0.9271
K2PIG	58	-1148.8790	2442.5143	0.0242	0.9513
HKYIG	61	-1144.3352	2442.8575	0.0204	0.9717
HKYI	60	-1146.3664	2443.7498	0.0131	0.9848
K2PI	57	-1151.4057	2444.4768	0.0091	0.9938
SYMG	61	-1145.8449	2445.8771	0.0045	0.9983
SYMIG	62	-1145.8267	2449.0380	0.0009	0.9993
SYMI	61	-1147.6541	2449.4955	0.0007	1.0000
JC69G	56	-1174.9241	2488.4482	0.0000	1.0000
JC69IG	57	-1174.0424	2489.7501	0.0000	1.0000
F81G	59	-1172.3304	2492.5343	0.0000	1.0000
JC69I	56	-1177.4959	2493.5919	0.0000	1.0000
F81IG	60	-1171.4960	2494.0090	0.0000	1.0000
F81I	59	-1174.7494	2497.3723	0.0000	1.0000
K2P	56	-1196.4890	2531.5780	0.0000	1.0000
HKY	59	-1193.2786	2534.4306	0.0000	1.0000
SYM	60	-1194.9377	2540.8924	0.0000	1.0000
GTR	63	-1190.1565	2540.9225	0.0000	1.0000
JC69	55	-1221.4929	2578.5460	0.0000	1.0000
F81	58	-1219.4301	2583.6165	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1136.4862	2400.9724	0.5567	0.5567
GTRIG	65	-1135.9609	2401.9217	0.3463	0.9031
GTRI	64	-1138.6471	2405.2942	0.0641	0.9672
HKYG	60	-1143.7738	2407.5475	0.0208	0.9880
HKYIG	61	-1144.3352	2410.6703	0.0044	0.9924
K2PG	57	-1148.5689	2411.1378	0.0035	0.9958
HKYI	60	-1146.3664	2412.7329	0.0016	0.9974
SYMG	61	-1145.8449	2413.6899	0.0010	0.9983
K2PIG	58	-1148.8790	2413.7580	0.0009	0.9993
SYMIG	62	-1145.8267	2415.6534	0.0004	0.9996
K2PI	57	-1151.4057	2416.8115	0.0002	0.9998
SYMI	61	-1147.6541	2417.3082	0.0002	1.0000
JC69G	56	-1174.9241	2461.8482	0.0000	1.0000
JC69IG	57	-1174.0424	2462.0848	0.0000	1.0000
F81G	59	-1172.3304	2462.6609	0.0000	1.0000
F81IG	60	-1171.4960	2462.9920	0.0000	1.0000
JC69I	56	-1177.4959	2466.9919	0.0000	1.0000
F81I	59	-1174.7494	2467.4989	0.0000	1.0000
HKY	59	-1193.2786	2504.5572	0.0000	1.0000
K2P	56	-1196.4890	2504.9780	0.0000	1.0000
GTR	63	-1190.1565	2506.3131	0.0000	1.0000
SYM	60	-1194.9377	2509.8755	0.0000	1.0000
JC69	55	-1221.4929	2552.9858	0.0000	1.0000
F81	58	-1219.4301	2554.8602	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1148.5689	2621.6806	0.8865	0.8865
K2PI	57	-1151.4057	2627.3542	0.0520	0.9385
K2PIG	58	-1148.8790	2627.9944	0.0377	0.9762
HKYG	60	-1143.7738	2629.1714	0.0209	0.9971
HKYI	60	-1146.3664	2634.3568	0.0016	0.9987
HKYIG	61	-1144.3352	2635.9880	0.0007	0.9994
GTRG	64	-1136.4862	2637.3713	0.0003	0.9997
SYMG	61	-1145.8449	2639.0076	0.0002	0.9999
GTRI	64	-1138.6471	2641.6931	0.0000	0.9999
GTRIG	65	-1135.9609	2642.0143	0.0000	1.0000
SYMI	61	-1147.6541	2642.6259	0.0000	1.0000
SYMIG	62	-1145.8267	2644.6648	0.0000	1.0000
JC69G	56	-1174.9241	2668.6972	0.0000	1.0000
JC69IG	57	-1174.0424	2672.6275	0.0000	1.0000
JC69I	56	-1177.4959	2673.8409	0.0000	1.0000
F81G	59	-1172.3304	2680.5911	0.0000	1.0000
F81IG	60	-1171.4960	2684.6159	0.0000	1.0000
F81I	59	-1174.7494	2685.4291	0.0000	1.0000
K2P	56	-1196.4890	2711.8270	0.0000	1.0000
HKY	59	-1193.2786	2722.4874	0.0000	1.0000
SYM	60	-1194.9377	2731.4994	0.0000	1.0000
GTR	63	-1190.1565	2739.0182	0.0000	1.0000
JC69	55	-1221.4929	2756.1411	0.0000	1.0000
F81	58	-1219.4301	2769.0967	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/chr7_4994.phylip = 1 - 297;
 Partition Dummy = 1:183-phylip/chr7_4994.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/chr7_4994.phylip = 1 - 297;
 Partition Dummy = 1:183-phylip/chr7_4994.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 183-phylip/chr7_4994.phylip = 1 - 297;
 Partition Dummy = 1:183-phylip/chr7_4994.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

