
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:04 02/24/2011

Input data from file "/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8051.phylip" (19 taxa, 524 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	44	-1683.3946	3463.0565	0.3070	0.3070
HKYI	40	-1688.2325	3463.2558	0.2779	0.5849
GTRG	44	-1683.9875	3464.2423	0.1697	0.7545
HKYG	40	-1688.8673	3464.5256	0.1473	0.9018
GTRIG	45	-1683.9875	3466.6360	0.0513	0.9531
HKYIG	41	-1688.8673	3466.8798	0.0454	0.9985
K2PI	37	-1697.7777	3475.3414	0.0007	0.9991
SYMI	41	-1693.8054	3476.7561	0.0003	0.9994
K2PG	37	-1698.5584	3476.9027	0.0003	0.9997
SYMG	41	-1694.7389	3478.6230	0.0001	0.9999
K2PIG	38	-1698.5584	3479.2281	0.0001	1.0000
SYMIG	42	-1694.7388	3480.9871	0.0000	1.0000
F81I	39	-1711.3583	3507.1629	0.0000	1.0000
F81G	39	-1712.1714	3508.7891	0.0000	1.0000
F81IG	40	-1712.1713	3511.1335	0.0000	1.0000
HKY	39	-1715.4722	3515.3906	0.0000	1.0000
GTR	43	-1710.8373	3515.5580	0.0000	1.0000
JC69I	36	-1719.8181	3517.1064	0.0000	1.0000
JC69G	36	-1720.7754	3519.0211	0.0000	1.0000
JC69IG	37	-1720.7754	3521.3368	0.0000	1.0000
K2P	36	-1724.9085	3527.2872	0.0000	1.0000
SYM	40	-1720.4708	3527.7325	0.0000	1.0000
F81	38	-1737.9142	3557.9398	0.0000	1.0000
JC69	35	-1746.1489	3567.4617	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	44	-1683.3946	3454.7892	0.3993	0.3993
GTRG	44	-1683.9875	3455.9751	0.2207	0.6200
HKYI	40	-1688.2325	3456.4649	0.1728	0.7928
HKYG	40	-1688.8673	3457.7347	0.0916	0.8843
GTRIG	45	-1683.9875	3457.9749	0.0812	0.9655
HKYIG	41	-1688.8673	3459.7345	0.0337	0.9992
K2PI	37	-1697.7777	3469.5554	0.0002	0.9995
SYMI	41	-1693.8054	3469.6108	0.0002	0.9997
K2PG	37	-1698.5584	3471.1167	0.0001	0.9998
SYMG	41	-1694.7389	3471.4778	0.0001	0.9999
K2PIG	38	-1698.5584	3473.1167	0.0000	1.0000
SYMIG	42	-1694.7388	3473.4777	0.0000	1.0000
F81I	39	-1711.3583	3500.7167	0.0000	1.0000
F81G	39	-1712.1714	3502.3428	0.0000	1.0000
F81IG	40	-1712.1713	3504.3426	0.0000	1.0000
GTR	43	-1710.8373	3507.6747	0.0000	1.0000
HKY	39	-1715.4722	3508.9444	0.0000	1.0000
JC69I	36	-1719.8181	3511.6362	0.0000	1.0000
JC69G	36	-1720.7754	3513.5509	0.0000	1.0000
JC69IG	37	-1720.7754	3515.5508	0.0000	1.0000
SYM	40	-1720.4708	3520.9416	0.0000	1.0000
K2P	36	-1724.9085	3521.8170	0.0000	1.0000
F81	38	-1737.9142	3551.8285	0.0000	1.0000
JC69	35	-1746.1489	3562.2977	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1688.2325	3626.9246	0.3543	0.3543
K2PI	37	-1697.7777	3627.2306	0.3040	0.6583
HKYG	40	-1688.8673	3628.1944	0.1878	0.8461
K2PG	37	-1698.5584	3628.7919	0.1393	0.9854
HKYIG	41	-1688.8673	3634.4557	0.0082	0.9936
K2PIG	38	-1698.5584	3635.0534	0.0061	0.9997
GTRI	44	-1683.3946	3642.2949	0.0002	0.9998
GTRG	44	-1683.9875	3643.4807	0.0001	0.9999
SYMI	41	-1693.8054	3644.3320	0.0001	1.0000
SYMG	41	-1694.7389	3646.1989	0.0000	1.0000
GTRIG	45	-1683.9875	3649.7421	0.0000	1.0000
SYMIG	42	-1694.7388	3652.4604	0.0000	1.0000
JC69I	36	-1719.8181	3665.0499	0.0000	1.0000
F81I	39	-1711.3583	3666.9148	0.0000	1.0000
JC69G	36	-1720.7754	3666.9646	0.0000	1.0000
F81G	39	-1712.1714	3668.5410	0.0000	1.0000
JC69IG	37	-1720.7754	3673.2260	0.0000	1.0000
F81IG	40	-1712.1713	3674.8023	0.0000	1.0000
HKY	39	-1715.4722	3675.1425	0.0000	1.0000
K2P	36	-1724.9085	3675.2307	0.0000	1.0000
GTR	43	-1710.8373	3690.9188	0.0000	1.0000
SYM	40	-1720.4708	3691.4013	0.0000	1.0000
JC69	35	-1746.1489	3711.4499	0.0000	1.0000
F81	38	-1737.9142	3713.7652	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 /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8051.phylip = 1 - 524;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8051.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 /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8051.phylip = 1 - 524;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8051.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 /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8051.phylip = 1 - 524;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_8051.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

