
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:56 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/chr6_5958.phylip" (19 taxa, 763 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	40	-1832.8611	3750.2652	0.3995	0.3995
HKYIG	41	-1832.0489	3750.8745	0.2946	0.6940
HKYI	40	-1833.1934	3750.9296	0.2865	0.9806
GTRG	44	-1832.3599	3758.2351	0.0074	0.9880
GTRIG	45	-1831.2684	3758.3108	0.0072	0.9951
GTRI	44	-1832.7859	3759.0870	0.0049	1.0000
K2PG	37	-1858.3093	3794.4972	0.0000	1.0000
K2PIG	38	-1857.2411	3794.5762	0.0000	1.0000
K2PI	37	-1858.4516	3794.7817	0.0000	1.0000
SYMIG	42	-1855.4284	3799.8734	0.0000	1.0000
SYMG	41	-1856.9523	3800.6813	0.0000	1.0000
SYMI	41	-1857.1216	3801.0199	0.0000	1.0000
HKY	39	-1860.3858	3803.0870	0.0000	1.0000
GTR	43	-1859.0959	3809.4547	0.0000	1.0000
F81G	39	-1867.0126	3816.3405	0.0000	1.0000
F81I	39	-1867.2077	3816.7308	0.0000	1.0000
F81IG	40	-1866.5485	3817.6399	0.0000	1.0000
K2P	36	-1888.2254	3852.1202	0.0000	1.0000
SYM	40	-1886.6863	3857.9155	0.0000	1.0000
JC69G	36	-1891.4751	3858.6196	0.0000	1.0000
JC69I	36	-1891.6287	3858.9267	0.0000	1.0000
JC69IG	37	-1890.8981	3859.6748	0.0000	1.0000
F81	38	-1894.6169	3869.3276	0.0000	1.0000
JC69	35	-1920.0388	3913.5439	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1832.8611	3745.7222	0.3803	0.3803
HKYIG	41	-1832.0489	3746.0978	0.3152	0.6956
HKYI	40	-1833.1934	3746.3867	0.2728	0.9684
GTRIG	45	-1831.2684	3752.5368	0.0126	0.9810
GTRG	44	-1832.3599	3752.7198	0.0115	0.9925
GTRI	44	-1832.7859	3753.5717	0.0075	1.0000
K2PIG	38	-1857.2411	3790.4823	0.0000	1.0000
K2PG	37	-1858.3093	3790.6186	0.0000	1.0000
K2PI	37	-1858.4516	3790.9031	0.0000	1.0000
SYMIG	42	-1855.4284	3794.8567	0.0000	1.0000
SYMG	41	-1856.9523	3795.9046	0.0000	1.0000
SYMI	41	-1857.1216	3796.2432	0.0000	1.0000
HKY	39	-1860.3858	3798.7717	0.0000	1.0000
GTR	43	-1859.0959	3804.1918	0.0000	1.0000
F81G	39	-1867.0126	3812.0251	0.0000	1.0000
F81I	39	-1867.2077	3812.4155	0.0000	1.0000
F81IG	40	-1866.5485	3813.0969	0.0000	1.0000
K2P	36	-1888.2254	3848.4508	0.0000	1.0000
SYM	40	-1886.6863	3853.3725	0.0000	1.0000
JC69G	36	-1891.4751	3854.9502	0.0000	1.0000
JC69I	36	-1891.6287	3855.2573	0.0000	1.0000
JC69IG	37	-1890.8981	3855.7962	0.0000	1.0000
F81	38	-1894.6169	3865.2337	0.0000	1.0000
JC69	35	-1920.0388	3910.0776	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1832.8611	3931.2125	0.5559	0.5559
HKYI	40	-1833.1934	3931.8770	0.3988	0.9547
HKYIG	41	-1832.0489	3936.2253	0.0453	1.0000
GTRG	44	-1832.3599	3956.7592	0.0000	1.0000
GTRI	44	-1832.7859	3957.6111	0.0000	1.0000
GTRIG	45	-1831.2684	3961.2134	0.0000	1.0000
K2PG	37	-1858.3093	3962.1971	0.0000	1.0000
K2PI	37	-1858.4516	3962.4817	0.0000	1.0000
K2PIG	38	-1857.2411	3966.6981	0.0000	1.0000
HKY	39	-1860.3858	3979.6247	0.0000	1.0000
SYMG	41	-1856.9523	3986.0322	0.0000	1.0000
SYMI	41	-1857.1216	3986.3707	0.0000	1.0000
SYMIG	42	-1855.4284	3989.6216	0.0000	1.0000
F81G	39	-1867.0126	3992.8782	0.0000	1.0000
F81I	39	-1867.2077	3993.2685	0.0000	1.0000
F81IG	40	-1866.5485	3998.5872	0.0000	1.0000
GTR	43	-1859.0959	4003.5939	0.0000	1.0000
K2P	36	-1888.2254	4015.3920	0.0000	1.0000
JC69G	36	-1891.4751	4021.8914	0.0000	1.0000
JC69I	36	-1891.6287	4022.1986	0.0000	1.0000
JC69IG	37	-1890.8981	4027.3747	0.0000	1.0000
SYM	40	-1886.6863	4038.8628	0.0000	1.0000
F81	38	-1894.6169	4041.4495	0.0000	1.0000
JC69	35	-1920.0388	4072.3816	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 /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_5958.phylip = 1 - 763;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_5958.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 /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_5958.phylip = 1 - 763;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_5958.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 /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_5958.phylip = 1 - 763;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_5958.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

