
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:15 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/chr1_12757.phylip" (19 taxa, 696 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	-2348.8482	4791.7793	0.5629	0.5629
GTRIG	45	-2348.6178	4793.6049	0.2260	0.7889
GTRI	44	-2349.8288	4793.7406	0.2111	1.0000
HKYG	40	-2365.3471	4815.7019	0.0000	1.0000
HKYIG	41	-2365.2640	4817.7941	0.0000	1.0000
HKYI	40	-2366.7224	4818.4525	0.0000	1.0000
SYMG	41	-2379.4167	4846.0995	0.0000	1.0000
SYMIG	42	-2379.2777	4848.0868	0.0000	1.0000
SYMI	41	-2380.7853	4848.8367	0.0000	1.0000
GTR	43	-2380.7176	4853.2388	0.0000	1.0000
K2PG	37	-2391.9344	4862.1423	0.0000	1.0000
K2PIG	38	-2391.8636	4864.2385	0.0000	1.0000
K2PI	37	-2393.4966	4865.2667	0.0000	1.0000
F81G	39	-2394.4221	4871.6003	0.0000	1.0000
HKY	39	-2394.8652	4872.4866	0.0000	1.0000
F81I	39	-2395.4527	4873.6616	0.0000	1.0000
F81IG	40	-2394.3558	4873.7193	0.0000	1.0000
SYM	40	-2406.8626	4898.7329	0.0000	1.0000
JC69G	36	-2417.6110	4911.2645	0.0000	1.0000
JC69IG	37	-2417.5573	4913.3882	0.0000	1.0000
K2P	36	-2418.9780	4913.9985	0.0000	1.0000
JC69I	36	-2418.9914	4914.0254	0.0000	1.0000
F81	38	-2421.8434	4924.1983	0.0000	1.0000
JC69	35	-2443.7952	4961.4086	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2348.8482	4785.6964	0.5440	0.5440
GTRIG	45	-2348.6178	4787.2356	0.2520	0.7960
GTRI	44	-2349.8288	4787.6577	0.2040	1.0000
HKYG	40	-2365.3471	4810.6942	0.0000	1.0000
HKYIG	41	-2365.2640	4812.5281	0.0000	1.0000
HKYI	40	-2366.7224	4813.4449	0.0000	1.0000
SYMG	41	-2379.4167	4840.8334	0.0000	1.0000
SYMIG	42	-2379.2777	4842.5554	0.0000	1.0000
SYMI	41	-2380.7853	4843.5706	0.0000	1.0000
GTR	43	-2380.7176	4847.4352	0.0000	1.0000
K2PG	37	-2391.9344	4857.8688	0.0000	1.0000
K2PIG	38	-2391.8636	4859.7271	0.0000	1.0000
K2PI	37	-2393.4966	4860.9932	0.0000	1.0000
F81G	39	-2394.4221	4866.8442	0.0000	1.0000
HKY	39	-2394.8652	4867.7305	0.0000	1.0000
F81IG	40	-2394.3558	4868.7117	0.0000	1.0000
F81I	39	-2395.4527	4868.9055	0.0000	1.0000
SYM	40	-2406.8626	4893.7253	0.0000	1.0000
JC69G	36	-2417.6110	4907.2220	0.0000	1.0000
JC69IG	37	-2417.5573	4909.1146	0.0000	1.0000
K2P	36	-2418.9780	4909.9560	0.0000	1.0000
JC69I	36	-2418.9914	4909.9829	0.0000	1.0000
F81	38	-2421.8434	4919.6868	0.0000	1.0000
JC69	35	-2443.7952	4957.5904	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-2348.8482	4985.6918	0.6823	0.6823
GTRI	44	-2349.8288	4987.6531	0.2559	0.9382
GTRIG	45	-2348.6178	4991.7764	0.0326	0.9708
HKYG	40	-2365.3471	4992.5082	0.0226	0.9934
HKYI	40	-2366.7224	4995.2589	0.0057	0.9991
HKYIG	41	-2365.2640	4998.8874	0.0009	1.0000
K2PG	37	-2391.9344	5026.0467	0.0000	1.0000
SYMG	41	-2379.4167	5027.1928	0.0000	1.0000
K2PI	37	-2393.4966	5029.1711	0.0000	1.0000
SYMI	41	-2380.7853	5029.9300	0.0000	1.0000
K2PIG	38	-2391.8636	5032.4504	0.0000	1.0000
SYMIG	42	-2379.2777	5033.4601	0.0000	1.0000
GTR	43	-2380.7176	5042.8852	0.0000	1.0000
F81G	39	-2394.4221	5044.1129	0.0000	1.0000
HKY	39	-2394.8652	5044.9991	0.0000	1.0000
F81I	39	-2395.4527	5046.1741	0.0000	1.0000
F81IG	40	-2394.3558	5050.5257	0.0000	1.0000
JC69G	36	-2417.6110	5070.8546	0.0000	1.0000
K2P	36	-2418.9780	5073.5886	0.0000	1.0000
JC69I	36	-2418.9914	5073.6155	0.0000	1.0000
SYM	40	-2406.8626	5075.5393	0.0000	1.0000
JC69IG	37	-2417.5573	5077.2926	0.0000	1.0000
F81	38	-2421.8434	5092.4101	0.0000	1.0000
JC69	35	-2443.7952	5116.6776	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/chr1_12757.phylip = 1 - 696;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_12757.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/chr1_12757.phylip = 1 - 696;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_12757.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/chr1_12757.phylip = 1 - 696;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_12757.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

