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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: K2PG


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	-1862.5848	3811.0059	0.3381	0.3381
GTRG	44	-1858.6130	3812.3228	0.1750	0.5131
HKYIG	41	-1862.3142	3812.7674	0.1401	0.6532
HKYI	40	-1863.5178	3812.8720	0.1330	0.7862
GTRI	44	-1859.5019	3814.1005	0.0720	0.8582
GTRIG	45	-1858.3372	3814.1071	0.0717	0.9299
K2PG	37	-1868.7242	3816.4253	0.0225	0.9524
SYMG	41	-1864.5553	3817.2496	0.0149	0.9673
K2PI	37	-1869.5420	3818.0610	0.0099	0.9772
K2PIG	38	-1868.4530	3818.1612	0.0094	0.9866
SYMI	41	-1865.2891	3818.7173	0.0072	0.9938
SYMIG	42	-1864.2770	3819.0040	0.0062	1.0000
HKY	39	-1882.8427	3849.2271	0.0000	1.0000
GTR	43	-1878.4036	3849.5764	0.0000	1.0000
SYM	40	-1883.3749	3852.5862	0.0000	1.0000
K2P	36	-1888.2596	3853.2259	0.0000	1.0000
F81G	39	-1910.0046	3903.5510	0.0000	1.0000
F81I	39	-1910.6238	3904.7893	0.0000	1.0000
F81IG	40	-1909.7522	3905.3406	0.0000	1.0000
JC69G	36	-1914.9302	3906.5672	0.0000	1.0000
JC69I	36	-1915.5547	3907.8161	0.0000	1.0000
JC69IG	37	-1914.6702	3908.3173	0.0000	1.0000
F81	38	-1929.0114	3939.2781	0.0000	1.0000
JC69	35	-1933.5507	3941.5458	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1862.5848	3805.1696	0.2562	0.2562
GTRG	44	-1858.6130	3805.2260	0.2491	0.5053
HKYIG	41	-1862.3142	3806.6283	0.1236	0.6289
GTRIG	45	-1858.3372	3806.6744	0.1207	0.7496
GTRI	44	-1859.5019	3807.0037	0.1024	0.8520
HKYI	40	-1863.5178	3807.0357	0.1008	0.9528
SYMG	41	-1864.5553	3811.1106	0.0131	0.9660
K2PG	37	-1868.7242	3811.4483	0.0111	0.9771
SYMIG	42	-1864.2770	3812.5540	0.0064	0.9834
SYMI	41	-1865.2891	3812.5782	0.0063	0.9897
K2PIG	38	-1868.4530	3812.9059	0.0054	0.9951
K2PI	37	-1869.5420	3813.0840	0.0049	1.0000
GTR	43	-1878.4036	3842.8071	0.0000	1.0000
HKY	39	-1882.8427	3843.6854	0.0000	1.0000
SYM	40	-1883.3749	3846.7499	0.0000	1.0000
K2P	36	-1888.2596	3848.5192	0.0000	1.0000
F81G	39	-1910.0046	3898.0093	0.0000	1.0000
F81I	39	-1910.6238	3899.2476	0.0000	1.0000
F81IG	40	-1909.7522	3899.5043	0.0000	1.0000
JC69G	36	-1914.9302	3901.8605	0.0000	1.0000
JC69I	36	-1915.5547	3903.1094	0.0000	1.0000
JC69IG	37	-1914.6702	3903.3403	0.0000	1.0000
F81	38	-1929.0114	3934.0228	0.0000	1.0000
JC69	35	-1933.5507	3937.1013	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1868.7242	3974.3193	0.6492	0.6492
K2PI	37	-1869.5420	3975.9550	0.2865	0.9357
K2PIG	38	-1868.4530	3980.1788	0.0347	0.9704
HKYG	40	-1862.5848	3981.2463	0.0203	0.9907
HKYI	40	-1863.5178	3983.1124	0.0080	0.9987
HKYIG	41	-1862.3142	3987.1069	0.0011	0.9998
SYMG	41	-1864.5553	3991.5892	0.0001	0.9999
SYMI	41	-1865.2891	3993.0568	0.0001	1.0000
SYMIG	42	-1864.2770	3997.4345	0.0000	1.0000
GTRG	44	-1858.6130	3998.9103	0.0000	1.0000
GTRI	44	-1859.5019	4000.6881	0.0000	1.0000
GTRIG	45	-1858.3372	4004.7607	0.0000	1.0000
K2P	36	-1888.2596	4006.9882	0.0000	1.0000
HKY	39	-1882.8427	4015.3601	0.0000	1.0000
SYM	40	-1883.3749	4022.8266	0.0000	1.0000
GTR	43	-1878.4036	4032.0896	0.0000	1.0000
JC69G	36	-1914.9302	4060.3295	0.0000	1.0000
JC69I	36	-1915.5547	4061.5784	0.0000	1.0000
JC69IG	37	-1914.6702	4066.2113	0.0000	1.0000
F81G	39	-1910.0046	4069.6840	0.0000	1.0000
F81I	39	-1910.6238	4070.9223	0.0000	1.0000
F81IG	40	-1909.7522	4075.5810	0.0000	1.0000
JC69	35	-1933.5507	4091.1684	0.0000	1.0000
F81	38	-1929.0114	4101.2956	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/chr2_23635.phylip = 1 - 603;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_23635.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/chr2_23635.phylip = 1 - 603;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_23635.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 (K2PG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_23635.phylip = 1 - 603;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_23635.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

