
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:17 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/chr14_3774.phylip" (19 taxa, 910 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	-4305.5071	8703.5922	0.7511	0.7511
GTRIG	45	-4305.5050	8705.8016	0.2489	1.0000
HKYG	40	-4323.3672	8730.5088	0.0000	1.0000
HKYIG	41	-4323.3527	8732.6731	0.0000	1.0000
SYMG	41	-4332.0573	8750.0824	0.0000	1.0000
GTRI	44	-4329.2781	8751.1343	0.0000	1.0000
SYMIG	42	-4332.0008	8752.1676	0.0000	1.0000
K2PG	37	-4339.5156	8756.2559	0.0000	1.0000
K2PIG	38	-4339.4600	8758.3230	0.0000	1.0000
HKYI	40	-4342.9265	8769.6275	0.0000	1.0000
SYMI	41	-4351.8309	8789.6296	0.0000	1.0000
K2PI	37	-4358.4345	8794.0938	0.0000	1.0000
GTR	43	-4417.8454	8926.0603	0.0000	1.0000
HKY	39	-4435.1086	8951.8035	0.0000	1.0000
SYM	40	-4446.1700	8976.1145	0.0000	1.0000
K2P	36	-4453.2813	8981.6142	0.0000	1.0000
F81G	39	-4483.3112	9048.2087	0.0000	1.0000
F81IG	40	-4483.3136	9050.4016	0.0000	1.0000
JC69G	36	-4492.0526	9059.1567	0.0000	1.0000
JC69IG	37	-4492.0524	9061.3295	0.0000	1.0000
F81I	39	-4498.7993	9079.1848	0.0000	1.0000
JC69I	36	-4507.2189	9089.4893	0.0000	1.0000
F81	38	-4586.1352	9251.6733	0.0000	1.0000
JC69	35	-4596.6388	9266.1610	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-4305.5071	8699.0142	0.7306	0.7306
GTRIG	45	-4305.5050	8701.0100	0.2694	1.0000
HKYG	40	-4323.3672	8726.7344	0.0000	1.0000
HKYIG	41	-4323.3527	8728.7053	0.0000	1.0000
SYMG	41	-4332.0573	8746.1146	0.0000	1.0000
GTRI	44	-4329.2781	8746.5562	0.0000	1.0000
SYMIG	42	-4332.0008	8748.0015	0.0000	1.0000
K2PG	37	-4339.5156	8753.0311	0.0000	1.0000
K2PIG	38	-4339.4600	8754.9200	0.0000	1.0000
HKYI	40	-4342.9265	8765.8531	0.0000	1.0000
SYMI	41	-4351.8309	8785.6618	0.0000	1.0000
K2PI	37	-4358.4345	8790.8690	0.0000	1.0000
GTR	43	-4417.8454	8921.6908	0.0000	1.0000
HKY	39	-4435.1086	8948.2173	0.0000	1.0000
SYM	40	-4446.1700	8972.3401	0.0000	1.0000
K2P	36	-4453.2813	8978.5626	0.0000	1.0000
F81G	39	-4483.3112	9044.6225	0.0000	1.0000
F81IG	40	-4483.3136	9046.6271	0.0000	1.0000
JC69G	36	-4492.0526	9056.1052	0.0000	1.0000
JC69IG	37	-4492.0524	9058.1048	0.0000	1.0000
F81I	39	-4498.7993	9075.5986	0.0000	1.0000
JC69I	36	-4507.2189	9086.4377	0.0000	1.0000
F81	38	-4586.1352	9248.2704	0.0000	1.0000
JC69	35	-4596.6388	9263.2777	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-4305.5071	8910.8057	0.9540	0.9540
GTRIG	45	-4305.5050	8917.6150	0.0317	0.9857
HKYG	40	-4323.3672	8919.2722	0.0138	0.9995
HKYIG	41	-4323.3527	8926.0566	0.0005	1.0000
K2PG	37	-4339.5156	8931.1286	0.0000	1.0000
K2PIG	38	-4339.4600	8937.8309	0.0000	1.0000
SYMG	41	-4332.0573	8943.4658	0.0000	1.0000
SYMIG	42	-4332.0008	8950.1662	0.0000	1.0000
GTRI	44	-4329.2781	8958.3478	0.0000	1.0000
HKYI	40	-4342.9265	8958.3908	0.0000	1.0000
K2PI	37	-4358.4345	8968.9664	0.0000	1.0000
SYMI	41	-4351.8309	8983.0131	0.0000	1.0000
GTR	43	-4417.8454	9128.6689	0.0000	1.0000
HKY	39	-4435.1086	9135.9416	0.0000	1.0000
K2P	36	-4453.2813	9151.8466	0.0000	1.0000
SYM	40	-4446.1700	9164.8779	0.0000	1.0000
JC69G	36	-4492.0526	9229.3892	0.0000	1.0000
F81G	39	-4483.3112	9232.3468	0.0000	1.0000
JC69IG	37	-4492.0524	9236.2022	0.0000	1.0000
F81IG	40	-4483.3136	9239.1649	0.0000	1.0000
JC69I	36	-4507.2189	9259.7217	0.0000	1.0000
F81I	39	-4498.7993	9263.3229	0.0000	1.0000
F81	38	-4586.1352	9431.1813	0.0000	1.0000
JC69	35	-4596.6388	9431.7482	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/chr14_3774.phylip = 1 - 910;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr14_3774.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/chr14_3774.phylip = 1 - 910;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr14_3774.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/chr14_3774.phylip = 1 - 910;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr14_3774.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

