
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:20 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/chr11_3129.phylip" (19 taxa, 1119 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	-5425.9654	10943.6180	0.5683	0.5683
GTRIG	45	-5425.1547	10944.1676	0.4317	1.0000
SYMG	41	-5447.7598	10980.7174	0.0000	1.0000
SYMIG	42	-5448.3201	10983.9970	0.0000	1.0000
HKYG	40	-5462.5707	11008.1840	0.0000	1.0000
HKYIG	41	-5461.6048	11008.4073	0.0000	1.0000
K2PG	37	-5467.6522	11011.9057	0.0000	1.0000
K2PIG	38	-5466.7114	11012.1672	0.0000	1.0000
GTRI	44	-5472.0187	11035.7245	0.0000	1.0000
SYMI	41	-5494.8883	11074.9744	0.0000	1.0000
HKYI	40	-5510.9391	11104.9208	0.0000	1.0000
K2PI	37	-5518.0529	11112.7070	0.0000	1.0000
F81G	39	-5589.0266	11258.9447	0.0000	1.0000
F81IG	40	-5588.9324	11260.9074	0.0000	1.0000
JC69G	36	-5593.4170	11261.2961	0.0000	1.0000
JC69IG	37	-5593.1314	11262.8642	0.0000	1.0000
F81I	39	-5631.4533	11343.7981	0.0000	1.0000
JC69I	36	-5636.2454	11346.9529	0.0000	1.0000
GTR	43	-5659.9803	11409.4805	0.0000	1.0000
SYM	40	-5688.9214	11460.8855	0.0000	1.0000
HKY	39	-5704.6605	11490.2125	0.0000	1.0000
K2P	36	-5715.8845	11506.2312	0.0000	1.0000
F81	38	-5820.5810	11719.9065	0.0000	1.0000
JC69	35	-5826.9691	11726.2650	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-5425.9654	10939.9309	0.5472	0.5472
GTRIG	45	-5425.1547	10940.3093	0.4528	1.0000
SYMG	41	-5447.7598	10977.5196	0.0000	1.0000
SYMIG	42	-5448.3201	10980.6401	0.0000	1.0000
HKYG	40	-5462.5707	11005.1414	0.0000	1.0000
HKYIG	41	-5461.6048	11005.2095	0.0000	1.0000
K2PG	37	-5467.6522	11009.3044	0.0000	1.0000
K2PIG	38	-5466.7114	11009.4228	0.0000	1.0000
GTRI	44	-5472.0187	11032.0373	0.0000	1.0000
SYMI	41	-5494.8883	11071.7766	0.0000	1.0000
HKYI	40	-5510.9391	11101.8781	0.0000	1.0000
K2PI	37	-5518.0529	11110.1057	0.0000	1.0000
F81G	39	-5589.0266	11256.0532	0.0000	1.0000
F81IG	40	-5588.9324	11257.8648	0.0000	1.0000
JC69G	36	-5593.4170	11258.8340	0.0000	1.0000
JC69IG	37	-5593.1314	11260.2629	0.0000	1.0000
F81I	39	-5631.4533	11340.9065	0.0000	1.0000
JC69I	36	-5636.2454	11344.4908	0.0000	1.0000
GTR	43	-5659.9803	11405.9605	0.0000	1.0000
SYM	40	-5688.9214	11457.8429	0.0000	1.0000
HKY	39	-5704.6605	11487.3209	0.0000	1.0000
K2P	36	-5715.8845	11503.7691	0.0000	1.0000
F81	38	-5820.5810	11717.1620	0.0000	1.0000
JC69	35	-5826.9691	11723.9382	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-5425.9654	11160.8193	0.9370	0.9370
GTRIG	45	-5425.1547	11166.2179	0.0630	1.0000
SYMG	41	-5447.7598	11183.3474	0.0000	1.0000
SYMIG	42	-5448.3201	11191.4881	0.0000	1.0000
K2PG	37	-5467.6522	11195.0514	0.0000	1.0000
K2PIG	38	-5466.7114	11200.1900	0.0000	1.0000
HKYG	40	-5462.5707	11205.9490	0.0000	1.0000
HKYIG	41	-5461.6048	11211.0373	0.0000	1.0000
GTRI	44	-5472.0187	11252.9257	0.0000	1.0000
SYMI	41	-5494.8883	11277.6044	0.0000	1.0000
K2PI	37	-5518.0529	11295.8528	0.0000	1.0000
HKYI	40	-5510.9391	11302.6858	0.0000	1.0000
JC69G	36	-5593.4170	11439.5609	0.0000	1.0000
JC69IG	37	-5593.1314	11446.0099	0.0000	1.0000
F81G	39	-5589.0266	11451.8406	0.0000	1.0000
F81IG	40	-5588.9324	11458.6724	0.0000	1.0000
JC69I	36	-5636.2454	11525.2176	0.0000	1.0000
F81I	39	-5631.4533	11536.6940	0.0000	1.0000
GTR	43	-5659.9803	11621.8287	0.0000	1.0000
SYM	40	-5688.9214	11658.6505	0.0000	1.0000
HKY	39	-5704.6605	11683.1084	0.0000	1.0000
K2P	36	-5715.8845	11684.4960	0.0000	1.0000
JC69	35	-5826.9691	11899.6449	0.0000	1.0000
F81	38	-5820.5810	11907.9293	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/chr11_3129.phylip = 1 - 1119;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_3129.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/chr11_3129.phylip = 1 - 1119;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_3129.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/chr11_3129.phylip = 1 - 1119;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_3129.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

