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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	44	-2882.5605	5858.3452	0.4423	0.4423
HKYG	40	-2887.5798	5859.4641	0.2528	0.6951
GTRIG	45	-2882.5001	5860.4692	0.1529	0.8480
HKYIG	41	-2887.5337	5861.5930	0.0872	0.9352
GTRI	44	-2885.0612	5863.3467	0.0363	0.9715
HKYI	40	-2889.8937	5864.0918	0.0250	0.9965
K2PG	37	-2895.7248	5869.1255	0.0020	0.9985
K2PIG	38	-2895.6724	5871.2243	0.0007	0.9992
SYMG	41	-2892.8340	5872.1935	0.0004	0.9996
K2PI	37	-2898.1356	5873.9470	0.0002	0.9998
SYMIG	42	-2892.7819	5874.3164	0.0002	1.0000
SYMI	41	-2895.3480	5877.2216	0.0000	1.0000
GTR	43	-2918.2231	5927.4317	0.0000	1.0000
HKY	39	-2923.5935	5929.2761	0.0000	1.0000
K2P	36	-2931.6074	5938.6926	0.0000	1.0000
SYM	40	-2928.0903	5940.4851	0.0000	1.0000
F81G	39	-2968.1797	6018.4486	0.0000	1.0000
F81IG	40	-2968.1734	6020.6514	0.0000	1.0000
F81I	39	-2970.1584	6022.4058	0.0000	1.0000
JC69G	36	-2973.9487	6023.3751	0.0000	1.0000
JC69IG	37	-2973.9411	6025.5581	0.0000	1.0000
JC69I	36	-2975.9894	6027.4567	0.0000	1.0000
F81	38	-3002.4619	6084.8034	0.0000	1.0000
JC69	35	-3008.3239	6089.9334	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2882.5605	5853.1209	0.4969	0.4969
GTRIG	45	-2882.5001	5855.0002	0.1942	0.6911
HKYG	40	-2887.5798	5855.1597	0.1793	0.8704
HKYIG	41	-2887.5337	5857.0674	0.0691	0.9395
GTRI	44	-2885.0612	5858.1224	0.0408	0.9802
HKYI	40	-2889.8937	5859.7874	0.0177	0.9979
K2PG	37	-2895.7248	5865.4497	0.0010	0.9990
K2PIG	38	-2895.6724	5867.3448	0.0004	0.9994
SYMG	41	-2892.8340	5867.6679	0.0003	0.9997
SYMIG	42	-2892.7819	5869.5638	0.0001	0.9999
K2PI	37	-2898.1356	5870.2711	0.0001	1.0000
SYMI	41	-2895.3480	5872.6959	0.0000	1.0000
GTR	43	-2918.2231	5922.4462	0.0000	1.0000
HKY	39	-2923.5935	5925.1870	0.0000	1.0000
K2P	36	-2931.6074	5935.2148	0.0000	1.0000
SYM	40	-2928.0903	5936.1806	0.0000	1.0000
F81G	39	-2968.1797	6014.3595	0.0000	1.0000
F81IG	40	-2968.1734	6016.3469	0.0000	1.0000
F81I	39	-2970.1584	6018.3167	0.0000	1.0000
JC69G	36	-2973.9487	6019.8973	0.0000	1.0000
JC69IG	37	-2973.9411	6021.8822	0.0000	1.0000
JC69I	36	-2975.9894	6023.9789	0.0000	1.0000
F81	38	-3002.4619	6080.9238	0.0000	1.0000
JC69	35	-3008.3239	6086.6478	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-2895.7248	6038.9188	0.7698	0.7698
HKYG	40	-2887.5798	6042.6939	0.1166	0.8864
K2PI	37	-2898.1356	6043.7403	0.0691	0.9555
K2PIG	38	-2895.6724	6045.5022	0.0286	0.9841
HKYI	40	-2889.8937	6047.3216	0.0115	0.9956
HKYIG	41	-2887.5337	6049.2899	0.0043	0.9999
GTRG	44	-2882.5605	6059.4085	0.0000	1.0000
SYMG	41	-2892.8340	6059.8904	0.0000	1.0000
GTRI	44	-2885.0612	6064.4100	0.0000	1.0000
SYMI	41	-2895.3480	6064.9185	0.0000	1.0000
GTRIG	45	-2882.5001	6065.9762	0.0000	1.0000
SYMIG	42	-2892.7819	6066.4747	0.0000	1.0000
K2P	36	-2931.6074	6103.9956	0.0000	1.0000
HKY	39	-2923.5935	6108.0328	0.0000	1.0000
SYM	40	-2928.0903	6123.7148	0.0000	1.0000
GTR	43	-2918.2231	6124.0454	0.0000	1.0000
JC69G	36	-2973.9487	6188.6781	0.0000	1.0000
JC69I	36	-2975.9894	6192.7596	0.0000	1.0000
JC69IG	37	-2973.9411	6195.3514	0.0000	1.0000
F81G	39	-2968.1797	6197.2053	0.0000	1.0000
F81I	39	-2970.1584	6201.1625	0.0000	1.0000
F81IG	40	-2968.1734	6203.8811	0.0000	1.0000
JC69	35	-3008.3239	6250.7403	0.0000	1.0000
F81	38	-3002.4619	6259.0813	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/chrZ_7928.phylip = 1 - 803;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_7928.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/chrZ_7928.phylip = 1 - 803;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_7928.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 (K2PG)]

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

