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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	-2321.3596	4737.4083	0.6308	0.6308
GTRIG	45	-2320.8423	4738.6898	0.3324	0.9632
GTRI	44	-2324.2021	4743.0933	0.0368	1.0000
HKYG	40	-2336.8028	4759.1090	0.0000	1.0000
HKYIG	41	-2336.5538	4760.8958	0.0000	1.0000
HKYI	40	-2338.4882	4762.4797	0.0000	1.0000
SYMG	41	-2370.8972	4829.5826	0.0000	1.0000
SYMIG	42	-2370.7250	4831.5308	0.0000	1.0000
SYMI	41	-2372.5161	4832.8205	0.0000	1.0000
K2PG	37	-2377.3244	4833.3432	0.0000	1.0000
K2PIG	38	-2377.2194	4835.3953	0.0000	1.0000
K2PI	37	-2378.7503	4836.1950	0.0000	1.0000
GTR	43	-2374.9160	4842.2130	0.0000	1.0000
HKY	39	-2390.7472	4864.7206	0.0000	1.0000
F81G	39	-2411.8498	4906.9257	0.0000	1.0000
F81IG	40	-2411.7875	4909.0784	0.0000	1.0000
F81I	39	-2413.2315	4909.6891	0.0000	1.0000
SYM	40	-2427.6489	4940.8012	0.0000	1.0000
K2P	36	-2435.2893	4947.0186	0.0000	1.0000
JC69G	36	-2443.3380	4963.1160	0.0000	1.0000
JC69I	36	-2443.3930	4963.2260	0.0000	1.0000
JC69IG	37	-2442.3172	4963.3289	0.0000	1.0000
F81	38	-2466.5770	5014.1106	0.0000	1.0000
JC69	35	-2498.0165	5070.2260	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2321.3596	4730.7191	0.5969	0.5969
GTRIG	45	-2320.8423	4731.6847	0.3683	0.9652
GTRI	44	-2324.2021	4736.4042	0.0348	1.0000
HKYG	40	-2336.8028	4753.6056	0.0000	1.0000
HKYIG	41	-2336.5538	4755.1075	0.0000	1.0000
HKYI	40	-2338.4882	4756.9764	0.0000	1.0000
SYMG	41	-2370.8972	4823.7943	0.0000	1.0000
SYMIG	42	-2370.7250	4825.4500	0.0000	1.0000
SYMI	41	-2372.5161	4827.0322	0.0000	1.0000
K2PG	37	-2377.3244	4828.6487	0.0000	1.0000
K2PIG	38	-2377.2194	4830.4388	0.0000	1.0000
K2PI	37	-2378.7503	4831.5005	0.0000	1.0000
GTR	43	-2374.9160	4835.8319	0.0000	1.0000
HKY	39	-2390.7472	4859.4944	0.0000	1.0000
F81G	39	-2411.8498	4901.6996	0.0000	1.0000
F81IG	40	-2411.7875	4903.5750	0.0000	1.0000
F81I	39	-2413.2315	4904.4630	0.0000	1.0000
SYM	40	-2427.6489	4935.2978	0.0000	1.0000
K2P	36	-2435.2893	4942.5786	0.0000	1.0000
JC69IG	37	-2442.3172	4958.6344	0.0000	1.0000
JC69G	36	-2443.3380	4958.6760	0.0000	1.0000
JC69I	36	-2443.3930	4958.7860	0.0000	1.0000
F81	38	-2466.5770	5009.1541	0.0000	1.0000
JC69	35	-2498.0165	5066.0330	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-2321.3596	4926.8170	0.8175	0.8175
HKYG	40	-2336.8028	4931.8764	0.0651	0.8826
GTRIG	45	-2320.8423	4932.2393	0.0543	0.9370
GTRI	44	-2324.2021	4932.5020	0.0476	0.9846
HKYI	40	-2338.4882	4935.2471	0.0121	0.9967
HKYIG	41	-2336.5538	4937.8351	0.0033	1.0000
K2PG	37	-2377.3244	4993.5492	0.0000	1.0000
K2PI	37	-2378.7503	4996.4010	0.0000	1.0000
K2PIG	38	-2377.2194	4999.7960	0.0000	1.0000
SYMG	41	-2370.8972	5006.5219	0.0000	1.0000
SYMI	41	-2372.5161	5009.7598	0.0000	1.0000
SYMIG	42	-2370.7250	5012.6343	0.0000	1.0000
GTR	43	-2374.9160	5027.4730	0.0000	1.0000
HKY	39	-2390.7472	5033.3085	0.0000	1.0000
F81G	39	-2411.8498	5075.5136	0.0000	1.0000
F81I	39	-2413.2315	5078.2770	0.0000	1.0000
F81IG	40	-2411.7875	5081.8458	0.0000	1.0000
K2P	36	-2435.2893	5103.0223	0.0000	1.0000
SYM	40	-2427.6489	5113.5686	0.0000	1.0000
JC69G	36	-2443.3380	5119.1197	0.0000	1.0000
JC69I	36	-2443.3930	5119.2297	0.0000	1.0000
JC69IG	37	-2442.3172	5123.5348	0.0000	1.0000
F81	38	-2466.5770	5178.5113	0.0000	1.0000
JC69	35	-2498.0165	5222.0199	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/chr5_14720.phylip = 1 - 637;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_14720.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/chr5_14720.phylip = 1 - 637;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_14720.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 (GTRG)]

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


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

