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

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
Minimum BIC  model: K2PIG


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
HKYIG	41	-1900.3888	3887.9721	0.8459	0.8459
GTRIG	45	-1898.2110	3892.7042	0.0794	0.9253
HKYI	40	-1904.6469	3894.2336	0.0370	0.9622
HKYG	40	-1904.9024	3894.7445	0.0286	0.9908
K2PIG	38	-1909.1927	3898.8359	0.0037	0.9945
GTRI	44	-1902.7711	3899.5422	0.0026	0.9971
GTRG	44	-1902.8549	3899.7097	0.0024	0.9995
K2PI	37	-1913.1327	3904.4814	0.0002	0.9997
K2PG	37	-1913.3934	3905.0027	0.0002	0.9999
SYMIG	42	-1908.4776	3906.4114	0.0001	1.0000
SYMI	41	-1912.5165	3912.2275	0.0000	1.0000
SYMG	41	-1912.6846	3912.5638	0.0000	1.0000
F81IG	40	-1923.8060	3932.5517	0.0000	1.0000
F81I	39	-1927.6088	3937.9094	0.0000	1.0000
F81G	39	-1927.7366	3938.1650	0.0000	1.0000
JC69IG	37	-1931.1368	3940.4896	0.0000	1.0000
JC69I	36	-1934.8212	3945.6303	0.0000	1.0000
JC69G	36	-1934.9510	3945.8900	0.0000	1.0000
HKY	39	-1957.8507	3998.3931	0.0000	1.0000
GTR	43	-1955.4931	4002.7108	0.0000	1.0000
K2P	36	-1966.0170	4008.0220	0.0000	1.0000
SYM	40	-1964.8010	4014.5418	0.0000	1.0000
F81	38	-1979.2466	4038.9436	0.0000	1.0000
JC69	35	-1986.3656	4046.4980	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	41	-1900.3888	3882.7775	0.8051	0.8051
GTRIG	45	-1898.2110	3886.4220	0.1301	0.9352
HKYI	40	-1904.6469	3889.2938	0.0310	0.9662
HKYG	40	-1904.9024	3889.8047	0.0240	0.9901
GTRI	44	-1902.7711	3893.5422	0.0037	0.9938
GTRG	44	-1902.8549	3893.7097	0.0034	0.9972
K2PIG	38	-1909.1927	3894.3855	0.0024	0.9997
K2PI	37	-1913.1327	3900.2655	0.0001	0.9998
K2PG	37	-1913.3934	3900.7868	0.0001	0.9999
SYMIG	42	-1908.4776	3900.9552	0.0001	1.0000
SYMI	41	-1912.5165	3907.0329	0.0000	1.0000
SYMG	41	-1912.6846	3907.3692	0.0000	1.0000
F81IG	40	-1923.8060	3927.6119	0.0000	1.0000
F81I	39	-1927.6088	3933.2177	0.0000	1.0000
F81G	39	-1927.7366	3933.4733	0.0000	1.0000
JC69IG	37	-1931.1368	3936.2737	0.0000	1.0000
JC69I	36	-1934.8212	3941.6423	0.0000	1.0000
JC69G	36	-1934.9510	3941.9020	0.0000	1.0000
HKY	39	-1957.8507	3993.7014	0.0000	1.0000
GTR	43	-1955.4931	3996.9861	0.0000	1.0000
K2P	36	-1966.0170	4004.0340	0.0000	1.0000
SYM	40	-1964.8010	4009.6020	0.0000	1.0000
F81	38	-1979.2466	4034.4931	0.0000	1.0000
JC69	35	-1986.3656	4042.7312	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PIG	38	-1909.1927	4067.5970	0.3988	0.3988
K2PI	37	-1913.1327	4068.9188	0.2060	0.6048
K2PG	37	-1913.3934	4069.4401	0.1587	0.7635
HKYIG	41	-1900.3888	4069.6636	0.1419	0.9054
HKYI	40	-1904.6469	4071.6218	0.0533	0.9587
HKYG	40	-1904.9024	4072.1327	0.0413	1.0000
GTRIG	45	-1898.2110	4091.5409	0.0000	1.0000
SYMIG	42	-1908.4776	4092.3995	0.0000	1.0000
SYMI	41	-1912.5165	4093.9190	0.0000	1.0000
GTRI	44	-1902.7711	4094.1029	0.0000	1.0000
SYMG	41	-1912.6846	4094.2553	0.0000	1.0000
GTRG	44	-1902.8549	4094.2704	0.0000	1.0000
JC69IG	37	-1931.1368	4104.9270	0.0000	1.0000
JC69I	36	-1934.8212	4105.7374	0.0000	1.0000
JC69G	36	-1934.9510	4105.9971	0.0000	1.0000
F81IG	40	-1923.8060	4109.9398	0.0000	1.0000
F81I	39	-1927.6088	4110.9874	0.0000	1.0000
F81G	39	-1927.7366	4111.2430	0.0000	1.0000
K2P	36	-1966.0170	4168.1291	0.0000	1.0000
HKY	39	-1957.8507	4171.4711	0.0000	1.0000
SYM	40	-1964.8010	4191.9300	0.0000	1.0000
GTR	43	-1955.4931	4192.9886	0.0000	1.0000
JC69	35	-1986.3656	4202.2681	0.0000	1.0000
F81	38	-1979.2466	4207.7046	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21324.phylip = 1 - 705;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21324.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21324.phylip = 1 - 705;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21324.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (K2PIG)]

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


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

