
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:43 04/04/2011

Input data from file "591-tomodel-phylip/chr6_7207.phylip" (7 taxa, 666 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	20	-2114.1388	4269.5799	0.5900	0.5900
GTRIG	21	-2114.1200	4271.6748	0.2070	0.7970
HKYG	16	-2119.9889	4272.8161	0.1170	0.9140
HKYIG	17	-2119.9702	4274.8849	0.0416	0.9556
GTRI	20	-2117.1076	4275.5174	0.0303	0.9859
HKYI	16	-2122.1667	4277.1716	0.0133	0.9991
K2PG	13	-2128.9618	4284.4818	0.0003	0.9995
SYMG	17	-2125.0725	4285.0895	0.0003	0.9997
K2PIG	14	-2128.9416	4286.5284	0.0001	0.9998
SYMIG	18	-2125.0557	4287.1685	0.0001	0.9999
K2PI	13	-2130.9954	4288.5492	0.0000	1.0000
SYMI	17	-2127.5754	4290.0953	0.0000	1.0000
GTR	19	-2135.7322	4310.6409	0.0000	1.0000
HKY	15	-2141.4977	4313.7339	0.0000	1.0000
F81G	15	-2145.9307	4322.5998	0.0000	1.0000
F81IG	16	-2145.9175	4324.6732	0.0000	1.0000
SYM	16	-2146.2873	4325.4128	0.0000	1.0000
K2P	12	-2150.7567	4325.9913	0.0000	1.0000
F81I	15	-2147.8211	4326.3807	0.0000	1.0000
JC69G	12	-2153.6042	4331.6862	0.0000	1.0000
JC69IG	13	-2153.5867	4333.7317	0.0000	1.0000
JC69I	12	-2155.4740	4335.4257	0.0000	1.0000
F81	14	-2166.4690	4361.5831	0.0000	1.0000
JC69	11	-2174.3444	4371.0925	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2114.1388	4268.2776	0.6019	0.6019
GTRIG	21	-2114.1200	4270.2400	0.2256	0.8276
HKYG	16	-2119.9889	4271.9779	0.0946	0.9222
HKYIG	17	-2119.9702	4273.9405	0.0355	0.9577
GTRI	20	-2117.1076	4274.2151	0.0309	0.9886
HKYI	16	-2122.1667	4276.3334	0.0107	0.9993
K2PG	13	-2128.9618	4283.9236	0.0002	0.9996
SYMG	17	-2125.0725	4284.1451	0.0002	0.9998
K2PIG	14	-2128.9416	4285.8832	0.0001	0.9999
SYMIG	18	-2125.0557	4286.1114	0.0001	1.0000
K2PI	13	-2130.9954	4287.9909	0.0000	1.0000
SYMI	17	-2127.5754	4289.1508	0.0000	1.0000
GTR	19	-2135.7322	4309.4644	0.0000	1.0000
HKY	15	-2141.4977	4312.9954	0.0000	1.0000
F81G	15	-2145.9307	4321.8613	0.0000	1.0000
F81IG	16	-2145.9175	4323.8350	0.0000	1.0000
SYM	16	-2146.2873	4324.5746	0.0000	1.0000
K2P	12	-2150.7567	4325.5135	0.0000	1.0000
F81I	15	-2147.8211	4325.6422	0.0000	1.0000
JC69G	12	-2153.6042	4331.2084	0.0000	1.0000
JC69IG	13	-2153.5867	4333.1734	0.0000	1.0000
JC69I	12	-2155.4740	4334.9479	0.0000	1.0000
F81	14	-2166.4690	4360.9380	0.0000	1.0000
JC69	11	-2174.3444	4370.6888	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2128.9618	4342.4403	0.5883	0.5883
HKYG	16	-2119.9889	4343.9985	0.2699	0.8582
K2PI	13	-2130.9954	4346.5076	0.0770	0.9352
HKYI	16	-2122.1667	4348.3540	0.0306	0.9658
K2PIG	14	-2128.9416	4348.9013	0.0233	0.9890
HKYIG	17	-2119.9702	4350.4624	0.0107	0.9997
GTRG	20	-2114.1388	4358.3034	0.0002	0.9999
SYMG	17	-2125.0725	4360.6670	0.0001	1.0000
GTRI	20	-2117.1076	4364.2409	0.0000	1.0000
GTRIG	21	-2114.1200	4364.7671	0.0000	1.0000
SYMI	17	-2127.5754	4365.6727	0.0000	1.0000
SYMIG	18	-2125.0557	4367.1346	0.0000	1.0000
K2P	12	-2150.7567	4379.5290	0.0000	1.0000
HKY	15	-2141.4977	4380.5148	0.0000	1.0000
JC69G	12	-2153.6042	4385.2239	0.0000	1.0000
JC69I	12	-2155.4740	4388.9634	0.0000	1.0000
F81G	15	-2145.9307	4389.3807	0.0000	1.0000
JC69IG	13	-2153.5867	4391.6902	0.0000	1.0000
F81I	15	-2147.8211	4393.1615	0.0000	1.0000
GTR	19	-2135.7322	4394.9889	0.0000	1.0000
F81IG	16	-2145.9175	4395.8556	0.0000	1.0000
SYM	16	-2146.2873	4396.5953	0.0000	1.0000
JC69	11	-2174.3444	4420.2030	0.0000	1.0000
F81	14	-2166.4690	4423.9560	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-phylip/chr6_7207.phylip = 1 - 666;
 Partition Dummy = 1:591-tomodel-phylip/chr6_7207.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-phylip/chr6_7207.phylip = 1 - 666;
 Partition Dummy = 1:591-tomodel-phylip/chr6_7207.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 591-tomodel-phylip/chr6_7207.phylip = 1 - 666;
 Partition Dummy = 1:591-tomodel-phylip/chr6_7207.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

