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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-1691.5385	3424.4341	0.4567	0.4567
GTRG	20	-1691.7786	3424.9142	0.3592	0.8159
GTRIG	21	-1691.7786	3427.0523	0.1233	0.9392
HKYI	16	-1698.5589	3429.9910	0.0284	0.9676
HKYG	16	-1698.7240	3430.3212	0.0241	0.9916
HKYIG	17	-1698.7241	3432.4320	0.0084	1.0000
GTR	19	-1707.6677	3454.5612	0.0000	1.0000
HKY	15	-1714.9604	3460.6900	0.0000	1.0000
K2PI	13	-1721.6759	3469.9333	0.0000	1.0000
K2PG	13	-1721.8635	3470.3085	0.0000	1.0000
SYMI	17	-1718.4697	3471.9234	0.0000	1.0000
K2PIG	14	-1721.8636	3472.3991	0.0000	1.0000
SYMG	17	-1718.7727	3472.5294	0.0000	1.0000
SYMIG	18	-1718.7728	3474.6470	0.0000	1.0000
F81I	15	-1726.9037	3484.5767	0.0000	1.0000
F81G	15	-1727.0259	3484.8211	0.0000	1.0000
F81IG	16	-1727.0259	3486.9251	0.0000	1.0000
SYM	16	-1734.5525	3501.9781	0.0000	1.0000
K2P	12	-1738.9973	3502.4921	0.0000	1.0000
F81	14	-1742.4483	3513.5687	0.0000	1.0000
JC69I	12	-1747.4267	3519.3511	0.0000	1.0000
JC69G	12	-1747.5672	3519.6320	0.0000	1.0000
JC69IG	13	-1747.5672	3521.7158	0.0000	1.0000
JC69	11	-1763.3012	3549.0227	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1691.5385	3423.0771	0.4584	0.4584
GTRG	20	-1691.7786	3423.5572	0.3606	0.8190
GTRIG	21	-1691.7786	3425.5572	0.1327	0.9517
HKYI	16	-1698.5589	3429.1178	0.0224	0.9741
HKYG	16	-1698.7240	3429.4480	0.0190	0.9930
HKYIG	17	-1698.7241	3431.4481	0.0070	1.0000
GTR	19	-1707.6677	3453.3354	0.0000	1.0000
HKY	15	-1714.9604	3459.9208	0.0000	1.0000
K2PI	13	-1721.6759	3469.3519	0.0000	1.0000
K2PG	13	-1721.8635	3469.7271	0.0000	1.0000
SYMI	17	-1718.4697	3470.9395	0.0000	1.0000
SYMG	17	-1718.7727	3471.5455	0.0000	1.0000
K2PIG	14	-1721.8636	3471.7271	0.0000	1.0000
SYMIG	18	-1718.7728	3473.5455	0.0000	1.0000
F81I	15	-1726.9037	3483.8075	0.0000	1.0000
F81G	15	-1727.0259	3484.0518	0.0000	1.0000
F81IG	16	-1727.0259	3486.0519	0.0000	1.0000
SYM	16	-1734.5525	3501.1049	0.0000	1.0000
K2P	12	-1738.9973	3501.9945	0.0000	1.0000
F81	14	-1742.4483	3512.8967	0.0000	1.0000
JC69I	12	-1747.4267	3518.8535	0.0000	1.0000
JC69G	12	-1747.5672	3519.1344	0.0000	1.0000
JC69IG	13	-1747.5672	3521.1343	0.0000	1.0000
JC69	11	-1763.3012	3548.6023	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1698.5589	3500.5013	0.5301	0.5301
HKYG	16	-1698.7240	3500.8315	0.4495	0.9796
HKYIG	17	-1698.7241	3507.2931	0.0178	0.9974
GTRI	20	-1691.5385	3512.3064	0.0014	0.9988
GTRG	20	-1691.7786	3512.7866	0.0011	1.0000
GTRIG	21	-1691.7786	3519.2480	0.0000	1.0000
HKY	15	-1714.9604	3526.8428	0.0000	1.0000
K2PI	13	-1721.6759	3527.3510	0.0000	1.0000
K2PG	13	-1721.8635	3527.7262	0.0000	1.0000
K2PIG	14	-1721.8636	3534.1877	0.0000	1.0000
GTR	19	-1707.6677	3538.1033	0.0000	1.0000
SYMI	17	-1718.4697	3546.7844	0.0000	1.0000
SYMG	17	-1718.7727	3547.3904	0.0000	1.0000
F81I	15	-1726.9037	3550.7295	0.0000	1.0000
F81G	15	-1727.0259	3550.9739	0.0000	1.0000
SYMIG	18	-1718.7728	3553.8520	0.0000	1.0000
K2P	12	-1738.9973	3555.5321	0.0000	1.0000
F81IG	16	-1727.0259	3557.4354	0.0000	1.0000
JC69I	12	-1747.4267	3572.3911	0.0000	1.0000
SYM	16	-1734.5525	3572.4884	0.0000	1.0000
JC69G	12	-1747.5672	3572.6720	0.0000	1.0000
F81	14	-1742.4483	3575.3572	0.0000	1.0000
JC69IG	13	-1747.5672	3579.1334	0.0000	1.0000
JC69	11	-1763.3012	3597.6784	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_25682.phylip = 1 - 640;
 Partition Dummy = 1:591-tomodel-phylip/chr1_25682.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_25682.phylip = 1 - 640;
 Partition Dummy = 1:591-tomodel-phylip/chr1_25682.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

