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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	45	-2807.6599	5711.2767	0.6550	0.6550
GTRG	44	-2809.7405	5713.1706	0.2541	0.9091
HKYIG	41	-2814.5881	5716.1032	0.0586	0.9677
HKYG	40	-2816.3139	5717.3135	0.0320	0.9998
GTRI	44	-2816.9204	5727.5305	0.0002	1.0000
HKYI	40	-2822.9987	5730.6832	0.0000	1.0000
SYMIG	42	-2856.7699	5802.7146	0.0000	1.0000
SYMG	41	-2858.1282	5803.1834	0.0000	1.0000
K2PIG	38	-2864.6703	5809.5629	0.0000	1.0000
K2PG	37	-2866.0902	5810.1805	0.0000	1.0000
SYMI	41	-2866.7553	5820.4377	0.0000	1.0000
K2PI	37	-2873.5592	5825.1183	0.0000	1.0000
F81IG	40	-2881.2000	5847.0857	0.0000	1.0000
F81G	39	-2883.3363	5849.1233	0.0000	1.0000
F81I	39	-2890.2426	5862.9360	0.0000	1.0000
JC69G	36	-2926.2157	5928.2156	0.0000	1.0000
JC69IG	37	-2925.1616	5928.3233	0.0000	1.0000
GTR	43	-2921.6740	5934.7769	0.0000	1.0000
HKY	39	-2926.9138	5936.2784	0.0000	1.0000
JC69I	36	-2932.5589	5940.9020	0.0000	1.0000
SYM	40	-2970.2260	6025.1377	0.0000	1.0000
K2P	36	-2975.5316	6026.8472	0.0000	1.0000
F81	38	-2987.8861	6055.9945	0.0000	1.0000
JC69	35	-3029.5145	6132.6034	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-2807.6599	5705.3199	0.7047	0.7047
GTRG	44	-2809.7405	5707.4810	0.2392	0.9439
HKYIG	41	-2814.5881	5711.1762	0.0377	0.9816
HKYG	40	-2816.3139	5712.6278	0.0182	0.9998
GTRI	44	-2816.9204	5721.8408	0.0002	1.0000
HKYI	40	-2822.9987	5725.9974	0.0000	1.0000
SYMIG	42	-2856.7699	5797.5398	0.0000	1.0000
SYMG	41	-2858.1282	5798.2563	0.0000	1.0000
K2PIG	38	-2864.6703	5805.3407	0.0000	1.0000
K2PG	37	-2866.0902	5806.1805	0.0000	1.0000
SYMI	41	-2866.7553	5815.5107	0.0000	1.0000
K2PI	37	-2873.5592	5821.1183	0.0000	1.0000
F81IG	40	-2881.2000	5842.4000	0.0000	1.0000
F81G	39	-2883.3363	5844.6725	0.0000	1.0000
F81I	39	-2890.2426	5858.4852	0.0000	1.0000
JC69IG	37	-2925.1616	5924.3233	0.0000	1.0000
JC69G	36	-2926.2157	5924.4315	0.0000	1.0000
GTR	43	-2921.6740	5929.3479	0.0000	1.0000
HKY	39	-2926.9138	5931.8277	0.0000	1.0000
JC69I	36	-2932.5589	5937.1179	0.0000	1.0000
SYM	40	-2970.2260	6020.4520	0.0000	1.0000
K2P	36	-2975.5316	6023.0631	0.0000	1.0000
F81	38	-2987.8861	6051.7723	0.0000	1.0000
JC69	35	-3029.5145	6129.0289	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2816.3139	5896.9478	0.8269	0.8269
HKYIG	41	-2814.5881	5900.1042	0.1706	0.9976
GTRG	44	-2809.7405	5910.2330	0.0011	0.9986
HKYI	40	-2822.9987	5910.3175	0.0010	0.9997
GTRIG	45	-2807.6599	5912.6799	0.0003	1.0000
GTRI	44	-2816.9204	5924.5928	0.0000	1.0000
K2PG	37	-2866.0902	5976.6765	0.0000	1.0000
K2PIG	38	-2864.6703	5980.4447	0.0000	1.0000
SYMG	41	-2858.1282	5987.1843	0.0000	1.0000
SYMIG	42	-2856.7699	5991.0758	0.0000	1.0000
K2PI	37	-2873.5592	5991.6143	0.0000	1.0000
SYMI	41	-2866.7553	6004.4387	0.0000	1.0000
F81G	39	-2883.3363	6024.3845	0.0000	1.0000
F81IG	40	-2881.2000	6026.7200	0.0000	1.0000
F81I	39	-2890.2426	6038.1973	0.0000	1.0000
JC69G	36	-2926.2157	6090.3195	0.0000	1.0000
JC69IG	37	-2925.1616	6094.8193	0.0000	1.0000
JC69I	36	-2932.5589	6103.0059	0.0000	1.0000
HKY	39	-2926.9138	6111.5397	0.0000	1.0000
GTR	43	-2921.6740	6127.4920	0.0000	1.0000
K2P	36	-2975.5316	6188.9511	0.0000	1.0000
SYM	40	-2970.2260	6204.7720	0.0000	1.0000
F81	38	-2987.8861	6226.8763	0.0000	1.0000
JC69	35	-3029.5145	6290.3090	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/chrZ_6047.phylip = 1 - 741;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6047.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 (GTRIG)]

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

