
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:18 04/04/2011

Input data from file "683-tomodel-phylip/chr7_4505.phylip" (7 taxa, 586 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-1702.7634	3438.4828	0.4282	0.4282
HKYG	16	-1702.9573	3438.8707	0.3527	0.7809
HKYIG	17	-1702.7411	3440.5596	0.1516	0.9325
GTRI	20	-1701.1029	3443.6925	0.0316	0.9641
GTRG	20	-1701.3934	3444.2734	0.0237	0.9878
GTRIG	21	-1701.0966	3445.8315	0.0109	0.9986
HKY	15	-1709.6196	3450.0812	0.0013	0.9999
GTR	19	-1708.2798	3455.9024	0.0001	1.0000
F81I	15	-1733.6245	3498.0910	0.0000	1.0000
F81G	15	-1733.7618	3498.3657	0.0000	1.0000
F81IG	16	-1733.6048	3500.1657	0.0000	1.0000
F81	14	-1740.0182	3508.7720	0.0000	1.0000
K2PI	13	-1745.4293	3517.4950	0.0000	1.0000
K2PG	13	-1745.6463	3517.9290	0.0000	1.0000
K2PIG	14	-1745.4190	3519.5735	0.0000	1.0000
SYMI	17	-1742.4600	3519.9975	0.0000	1.0000
SYMG	17	-1742.7920	3520.6615	0.0000	1.0000
SYMIG	18	-1742.4624	3522.1311	0.0000	1.0000
K2P	12	-1752.4693	3529.4832	0.0000	1.0000
SYM	16	-1749.9927	3532.9415	0.0000	1.0000
JC69I	12	-1771.5522	3567.6489	0.0000	1.0000
JC69G	12	-1771.6487	3567.8419	0.0000	1.0000
JC69IG	13	-1771.5396	3569.7156	0.0000	1.0000
JC69	11	-1778.0949	3578.6498	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1702.7634	3437.5267	0.4154	0.4154
HKYG	16	-1702.9573	3437.9146	0.3422	0.7576
HKYIG	17	-1702.7411	3439.4822	0.1563	0.9139
GTRI	20	-1701.1029	3442.2058	0.0400	0.9540
GTRG	20	-1701.3934	3442.7867	0.0299	0.9839
GTRIG	21	-1701.0966	3444.1932	0.0148	0.9987
HKY	15	-1709.6196	3449.2391	0.0012	0.9999
GTR	19	-1708.2798	3454.5596	0.0001	1.0000
F81I	15	-1733.6245	3497.2489	0.0000	1.0000
F81G	15	-1733.7618	3497.5236	0.0000	1.0000
F81IG	16	-1733.6048	3499.2097	0.0000	1.0000
F81	14	-1740.0182	3508.0364	0.0000	1.0000
K2PI	13	-1745.4293	3516.8587	0.0000	1.0000
K2PG	13	-1745.6463	3517.2926	0.0000	1.0000
K2PIG	14	-1745.4190	3518.8380	0.0000	1.0000
SYMI	17	-1742.4600	3518.9201	0.0000	1.0000
SYMG	17	-1742.7920	3519.5840	0.0000	1.0000
SYMIG	18	-1742.4624	3520.9248	0.0000	1.0000
K2P	12	-1752.4693	3528.9387	0.0000	1.0000
SYM	16	-1749.9927	3531.9854	0.0000	1.0000
JC69I	12	-1771.5522	3567.1044	0.0000	1.0000
JC69G	12	-1771.6487	3567.2974	0.0000	1.0000
JC69IG	13	-1771.5396	3569.0792	0.0000	1.0000
JC69	11	-1778.0949	3578.1898	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1702.7634	3507.4998	0.5287	0.5287
HKYG	16	-1702.9573	3507.8877	0.4355	0.9642
HKYIG	17	-1702.7411	3513.8286	0.0223	0.9865
HKY	15	-1709.6196	3514.8389	0.0135	1.0000
GTRI	20	-1701.1029	3529.6722	0.0000	1.0000
GTRG	20	-1701.3934	3530.2531	0.0000	1.0000
GTRIG	21	-1701.0966	3536.0329	0.0000	1.0000
GTR	19	-1708.2798	3537.6527	0.0000	1.0000
F81I	15	-1733.6245	3562.8487	0.0000	1.0000
F81G	15	-1733.7618	3563.1234	0.0000	1.0000
F81IG	16	-1733.6048	3569.1828	0.0000	1.0000
F81	14	-1740.0182	3569.2629	0.0000	1.0000
K2PI	13	-1745.4293	3573.7118	0.0000	1.0000
K2PG	13	-1745.6463	3574.1458	0.0000	1.0000
K2PIG	14	-1745.4190	3580.0645	0.0000	1.0000
K2P	12	-1752.4693	3581.4185	0.0000	1.0000
SYMI	17	-1742.4600	3593.2665	0.0000	1.0000
SYMG	17	-1742.7920	3593.9305	0.0000	1.0000
SYMIG	18	-1742.4624	3599.6446	0.0000	1.0000
SYM	16	-1749.9927	3601.9586	0.0000	1.0000
JC69I	12	-1771.5522	3619.5842	0.0000	1.0000
JC69G	12	-1771.6487	3619.7772	0.0000	1.0000
JC69IG	13	-1771.5396	3625.9324	0.0000	1.0000
JC69	11	-1778.0949	3626.2964	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_4505.phylip = 1 - 586;
 Partition Dummy = 1:683-tomodel-phylip/chr7_4505.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_4505.phylip = 1 - 586;
 Partition Dummy = 1:683-tomodel-phylip/chr7_4505.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_4505.phylip = 1 - 586;
 Partition Dummy = 1:683-tomodel-phylip/chr7_4505.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

