
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:35 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_17003.phylip" (19 taxa, 599 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	-2454.4326	5004.0131	0.6737	0.6737
GTRIG	45	-2454.0646	5005.6157	0.3023	0.9761
GTRI	44	-2457.7785	5010.7050	0.0237	0.9998
HKYG	40	-2467.6386	5021.1554	0.0001	0.9999
HKYIG	41	-2467.3426	5022.8683	0.0001	1.0000
HKYI	40	-2469.3060	5024.4901	0.0000	1.0000
K2PG	37	-2488.6197	5056.2518	0.0000	1.0000
K2PIG	38	-2488.2887	5057.8703	0.0000	1.0000
K2PI	37	-2489.9125	5058.8374	0.0000	1.0000
SYMG	41	-2485.6981	5059.5794	0.0000	1.0000
SYMIG	42	-2485.2456	5060.9875	0.0000	1.0000
SYMI	41	-2487.3111	5062.8053	0.0000	1.0000
GTR	43	-2493.7541	5080.3263	0.0000	1.0000
HKY	39	-2507.5652	5098.7118	0.0000	1.0000
F81G	39	-2525.1132	5133.8078	0.0000	1.0000
F81IG	40	-2524.6501	5135.1783	0.0000	1.0000
F81I	39	-2526.2989	5136.1792	0.0000	1.0000
SYM	40	-2527.4905	5140.8592	0.0000	1.0000
K2P	36	-2532.1821	5141.1044	0.0000	1.0000
JC69G	36	-2541.3132	5159.3666	0.0000	1.0000
JC69IG	37	-2540.8730	5160.7585	0.0000	1.0000
JC69I	36	-2542.5701	5161.8805	0.0000	1.0000
F81	38	-2566.3169	5213.9267	0.0000	1.0000
JC69	35	-2584.6856	5243.8472	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2454.4326	4996.8651	0.6382	0.6382
GTRIG	45	-2454.0646	4998.1292	0.3392	0.9774
GTRI	44	-2457.7785	5003.5570	0.0225	0.9999
HKYG	40	-2467.6386	5015.2772	0.0001	1.0000
HKYIG	41	-2467.3426	5016.6852	0.0000	1.0000
HKYI	40	-2469.3060	5018.6119	0.0000	1.0000
K2PG	37	-2488.6197	5051.2393	0.0000	1.0000
K2PIG	38	-2488.2887	5052.5774	0.0000	1.0000
SYMG	41	-2485.6981	5053.3963	0.0000	1.0000
K2PI	37	-2489.9125	5053.8249	0.0000	1.0000
SYMIG	42	-2485.2456	5054.4911	0.0000	1.0000
SYMI	41	-2487.3111	5056.6222	0.0000	1.0000
GTR	43	-2493.7541	5073.5083	0.0000	1.0000
HKY	39	-2507.5652	5093.1304	0.0000	1.0000
F81G	39	-2525.1132	5128.2264	0.0000	1.0000
F81IG	40	-2524.6501	5129.3002	0.0000	1.0000
F81I	39	-2526.2989	5130.5978	0.0000	1.0000
SYM	40	-2527.4905	5134.9810	0.0000	1.0000
K2P	36	-2532.1821	5136.3642	0.0000	1.0000
JC69G	36	-2541.3132	5154.6264	0.0000	1.0000
JC69IG	37	-2540.8730	5155.7460	0.0000	1.0000
JC69I	36	-2542.5701	5157.1403	0.0000	1.0000
F81	38	-2566.3169	5208.6338	0.0000	1.0000
JC69	35	-2584.6856	5239.3712	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-2454.4326	5190.2566	0.5222	0.5222
HKYG	40	-2467.6386	5191.0877	0.3446	0.8668
HKYI	40	-2469.3060	5194.4224	0.0650	0.9318
GTRIG	45	-2454.0646	5195.9160	0.0308	0.9627
HKYIG	41	-2467.3426	5196.8909	0.0189	0.9816
GTRI	44	-2457.7785	5196.9485	0.0184	1.0000
K2PG	37	-2488.6197	5213.8640	0.0000	1.0000
K2PI	37	-2489.9125	5216.4496	0.0000	1.0000
K2PIG	38	-2488.2887	5219.5974	0.0000	1.0000
SYMG	41	-2485.6981	5233.6020	0.0000	1.0000
SYMI	41	-2487.3111	5236.8279	0.0000	1.0000
SYMIG	42	-2485.2456	5239.0921	0.0000	1.0000
GTR	43	-2493.7541	5262.5045	0.0000	1.0000
HKY	39	-2507.5652	5264.5456	0.0000	1.0000
K2P	36	-2532.1821	5294.5936	0.0000	1.0000
F81G	39	-2525.1132	5299.6416	0.0000	1.0000
F81I	39	-2526.2989	5302.0130	0.0000	1.0000
F81IG	40	-2524.6501	5305.1106	0.0000	1.0000
SYM	40	-2527.4905	5310.7915	0.0000	1.0000
JC69G	36	-2541.3132	5312.8558	0.0000	1.0000
JC69I	36	-2542.5701	5315.3697	0.0000	1.0000
JC69IG	37	-2540.8730	5318.3707	0.0000	1.0000
F81	38	-2566.3169	5375.6538	0.0000	1.0000
JC69	35	-2584.6856	5393.2053	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/chr2_17003.phylip = 1 - 599;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_17003.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/chr2_17003.phylip = 1 - 599;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_17003.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/chr2_17003.phylip = 1 - 599;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_17003.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

