
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:46 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_16833.phylip" (19 taxa, 529 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: GTRI


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	-1845.8559	3787.8936	0.5121	0.5121
GTRG	44	-1846.3261	3788.8340	0.3200	0.8321
GTRIG	45	-1845.7767	3790.1247	0.1678	0.9999
HKYI	40	-1859.6267	3805.9748	0.0001	1.0000
HKYG	40	-1860.4548	3807.6309	0.0000	1.0000
HKYIG	41	-1859.6342	3808.3403	0.0000	1.0000
GTR	43	-1871.5987	3836.9994	0.0000	1.0000
F81I	39	-1877.8254	3840.0311	0.0000	1.0000
F81G	39	-1878.4938	3841.3680	0.0000	1.0000
F81IG	40	-1877.8277	3842.3767	0.0000	1.0000
HKY	39	-1882.0402	3848.4608	0.0000	1.0000
K2PI	37	-1888.6235	3856.9741	0.0000	1.0000
SYMI	41	-1884.3098	3857.6914	0.0000	1.0000
K2PG	37	-1889.3561	3858.4392	0.0000	1.0000
SYMG	41	-1884.9383	3858.9485	0.0000	1.0000
K2PIG	38	-1888.6303	3859.3095	0.0000	1.0000
SYMIG	42	-1884.3109	3860.0539	0.0000	1.0000
F81	38	-1898.6219	3879.2927	0.0000	1.0000
JC69I	36	-1901.8814	3881.1774	0.0000	1.0000
JC69G	36	-1902.4823	3882.3791	0.0000	1.0000
JC69IG	37	-1901.8825	3883.4922	0.0000	1.0000
SYM	40	-1903.4524	3893.6261	0.0000	1.0000
K2P	36	-1909.5922	3896.5991	0.0000	1.0000
JC69	35	-1922.3016	3919.7148	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	44	-1845.8559	3779.7118	0.4943	0.4943
GTRG	44	-1846.3261	3780.6522	0.3089	0.8031
GTRIG	45	-1845.7767	3781.5533	0.1968	0.9999
HKYI	40	-1859.6267	3799.2535	0.0000	1.0000
HKYG	40	-1860.4548	3800.9096	0.0000	1.0000
HKYIG	41	-1859.6342	3801.2685	0.0000	1.0000
GTR	43	-1871.5987	3829.1973	0.0000	1.0000
F81I	39	-1877.8254	3833.6507	0.0000	1.0000
F81G	39	-1878.4938	3834.9876	0.0000	1.0000
F81IG	40	-1877.8277	3835.6554	0.0000	1.0000
HKY	39	-1882.0402	3842.0805	0.0000	1.0000
SYMI	41	-1884.3098	3850.6195	0.0000	1.0000
K2PI	37	-1888.6235	3851.2470	0.0000	1.0000
SYMG	41	-1884.9383	3851.8767	0.0000	1.0000
SYMIG	42	-1884.3109	3852.6218	0.0000	1.0000
K2PG	37	-1889.3561	3852.7121	0.0000	1.0000
K2PIG	38	-1888.6303	3853.2606	0.0000	1.0000
F81	38	-1898.6219	3873.2437	0.0000	1.0000
JC69I	36	-1901.8814	3875.7628	0.0000	1.0000
JC69G	36	-1902.4823	3876.9645	0.0000	1.0000
JC69IG	37	-1901.8825	3877.7651	0.0000	1.0000
SYM	40	-1903.4524	3886.9048	0.0000	1.0000
K2P	36	-1909.5922	3891.1845	0.0000	1.0000
JC69	35	-1922.3016	3914.6033	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	44	-1845.8559	3967.6353	0.4751	0.4751
GTRG	44	-1846.3261	3968.5757	0.2968	0.7719
HKYI	40	-1859.6267	3970.0930	0.1390	0.9109
HKYG	40	-1860.4548	3971.7491	0.0607	0.9716
GTRIG	45	-1845.7767	3973.7478	0.0224	0.9940
HKYIG	41	-1859.6342	3976.3790	0.0060	1.0000
F81I	39	-1877.8254	4000.2192	0.0000	1.0000
F81G	39	-1878.4938	4001.5561	0.0000	1.0000
F81IG	40	-1877.8277	4006.4950	0.0000	1.0000
HKY	39	-1882.0402	4008.6490	0.0000	1.0000
K2PI	37	-1888.6235	4009.2736	0.0000	1.0000
K2PG	37	-1889.3561	4010.7387	0.0000	1.0000
GTR	43	-1871.5987	4012.8498	0.0000	1.0000
K2PIG	38	-1888.6303	4015.5581	0.0000	1.0000
SYMI	41	-1884.3098	4025.7301	0.0000	1.0000
SYMG	41	-1884.9383	4026.9872	0.0000	1.0000
JC69I	36	-1901.8814	4029.5184	0.0000	1.0000
JC69G	36	-1902.4823	4030.7201	0.0000	1.0000
SYMIG	42	-1884.3109	4032.0033	0.0000	1.0000
F81	38	-1898.6219	4035.5413	0.0000	1.0000
JC69IG	37	-1901.8825	4035.7917	0.0000	1.0000
K2P	36	-1909.5922	4044.9401	0.0000	1.0000
SYM	40	-1903.4524	4057.7444	0.0000	1.0000
JC69	35	-1922.3016	4064.0879	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/chr3_16833.phylip = 1 - 529;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_16833.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/chr3_16833.phylip = 1 - 529;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_16833.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 (GTRI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_16833.phylip = 1 - 529;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_16833.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;


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

