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

Minimum AIC  model: GTRIG
Minimum AICc model: HKYIG
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
HKYIG	41	-2826.2257	5738.1467	0.3331	0.3331
GTRIG	45	-2822.0184	5738.4980	0.2794	0.6125
HKYG	40	-2827.7963	5739.1081	0.2060	0.8185
GTRG	44	-2823.6985	5739.6597	0.1563	0.9748
HKYI	40	-2830.3707	5744.2570	0.0157	0.9905
GTRI	44	-2826.5068	5745.2763	0.0094	0.9999
K2PIG	38	-2838.8670	5756.9040	0.0000	1.0000
K2PG	37	-2840.4030	5757.8103	0.0000	1.0000
SYMIG	42	-2835.3760	5758.6316	0.0000	1.0000
SYMG	41	-2837.0261	5759.7475	0.0000	1.0000
K2PI	37	-2842.9160	5762.8363	0.0000	1.0000
SYMI	41	-2839.6694	5765.0341	0.0000	1.0000
HKY	39	-2869.2115	5819.7634	0.0000	1.0000
GTR	43	-2865.6913	5821.4514	0.0000	1.0000
K2P	36	-2881.2754	5837.3938	0.0000	1.0000
SYM	40	-2878.3331	5840.1817	0.0000	1.0000
F81IG	40	-2892.0296	5867.5748	0.0000	1.0000
F81G	39	-2893.2175	5867.7754	0.0000	1.0000
F81I	39	-2895.7348	5872.8101	0.0000	1.0000
JC69IG	37	-2902.3503	5881.7050	0.0000	1.0000
JC69G	36	-2903.5310	5881.9052	0.0000	1.0000
JC69I	36	-2905.9920	5886.8270	0.0000	1.0000
F81	38	-2931.2517	5941.6734	0.0000	1.0000
JC69	35	-2941.6885	5956.0636	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-2822.0184	5734.0368	0.3516	0.3516
HKYIG	41	-2826.2257	5734.4515	0.2857	0.6373
GTRG	44	-2823.6985	5735.3970	0.1781	0.8154
HKYG	40	-2827.7963	5735.5926	0.1615	0.9769
HKYI	40	-2830.3707	5740.7415	0.0123	0.9892
GTRI	44	-2826.5068	5741.0136	0.0107	1.0000
K2PIG	38	-2838.8670	5753.7339	0.0000	1.0000
SYMIG	42	-2835.3760	5754.7519	0.0000	1.0000
K2PG	37	-2840.4030	5754.8060	0.0000	1.0000
SYMG	41	-2837.0261	5756.0522	0.0000	1.0000
K2PI	37	-2842.9160	5759.8320	0.0000	1.0000
SYMI	41	-2839.6694	5761.3388	0.0000	1.0000
HKY	39	-2869.2115	5816.4229	0.0000	1.0000
GTR	43	-2865.6913	5817.3826	0.0000	1.0000
K2P	36	-2881.2754	5834.5507	0.0000	1.0000
SYM	40	-2878.3331	5836.6661	0.0000	1.0000
F81IG	40	-2892.0296	5864.0592	0.0000	1.0000
F81G	39	-2893.2175	5864.4349	0.0000	1.0000
F81I	39	-2895.7348	5869.4696	0.0000	1.0000
JC69IG	37	-2902.3503	5878.7007	0.0000	1.0000
JC69G	36	-2903.5310	5879.0621	0.0000	1.0000
JC69I	36	-2905.9920	5883.9839	0.0000	1.0000
F81	38	-2931.2517	5938.5033	0.0000	1.0000
JC69	35	-2941.6885	5953.3770	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2827.7963	5930.8490	0.7377	0.7377
HKYIG	41	-2826.2257	5934.5893	0.1137	0.8514
K2PG	37	-2840.4030	5935.4182	0.0751	0.9265
HKYI	40	-2830.3707	5935.9979	0.0562	0.9827
K2PIG	38	-2838.8670	5939.2276	0.0112	0.9939
K2PI	37	-2842.9160	5940.4442	0.0061	0.9999
GTRG	44	-2823.6985	5950.1791	0.0000	1.0000
GTRIG	45	-2822.0184	5953.7003	0.0000	1.0000
GTRI	44	-2826.5068	5955.7957	0.0000	1.0000
SYMG	41	-2837.0261	5956.1901	0.0000	1.0000
SYMIG	42	-2835.3760	5959.7712	0.0000	1.0000
SYMI	41	-2839.6694	5961.4766	0.0000	1.0000
HKY	39	-2869.2115	6006.7980	0.0000	1.0000
K2P	36	-2881.2754	6010.2815	0.0000	1.0000
GTR	43	-2865.6913	6027.2833	0.0000	1.0000
SYM	40	-2878.3331	6031.9226	0.0000	1.0000
JC69G	36	-2903.5310	6054.7929	0.0000	1.0000
F81G	39	-2893.2175	6054.8100	0.0000	1.0000
JC69IG	37	-2902.3503	6059.3129	0.0000	1.0000
F81IG	40	-2892.0296	6059.3157	0.0000	1.0000
JC69I	36	-2905.9920	6059.7147	0.0000	1.0000
F81I	39	-2895.7348	6059.8446	0.0000	1.0000
F81	38	-2931.2517	6123.9969	0.0000	1.0000
JC69	35	-2941.6885	6124.2264	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRIG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_9755.phylip = 1 - 974;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_9755.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYIG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_9755.phylip = 1 - 974;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_9755.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.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/chr3_9755.phylip = 1 - 974;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_9755.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

