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

Input data from file "591-tomodel-phylip/chr5_3203.phylip" (7 taxa, 589 characters)

Minimum AIC  model: K2PG
Minimum AICc model: K2PG
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
K2PG	13	-2061.4701	4149.5732	0.3247	0.3247
K2PI	13	-2061.5951	4149.8233	0.2865	0.6112
K2PIG	14	-2061.3814	4151.4944	0.1242	0.7354
HKYG	16	-2059.4547	4151.8605	0.1035	0.8389
HKYI	16	-2059.4912	4151.9334	0.0998	0.9386
HKYIG	17	-2059.3479	4153.7675	0.0399	0.9785
SYMG	17	-2061.1278	4157.3273	0.0067	0.9852
SYMI	17	-2061.2728	4157.6173	0.0058	0.9910
GTRI	20	-2058.8118	4159.1024	0.0028	0.9938
GTRG	20	-2058.8178	4159.1146	0.0028	0.9965
SYMIG	18	-2061.0912	4159.3823	0.0024	0.9989
GTRIG	21	-2058.7033	4161.0363	0.0011	1.0000
K2P	12	-2075.6091	4175.7599	0.0000	1.0000
HKY	15	-2073.6062	4178.0501	0.0000	1.0000
SYM	16	-2075.0451	4183.0412	0.0000	1.0000
GTR	19	-2072.6980	4184.7316	0.0000	1.0000
JC69I	12	-2119.4643	4263.4703	0.0000	1.0000
JC69G	12	-2119.8747	4264.2910	0.0000	1.0000
JC69IG	13	-2119.4595	4265.5520	0.0000	1.0000
F81I	15	-2117.5012	4265.8402	0.0000	1.0000
F81G	15	-2118.0228	4266.8834	0.0000	1.0000
F81IG	16	-2117.5033	4267.9576	0.0000	1.0000
JC69	11	-2131.2936	4285.0447	0.0000	1.0000
F81	14	-2129.4373	4287.6064	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	13	-2061.4701	4148.9402	0.3067	0.3067
K2PI	13	-2061.5951	4149.1903	0.2707	0.5774
K2PIG	14	-2061.3814	4150.7627	0.1233	0.7007
HKYG	16	-2059.4547	4150.9095	0.1146	0.8152
HKYI	16	-2059.4912	4150.9824	0.1105	0.9257
HKYIG	17	-2059.3479	4152.6957	0.0469	0.9726
SYMG	17	-2061.1278	4156.2555	0.0079	0.9805
SYMI	17	-2061.2728	4156.5455	0.0068	0.9874
GTRI	20	-2058.8118	4157.6235	0.0040	0.9914
GTRG	20	-2058.8178	4157.6357	0.0040	0.9953
SYMIG	18	-2061.0912	4158.1823	0.0030	0.9984
GTRIG	21	-2058.7033	4159.4067	0.0016	1.0000
K2P	12	-2075.6091	4175.2182	0.0000	1.0000
HKY	15	-2073.6062	4177.2124	0.0000	1.0000
SYM	16	-2075.0451	4182.0902	0.0000	1.0000
GTR	19	-2072.6980	4183.3959	0.0000	1.0000
JC69I	12	-2119.4643	4262.9286	0.0000	1.0000
JC69G	12	-2119.8747	4263.7493	0.0000	1.0000
JC69IG	13	-2119.4595	4264.9190	0.0000	1.0000
F81I	15	-2117.5012	4265.0025	0.0000	1.0000
F81G	15	-2118.0228	4266.0457	0.0000	1.0000
F81IG	16	-2117.5033	4267.0066	0.0000	1.0000
JC69	11	-2131.2936	4284.5871	0.0000	1.0000
F81	14	-2129.4373	4286.8747	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2061.4701	4205.8597	0.5185	0.5185
K2PI	13	-2061.5951	4206.1098	0.4576	0.9761
K2PIG	14	-2061.3814	4212.0607	0.0233	0.9994
HKYG	16	-2059.4547	4220.9643	0.0003	0.9997
HKYI	16	-2059.4912	4221.0372	0.0003	1.0000
HKYIG	17	-2059.3479	4227.1290	0.0000	1.0000
K2P	12	-2075.6091	4227.7594	0.0000	1.0000
SYMG	17	-2061.1278	4230.6888	0.0000	1.0000
SYMI	17	-2061.2728	4230.9787	0.0000	1.0000
SYMIG	18	-2061.0912	4236.9940	0.0000	1.0000
HKY	15	-2073.6062	4242.8888	0.0000	1.0000
GTRI	20	-2058.8118	4245.1920	0.0000	1.0000
GTRG	20	-2058.8178	4245.2042	0.0000	1.0000
GTRIG	21	-2058.7033	4251.3536	0.0000	1.0000
SYM	16	-2075.0451	4252.1450	0.0000	1.0000
GTR	19	-2072.6980	4266.5860	0.0000	1.0000
JC69I	12	-2119.4643	4315.4698	0.0000	1.0000
JC69G	12	-2119.8747	4316.2905	0.0000	1.0000
JC69IG	13	-2119.4595	4321.8385	0.0000	1.0000
F81I	15	-2117.5012	4330.6789	0.0000	1.0000
F81G	15	-2118.0228	4331.7221	0.0000	1.0000
JC69	11	-2131.2936	4332.7498	0.0000	1.0000
F81IG	16	-2117.5033	4337.0614	0.0000	1.0000
F81	14	-2129.4373	4348.1726	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_3203.phylip = 1 - 589;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3203.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;


[Mrbayes block for the best AICc model (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_3203.phylip = 1 - 589;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3203.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;


[Mrbayes block for the best BIC model (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_3203.phylip = 1 - 589;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3203.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

