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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	44	-1822.4132	3744.1407	0.5773	0.5773
GTRIG	45	-1821.9274	3745.7172	0.2625	0.8398
HKYG	40	-1829.3748	3748.0152	0.0832	0.9230
HKYIG	41	-1829.1336	3750.0235	0.0305	0.9534
GTRI	44	-1825.4447	3750.2037	0.0279	0.9813
HKYI	40	-1831.9082	3753.0819	0.0066	0.9879
K2PG	37	-1835.8055	3753.4877	0.0054	0.9933
SYMG	41	-1831.6875	3755.1313	0.0024	0.9957
K2PIG	38	-1835.6050	3755.5359	0.0019	0.9976
SYMIG	42	-1831.0745	3756.4103	0.0013	0.9988
SYMI	41	-1832.9677	3757.6917	0.0007	0.9995
K2PI	37	-1838.2011	3758.2790	0.0005	1.0000
GTR	43	-1875.2055	3847.1916	0.0000	1.0000
F81G	39	-1880.3700	3847.5288	0.0000	1.0000
HKY	39	-1880.7881	3848.3649	0.0000	1.0000
JC69G	36	-1885.0883	3849.6180	0.0000	1.0000
F81IG	40	-1880.3020	3849.8694	0.0000	1.0000
JC69IG	37	-1884.9256	3851.7279	0.0000	1.0000
F81I	39	-1882.6998	3852.1884	0.0000	1.0000
JC69I	36	-1886.9388	3853.3189	0.0000	1.0000
K2P	36	-1887.5396	3854.5206	0.0000	1.0000
SYM	40	-1883.9717	3857.2089	0.0000	1.0000
F81	38	-1928.6296	3941.5851	0.0000	1.0000
JC69	35	-1933.8255	3944.6705	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1822.4132	3732.8265	0.5772	0.5772
GTRIG	45	-1821.9274	3733.8548	0.3452	0.9223
HKYG	40	-1829.3748	3738.7496	0.0299	0.9522
GTRI	44	-1825.4447	3738.8895	0.0278	0.9800
HKYIG	41	-1829.1336	3740.2672	0.0140	0.9940
HKYI	40	-1831.9082	3743.8164	0.0024	0.9964
SYMG	41	-1831.6875	3745.3750	0.0011	0.9975
K2PG	37	-1835.8055	3745.6109	0.0010	0.9984
SYMIG	42	-1831.0745	3746.1489	0.0007	0.9992
K2PIG	38	-1835.6050	3747.2101	0.0004	0.9996
SYMI	41	-1832.9677	3747.9353	0.0003	0.9999
K2PI	37	-1838.2011	3750.4023	0.0001	1.0000
GTR	43	-1875.2055	3836.4110	0.0000	1.0000
F81G	39	-1880.3700	3838.7400	0.0000	1.0000
HKY	39	-1880.7881	3839.5761	0.0000	1.0000
F81IG	40	-1880.3020	3840.6039	0.0000	1.0000
JC69G	36	-1885.0883	3842.1767	0.0000	1.0000
F81I	39	-1882.6998	3843.3997	0.0000	1.0000
JC69IG	37	-1884.9256	3843.8512	0.0000	1.0000
JC69I	36	-1886.9388	3845.8776	0.0000	1.0000
K2P	36	-1887.5396	3847.0792	0.0000	1.0000
SYM	40	-1883.9717	3847.9434	0.0000	1.0000
F81	38	-1928.6296	3933.2592	0.0000	1.0000
JC69	35	-1933.8255	3937.6510	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1835.8055	3892.8297	0.8038	0.8038
K2PI	37	-1838.2011	3897.6210	0.0732	0.8770
HKYG	40	-1829.3748	3897.9051	0.0635	0.9405
K2PIG	38	-1835.6050	3898.4077	0.0494	0.9900
HKYI	40	-1831.9082	3902.9718	0.0050	0.9950
HKYIG	41	-1829.1336	3903.4015	0.0041	0.9991
GTRG	44	-1822.4132	3907.8974	0.0004	0.9995
SYMG	41	-1831.6875	3908.5093	0.0003	0.9998
SYMI	41	-1832.9677	3911.0696	0.0001	0.9999
GTRIG	45	-1821.9274	3912.9046	0.0000	0.9999
SYMIG	42	-1831.0745	3913.2621	0.0000	1.0000
GTRI	44	-1825.4447	3913.9604	0.0000	1.0000
JC69G	36	-1885.0883	3985.4166	0.0000	1.0000
JC69I	36	-1886.9388	3989.1174	0.0000	1.0000
K2P	36	-1887.5396	3990.3191	0.0000	1.0000
JC69IG	37	-1884.9256	3991.0700	0.0000	1.0000
F81G	39	-1880.3700	3993.9166	0.0000	1.0000
HKY	39	-1880.7881	3994.7527	0.0000	1.0000
F81I	39	-1882.6998	3998.5762	0.0000	1.0000
F81IG	40	-1880.3020	3999.7593	0.0000	1.0000
SYM	40	-1883.9717	4007.0988	0.0000	1.0000
GTR	43	-1875.2055	4007.5031	0.0000	1.0000
JC69	35	-1933.8255	4076.9120	0.0000	1.0000
F81	38	-1928.6296	4084.4569	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_5924.phylip = 1 - 395;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_5924.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_5924.phylip = 1 - 395;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_5924.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) 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/chr3_5924.phylip = 1 - 395;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_5924.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

