
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:56 04/30/2011

Input data from file "591-tomodel-temp-phylip/chr10_4411.phylip" (7 taxa, 715 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	20	-2035.5811	4112.3725	0.5766	0.5766
GTRIG	21	-2035.5128	4114.3590	0.2136	0.7901
GTRI	20	-2036.6095	4114.4294	0.2062	0.9963
HKYG	16	-2045.3703	4123.5200	0.0022	0.9985
HKYIG	17	-2045.3097	4125.4974	0.0008	0.9993
HKYI	16	-2046.5005	4125.7804	0.0007	1.0000
SYMG	17	-2050.4887	4135.8555	0.0000	1.0000
SYMI	17	-2051.4716	4137.8212	0.0000	1.0000
SYMIG	18	-2050.4345	4137.8517	0.0000	1.0000
GTR	19	-2050.0678	4139.2290	0.0000	1.0000
K2PG	13	-2057.5696	4141.6584	0.0000	1.0000
K2PIG	14	-2057.5157	4143.6314	0.0000	1.0000
K2PI	13	-2058.7332	4143.9857	0.0000	1.0000
HKY	15	-2059.2319	4149.1505	0.0000	1.0000
SYM	16	-2063.9624	4160.7042	0.0000	1.0000
K2P	12	-2071.2053	4166.8550	0.0000	1.0000
F81G	15	-2092.3347	4215.3561	0.0000	1.0000
F81IG	16	-2092.2825	4217.3443	0.0000	1.0000
F81I	15	-2094.4119	4219.5105	0.0000	1.0000
JC69G	12	-2103.2657	4230.9759	0.0000	1.0000
JC69I	12	-2104.0431	4232.5307	0.0000	1.0000
JC69IG	13	-2103.2175	4232.9542	0.0000	1.0000
F81	14	-2106.6483	4241.8966	0.0000	1.0000
JC69	11	-2116.1740	4254.7235	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2035.5811	4111.1621	0.5693	0.5693
GTRIG	21	-2035.5128	4113.0256	0.2242	0.7935
GTRI	20	-2036.6095	4113.2191	0.2035	0.9970
HKYG	16	-2045.3703	4122.7406	0.0017	0.9987
HKYIG	17	-2045.3097	4124.6194	0.0007	0.9994
HKYI	16	-2046.5005	4125.0010	0.0006	1.0000
SYMG	17	-2050.4887	4134.9775	0.0000	1.0000
SYMIG	18	-2050.4345	4136.8689	0.0000	1.0000
SYMI	17	-2051.4716	4136.9431	0.0000	1.0000
GTR	19	-2050.0678	4138.1355	0.0000	1.0000
K2PG	13	-2057.5696	4141.1391	0.0000	1.0000
K2PIG	14	-2057.5157	4143.0314	0.0000	1.0000
K2PI	13	-2058.7332	4143.4665	0.0000	1.0000
HKY	15	-2059.2319	4148.4638	0.0000	1.0000
SYM	16	-2063.9624	4159.9249	0.0000	1.0000
K2P	12	-2071.2053	4166.4105	0.0000	1.0000
F81G	15	-2092.3347	4214.6694	0.0000	1.0000
F81IG	16	-2092.2825	4216.5649	0.0000	1.0000
F81I	15	-2094.4119	4218.8238	0.0000	1.0000
JC69G	12	-2103.2657	4230.5314	0.0000	1.0000
JC69I	12	-2104.0431	4232.0862	0.0000	1.0000
JC69IG	13	-2103.2175	4232.4350	0.0000	1.0000
F81	14	-2106.6483	4241.2966	0.0000	1.0000
JC69	11	-2116.1740	4254.3480	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2045.3703	4195.8971	0.6485	0.6485
HKYI	16	-2046.5005	4198.1575	0.2095	0.8580
K2PG	13	-2057.5696	4200.5788	0.0624	0.9204
HKYIG	17	-2045.3097	4202.3482	0.0258	0.9462
GTRG	20	-2035.5811	4202.6078	0.0226	0.9688
K2PI	13	-2058.7332	4202.9062	0.0195	0.9883
GTRI	20	-2036.6095	4204.6647	0.0081	0.9964
K2PIG	14	-2057.5157	4207.0434	0.0025	0.9989
GTRIG	21	-2035.5128	4209.0436	0.0009	0.9998
SYMG	17	-2050.4887	4212.7063	0.0001	0.9999
SYMI	17	-2051.4716	4214.6719	0.0001	1.0000
HKY	15	-2059.2319	4217.0480	0.0000	1.0000
SYMIG	18	-2050.4345	4219.1700	0.0000	1.0000
K2P	12	-2071.2053	4221.2779	0.0000	1.0000
GTR	19	-2050.0678	4225.0089	0.0000	1.0000
SYM	16	-2063.9624	4233.0814	0.0000	1.0000
F81G	15	-2092.3347	4283.2536	0.0000	1.0000
JC69G	12	-2103.2657	4285.3988	0.0000	1.0000
JC69I	12	-2104.0431	4286.9536	0.0000	1.0000
F81I	15	-2094.4119	4287.4081	0.0000	1.0000
F81IG	16	-2092.2825	4289.7214	0.0000	1.0000
JC69IG	13	-2103.2175	4291.8746	0.0000	1.0000
JC69	11	-2116.1740	4304.6431	0.0000	1.0000
F81	14	-2106.6483	4305.3085	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 591-tomodel-temp-phylip/chr10_4411.phylip = 1 - 715;
 Partition Dummy = 1:591-tomodel-temp-phylip/chr10_4411.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 591-tomodel-temp-phylip/chr10_4411.phylip = 1 - 715;
 Partition Dummy = 1:591-tomodel-temp-phylip/chr10_4411.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 591-tomodel-temp-phylip/chr10_4411.phylip = 1 - 715;
 Partition Dummy = 1:591-tomodel-temp-phylip/chr10_4411.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

