
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:25 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/chr6_7207.phylip" (19 taxa, 595 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	-2511.4174	5118.0349	0.9769	0.9769
GTRIG	45	-2513.9945	5125.5300	0.0230	0.9999
HKYG	40	-2525.4622	5136.8450	0.0001	1.0000
GTRI	44	-2522.1841	5139.5682	0.0000	1.0000
HKYIG	41	-2526.1562	5140.5403	0.0000	1.0000
HKYI	40	-2533.2215	5152.3635	0.0000	1.0000
SYMG	41	-2540.8295	5169.8869	0.0000	1.0000
K2PG	37	-2545.6241	5170.2966	0.0000	1.0000
GTR	43	-2541.1410	5175.1495	0.0000	1.0000
K2PIG	38	-2547.2121	5175.7550	0.0000	1.0000
SYMIG	42	-2544.0879	5178.7193	0.0000	1.0000
SYMI	41	-2548.0392	5184.3063	0.0000	1.0000
K2PI	37	-2553.4458	5185.9400	0.0000	1.0000
HKY	39	-2557.8150	5199.2516	0.0000	1.0000
SYM	40	-2569.1142	5224.1491	0.0000	1.0000
K2P	36	-2578.6950	5234.1642	0.0000	1.0000
F81G	39	-2586.2421	5256.1059	0.0000	1.0000
F81IG	40	-2586.2371	5258.3948	0.0000	1.0000
F81I	39	-2591.2705	5266.1626	0.0000	1.0000
JC69G	36	-2602.5905	5281.9552	0.0000	1.0000
JC69IG	37	-2605.6753	5290.3991	0.0000	1.0000
JC69I	36	-2611.1152	5299.0046	0.0000	1.0000
F81	38	-2612.0888	5305.5086	0.0000	1.0000
JC69	35	-2629.0484	5332.6048	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2511.4174	5110.8349	0.9727	0.9727
GTRIG	45	-2513.9945	5117.9890	0.0272	0.9999
HKYG	40	-2525.4622	5130.9245	0.0000	1.0000
GTRI	44	-2522.1841	5132.3682	0.0000	1.0000
HKYIG	41	-2526.1562	5134.3125	0.0000	1.0000
HKYI	40	-2533.2215	5146.4429	0.0000	1.0000
SYMG	41	-2540.8295	5163.6591	0.0000	1.0000
K2PG	37	-2545.6241	5165.2481	0.0000	1.0000
GTR	43	-2541.1410	5168.2820	0.0000	1.0000
K2PIG	38	-2547.2121	5170.4241	0.0000	1.0000
SYMIG	42	-2544.0879	5172.1759	0.0000	1.0000
SYMI	41	-2548.0392	5178.0785	0.0000	1.0000
K2PI	37	-2553.4458	5180.8915	0.0000	1.0000
HKY	39	-2557.8150	5193.6300	0.0000	1.0000
SYM	40	-2569.1142	5218.2285	0.0000	1.0000
K2P	36	-2578.6950	5229.3900	0.0000	1.0000
F81G	39	-2586.2421	5250.4843	0.0000	1.0000
F81IG	40	-2586.2371	5252.4742	0.0000	1.0000
F81I	39	-2591.2705	5260.5409	0.0000	1.0000
JC69G	36	-2602.5905	5277.1810	0.0000	1.0000
JC69IG	37	-2605.6753	5285.3506	0.0000	1.0000
JC69I	36	-2611.1152	5294.2304	0.0000	1.0000
F81	38	-2612.0888	5300.1776	0.0000	1.0000
JC69	35	-2629.0484	5328.0967	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-2511.4174	5303.9316	0.7749	0.7749
HKYG	40	-2525.4622	5306.4669	0.2181	0.9930
HKYIG	41	-2526.1562	5314.2435	0.0045	0.9975
GTRIG	45	-2513.9945	5315.4742	0.0024	0.9999
HKYI	40	-2533.2215	5321.9854	0.0001	1.0000
GTRI	44	-2522.1841	5325.4649	0.0000	1.0000
K2PG	37	-2545.6241	5327.6249	0.0000	1.0000
K2PIG	38	-2547.2121	5337.1894	0.0000	1.0000
K2PI	37	-2553.4458	5343.2683	0.0000	1.0000
SYMG	41	-2540.8295	5343.5901	0.0000	1.0000
SYMIG	42	-2544.0879	5356.4954	0.0000	1.0000
GTR	43	-2541.1410	5356.9901	0.0000	1.0000
SYMI	41	-2548.0392	5358.0095	0.0000	1.0000
HKY	39	-2557.8150	5364.7839	0.0000	1.0000
K2P	36	-2578.6950	5387.3782	0.0000	1.0000
SYM	40	-2569.1142	5393.7709	0.0000	1.0000
F81G	39	-2586.2421	5421.6382	0.0000	1.0000
F81IG	40	-2586.2371	5428.0167	0.0000	1.0000
F81I	39	-2591.2705	5431.6948	0.0000	1.0000
JC69G	36	-2602.5905	5435.1692	0.0000	1.0000
JC69IG	37	-2605.6753	5447.7274	0.0000	1.0000
JC69I	36	-2611.1152	5452.2186	0.0000	1.0000
F81	38	-2612.0888	5466.9430	0.0000	1.0000
JC69	35	-2629.0484	5481.6964	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/chr6_7207.phylip = 1 - 595;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7207.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/chr6_7207.phylip = 1 - 595;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7207.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/chr6_7207.phylip = 1 - 595;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7207.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

