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

Input data from file "591-tomodel-phylip/chr5_10375.phylip" (7 taxa, 597 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	-2028.5767	4090.0913	0.6039	0.6039
HKYIG	17	-2028.5774	4092.2118	0.2092	0.8130
HKYI	16	-2029.9577	4092.8533	0.1518	0.9648
GTRG	20	-2027.6314	4096.7211	0.0219	0.9868
GTRIG	21	-2027.6334	4098.8737	0.0075	0.9942
GTRI	20	-2028.9711	4099.4006	0.0057	1.0000
HKY	15	-2041.0596	4112.9453	0.0000	1.0000
GTR	19	-2040.3674	4120.0520	0.0000	1.0000
K2PG	13	-2061.5496	4149.7236	0.0000	1.0000
K2PIG	14	-2061.5497	4151.8211	0.0000	1.0000
SYMG	17	-2058.7300	4152.5171	0.0000	1.0000
K2PI	13	-2063.1622	4152.9487	0.0000	1.0000
SYMIG	18	-2058.7321	4154.6477	0.0000	1.0000
SYMI	17	-2060.8771	4156.8112	0.0000	1.0000
K2P	12	-2074.3940	4173.3222	0.0000	1.0000
SYM	16	-2072.1121	4177.1622	0.0000	1.0000
F81G	15	-2085.8403	4202.5068	0.0000	1.0000
F81I	15	-2086.3213	4203.4688	0.0000	1.0000
F81IG	16	-2085.8415	4204.6209	0.0000	1.0000
F81	14	-2097.4613	4223.6443	0.0000	1.0000
JC69G	12	-2125.2526	4275.0395	0.0000	1.0000
JC69I	12	-2126.0333	4276.6008	0.0000	1.0000
JC69IG	13	-2125.2531	4277.1306	0.0000	1.0000
JC69	11	-2136.3133	4295.0780	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2028.5767	4089.1534	0.5897	0.5897
HKYIG	17	-2028.5774	4091.1548	0.2168	0.8065
HKYI	16	-2029.9577	4091.9153	0.1482	0.9547
GTRG	20	-2027.6314	4095.2628	0.0278	0.9825
GTRIG	21	-2027.6334	4097.2668	0.0102	0.9927
GTRI	20	-2028.9711	4097.9423	0.0073	1.0000
HKY	15	-2041.0596	4112.1192	0.0000	1.0000
GTR	19	-2040.3674	4118.7348	0.0000	1.0000
K2PG	13	-2061.5496	4149.0992	0.0000	1.0000
K2PIG	14	-2061.5497	4151.0995	0.0000	1.0000
SYMG	17	-2058.7300	4151.4601	0.0000	1.0000
K2PI	13	-2063.1622	4152.3244	0.0000	1.0000
SYMIG	18	-2058.7321	4153.4643	0.0000	1.0000
SYMI	17	-2060.8771	4155.7542	0.0000	1.0000
K2P	12	-2074.3940	4172.7880	0.0000	1.0000
SYM	16	-2072.1121	4176.2243	0.0000	1.0000
F81G	15	-2085.8403	4201.6807	0.0000	1.0000
F81I	15	-2086.3213	4202.6426	0.0000	1.0000
F81IG	16	-2085.8415	4203.6830	0.0000	1.0000
F81	14	-2097.4613	4222.9226	0.0000	1.0000
JC69G	12	-2125.2526	4274.5052	0.0000	1.0000
JC69I	12	-2126.0333	4276.0665	0.0000	1.0000
JC69IG	13	-2125.2531	4276.5062	0.0000	1.0000
JC69	11	-2136.3133	4294.6267	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2028.5767	4159.4241	0.7738	0.7738
HKYI	16	-2029.9577	4162.1860	0.1945	0.9683
HKYIG	17	-2028.5774	4165.8174	0.0316	0.9999
HKY	15	-2041.0596	4177.9979	0.0001	1.0000
GTRG	20	-2027.6314	4183.1011	0.0000	1.0000
GTRI	20	-2028.9711	4185.7806	0.0000	1.0000
GTRIG	21	-2027.6334	4189.4970	0.0000	1.0000
GTR	19	-2040.3674	4202.1812	0.0000	1.0000
K2PG	13	-2061.5496	4206.1941	0.0000	1.0000
K2PI	13	-2063.1622	4209.4193	0.0000	1.0000
K2PIG	14	-2061.5497	4212.5863	0.0000	1.0000
K2P	12	-2074.3940	4225.4910	0.0000	1.0000
SYMG	17	-2058.7300	4226.1227	0.0000	1.0000
SYMI	17	-2060.8771	4230.4168	0.0000	1.0000
SYMIG	18	-2058.7321	4232.5188	0.0000	1.0000
SYM	16	-2072.1121	4246.4950	0.0000	1.0000
F81G	15	-2085.8403	4267.5594	0.0000	1.0000
F81I	15	-2086.3213	4268.5214	0.0000	1.0000
F81IG	16	-2085.8415	4273.9537	0.0000	1.0000
F81	14	-2097.4613	4284.4095	0.0000	1.0000
JC69G	12	-2125.2526	4327.2082	0.0000	1.0000
JC69I	12	-2126.0333	4328.7695	0.0000	1.0000
JC69IG	13	-2125.2531	4333.6012	0.0000	1.0000
JC69	11	-2136.3133	4342.9378	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/chr5_10375.phylip = 1 - 597;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10375.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/chr5_10375.phylip = 1 - 597;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10375.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/chr5_10375.phylip = 1 - 597;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10375.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

