
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:34 02/17/2011

Input data from file "183-phylip/chr13_2870.phylip" (29 taxa, 462 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-2022.5591	4183.3726	0.6764	0.6764
K2PG	57	-2027.6937	4185.7537	0.2057	0.8820
GTRG	64	-2019.0485	4187.0542	0.1073	0.9893
SYMG	61	-2026.0706	4193.0512	0.0054	0.9947
HKYIG	61	-2026.4805	4193.8710	0.0036	0.9983
K2PIG	58	-2031.5123	4196.0072	0.0012	0.9995
GTRIG	65	-2023.0701	4197.8068	0.0005	1.0000
SYMIG	62	-2029.8816	4203.3421	0.0000	1.0000
HKYI	60	-2040.6091	4219.4726	0.0000	1.0000
K2PI	57	-2045.1876	4220.7416	0.0000	1.0000
GTRI	64	-2038.1368	4225.2308	0.0000	1.0000
SYMI	61	-2044.2423	4229.3947	0.0000	1.0000
F81IG	60	-2104.1346	4346.5235	0.0000	1.0000
JC69IG	57	-2108.2816	4346.9295	0.0000	1.0000
F81G	59	-2107.2263	4350.0645	0.0000	1.0000
JC69G	56	-2111.4789	4350.7208	0.0000	1.0000
HKY	59	-2116.6050	4368.8219	0.0000	1.0000
F81I	59	-2116.6482	4368.9083	0.0000	1.0000
JC69I	56	-2121.0829	4369.9289	0.0000	1.0000
GTR	63	-2112.2186	4370.6985	0.0000	1.0000
K2P	56	-2123.8515	4375.4660	0.0000	1.0000
SYM	60	-2121.1555	4380.5653	0.0000	1.0000
F81	58	-2191.9916	4516.9658	0.0000	1.0000
JC69	55	-2197.0548	4519.2820	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-2022.5591	4165.1183	0.5699	0.5699
GTRG	64	-2019.0485	4166.0970	0.3494	0.9193
K2PG	57	-2027.6937	4169.3873	0.0674	0.9867
SYMG	61	-2026.0706	4174.1412	0.0063	0.9929
HKYIG	61	-2026.4805	4174.9610	0.0042	0.9971
GTRIG	65	-2023.0701	4176.1401	0.0023	0.9994
K2PIG	58	-2031.5123	4179.0245	0.0005	0.9999
SYMIG	62	-2029.8816	4183.7632	0.0001	1.0000
HKYI	60	-2040.6091	4201.2182	0.0000	1.0000
GTRI	64	-2038.1368	4204.2736	0.0000	1.0000
K2PI	57	-2045.1876	4204.3753	0.0000	1.0000
SYMI	61	-2044.2423	4210.4847	0.0000	1.0000
F81IG	60	-2104.1346	4328.2691	0.0000	1.0000
JC69IG	57	-2108.2816	4330.5632	0.0000	1.0000
F81G	59	-2107.2263	4332.4526	0.0000	1.0000
JC69G	56	-2111.4789	4334.9578	0.0000	1.0000
GTR	63	-2112.2186	4350.4372	0.0000	1.0000
HKY	59	-2116.6050	4351.2100	0.0000	1.0000
F81I	59	-2116.6482	4351.2964	0.0000	1.0000
JC69I	56	-2121.0829	4354.1659	0.0000	1.0000
K2P	56	-2123.8515	4359.7030	0.0000	1.0000
SYM	60	-2121.1555	4362.3109	0.0000	1.0000
F81	58	-2191.9916	4499.9832	0.0000	1.0000
JC69	55	-2197.0548	4504.1095	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-2027.6937	4405.1145	0.9822	0.9822
HKYG	60	-2022.5591	4413.2522	0.0168	0.9990
K2PIG	58	-2031.5123	4418.8873	0.0010	1.0000
SYMG	61	-2026.0706	4426.4107	0.0000	1.0000
HKYIG	61	-2026.4805	4427.2305	0.0000	1.0000
GTRG	64	-2019.0485	4430.7732	0.0000	1.0000
K2PI	57	-2045.1876	4440.1025	0.0000	1.0000
SYMIG	62	-2029.8816	4440.1682	0.0000	1.0000
GTRIG	65	-2023.0701	4444.9518	0.0000	1.0000
HKYI	60	-2040.6091	4449.3521	0.0000	1.0000
SYMI	61	-2044.2423	4462.7541	0.0000	1.0000
GTRI	64	-2038.1368	4468.9497	0.0000	1.0000
JC69IG	57	-2108.2816	4566.2904	0.0000	1.0000
JC69G	56	-2111.4789	4566.5495	0.0000	1.0000
F81IG	60	-2104.1346	4576.4030	0.0000	1.0000
F81G	59	-2107.2263	4576.4509	0.0000	1.0000
JC69I	56	-2121.0829	4585.7575	0.0000	1.0000
K2P	56	-2123.8515	4591.2947	0.0000	1.0000
HKY	59	-2116.6050	4595.2083	0.0000	1.0000
F81I	59	-2116.6482	4595.2947	0.0000	1.0000
SYM	60	-2121.1555	4610.4448	0.0000	1.0000
GTR	63	-2112.2186	4610.9778	0.0000	1.0000
JC69	55	-2197.0548	4731.5656	0.0000	1.0000
F81	58	-2191.9916	4739.8459	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 183-phylip/chr13_2870.phylip = 1 - 462;
 Partition Dummy = 1:183-phylip/chr13_2870.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 183-phylip/chr13_2870.phylip = 1 - 462;
 Partition Dummy = 1:183-phylip/chr13_2870.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 (K2PG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr13_2870.phylip = 1 - 462;
 Partition Dummy = 1:183-phylip/chr13_2870.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

