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

Input data from file "591-tomodel-phylip/chr1_32424.phylip" (7 taxa, 717 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	-2214.9965	4462.7702	0.4336	0.4336
HKYIG	17	-2214.9648	4464.8050	0.1568	0.5904
GTRG	20	-2211.8132	4464.8333	0.1546	0.7450
HKYI	16	-2216.0640	4464.9051	0.1491	0.8941
GTRIG	21	-2211.7816	4466.8927	0.0552	0.9493
GTRI	20	-2212.9280	4467.0630	0.0507	1.0000
HKY	15	-2233.5336	4497.7519	0.0000	1.0000
GTR	19	-2231.4962	4502.0828	0.0000	1.0000
F81G	15	-2269.8285	4570.3417	0.0000	1.0000
F81I	15	-2270.4863	4571.6573	0.0000	1.0000
F81IG	16	-2269.8183	4572.4137	0.0000	1.0000
K2PG	13	-2280.8768	4588.2714	0.0000	1.0000
K2PIG	14	-2280.8451	4590.2884	0.0000	1.0000
K2PI	13	-2281.9558	4590.4293	0.0000	1.0000
SYMG	17	-2279.9309	4594.7373	0.0000	1.0000
SYMIG	18	-2279.9020	4596.7840	0.0000	1.0000
SYMI	17	-2281.1502	4597.1759	0.0000	1.0000
F81	14	-2284.9553	4598.5088	0.0000	1.0000
K2P	12	-2300.8265	4626.0961	0.0000	1.0000
SYM	16	-2299.1565	4631.0902	0.0000	1.0000
JC69G	12	-2323.8124	4672.0680	0.0000	1.0000
JC69I	12	-2324.5461	4673.5354	0.0000	1.0000
JC69IG	13	-2323.7964	4674.1105	0.0000	1.0000
JC69	11	-2340.1618	4702.6981	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2214.9965	4461.9930	0.4035	0.4035
GTRG	20	-2211.8132	4463.6264	0.1783	0.5818
HKYIG	17	-2214.9648	4463.9295	0.1532	0.7351
HKYI	16	-2216.0640	4464.1280	0.1388	0.8738
GTRIG	21	-2211.7816	4465.5632	0.0677	0.9415
GTRI	20	-2212.9280	4465.8561	0.0585	1.0000
HKY	15	-2233.5336	4497.0672	0.0000	1.0000
GTR	19	-2231.4962	4500.9925	0.0000	1.0000
F81G	15	-2269.8285	4569.6570	0.0000	1.0000
F81I	15	-2270.4863	4570.9726	0.0000	1.0000
F81IG	16	-2269.8183	4571.6366	0.0000	1.0000
K2PG	13	-2280.8768	4587.7536	0.0000	1.0000
K2PIG	14	-2280.8451	4589.6901	0.0000	1.0000
K2PI	13	-2281.9558	4589.9116	0.0000	1.0000
SYMG	17	-2279.9309	4593.8617	0.0000	1.0000
SYMIG	18	-2279.9020	4595.8041	0.0000	1.0000
SYMI	17	-2281.1502	4596.3004	0.0000	1.0000
F81	14	-2284.9553	4597.9105	0.0000	1.0000
K2P	12	-2300.8265	4625.6529	0.0000	1.0000
SYM	16	-2299.1565	4630.3130	0.0000	1.0000
JC69G	12	-2323.8124	4671.6248	0.0000	1.0000
JC69I	12	-2324.5461	4673.0922	0.0000	1.0000
JC69IG	13	-2323.7964	4673.5927	0.0000	1.0000
JC69	11	-2340.1618	4702.3236	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2214.9965	4535.1942	0.7233	0.7233
HKYI	16	-2216.0640	4537.3292	0.2487	0.9721
HKYIG	17	-2214.9648	4541.7058	0.0279	1.0000
GTRG	20	-2211.8132	4555.1280	0.0000	1.0000
GTRI	20	-2212.9280	4557.3576	0.0000	1.0000
GTRIG	21	-2211.7816	4561.6398	0.0000	1.0000
HKY	15	-2233.5336	4565.6933	0.0000	1.0000
GTR	19	-2231.4962	4587.9189	0.0000	1.0000
F81G	15	-2269.8285	4638.2831	0.0000	1.0000
F81I	15	-2270.4863	4639.5987	0.0000	1.0000
F81IG	16	-2269.8183	4644.8378	0.0000	1.0000
K2PG	13	-2280.8768	4647.2296	0.0000	1.0000
K2PI	13	-2281.9558	4649.3875	0.0000	1.0000
K2PIG	14	-2280.8451	4653.7412	0.0000	1.0000
F81	14	-2284.9553	4661.9616	0.0000	1.0000
SYMG	17	-2279.9309	4671.6380	0.0000	1.0000
SYMI	17	-2281.1502	4674.0767	0.0000	1.0000
SYMIG	18	-2279.9020	4678.1554	0.0000	1.0000
K2P	12	-2300.8265	4680.5538	0.0000	1.0000
SYM	16	-2299.1565	4703.5142	0.0000	1.0000
JC69G	12	-2323.8124	4726.5257	0.0000	1.0000
JC69I	12	-2324.5461	4727.9931	0.0000	1.0000
JC69IG	13	-2323.7964	4733.0687	0.0000	1.0000
JC69	11	-2340.1618	4752.6495	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_32424.phylip = 1 - 717;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32424.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_32424.phylip = 1 - 717;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32424.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_32424.phylip = 1 - 717;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32424.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

