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

Input data from file "591-tomodel-phylip/chr8_5150.phylip" (7 taxa, 708 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	21	-2234.9274	4513.2017	0.3995	0.3995
HKYG	16	-2240.3087	4513.4047	0.3610	0.7605
HKYIG	17	-2240.3088	4515.5046	0.1263	0.8868
GTRG	20	-2237.6669	4516.5565	0.0747	0.9615
HKYI	16	-2242.7080	4518.2032	0.0328	0.9943
GTRI	20	-2240.5422	4522.3070	0.0042	0.9985
HKY	15	-2246.8746	4524.4428	0.0014	0.9999
GTR	19	-2245.4459	4529.9965	0.0001	1.0000
K2PG	13	-2264.1885	4554.9015	0.0000	1.0000
K2PIG	14	-2264.1886	4556.9833	0.0000	1.0000
SYMIG	18	-2260.8612	4558.7151	0.0000	1.0000
K2PI	13	-2266.3817	4559.2879	0.0000	1.0000
SYMG	17	-2262.8927	4560.6724	0.0000	1.0000
SYMI	17	-2265.2221	4565.3312	0.0000	1.0000
K2P	12	-2270.9272	4566.3033	0.0000	1.0000
SYM	16	-2270.5559	4573.8991	0.0000	1.0000
F81G	15	-2281.5498	4593.7931	0.0000	1.0000
F81IG	16	-2281.5499	4595.8870	0.0000	1.0000
F81I	15	-2283.6455	4597.9846	0.0000	1.0000
F81	14	-2287.8252	4604.2564	0.0000	1.0000
JC69G	12	-2301.1140	4626.6769	0.0000	1.0000
JC69IG	13	-2301.1141	4628.7527	0.0000	1.0000
JC69I	12	-2303.1514	4630.7517	0.0000	1.0000
JC69	11	-2307.5532	4637.4857	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	21	-2234.9274	4511.8548	0.4578	0.4578
HKYG	16	-2240.3087	4512.6174	0.3126	0.7704
HKYIG	17	-2240.3088	4514.6177	0.1150	0.8854
GTRG	20	-2237.6669	4515.3338	0.0804	0.9658
HKYI	16	-2242.7080	4517.4160	0.0284	0.9942
GTRI	20	-2240.5422	4521.0843	0.0045	0.9987
HKY	15	-2246.8746	4523.7491	0.0012	0.9999
GTR	19	-2245.4459	4528.8918	0.0001	1.0000
K2PG	13	-2264.1885	4554.3770	0.0000	1.0000
K2PIG	14	-2264.1886	4556.3772	0.0000	1.0000
SYMIG	18	-2260.8612	4557.7224	0.0000	1.0000
K2PI	13	-2266.3817	4558.7634	0.0000	1.0000
SYMG	17	-2262.8927	4559.7855	0.0000	1.0000
SYMI	17	-2265.2221	4564.4443	0.0000	1.0000
K2P	12	-2270.9272	4565.8544	0.0000	1.0000
SYM	16	-2270.5559	4573.1118	0.0000	1.0000
F81G	15	-2281.5498	4593.0995	0.0000	1.0000
F81IG	16	-2281.5499	4595.0998	0.0000	1.0000
F81I	15	-2283.6455	4597.2910	0.0000	1.0000
F81	14	-2287.8252	4603.6503	0.0000	1.0000
JC69G	12	-2301.1140	4626.2279	0.0000	1.0000
JC69IG	13	-2301.1141	4628.2282	0.0000	1.0000
JC69I	12	-2303.1514	4630.3028	0.0000	1.0000
JC69	11	-2307.5532	4637.1064	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2240.3087	4585.6165	0.8577	0.8577
HKYI	16	-2242.7080	4590.4151	0.0779	0.9356
HKYIG	17	-2240.3088	4592.1792	0.0322	0.9678
HKY	15	-2246.8746	4592.1858	0.0321	1.0000
GTRG	20	-2237.6669	4606.5827	0.0000	1.0000
GTRIG	21	-2234.9274	4607.6661	0.0000	1.0000
GTRI	20	-2240.5422	4612.3332	0.0000	1.0000
K2PG	13	-2264.1885	4613.6888	0.0000	1.0000
GTR	19	-2245.4459	4615.5783	0.0000	1.0000
K2PI	13	-2266.3817	4618.0752	0.0000	1.0000
K2PIG	14	-2264.1886	4620.2514	0.0000	1.0000
K2P	12	-2270.9272	4620.6037	0.0000	1.0000
SYMG	17	-2262.8927	4637.3470	0.0000	1.0000
SYMIG	18	-2260.8612	4639.8464	0.0000	1.0000
SYMI	17	-2265.2221	4642.0058	0.0000	1.0000
SYM	16	-2270.5559	4646.1109	0.0000	1.0000
F81G	15	-2281.5498	4661.5362	0.0000	1.0000
F81I	15	-2283.6455	4665.7276	0.0000	1.0000
F81	14	-2287.8252	4667.5245	0.0000	1.0000
F81IG	16	-2281.5499	4668.0989	0.0000	1.0000
JC69G	12	-2301.1140	4680.9773	0.0000	1.0000
JC69I	12	-2303.1514	4685.0521	0.0000	1.0000
JC69	11	-2307.5532	4687.2933	0.0000	1.0000
JC69IG	13	-2301.1141	4687.5400	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr8_5150.phylip = 1 - 708;
 Partition Dummy = 1:591-tomodel-phylip/chr8_5150.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr8_5150.phylip = 1 - 708;
 Partition Dummy = 1:591-tomodel-phylip/chr8_5150.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr8_5150.phylip = 1 - 708;
 Partition Dummy = 1:591-tomodel-phylip/chr8_5150.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

