
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:32 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/chrZ_8694.phylip" (19 taxa, 777 characters)

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


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	-2284.0581	4661.5260	0.3413	0.3413
HKYG	40	-2288.9562	4662.3689	0.2239	0.5652
GTRIG	45	-2283.5106	4662.6847	0.1912	0.7564
HKYIG	41	-2288.3461	4663.3779	0.1352	0.8916
GTRI	44	-2285.5869	4664.5835	0.0740	0.9656
HKYI	40	-2290.8292	4666.1148	0.0344	1.0000
SYMG	41	-2304.7166	4696.1190	0.0000	1.0000
SYMIG	42	-2304.0244	4696.9698	0.0000	1.0000
K2PG	37	-2311.1024	4700.0099	0.0000	1.0000
SYMI	41	-2306.9464	4700.5785	0.0000	1.0000
K2PIG	38	-2310.3883	4700.7928	0.0000	1.0000
K2PI	37	-2313.3956	4704.5964	0.0000	1.0000
HKY	39	-2321.2249	4724.6832	0.0000	1.0000
GTR	43	-2317.0598	4725.2819	0.0000	1.0000
SYM	40	-2336.6847	4757.8259	0.0000	1.0000
K2P	36	-2343.1523	4761.9047	0.0000	1.0000
F81G	39	-2353.4407	4789.1148	0.0000	1.0000
F81IG	40	-2352.8371	4790.1307	0.0000	1.0000
F81I	39	-2355.5952	4793.4238	0.0000	1.0000
JC69G	36	-2374.1938	4823.9876	0.0000	1.0000
JC69IG	37	-2374.9528	4827.7108	0.0000	1.0000
JC69I	36	-2377.6447	4830.8895	0.0000	1.0000
F81	38	-2385.2364	4850.4890	0.0000	1.0000
JC69	35	-2406.4088	4886.2184	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2284.0581	4656.1162	0.3848	0.3848
GTRIG	45	-2283.5106	4657.0212	0.2448	0.6296
HKYG	40	-2288.9562	4657.9124	0.1568	0.7863
HKYIG	41	-2288.3461	4658.6922	0.1061	0.8925
GTRI	44	-2285.5869	4659.1737	0.0834	0.9759
HKYI	40	-2290.8292	4661.6583	0.0241	1.0000
SYMG	41	-2304.7166	4691.4332	0.0000	1.0000
SYMIG	42	-2304.0244	4692.0488	0.0000	1.0000
SYMI	41	-2306.9464	4695.8928	0.0000	1.0000
K2PG	37	-2311.1024	4696.2048	0.0000	1.0000
K2PIG	38	-2310.3883	4696.7765	0.0000	1.0000
K2PI	37	-2313.3956	4700.7913	0.0000	1.0000
GTR	43	-2317.0598	4720.1195	0.0000	1.0000
HKY	39	-2321.2249	4720.4498	0.0000	1.0000
SYM	40	-2336.6847	4753.3693	0.0000	1.0000
K2P	36	-2343.1523	4758.3047	0.0000	1.0000
F81G	39	-2353.4407	4784.8815	0.0000	1.0000
F81IG	40	-2352.8371	4785.6742	0.0000	1.0000
F81I	39	-2355.5952	4789.1904	0.0000	1.0000
JC69G	36	-2374.1938	4820.3876	0.0000	1.0000
JC69IG	37	-2374.9528	4823.9056	0.0000	1.0000
JC69I	36	-2377.6447	4827.2895	0.0000	1.0000
F81	38	-2385.2364	4846.4728	0.0000	1.0000
JC69	35	-2406.4088	4882.8176	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2288.9562	4844.1300	0.8197	0.8197
HKYI	40	-2290.8292	4847.8759	0.1260	0.9456
HKYIG	41	-2288.3461	4849.5652	0.0541	0.9998
GTRG	44	-2284.0581	4860.9556	0.0002	0.9999
GTRI	44	-2285.5869	4864.0131	0.0000	1.0000
GTRIG	45	-2283.5106	4866.5161	0.0000	1.0000
K2PG	37	-2311.1024	4868.4561	0.0000	1.0000
K2PI	37	-2313.3956	4873.0426	0.0000	1.0000
K2PIG	38	-2310.3883	4873.6833	0.0000	1.0000
SYMG	41	-2304.7166	4882.3063	0.0000	1.0000
SYMI	41	-2306.9464	4886.7658	0.0000	1.0000
SYMIG	42	-2304.0244	4887.5773	0.0000	1.0000
HKY	39	-2321.2249	4902.0120	0.0000	1.0000
GTR	43	-2317.0598	4920.3035	0.0000	1.0000
K2P	36	-2343.1523	4925.9005	0.0000	1.0000
SYM	40	-2336.6847	4939.5870	0.0000	1.0000
F81G	39	-2353.4407	4966.4436	0.0000	1.0000
F81I	39	-2355.5952	4970.7526	0.0000	1.0000
F81IG	40	-2352.8371	4971.8918	0.0000	1.0000
JC69G	36	-2374.1938	4987.9834	0.0000	1.0000
JC69I	36	-2377.6447	4994.8853	0.0000	1.0000
JC69IG	37	-2374.9528	4996.1569	0.0000	1.0000
F81	38	-2385.2364	5023.3795	0.0000	1.0000
JC69	35	-2406.4088	5045.7581	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/chrZ_8694.phylip = 1 - 777;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_8694.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/chrZ_8694.phylip = 1 - 777;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_8694.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 (HKYG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_8694.phylip = 1 - 777;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_8694.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

