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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-2126.1681	4285.1705	0.4141	0.4141
HKYI	16	-2126.4595	4285.7533	0.3094	0.7235
HKYIG	17	-2126.1775	4287.2951	0.1431	0.8666
GTRG	20	-2124.1255	4289.5473	0.0464	0.9130
GTRI	20	-2124.4039	4290.1040	0.0351	0.9481
HKY	15	-2129.7812	4290.2976	0.0319	0.9800
GTRIG	21	-2124.1385	4291.7052	0.0158	0.9958
GTR	19	-2127.5949	4294.3609	0.0042	1.0000
K2PG	13	-2169.0646	4364.6850	0.0000	1.0000
K2PI	13	-2169.2585	4365.0727	0.0000	1.0000
K2PIG	14	-2169.0710	4366.7841	0.0000	1.0000
SYMG	17	-2168.5785	4372.0971	0.0000	1.0000
K2P	12	-2173.9791	4372.4338	0.0000	1.0000
SYMI	17	-2168.7735	4372.4871	0.0000	1.0000
SYMIG	18	-2168.5884	4374.2292	0.0000	1.0000
SYM	16	-2173.6431	4380.1205	0.0000	1.0000
F81G	15	-2178.4626	4387.6603	0.0000	1.0000
F81I	15	-2178.4977	4387.7305	0.0000	1.0000
F81IG	16	-2178.4555	4389.7454	0.0000	1.0000
F81	14	-2181.1565	4390.9551	0.0000	1.0000
JC69I	12	-2214.1838	4452.8431	0.0000	1.0000
JC69G	12	-2214.1888	4452.8531	0.0000	1.0000
JC69IG	13	-2214.1694	4454.8945	0.0000	1.0000
JC69	11	-2217.4988	4457.3993	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2126.1681	4284.3361	0.4006	0.4006
HKYI	16	-2126.4595	4284.9189	0.2993	0.6999
HKYIG	17	-2126.1775	4286.3550	0.1460	0.8459
GTRG	20	-2124.1255	4288.2510	0.0566	0.9025
GTRI	20	-2124.4039	4288.8077	0.0428	0.9453
HKY	15	-2129.7812	4289.5625	0.0294	0.9747
GTRIG	21	-2124.1385	4290.2771	0.0205	0.9952
GTR	19	-2127.5949	4293.1899	0.0048	1.0000
K2PG	13	-2169.0646	4364.1293	0.0000	1.0000
K2PI	13	-2169.2585	4364.5170	0.0000	1.0000
K2PIG	14	-2169.0710	4366.1419	0.0000	1.0000
SYMG	17	-2168.5785	4371.1570	0.0000	1.0000
SYMI	17	-2168.7735	4371.5470	0.0000	1.0000
K2P	12	-2173.9791	4371.9582	0.0000	1.0000
SYMIG	18	-2168.5884	4373.1769	0.0000	1.0000
SYM	16	-2173.6431	4379.2861	0.0000	1.0000
F81G	15	-2178.4626	4386.9253	0.0000	1.0000
F81I	15	-2178.4977	4386.9954	0.0000	1.0000
F81IG	16	-2178.4555	4388.9110	0.0000	1.0000
F81	14	-2181.1565	4390.3129	0.0000	1.0000
JC69I	12	-2214.1838	4452.3675	0.0000	1.0000
JC69G	12	-2214.1888	4452.3775	0.0000	1.0000
JC69IG	13	-2214.1694	4454.3387	0.0000	1.0000
JC69	11	-2217.4988	4456.9975	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2126.1681	4356.4286	0.4027	0.4027
HKYI	16	-2126.4595	4357.0115	0.3009	0.7037
HKY	15	-2129.7812	4357.1493	0.2809	0.9846
HKYIG	17	-2126.1775	4362.9533	0.0154	1.0000
GTRG	20	-2124.1255	4378.3667	0.0000	1.0000
GTR	19	-2127.5949	4378.7998	0.0000	1.0000
GTRI	20	-2124.4039	4378.9234	0.0000	1.0000
GTRIG	21	-2124.1385	4384.8985	0.0000	1.0000
K2PG	13	-2169.0646	4422.7045	0.0000	1.0000
K2PI	13	-2169.2585	4423.0922	0.0000	1.0000
K2P	12	-2173.9791	4426.0276	0.0000	1.0000
K2PIG	14	-2169.0710	4429.2229	0.0000	1.0000
SYMG	17	-2168.5785	4447.7554	0.0000	1.0000
SYMI	17	-2168.7735	4448.1453	0.0000	1.0000
SYM	16	-2173.6431	4451.3787	0.0000	1.0000
F81	14	-2181.1565	4453.3939	0.0000	1.0000
SYMIG	18	-2168.5884	4454.2810	0.0000	1.0000
F81G	15	-2178.4626	4454.5120	0.0000	1.0000
F81I	15	-2178.4977	4454.5822	0.0000	1.0000
F81IG	16	-2178.4555	4461.0036	0.0000	1.0000
JC69I	12	-2214.1838	4506.4369	0.0000	1.0000
JC69G	12	-2214.1888	4506.4469	0.0000	1.0000
JC69	11	-2217.4988	4506.5611	0.0000	1.0000
JC69IG	13	-2214.1694	4512.9139	0.0000	1.0000



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

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

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


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

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


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

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

