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

Input data from file "591-tomodel-phylip/chr1_14908.phylip" (7 taxa, 695 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	-2246.7905	4526.3833	0.4312	0.4312
HKYI	16	-2247.3339	4527.4701	0.2504	0.6816
HKYIG	17	-2246.7330	4528.3700	0.1597	0.8413
GTRG	20	-2244.2052	4529.6567	0.0839	0.9252
GTRI	20	-2244.8453	4530.9368	0.0442	0.9695
GTRIG	21	-2244.1535	4531.6799	0.0305	1.0000
HKY	15	-2260.7144	4552.1358	0.0000	1.0000
GTR	19	-2258.2817	4555.6893	0.0000	1.0000
SYMG	17	-2282.9322	4600.7683	0.0000	1.0000
SYMIG	18	-2282.9256	4602.8631	0.0000	1.0000
SYMI	17	-2284.2087	4603.3214	0.0000	1.0000
K2PG	13	-2290.5399	4607.6142	0.0000	1.0000
K2PI	13	-2291.1467	4608.8280	0.0000	1.0000
K2PIG	14	-2290.4976	4609.6128	0.0000	1.0000
SYM	16	-2298.3449	4629.4922	0.0000	1.0000
K2P	12	-2304.8502	4634.1580	0.0000	1.0000
F81G	15	-2313.7597	4658.2263	0.0000	1.0000
F81I	15	-2313.8526	4658.4121	0.0000	1.0000
F81IG	16	-2313.5866	4659.9755	0.0000	1.0000
F81	14	-2325.6759	4679.9694	0.0000	1.0000
JC69G	12	-2352.6950	4729.8475	0.0000	1.0000
JC69I	12	-2353.0095	4730.4764	0.0000	1.0000
JC69IG	13	-2352.6130	4731.7604	0.0000	1.0000
JC69	11	-2364.8149	4752.0164	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2246.7905	4525.5810	0.4106	0.4106
HKYI	16	-2247.3339	4526.6678	0.2384	0.6490
HKYIG	17	-2246.7330	4527.4661	0.1600	0.8090
GTRG	20	-2244.2052	4528.4104	0.0998	0.9087
GTRI	20	-2244.8453	4529.6905	0.0526	0.9613
GTRIG	21	-2244.1535	4530.3070	0.0386	1.0000
HKY	15	-2260.7144	4551.4288	0.0000	1.0000
GTR	19	-2258.2817	4554.5634	0.0000	1.0000
SYMG	17	-2282.9322	4599.8643	0.0000	1.0000
SYMIG	18	-2282.9256	4601.8512	0.0000	1.0000
SYMI	17	-2284.2087	4602.4174	0.0000	1.0000
K2PG	13	-2290.5399	4607.0797	0.0000	1.0000
K2PI	13	-2291.1467	4608.2935	0.0000	1.0000
K2PIG	14	-2290.4976	4608.9951	0.0000	1.0000
SYM	16	-2298.3449	4628.6898	0.0000	1.0000
K2P	12	-2304.8502	4633.7005	0.0000	1.0000
F81G	15	-2313.7597	4657.5193	0.0000	1.0000
F81I	15	-2313.8526	4657.7051	0.0000	1.0000
F81IG	16	-2313.5866	4659.1731	0.0000	1.0000
F81	14	-2325.6759	4679.3517	0.0000	1.0000
JC69G	12	-2352.6950	4729.3900	0.0000	1.0000
JC69I	12	-2353.0095	4730.0189	0.0000	1.0000
JC69IG	13	-2352.6130	4731.2259	0.0000	1.0000
JC69	11	-2364.8149	4751.6299	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2246.7905	4598.2836	0.6169	0.6169
HKYI	16	-2247.3339	4599.3703	0.3583	0.9752
HKYIG	17	-2246.7330	4604.7126	0.0248	1.0000
GTRG	20	-2244.2052	4619.2887	0.0000	1.0000
HKY	15	-2260.7144	4619.5875	0.0000	1.0000
GTRI	20	-2244.8453	4620.5688	0.0000	1.0000
GTRIG	21	-2244.1535	4625.7291	0.0000	1.0000
GTR	19	-2258.2817	4640.8977	0.0000	1.0000
K2PG	13	-2290.5399	4666.1506	0.0000	1.0000
K2PI	13	-2291.1467	4667.3643	0.0000	1.0000
K2PIG	14	-2290.4976	4672.6099	0.0000	1.0000
SYMG	17	-2282.9322	4677.1108	0.0000	1.0000
SYMI	17	-2284.2087	4679.6639	0.0000	1.0000
SYMIG	18	-2282.9256	4683.6417	0.0000	1.0000
K2P	12	-2304.8502	4688.2274	0.0000	1.0000
SYM	16	-2298.3449	4701.3924	0.0000	1.0000
F81G	15	-2313.7597	4725.6780	0.0000	1.0000
F81I	15	-2313.8526	4725.8638	0.0000	1.0000
F81IG	16	-2313.5866	4731.8757	0.0000	1.0000
F81	14	-2325.6759	4742.9665	0.0000	1.0000
JC69G	12	-2352.6950	4783.9169	0.0000	1.0000
JC69I	12	-2353.0095	4784.5459	0.0000	1.0000
JC69IG	13	-2352.6130	4790.2968	0.0000	1.0000
JC69	11	-2364.8149	4801.6129	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/chr1_14908.phylip = 1 - 695;
 Partition Dummy = 1:591-tomodel-phylip/chr1_14908.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/chr1_14908.phylip = 1 - 695;
 Partition Dummy = 1:591-tomodel-phylip/chr1_14908.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/chr1_14908.phylip = 1 - 695;
 Partition Dummy = 1:591-tomodel-phylip/chr1_14908.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

