
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:58 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/chr6_5980.phylip" (19 taxa, 336 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	40	-1005.6652	2102.4490	0.3824	0.3824
HKYG	40	-1005.9135	2102.9456	0.2983	0.6807
HKYIG	41	-1005.6004	2104.9151	0.1114	0.7922
GTRI	44	-1001.7002	2105.0086	0.1063	0.8985
GTRG	44	-1002.1696	2105.9474	0.0665	0.9650
GTRIG	45	-1001.6926	2107.6610	0.0282	0.9933
K2PI	37	-1014.4220	2112.2801	0.0028	0.9961
K2PG	37	-1014.4487	2112.3337	0.0027	0.9988
K2PIG	38	-1014.2487	2114.4772	0.0009	0.9997
SYMI	41	-1012.3425	2118.3992	0.0001	0.9999
SYMG	41	-1012.5659	2118.8460	0.0001	1.0000
SYMIG	42	-1012.2889	2120.9054	0.0000	1.0000
HKY	39	-1020.7075	2129.9555	0.0000	1.0000
GTR	43	-1016.5793	2132.1174	0.0000	1.0000
K2P	36	-1029.1907	2139.2910	0.0000	1.0000
F81I	39	-1028.9031	2146.3468	0.0000	1.0000
SYM	40	-1027.7594	2146.6374	0.0000	1.0000
F81G	39	-1029.0692	2146.6790	0.0000	1.0000
F81IG	40	-1028.8488	2148.8163	0.0000	1.0000
JC69I	36	-1035.9479	2152.8054	0.0000	1.0000
JC69G	36	-1036.0207	2152.9511	0.0000	1.0000
JC69IG	37	-1036.0315	2155.4992	0.0000	1.0000
F81	38	-1043.3169	2172.6137	0.0000	1.0000
JC69	35	-1050.3650	2179.1300	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1005.6652	2091.3304	0.2434	0.2434
GTRI	44	-1001.7002	2091.4003	0.2350	0.4784
HKYG	40	-1005.9135	2091.8270	0.1899	0.6683
GTRG	44	-1002.1696	2092.3392	0.1470	0.8152
HKYIG	41	-1005.6004	2093.2008	0.0955	0.9108
GTRIG	45	-1001.6926	2093.3851	0.0871	0.9979
K2PI	37	-1014.4220	2102.8439	0.0008	0.9987
K2PG	37	-1014.4487	2102.8974	0.0007	0.9994
K2PIG	38	-1014.2487	2104.4974	0.0003	0.9998
SYMI	41	-1012.3425	2106.6849	0.0001	0.9999
SYMG	41	-1012.5659	2107.1317	0.0001	1.0000
SYMIG	42	-1012.2889	2108.5778	0.0000	1.0000
GTR	43	-1016.5793	2119.1585	0.0000	1.0000
HKY	39	-1020.7075	2119.4149	0.0000	1.0000
K2P	36	-1029.1907	2130.3813	0.0000	1.0000
SYM	40	-1027.7594	2135.5188	0.0000	1.0000
F81I	39	-1028.9031	2135.8062	0.0000	1.0000
F81G	39	-1029.0692	2136.1385	0.0000	1.0000
F81IG	40	-1028.8488	2137.6976	0.0000	1.0000
JC69I	36	-1035.9479	2143.8957	0.0000	1.0000
JC69G	36	-1036.0207	2144.0414	0.0000	1.0000
JC69IG	37	-1036.0315	2146.0630	0.0000	1.0000
F81	38	-1043.3169	2162.6339	0.0000	1.0000
JC69	35	-1050.3650	2170.7300	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1005.6652	2244.0148	0.2621	0.2621
K2PI	37	-1014.4220	2244.0770	0.2541	0.5162
K2PG	37	-1014.4487	2244.1305	0.2474	0.7635
HKYG	40	-1005.9135	2244.5114	0.2045	0.9680
K2PIG	38	-1014.2487	2249.5476	0.0165	0.9845
HKYIG	41	-1005.6004	2249.7024	0.0153	0.9998
GTRI	44	-1001.7002	2259.3532	0.0001	0.9999
GTRG	44	-1002.1696	2260.2921	0.0001	1.0000
SYMI	41	-1012.3425	2263.1865	0.0000	1.0000
SYMG	41	-1012.5659	2263.6333	0.0000	1.0000
GTRIG	45	-1001.6926	2265.1551	0.0000	1.0000
K2P	36	-1029.1907	2267.7973	0.0000	1.0000
HKY	39	-1020.7075	2268.2823	0.0000	1.0000
SYMIG	42	-1012.2889	2268.8965	0.0000	1.0000
JC69I	36	-1035.9479	2281.3117	0.0000	1.0000
JC69G	36	-1036.0207	2281.4574	0.0000	1.0000
GTR	43	-1016.5793	2283.2943	0.0000	1.0000
F81I	39	-1028.9031	2284.6736	0.0000	1.0000
F81G	39	-1029.0692	2285.0058	0.0000	1.0000
JC69IG	37	-1036.0315	2287.2961	0.0000	1.0000
SYM	40	-1027.7594	2288.2032	0.0000	1.0000
F81IG	40	-1028.8488	2290.3821	0.0000	1.0000
JC69	35	-1050.3650	2304.3289	0.0000	1.0000
F81	38	-1043.3169	2307.6841	0.0000	1.0000



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

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

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


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

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


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

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


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

