
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:53 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/chr7_10409.phylip" (19 taxa, 273 characters)

Minimum AIC  model: K2PG
Minimum AICc model: K2PG
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
K2PG	37	-845.4059	1776.7777	0.4849	0.4849
K2PI	37	-846.3414	1778.6487	0.1903	0.6752
K2PIG	38	-845.3090	1779.2846	0.1384	0.8136
HKYG	40	-843.4513	1781.0405	0.0575	0.8712
HKYI	40	-844.1295	1782.3970	0.0292	0.9004
SYMG	41	-842.7880	1782.4851	0.0279	0.9283
HKYIG	41	-843.3509	1783.6108	0.0159	0.9442
JC69G	36	-850.3871	1784.0624	0.0127	0.9569
SYMI	41	-843.7913	1784.4917	0.0102	0.9672
SYMIG	42	-842.6935	1785.0913	0.0076	0.9748
JC69I	36	-850.9394	1785.1669	0.0073	0.9821
GTRG	44	-840.0188	1785.4060	0.0065	0.9886
JC69IG	37	-850.3504	1786.6668	0.0035	0.9920
GTRI	44	-840.8585	1787.0854	0.0028	0.9948
GTRIG	45	-839.9125	1788.0629	0.0017	0.9965
F81G	39	-848.5105	1788.4115	0.0014	0.9980
F81I	39	-849.0304	1789.4514	0.0009	0.9988
K2P	36	-853.3895	1790.0672	0.0006	0.9995
F81IG	40	-848.4761	1791.0902	0.0004	0.9998
HKY	39	-851.3985	1794.1875	0.0001	0.9999
JC69	35	-857.5993	1795.8315	0.0000	1.0000
SYM	40	-850.8701	1795.8782	0.0000	1.0000
GTR	43	-848.4769	1799.4777	0.0000	1.0000
F81	38	-855.8176	1800.3018	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	37	-845.4059	1764.8118	0.3001	0.3001
K2PIG	38	-845.3090	1766.6180	0.1216	0.4217
K2PI	37	-846.3414	1766.6828	0.1178	0.5395
HKYG	40	-843.4513	1766.9026	0.1055	0.6450
SYMG	41	-842.7880	1767.5760	0.0753	0.7203
GTRG	44	-840.0188	1768.0376	0.0598	0.7801
HKYI	40	-844.1295	1768.2590	0.0535	0.8337
HKYIG	41	-843.3509	1768.7017	0.0429	0.8766
SYMIG	42	-842.6935	1769.3870	0.0305	0.9070
SYMI	41	-843.7913	1769.5826	0.0276	0.9347
GTRI	44	-840.8585	1769.7170	0.0258	0.9605
GTRIG	45	-839.9125	1769.8250	0.0245	0.9849
JC69G	36	-850.3871	1772.7742	0.0056	0.9905
JC69I	36	-850.9394	1773.8788	0.0032	0.9938
JC69IG	37	-850.3504	1774.7008	0.0021	0.9959
F81G	39	-848.5105	1775.0209	0.0018	0.9977
F81I	39	-849.0304	1776.0608	0.0011	0.9988
F81IG	40	-848.4761	1776.9522	0.0007	0.9995
K2P	36	-853.3895	1778.7790	0.0003	0.9998
HKY	39	-851.3985	1780.7970	0.0001	0.9999
SYM	40	-850.8701	1781.7403	0.0001	1.0000
GTR	43	-848.4769	1782.9537	0.0000	1.0000
JC69	35	-857.5993	1785.1986	0.0000	1.0000
F81	38	-855.8176	1787.6351	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-845.4059	1898.3622	0.6042	0.6042
K2PI	37	-846.3414	1900.2332	0.2371	0.8413
JC69G	36	-850.3871	1902.7152	0.0685	0.9099
K2PIG	38	-845.3090	1903.7779	0.0403	0.9502
JC69I	36	-850.9394	1903.8198	0.0395	0.9896
JC69IG	37	-850.3504	1908.2513	0.0043	0.9939
K2P	36	-853.3895	1908.7200	0.0034	0.9973
HKYG	40	-843.4513	1911.2814	0.0009	0.9983
JC69	35	-857.5993	1911.5301	0.0008	0.9991
HKYI	40	-844.1295	1912.6379	0.0005	0.9996
SYMG	41	-842.7880	1915.5644	0.0001	0.9997
F81G	39	-848.5105	1915.7903	0.0001	0.9998
HKYIG	41	-843.3509	1916.6900	0.0001	0.9999
F81I	39	-849.0304	1916.8302	0.0001	0.9999
SYMI	41	-843.7913	1917.5710	0.0000	1.0000
SYMIG	42	-842.6935	1920.9848	0.0000	1.0000
F81IG	40	-848.4761	1921.3311	0.0000	1.0000
HKY	39	-851.3985	1921.5664	0.0000	1.0000
F81	38	-855.8176	1924.7951	0.0000	1.0000
SYM	40	-850.8701	1926.1192	0.0000	1.0000
GTRG	44	-840.0188	1926.8543	0.0000	1.0000
GTRI	44	-840.8585	1928.5337	0.0000	1.0000
GTRIG	45	-839.9125	1932.2513	0.0000	1.0000
GTR	43	-848.4769	1938.1610	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (K2PG)]

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


[Mrbayes block for the best AICc model (K2PG)]

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


[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/chr7_10409.phylip = 1 - 273;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10409.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

