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

Input data from file "591-tomodel-phylip/chr7_2589.phylip" (7 taxa, 834 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	-2363.0665	4758.7989	0.4351	0.4351
HKYI	16	-2363.1790	4759.0239	0.3888	0.8240
HKYIG	17	-2363.0748	4760.8997	0.1522	0.9762
GTRG	20	-2362.6544	4766.3421	0.0100	0.9862
GTRI	20	-2362.6698	4766.3728	0.0099	0.9961
GTRIG	21	-2362.5372	4768.2124	0.0039	1.0000
HKY	15	-2383.8759	4798.3386	0.0000	1.0000
GTR	19	-2383.6124	4806.1585	0.0000	1.0000
F81I	15	-2402.2177	4835.0223	0.0000	1.0000
F81G	15	-2402.2510	4835.0887	0.0000	1.0000
F81IG	16	-2402.2610	4837.1878	0.0000	1.0000
F81	14	-2421.0236	4870.5601	0.0000	1.0000
SYMI	17	-2432.4336	4899.6173	0.0000	1.0000
SYMG	17	-2432.9535	4900.6570	0.0000	1.0000
K2PG	13	-2437.5259	4901.4958	0.0000	1.0000
K2PI	13	-2437.7265	4901.8970	0.0000	1.0000
SYMIG	18	-2432.9579	4902.7551	0.0000	1.0000
K2PIG	14	-2437.5403	4903.5934	0.0000	1.0000
K2P	12	-2455.5465	4935.4731	0.0000	1.0000
SYM	16	-2451.6456	4935.9571	0.0000	1.0000
JC69G	12	-2466.3209	4957.0219	0.0000	1.0000
JC69I	12	-2466.5407	4957.4615	0.0000	1.0000
JC69IG	13	-2466.3398	4959.1234	0.0000	1.0000
JC69	11	-2482.7927	4987.9065	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2363.0665	4758.1331	0.4301	0.4301
HKYI	16	-2363.1790	4758.3580	0.3844	0.8145
HKYIG	17	-2363.0748	4760.1497	0.1569	0.9715
GTRG	20	-2362.6544	4765.3088	0.0119	0.9834
GTRI	20	-2362.6698	4765.3396	0.0117	0.9951
GTRIG	21	-2362.5372	4767.0744	0.0049	1.0000
HKY	15	-2383.8759	4797.7518	0.0000	1.0000
GTR	19	-2383.6124	4805.2249	0.0000	1.0000
F81I	15	-2402.2177	4834.4355	0.0000	1.0000
F81G	15	-2402.2510	4834.5019	0.0000	1.0000
F81IG	16	-2402.2610	4836.5219	0.0000	1.0000
F81	14	-2421.0236	4870.0472	0.0000	1.0000
SYMI	17	-2432.4336	4898.8673	0.0000	1.0000
SYMG	17	-2432.9535	4899.9070	0.0000	1.0000
K2PG	13	-2437.5259	4901.0519	0.0000	1.0000
K2PI	13	-2437.7265	4901.4531	0.0000	1.0000
SYMIG	18	-2432.9579	4901.9158	0.0000	1.0000
K2PIG	14	-2437.5403	4903.0806	0.0000	1.0000
K2P	12	-2455.5465	4935.0931	0.0000	1.0000
SYM	16	-2451.6456	4935.2912	0.0000	1.0000
JC69G	12	-2466.3209	4956.6419	0.0000	1.0000
JC69I	12	-2466.5407	4957.0815	0.0000	1.0000
JC69IG	13	-2466.3398	4958.6795	0.0000	1.0000
JC69	11	-2482.7927	4987.5853	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2363.0665	4833.7528	0.5187	0.5187
HKYI	16	-2363.1790	4833.9777	0.4635	0.9822
HKYIG	17	-2363.0748	4840.4957	0.0178	1.0000
GTRG	20	-2362.6544	4859.8335	0.0000	1.0000
GTRI	20	-2362.6698	4859.8642	0.0000	1.0000
GTRIG	21	-2362.5372	4866.3253	0.0000	1.0000
HKY	15	-2383.8759	4868.6453	0.0000	1.0000
GTR	19	-2383.6124	4895.0233	0.0000	1.0000
F81I	15	-2402.2177	4905.3290	0.0000	1.0000
F81G	15	-2402.2510	4905.3954	0.0000	1.0000
F81IG	16	-2402.2610	4912.1417	0.0000	1.0000
F81	14	-2421.0236	4936.2145	0.0000	1.0000
K2PG	13	-2437.5259	4962.4929	0.0000	1.0000
K2PI	13	-2437.7265	4962.8941	0.0000	1.0000
K2PIG	14	-2437.5403	4969.2478	0.0000	1.0000
SYMI	17	-2432.4336	4979.2132	0.0000	1.0000
SYMG	17	-2432.9535	4980.2530	0.0000	1.0000
SYMIG	18	-2432.9579	4986.9880	0.0000	1.0000
K2P	12	-2455.5465	4991.8079	0.0000	1.0000
SYM	16	-2451.6456	5010.9110	0.0000	1.0000
JC69G	12	-2466.3209	5013.3567	0.0000	1.0000
JC69I	12	-2466.5407	5013.7963	0.0000	1.0000
JC69IG	13	-2466.3398	5020.1206	0.0000	1.0000
JC69	11	-2482.7927	5039.5739	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/chr7_2589.phylip = 1 - 834;
 Partition Dummy = 1:591-tomodel-phylip/chr7_2589.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/chr7_2589.phylip = 1 - 834;
 Partition Dummy = 1:591-tomodel-phylip/chr7_2589.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/chr7_2589.phylip = 1 - 834;
 Partition Dummy = 1:591-tomodel-phylip/chr7_2589.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

