
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/chr3_3073.phylip" (19 taxa, 566 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	-1941.3377	3962.0011	0.3914	0.3914
K2PIG	38	-1940.8908	3963.4058	0.1939	0.5853
K2PI	37	-1942.3322	3963.9901	0.1448	0.7301
HKYG	40	-1939.1306	3964.5087	0.1117	0.8418
HKYIG	41	-1938.6308	3965.8341	0.0576	0.8994
HKYI	40	-1940.1381	3966.5238	0.0408	0.9401
GTRG	44	-1936.3269	3968.2545	0.0172	0.9573
SYMG	41	-1940.0600	3968.6926	0.0138	0.9711
GTRIG	45	-1935.6389	3969.2393	0.0105	0.9816
SYMIG	42	-1939.5042	3969.9148	0.0075	0.9891
GTRI	44	-1937.3744	3970.3495	0.0060	0.9951
SYMI	41	-1941.0959	3970.7644	0.0049	1.0000
JC69G	36	-1974.0711	4025.1780	0.0000	1.0000
JC69IG	37	-1973.6992	4026.7242	0.0000	1.0000
JC69I	36	-1975.0905	4027.2169	0.0000	1.0000
F81G	39	-1972.6393	4029.2101	0.0000	1.0000
F81IG	40	-1972.2543	4030.7562	0.0000	1.0000
F81I	39	-1973.6537	4031.2389	0.0000	1.0000
K2P	36	-2017.3132	4111.6623	0.0000	1.0000
HKY	39	-2015.6840	4115.2996	0.0000	1.0000
SYM	40	-2016.4346	4119.1169	0.0000	1.0000
GTR	43	-2014.2086	4121.6662	0.0000	1.0000
JC69	35	-2048.4841	4171.7228	0.0000	1.0000
F81	38	-2047.2645	4176.1533	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	37	-1941.3377	3956.6753	0.3074	0.3074
K2PIG	38	-1940.8908	3957.7815	0.1768	0.4842
HKYG	40	-1939.1306	3958.2611	0.1391	0.6233
K2PI	37	-1942.3322	3958.6643	0.1137	0.7370
HKYIG	41	-1938.6308	3959.2616	0.0843	0.8213
HKYI	40	-1940.1381	3960.2762	0.0508	0.8721
GTRG	44	-1936.3269	3960.6537	0.0420	0.9141
GTRIG	45	-1935.6389	3961.2778	0.0308	0.9449
SYMG	41	-1940.0600	3962.1200	0.0202	0.9651
GTRI	44	-1937.3744	3962.7487	0.0148	0.9799
SYMIG	42	-1939.5042	3963.0085	0.0130	0.9928
SYMI	41	-1941.0959	3964.1919	0.0072	1.0000
JC69G	36	-1974.0711	4020.1421	0.0000	1.0000
JC69IG	37	-1973.6992	4021.3984	0.0000	1.0000
JC69I	36	-1975.0905	4022.1809	0.0000	1.0000
F81G	39	-1972.6393	4023.2785	0.0000	1.0000
F81IG	40	-1972.2543	4024.5086	0.0000	1.0000
F81I	39	-1973.6537	4025.3073	0.0000	1.0000
K2P	36	-2017.3132	4106.6264	0.0000	1.0000
HKY	39	-2015.6840	4109.3681	0.0000	1.0000
SYM	40	-2016.4346	4112.8693	0.0000	1.0000
GTR	43	-2014.2086	4114.4172	0.0000	1.0000
JC69	35	-2048.4841	4166.9681	0.0000	1.0000
F81	38	-2047.2645	4170.5290	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1941.3377	4117.2033	0.6961	0.6961
K2PI	37	-1942.3322	4119.1923	0.2575	0.9536
K2PIG	38	-1940.8908	4122.6481	0.0457	0.9993
HKYG	40	-1939.1306	4131.8049	0.0005	0.9998
HKYI	40	-1940.1381	4133.8199	0.0002	1.0000
HKYIG	41	-1938.6308	4137.1440	0.0000	1.0000
SYMG	41	-1940.0600	4140.0024	0.0000	1.0000
SYMI	41	-1941.0959	4142.0743	0.0000	1.0000
SYMIG	42	-1939.5042	4145.2294	0.0000	1.0000
GTRG	44	-1936.3269	4151.5519	0.0000	1.0000
GTRI	44	-1937.3744	4153.6468	0.0000	1.0000
GTRIG	45	-1935.6389	4156.5145	0.0000	1.0000
JC69G	36	-1974.0711	4176.3315	0.0000	1.0000
JC69I	36	-1975.0905	4178.3703	0.0000	1.0000
JC69IG	37	-1973.6992	4181.9264	0.0000	1.0000
F81G	39	-1972.6393	4192.4837	0.0000	1.0000
F81I	39	-1973.6537	4194.5125	0.0000	1.0000
F81IG	40	-1972.2543	4198.0524	0.0000	1.0000
K2P	36	-2017.3132	4262.8158	0.0000	1.0000
HKY	39	-2015.6840	4278.5732	0.0000	1.0000
SYM	40	-2016.4346	4286.4131	0.0000	1.0000
GTR	43	-2014.2086	4300.9767	0.0000	1.0000
JC69	35	-2048.4841	4318.8189	0.0000	1.0000
F81	38	-2047.2645	4335.3956	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/chr3_3073.phylip = 1 - 566;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_3073.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/chr3_3073.phylip = 1 - 566;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_3073.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/chr3_3073.phylip = 1 - 566;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_3073.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

