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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKY


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
HKYG	40	-1023.2704	2136.0756	0.3541	0.3541
HKYI	40	-1023.4084	2136.3517	0.3084	0.6625
HKYIG	41	-1023.2827	2138.6062	0.0999	0.7624
HKY	39	-1026.0676	2139.1787	0.0750	0.8375
GTRG	44	-1019.7812	2139.2094	0.0739	0.9114
GTRI	44	-1020.0815	2139.8100	0.0547	0.9661
GTRIG	45	-1019.7875	2141.7874	0.0204	0.9865
GTR	43	-1022.7564	2142.6097	0.0135	1.0000
F81G	39	-1034.9413	2156.9261	0.0000	1.0000
F81I	39	-1035.0626	2157.1686	0.0000	1.0000
F81IG	40	-1034.9510	2159.4368	0.0000	1.0000
F81	38	-1037.6905	2159.9475	0.0000	1.0000
K2PG	37	-1043.1380	2168.3797	0.0000	1.0000
K2PI	37	-1043.2089	2168.5215	0.0000	1.0000
K2PIG	38	-1043.1448	2170.8560	0.0000	1.0000
K2P	36	-1046.1666	2171.9884	0.0000	1.0000
SYMG	41	-1042.5167	2177.0743	0.0000	1.0000
SYMI	41	-1042.6330	2177.3068	0.0000	1.0000
SYMIG	42	-1042.5286	2179.6186	0.0000	1.0000
SYM	40	-1045.4473	2180.4295	0.0000	1.0000
JC69G	36	-1052.5754	2184.8060	0.0000	1.0000
JC69I	36	-1052.6729	2185.0010	0.0000	1.0000
JC69IG	37	-1052.5829	2187.2695	0.0000	1.0000
JC69	35	-1055.5088	2188.2383	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1023.2704	2126.5408	0.2669	0.2669
HKYI	40	-1023.4084	2126.8169	0.2325	0.4993
GTRG	44	-1019.7812	2127.5624	0.1601	0.6594
GTRI	44	-1020.0815	2128.1629	0.1186	0.7780
HKYIG	41	-1023.2827	2128.5654	0.0970	0.8750
GTRIG	45	-1019.7875	2129.5750	0.0585	0.9335
HKY	39	-1026.0676	2130.1352	0.0442	0.9778
GTR	43	-1022.7564	2131.5129	0.0222	1.0000
F81G	39	-1034.9413	2147.8826	0.0000	1.0000
F81I	39	-1035.0626	2148.1252	0.0000	1.0000
F81IG	40	-1034.9510	2149.9019	0.0000	1.0000
F81	38	-1037.6905	2151.3810	0.0000	1.0000
K2PG	37	-1043.1380	2160.2759	0.0000	1.0000
K2PI	37	-1043.2089	2160.4177	0.0000	1.0000
K2PIG	38	-1043.1448	2162.2896	0.0000	1.0000
K2P	36	-1046.1666	2164.3333	0.0000	1.0000
SYMG	41	-1042.5167	2167.0335	0.0000	1.0000
SYMI	41	-1042.6330	2167.2659	0.0000	1.0000
SYMIG	42	-1042.5286	2169.0572	0.0000	1.0000
SYM	40	-1045.4473	2170.8946	0.0000	1.0000
JC69G	36	-1052.5754	2177.1508	0.0000	1.0000
JC69I	36	-1052.6729	2177.3458	0.0000	1.0000
JC69IG	37	-1052.5829	2179.1657	0.0000	1.0000
JC69	35	-1055.5088	2181.0176	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	39	-1026.0676	2284.3117	0.3836	0.3836
HKYG	40	-1023.2704	2284.6705	0.3206	0.7042
HKYI	40	-1023.4084	2284.9466	0.2793	0.9835
HKYIG	41	-1023.2827	2290.6483	0.0161	0.9996
GTR	43	-1022.7564	2301.5023	0.0001	0.9997
GTRG	44	-1019.7812	2301.5051	0.0001	0.9998
F81	38	-1037.6905	2301.6042	0.0001	0.9998
F81G	39	-1034.9413	2302.0591	0.0001	0.9999
GTRI	44	-1020.0815	2302.1056	0.0001	0.9999
F81I	39	-1035.0626	2302.3017	0.0000	1.0000
K2PG	37	-1043.1380	2306.5459	0.0000	1.0000
K2P	36	-1046.1666	2306.6500	0.0000	1.0000
K2PI	37	-1043.2089	2306.6877	0.0000	1.0000
GTRIG	45	-1019.7875	2307.4709	0.0000	1.0000
F81IG	40	-1034.9510	2308.0317	0.0000	1.0000
K2PIG	38	-1043.1448	2312.5128	0.0000	1.0000
JC69	35	-1055.5088	2319.3812	0.0000	1.0000
JC69G	36	-1052.5754	2319.4676	0.0000	1.0000
JC69I	36	-1052.6729	2319.6626	0.0000	1.0000
JC69IG	37	-1052.5829	2325.4357	0.0000	1.0000
SYM	40	-1045.4473	2329.0244	0.0000	1.0000
SYMG	41	-1042.5167	2329.1165	0.0000	1.0000
SYMI	41	-1042.6330	2329.3489	0.0000	1.0000
SYMIG	42	-1042.5286	2335.0934	0.0000	1.0000



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

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

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


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

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


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_14921.phylip = 1 - 385;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_14921.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

