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

Minimum AIC  model: SYMG
Minimum AICc model: SYMG
Minimum BIC  model: K2PG


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
SYMG	41	-4050.4690	8186.0632	0.4846	0.4846
SYMIG	42	-4049.4811	8186.2429	0.4430	0.9276
GTRG	44	-4050.0743	8191.7519	0.0282	0.9558
GTRIG	45	-4049.1023	8191.9750	0.0252	0.9810
SYMI	41	-4054.4400	8194.0052	0.0091	0.9901
K2PG	37	-4059.4479	8195.4383	0.0045	0.9946
GTRI	44	-4052.2647	8196.1327	0.0032	0.9978
K2PIG	38	-4059.4252	8197.5327	0.0016	0.9993
HKYG	40	-4058.4766	8199.9269	0.0005	0.9998
HKYIG	41	-4058.4505	8202.0262	0.0002	1.0000
K2PI	37	-4064.1715	8204.8855	0.0000	1.0000
HKYI	40	-4063.1812	8209.3361	0.0000	1.0000
SYM	40	-4092.5826	8268.1390	0.0000	1.0000
GTR	43	-4091.8384	8273.1168	0.0000	1.0000
K2P	36	-4103.4053	8281.2171	0.0000	1.0000
JC69G	36	-4104.6602	8283.7270	0.0000	1.0000
HKY	39	-4102.4748	8285.7757	0.0000	1.0000
F81G	39	-4103.5945	8288.0151	0.0000	1.0000
JC69IG	37	-4105.8181	8288.1787	0.0000	1.0000
JC69I	36	-4109.0096	8292.4257	0.0000	1.0000
F81IG	40	-4104.7677	8292.5091	0.0000	1.0000
F81I	39	-4107.9508	8296.7277	0.0000	1.0000
JC69	35	-4145.8897	8364.0537	0.0000	1.0000
F81	38	-4144.7878	8368.2580	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	41	-4050.4690	8182.9380	0.4605	0.4605
SYMIG	42	-4049.4811	8182.9622	0.4549	0.9154
GTRG	44	-4050.0743	8188.1486	0.0340	0.9495
GTRIG	45	-4049.1023	8188.2045	0.0331	0.9825
SYMI	41	-4054.4400	8190.8800	0.0087	0.9912
GTRI	44	-4052.2647	8192.5294	0.0038	0.9950
K2PG	37	-4059.4479	8192.8958	0.0032	0.9982
K2PIG	38	-4059.4252	8194.8504	0.0012	0.9994
HKYG	40	-4058.4766	8196.9532	0.0004	0.9998
HKYIG	41	-4058.4505	8198.9010	0.0002	1.0000
K2PI	37	-4064.1715	8202.3430	0.0000	1.0000
HKYI	40	-4063.1812	8206.3624	0.0000	1.0000
SYM	40	-4092.5826	8265.1653	0.0000	1.0000
GTR	43	-4091.8384	8269.6768	0.0000	1.0000
K2P	36	-4103.4053	8278.8106	0.0000	1.0000
JC69G	36	-4104.6602	8281.3205	0.0000	1.0000
HKY	39	-4102.4748	8282.9497	0.0000	1.0000
F81G	39	-4103.5945	8285.1890	0.0000	1.0000
JC69IG	37	-4105.8181	8285.6362	0.0000	1.0000
F81IG	40	-4104.7677	8289.5354	0.0000	1.0000
JC69I	36	-4109.0096	8290.0192	0.0000	1.0000
F81I	39	-4107.9508	8293.9016	0.0000	1.0000
JC69	35	-4145.8897	8361.7794	0.0000	1.0000
F81	38	-4144.7878	8365.5757	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-4059.4479	8379.4604	0.9562	0.9562
K2PIG	38	-4059.4252	8386.4573	0.0289	0.9851
K2PI	37	-4064.1715	8388.9076	0.0085	0.9936
SYMG	41	-4050.4690	8389.6717	0.0058	0.9994
SYMIG	42	-4049.4811	8394.7382	0.0005	0.9998
SYMI	41	-4054.4400	8397.6137	0.0001	0.9999
HKYG	40	-4058.4766	8398.6446	0.0001	1.0000
HKYIG	41	-4058.4505	8405.6347	0.0000	1.0000
HKYI	40	-4063.1812	8408.0538	0.0000	1.0000
GTRG	44	-4050.0743	8410.0092	0.0000	1.0000
GTRI	44	-4052.2647	8414.3900	0.0000	1.0000
GTRIG	45	-4049.1023	8415.1074	0.0000	1.0000
K2P	36	-4103.4053	8460.3329	0.0000	1.0000
JC69G	36	-4104.6602	8462.8428	0.0000	1.0000
SYM	40	-4092.5826	8466.8567	0.0000	1.0000
JC69I	36	-4109.0096	8471.5415	0.0000	1.0000
JC69IG	37	-4105.8181	8472.2008	0.0000	1.0000
HKY	39	-4102.4748	8479.5988	0.0000	1.0000
F81G	39	-4103.5945	8481.8382	0.0000	1.0000
GTR	43	-4091.8384	8486.4951	0.0000	1.0000
F81I	39	-4107.9508	8490.5508	0.0000	1.0000
F81IG	40	-4104.7677	8491.2269	0.0000	1.0000
JC69	35	-4145.8897	8538.2594	0.0000	1.0000
F81	38	-4144.7878	8557.1825	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_7393.phylip = 1 - 1144;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_7393.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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_7393.phylip = 1 - 1144;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_7393.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=Fixed(Equal);
END;


[Mrbayes block for the best BIC model (K2PG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_7393.phylip = 1 - 1144;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_7393.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

