
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:59 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_27825.phylip" (19 taxa, 806 characters)

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
Minimum BIC  model: HKYIG


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	-2509.3995	5105.3069	0.6910	0.6910
GTRG	44	-2507.1502	5107.5041	0.2303	0.9213
GTRIG	45	-2507.3110	5110.0695	0.0639	0.9852
HKYG	40	-2514.4916	5113.2708	0.0129	0.9981
HKYI	40	-2517.1225	5118.5326	0.0009	0.9990
SYMIG	42	-2515.6897	5120.1134	0.0004	0.9994
GTRI	44	-2513.8008	5120.8054	0.0003	0.9997
K2PIG	38	-2520.7273	5121.3190	0.0002	1.0000
SYMG	41	-2520.1884	5126.8846	0.0000	1.0000
K2PG	37	-2525.4579	5128.5772	0.0000	1.0000
SYMI	41	-2522.2732	5131.0542	0.0000	1.0000
K2PI	37	-2527.9088	5133.4791	0.0000	1.0000
F81G	39	-2580.8619	5243.7969	0.0000	1.0000
JC69G	36	-2590.3194	5256.1030	0.0000	1.0000
F81IG	40	-2587.0012	5258.2899	0.0000	1.0000
F81I	39	-2592.7518	5267.5766	0.0000	1.0000
JC69IG	37	-2596.4811	5270.6236	0.0000	1.0000
JC69I	36	-2602.0110	5279.4862	0.0000	1.0000
HKY	39	-2600.0507	5282.1744	0.0000	1.0000
K2P	36	-2609.7720	5295.0083	0.0000	1.0000
GTR	43	-2604.8787	5300.7234	0.0000	1.0000
SYM	40	-2612.9529	5310.1933	0.0000	1.0000
F81	38	-2663.2214	5406.3072	0.0000	1.0000
JC69	35	-2684.9448	5443.1623	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	41	-2509.3995	5100.7990	0.6100	0.6100
GTRG	44	-2507.1502	5102.3004	0.2879	0.8979
GTRIG	45	-2507.3110	5104.6221	0.0902	0.9881
HKYG	40	-2514.4916	5108.9832	0.0102	0.9983
HKYI	40	-2517.1225	5114.2450	0.0007	0.9990
SYMIG	42	-2515.6897	5115.3794	0.0004	0.9995
GTRI	44	-2513.8008	5115.6017	0.0004	0.9998
K2PIG	38	-2520.7273	5117.4546	0.0001	1.0000
SYMG	41	-2520.1884	5122.3767	0.0000	1.0000
K2PG	37	-2525.4579	5124.9157	0.0000	1.0000
SYMI	41	-2522.2732	5126.5464	0.0000	1.0000
K2PI	37	-2527.9088	5129.8176	0.0000	1.0000
F81G	39	-2580.8619	5239.7238	0.0000	1.0000
JC69G	36	-2590.3194	5252.6388	0.0000	1.0000
F81IG	40	-2587.0012	5254.0023	0.0000	1.0000
F81I	39	-2592.7518	5263.5035	0.0000	1.0000
JC69IG	37	-2596.4811	5266.9622	0.0000	1.0000
JC69I	36	-2602.0110	5276.0219	0.0000	1.0000
HKY	39	-2600.0507	5278.1013	0.0000	1.0000
K2P	36	-2609.7720	5291.5441	0.0000	1.0000
GTR	43	-2604.8787	5295.7575	0.0000	1.0000
SYM	40	-2612.9529	5305.9057	0.0000	1.0000
F81	38	-2663.2214	5402.4428	0.0000	1.0000
JC69	35	-2684.9448	5439.8895	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYIG	41	-2509.3995	5293.1744	0.6503	0.6503
K2PIG	38	-2520.7273	5295.7537	0.1791	0.8293
HKYG	40	-2514.4916	5296.6666	0.1134	0.9428
K2PG	37	-2525.4579	5298.5228	0.0448	0.9876
HKYI	40	-2517.1225	5301.9283	0.0082	0.9958
K2PI	37	-2527.9088	5303.4247	0.0039	0.9997
GTRG	44	-2507.1502	5308.7521	0.0003	0.9999
SYMIG	42	-2515.6897	5312.4469	0.0000	1.0000
SYMG	41	-2520.1884	5314.7522	0.0000	1.0000
GTRIG	45	-2507.3110	5315.7659	0.0000	1.0000
SYMI	41	-2522.2732	5318.9218	0.0000	1.0000
GTRI	44	-2513.8008	5322.0534	0.0000	1.0000
JC69G	36	-2590.3194	5421.5538	0.0000	1.0000
F81G	39	-2580.8619	5422.7150	0.0000	1.0000
JC69IG	37	-2596.4811	5440.5693	0.0000	1.0000
F81IG	40	-2587.0012	5441.6857	0.0000	1.0000
JC69I	36	-2602.0110	5444.9370	0.0000	1.0000
F81I	39	-2592.7518	5446.4948	0.0000	1.0000
K2P	36	-2609.7720	5460.4591	0.0000	1.0000
HKY	39	-2600.0507	5461.0926	0.0000	1.0000
SYM	40	-2612.9529	5493.5891	0.0000	1.0000
GTR	43	-2604.8787	5497.5171	0.0000	1.0000
F81	38	-2663.2214	5580.7420	0.0000	1.0000
JC69	35	-2684.9448	5604.1125	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_27825.phylip = 1 - 806;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_27825.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 AICc model (HKYIG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_27825.phylip = 1 - 806;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_27825.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 (HKYIG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_27825.phylip = 1 - 806;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_27825.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;


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

