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

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

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


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	-2091.5785	4327.8303	0.6938	0.6938
GTRIG	65	-2091.2772	4329.7834	0.2613	0.9551
GTRI	64	-2094.6900	4334.0533	0.0309	0.9860
HKYG	60	-2100.9041	4336.3609	0.0097	0.9957
HKYIG	61	-2100.6050	4338.2778	0.0037	0.9994
HKYI	60	-2103.7649	4342.0825	0.0006	1.0000
SYMG	61	-2109.3408	4355.7494	0.0000	1.0000
K2PG	57	-2115.3687	4357.8046	0.0000	1.0000
SYMIG	62	-2109.1207	4357.8342	0.0000	1.0000
K2PIG	58	-2115.1234	4359.7993	0.0000	1.0000
SYMI	61	-2111.7291	4360.5259	0.0000	1.0000
K2PI	57	-2117.8172	4362.7016	0.0000	1.0000
GTR	63	-2126.7268	4395.5816	0.0000	1.0000
HKY	59	-2136.3364	4404.7204	0.0000	1.0000
SYM	60	-2143.2165	4420.9856	0.0000	1.0000
K2P	56	-2150.3143	4425.2203	0.0000	1.0000
F81G	59	-2165.6429	4463.3335	0.0000	1.0000
F81IG	60	-2165.4464	4465.4456	0.0000	1.0000
F81I	59	-2168.0479	4468.1434	0.0000	1.0000
JC69G	56	-2175.4734	4475.5385	0.0000	1.0000
JC69IG	57	-2175.2873	4477.6419	0.0000	1.0000
JC69I	56	-2177.7641	4480.1200	0.0000	1.0000
F81	58	-2199.6437	4528.8399	0.0000	1.0000
JC69	55	-2209.5548	4541.2357	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2091.5785	4311.1570	0.6454	0.6454
GTRIG	65	-2091.2772	4312.5544	0.3209	0.9664
GTRI	64	-2094.6900	4317.3800	0.0287	0.9951
HKYG	60	-2100.9041	4321.8082	0.0031	0.9983
HKYIG	61	-2100.6050	4323.2100	0.0016	0.9998
HKYI	60	-2103.7649	4327.5298	0.0002	1.0000
SYMG	61	-2109.3408	4340.6817	0.0000	1.0000
SYMIG	62	-2109.1207	4342.2414	0.0000	1.0000
K2PG	57	-2115.3687	4344.7374	0.0000	1.0000
SYMI	61	-2111.7291	4345.4582	0.0000	1.0000
K2PIG	58	-2115.1234	4346.2469	0.0000	1.0000
K2PI	57	-2117.8172	4349.6344	0.0000	1.0000
GTR	63	-2126.7268	4379.4536	0.0000	1.0000
HKY	59	-2136.3364	4390.6728	0.0000	1.0000
SYM	60	-2143.2165	4406.4330	0.0000	1.0000
K2P	56	-2150.3143	4412.6286	0.0000	1.0000
F81G	59	-2165.6429	4449.2859	0.0000	1.0000
F81IG	60	-2165.4464	4450.8929	0.0000	1.0000
F81I	59	-2168.0479	4454.0958	0.0000	1.0000
JC69G	56	-2175.4734	4462.9468	0.0000	1.0000
JC69IG	57	-2175.2873	4464.5747	0.0000	1.0000
JC69I	56	-2177.7641	4467.5283	0.0000	1.0000
F81	58	-2199.6437	4515.2874	0.0000	1.0000
JC69	55	-2209.5548	4529.1097	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2100.9041	4581.9114	0.8614	0.8614
HKYI	60	-2103.7649	4587.6331	0.0493	0.9107
HKYIG	61	-2100.6050	4587.6483	0.0489	0.9596
GTRG	64	-2091.5785	4588.6005	0.0304	0.9900
K2PG	57	-2115.3687	4591.8355	0.0060	0.9961
GTRIG	65	-2091.2772	4594.3330	0.0017	0.9978
GTRI	64	-2094.6900	4594.8234	0.0014	0.9991
K2PI	57	-2117.8172	4596.7325	0.0005	0.9997
K2PIG	58	-2115.1234	4597.6800	0.0003	1.0000
SYMG	61	-2109.3408	4605.1200	0.0000	1.0000
SYMI	61	-2111.7291	4609.8965	0.0000	1.0000
SYMIG	62	-2109.1207	4611.0148	0.0000	1.0000
HKY	59	-2136.3364	4646.4410	0.0000	1.0000
GTR	63	-2126.7268	4652.5620	0.0000	1.0000
K2P	56	-2150.3143	4655.3916	0.0000	1.0000
SYM	60	-2143.2165	4666.5362	0.0000	1.0000
F81G	59	-2165.6429	4705.0541	0.0000	1.0000
JC69G	56	-2175.4734	4705.7098	0.0000	1.0000
F81I	59	-2168.0479	4709.8640	0.0000	1.0000
JC69I	56	-2177.7641	4710.2913	0.0000	1.0000
F81IG	60	-2165.4464	4710.9961	0.0000	1.0000
JC69IG	57	-2175.2873	4711.6728	0.0000	1.0000
F81	58	-2199.6437	4766.7206	0.0000	1.0000
JC69	55	-2209.5548	4767.5377	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_10445.phylip = 1 - 564;
 Partition Dummy = 1:183-phylip/chr7_10445.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_10445.phylip = 1 - 564;
 Partition Dummy = 1:183-phylip/chr7_10445.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 (HKYG)]

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


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

