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

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

Minimum AIC  model: SYMG
Minimum AICc model: SYMG
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
SYMG	17	-2165.8560	4366.6255	0.4433	0.4433
GTRG	20	-2163.6961	4368.6515	0.1610	0.6043
SYMIG	18	-2165.8525	4368.7275	0.1550	0.7593
SYMI	17	-2167.0427	4368.9988	0.1353	0.8946
GTRIG	21	-2163.6950	4370.7775	0.0556	0.9502
GTRI	20	-2164.8732	4371.0058	0.0496	0.9998
K2PG	13	-2178.3057	4383.1514	0.0001	0.9999
K2PIG	14	-2178.3005	4385.2251	0.0000	0.9999
K2PI	13	-2179.5609	4385.6619	0.0000	1.0000
HKYG	16	-2177.0309	4386.8726	0.0000	1.0000
HKYIG	17	-2177.0199	4388.9532	0.0000	1.0000
HKYI	16	-2178.2272	4389.2652	0.0000	1.0000
SYM	16	-2179.4826	4391.7760	0.0000	1.0000
GTR	19	-2176.8892	4392.9161	0.0000	1.0000
JC69G	12	-2187.8239	4400.1101	0.0000	1.0000
JC69I	12	-2188.7047	4401.8716	0.0000	1.0000
JC69IG	13	-2187.8287	4402.1975	0.0000	1.0000
F81G	15	-2186.6315	4403.9773	0.0000	1.0000
F81I	15	-2187.5228	4405.7599	0.0000	1.0000
F81IG	16	-2186.6267	4406.0641	0.0000	1.0000
K2P	12	-2194.6097	4413.6816	0.0000	1.0000
HKY	15	-2193.2764	4417.2670	0.0000	1.0000
JC69	11	-2203.4543	4429.2991	0.0000	1.0000
F81	14	-2202.3217	4433.2675	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	17	-2165.8560	4365.7120	0.4169	0.4169
GTRG	20	-2163.6961	4367.3921	0.1800	0.5969
SYMIG	18	-2165.8525	4367.7051	0.1539	0.7508
SYMI	17	-2167.0427	4368.0853	0.1273	0.8781
GTRIG	21	-2163.6950	4369.3901	0.0663	0.9444
GTRI	20	-2164.8732	4369.7464	0.0555	0.9998
K2PG	13	-2178.3057	4382.6114	0.0001	0.9999
K2PIG	14	-2178.3005	4384.6010	0.0000	0.9999
K2PI	13	-2179.5609	4385.1218	0.0000	1.0000
HKYG	16	-2177.0309	4386.0619	0.0000	1.0000
HKYIG	17	-2177.0199	4388.0398	0.0000	1.0000
HKYI	16	-2178.2272	4388.4545	0.0000	1.0000
SYM	16	-2179.4826	4390.9653	0.0000	1.0000
GTR	19	-2176.8892	4391.7784	0.0000	1.0000
JC69G	12	-2187.8239	4399.6479	0.0000	1.0000
JC69I	12	-2188.7047	4401.4094	0.0000	1.0000
JC69IG	13	-2187.8287	4401.6574	0.0000	1.0000
F81G	15	-2186.6315	4403.2631	0.0000	1.0000
F81I	15	-2187.5228	4405.0456	0.0000	1.0000
F81IG	16	-2186.6267	4405.2534	0.0000	1.0000
K2P	12	-2194.6097	4413.2194	0.0000	1.0000
HKY	15	-2193.2764	4416.5527	0.0000	1.0000
JC69	11	-2203.4543	4428.9086	0.0000	1.0000
F81	14	-2202.3217	4432.6434	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2178.3057	4441.5506	0.4876	0.4876
SYMG	17	-2165.8560	4442.7865	0.2628	0.7504
K2PI	13	-2179.5609	4444.0611	0.1390	0.8894
SYMI	17	-2167.0427	4445.1597	0.0802	0.9696
K2PIG	14	-2178.3005	4448.0741	0.0187	0.9883
SYMIG	18	-2165.8525	4449.3133	0.0101	0.9983
JC69G	12	-2187.8239	4454.0533	0.0009	0.9993
JC69I	12	-2188.7047	4455.8149	0.0004	0.9997
GTRG	20	-2163.6961	4458.0679	0.0001	0.9998
HKYG	16	-2177.0309	4458.6025	0.0001	0.9999
GTRI	20	-2164.8732	4460.4222	0.0000	0.9999
JC69IG	13	-2187.8287	4460.5967	0.0000	1.0000
HKYI	16	-2178.2272	4460.9951	0.0000	1.0000
SYM	16	-2179.4826	4463.5059	0.0000	1.0000
GTRIG	21	-2163.6950	4464.5996	0.0000	1.0000
HKYIG	17	-2177.0199	4465.1142	0.0000	1.0000
K2P	12	-2194.6097	4467.6248	0.0000	1.0000
F81G	15	-2186.6315	4471.2699	0.0000	1.0000
F81I	15	-2187.5228	4473.0525	0.0000	1.0000
F81IG	16	-2186.6267	4477.7940	0.0000	1.0000
GTR	19	-2176.8892	4477.9204	0.0000	1.0000
JC69	11	-2203.4543	4478.7802	0.0000	1.0000
HKY	15	-2193.2764	4484.5596	0.0000	1.0000
F81	14	-2202.3217	4496.1164	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_9464.phylip = 1 - 688;
 Partition Dummy = 1:591-tomodel-phylip/chr6_9464.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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_9464.phylip = 1 - 688;
 Partition Dummy = 1:591-tomodel-phylip/chr6_9464.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=Fixed(Equal);
END;


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

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

