
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:50 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/chr3_2733.phylip" (19 taxa, 758 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	40	-1787.9535	3660.4815	0.4231	0.4231
HKYG	40	-1788.1047	3660.7840	0.3637	0.7867
HKYIG	41	-1788.1044	3663.0189	0.1190	0.9057
GTRI	44	-1785.7182	3664.9905	0.0444	0.9501
GTRG	44	-1785.8820	3665.3180	0.0377	0.9878
GTRIG	45	-1785.8818	3667.5782	0.0122	1.0000
HKY	39	-1798.5526	3679.4506	0.0000	1.0000
GTR	43	-1796.9007	3685.1012	0.0000	1.0000
SYMI	41	-1800.8210	3688.4520	0.0000	1.0000
SYMG	41	-1800.9385	3688.6870	0.0000	1.0000
SYMIG	42	-1800.9382	3690.9281	0.0000	1.0000
K2PI	37	-1806.5831	3691.0718	0.0000	1.0000
K2PG	37	-1806.7112	3691.3280	0.0000	1.0000
K2PIG	38	-1806.7111	3693.5445	0.0000	1.0000
SYM	40	-1811.3861	3707.3468	0.0000	1.0000
K2P	36	-1816.6258	3708.9464	0.0000	1.0000
F81I	39	-1840.4189	3763.1831	0.0000	1.0000
F81G	39	-1840.5095	3763.3645	0.0000	1.0000
F81IG	40	-1840.5139	3765.6024	0.0000	1.0000
F81	38	-1850.3892	3780.9008	0.0000	1.0000
JC69I	36	-1858.6244	3792.9436	0.0000	1.0000
JC69G	36	-1858.7387	3793.1722	0.0000	1.0000
JC69IG	37	-1858.7538	3795.4131	0.0000	1.0000
JC69	35	-1868.3322	3810.1546	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1787.9535	3655.9069	0.3928	0.3928
HKYG	40	-1788.1047	3656.2094	0.3376	0.7304
HKYIG	41	-1788.1044	3658.2089	0.1242	0.8546
GTRI	44	-1785.7182	3659.4365	0.0672	0.9219
GTRG	44	-1785.8820	3659.7640	0.0571	0.9790
GTRIG	45	-1785.8818	3661.7636	0.0210	1.0000
HKY	39	-1798.5526	3675.1052	0.0000	1.0000
GTR	43	-1796.9007	3679.8015	0.0000	1.0000
SYMI	41	-1800.8210	3683.6419	0.0000	1.0000
SYMG	41	-1800.9385	3683.8769	0.0000	1.0000
SYMIG	42	-1800.9382	3685.8764	0.0000	1.0000
K2PI	37	-1806.5831	3687.1662	0.0000	1.0000
K2PG	37	-1806.7112	3687.4225	0.0000	1.0000
K2PIG	38	-1806.7111	3689.4221	0.0000	1.0000
SYM	40	-1811.3861	3702.7722	0.0000	1.0000
K2P	36	-1816.6258	3705.2516	0.0000	1.0000
F81I	39	-1840.4189	3758.8377	0.0000	1.0000
F81G	39	-1840.5095	3759.0191	0.0000	1.0000
F81IG	40	-1840.5139	3761.0278	0.0000	1.0000
F81	38	-1850.3892	3776.7784	0.0000	1.0000
JC69I	36	-1858.6244	3789.2487	0.0000	1.0000
JC69G	36	-1858.7387	3789.4773	0.0000	1.0000
JC69IG	37	-1858.7538	3791.5075	0.0000	1.0000
JC69	35	-1868.3322	3806.6643	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1787.9535	3841.1343	0.5286	0.5286
HKYG	40	-1788.1047	3841.4367	0.4544	0.9829
HKYIG	41	-1788.1044	3848.0669	0.0165	0.9994
HKY	39	-1798.5526	3855.7019	0.0004	0.9998
K2PI	37	-1806.5831	3858.5015	0.0001	0.9999
K2PG	37	-1806.7112	3858.7577	0.0001	1.0000
GTRI	44	-1785.7182	3863.1866	0.0000	1.0000
GTRG	44	-1785.8820	3863.5141	0.0000	1.0000
K2PIG	38	-1806.7111	3865.3881	0.0000	1.0000
GTRIG	45	-1785.8818	3870.1443	0.0000	1.0000
K2P	36	-1816.6258	3871.9562	0.0000	1.0000
SYMI	41	-1800.8210	3873.4999	0.0000	1.0000
SYMG	41	-1800.9385	3873.7349	0.0000	1.0000
GTR	43	-1796.9007	3878.9209	0.0000	1.0000
SYMIG	42	-1800.9382	3880.3651	0.0000	1.0000
SYM	40	-1811.3861	3887.9995	0.0000	1.0000
F81I	39	-1840.4189	3939.4344	0.0000	1.0000
F81G	39	-1840.5095	3939.6157	0.0000	1.0000
F81IG	40	-1840.5139	3946.2551	0.0000	1.0000
F81	38	-1850.3892	3952.7444	0.0000	1.0000
JC69I	36	-1858.6244	3955.9533	0.0000	1.0000
JC69G	36	-1858.7387	3956.1819	0.0000	1.0000
JC69IG	37	-1858.7538	3962.8428	0.0000	1.0000
JC69	35	-1868.3322	3968.7383	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 /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_2733.phylip = 1 - 758;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_2733.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 /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_2733.phylip = 1 - 758;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_2733.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 /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_2733.phylip = 1 - 758;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_2733.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

