
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:02 08/01/2011

Input data from file "261-aggressive-phylip/12486_hg18.phylip" (29 taxa, 269 characters)

Minimum AIC  model: K2PG
Minimum AICc model: K2PG
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
K2PG	57	-968.6647	2082.6660	0.5720	0.5720
K2PI	57	-969.4498	2084.2361	0.2609	0.8328
K2PIG	58	-968.6535	2085.8976	0.1137	0.9465
HKYG	60	-966.8739	2088.9401	0.0248	0.9713
HKYI	60	-967.6377	2090.4677	0.0116	0.9829
SYMG	61	-966.4997	2091.5406	0.0068	0.9897
HKYIG	61	-966.8533	2092.2477	0.0048	0.9944
SYMI	61	-967.4086	2093.3582	0.0027	0.9971
SYMIG	62	-966.4874	2094.8971	0.0013	0.9984
K2P	56	-977.0926	2096.2984	0.0006	0.9990
GTRG	64	-963.8637	2096.5118	0.0006	0.9996
GTRI	64	-964.6628	2098.1100	0.0003	0.9999
GTRIG	65	-963.8505	2099.9669	0.0001	1.0000
HKY	59	-975.1163	2102.1082	0.0000	1.0000
SYM	60	-974.7229	2104.6382	0.0000	1.0000
GTR	63	-972.2058	2109.7482	0.0000	1.0000
JC69G	56	-988.3321	2118.7773	0.0000	1.0000
JC69I	56	-989.0867	2120.2866	0.0000	1.0000
JC69IG	57	-988.2956	2121.9276	0.0000	1.0000
F81G	59	-985.8116	2123.4987	0.0000	1.0000
F81I	59	-986.4368	2124.7492	0.0000	1.0000
F81IG	60	-985.7719	2126.7361	0.0000	1.0000
JC69	55	-996.6643	2132.2489	0.0000	1.0000
F81	58	-994.2140	2137.0184	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	57	-968.6647	2051.3295	0.3491	0.3491
K2PI	57	-969.4498	2052.8996	0.1592	0.5084
K2PIG	58	-968.6535	2053.3071	0.1299	0.6383
HKYG	60	-966.8739	2053.7478	0.1042	0.7425
SYMG	61	-966.4997	2054.9995	0.0557	0.7982
HKYI	60	-967.6377	2055.2754	0.0485	0.8467
HKYIG	61	-966.8533	2055.7066	0.0391	0.8859
GTRG	64	-963.8637	2055.7275	0.0387	0.9246
SYMI	61	-967.4086	2056.8171	0.0225	0.9470
SYMIG	62	-966.4874	2056.9747	0.0208	0.9678
GTRI	64	-964.6628	2057.3257	0.0174	0.9852
GTRIG	65	-963.8505	2057.7009	0.0144	0.9997
K2P	56	-977.0926	2066.1852	0.0002	0.9999
HKY	59	-975.1163	2068.2326	0.0001	0.9999
SYM	60	-974.7229	2069.4459	0.0000	1.0000
GTR	63	-972.2058	2070.4116	0.0000	1.0000
JC69G	56	-988.3321	2088.6641	0.0000	1.0000
F81G	59	-985.8116	2089.6231	0.0000	1.0000
JC69I	56	-989.0867	2090.1734	0.0000	1.0000
JC69IG	57	-988.2956	2090.5911	0.0000	1.0000
F81I	59	-986.4368	2090.8736	0.0000	1.0000
F81IG	60	-985.7719	2091.5438	0.0000	1.0000
JC69	55	-996.6643	2103.3287	0.0000	1.0000
F81	58	-994.2140	2104.4280	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-968.6647	2256.2280	0.6563	0.6563
K2PI	57	-969.4498	2257.7982	0.2994	0.9557
K2PIG	58	-968.6535	2261.8003	0.0405	0.9962
K2P	56	-977.0926	2267.4891	0.0024	0.9985
HKYG	60	-966.8739	2269.4305	0.0009	0.9994
HKYI	60	-967.6377	2270.9581	0.0004	0.9998
SYMG	61	-966.4997	2274.2769	0.0001	0.9999
HKYIG	61	-966.8533	2274.9840	0.0001	1.0000
SYMI	61	-967.4086	2276.0945	0.0000	1.0000
SYMIG	62	-966.4874	2279.8468	0.0000	1.0000
HKY	59	-975.1163	2280.3206	0.0000	1.0000
SYM	60	-974.7229	2285.1286	0.0000	1.0000
GTRG	64	-963.8637	2285.7890	0.0000	1.0000
GTRI	64	-964.6628	2287.3872	0.0000	1.0000
JC69G	56	-988.3321	2289.9680	0.0000	1.0000
GTRIG	65	-963.8505	2291.3572	0.0000	1.0000
JC69I	56	-989.0867	2291.4772	0.0000	1.0000
JC69IG	57	-988.2956	2295.4897	0.0000	1.0000
GTR	63	-972.2058	2296.8784	0.0000	1.0000
JC69	55	-996.6643	2301.0378	0.0000	1.0000
F81G	59	-985.8116	2301.7111	0.0000	1.0000
F81I	59	-986.4368	2302.9616	0.0000	1.0000
F81IG	60	-985.7719	2307.2265	0.0000	1.0000
F81	58	-994.2140	2312.9212	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/12486_hg18.phylip = 1 - 269;
 Partition Dummy = 1:261-aggressive-phylip/12486_hg18.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;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/12486_hg18.phylip = 1 - 269;
 Partition Dummy = 1:261-aggressive-phylip/12486_hg18.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;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/12486_hg18.phylip = 1 - 269;
 Partition Dummy = 1:261-aggressive-phylip/12486_hg18.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

