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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	40	-1071.2885	2233.6582	0.5687	0.5687
HKYI	40	-1072.4064	2235.8939	0.1859	0.7546
HKYIG	41	-1071.2062	2236.0871	0.1688	0.9235
GTRG	44	-1068.5180	2238.5977	0.0481	0.9716
GTRI	44	-1069.6865	2240.9347	0.0150	0.9865
GTRIG	45	-1068.4600	2241.1468	0.0135	1.0000
K2PG	37	-1085.8491	2255.1028	0.0000	1.0000
HKY	39	-1085.0901	2258.6854	0.0000	1.0000
K2PI	37	-1089.2696	2261.9438	0.0000	1.0000
K2PIG	38	-1088.0228	2261.9918	0.0000	1.0000
GTR	43	-1082.4304	2263.7756	0.0000	1.0000
SYMG	41	-1085.2315	2264.1376	0.0000	1.0000
SYMI	41	-1088.5707	2270.8160	0.0000	1.0000
SYMIG	42	-1087.4196	2271.1250	0.0000	1.0000
K2P	36	-1099.9355	2280.7510	0.0000	1.0000
SYM	40	-1098.9950	2289.0711	0.0000	1.0000
F81G	39	-1100.7968	2290.0986	0.0000	1.0000
F81IG	40	-1100.6994	2292.4798	0.0000	1.0000
F81I	39	-1102.0413	2292.5878	0.0000	1.0000
JC69G	36	-1112.8591	2306.5982	0.0000	1.0000
F81	38	-1114.2890	2314.5243	0.0000	1.0000
JC69IG	37	-1115.6478	2314.7002	0.0000	1.0000
JC69I	36	-1116.9278	2314.7355	0.0000	1.0000
JC69	35	-1126.5055	2331.3831	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1071.2885	2222.5771	0.4497	0.4497
HKYIG	41	-1071.2062	2224.4125	0.1796	0.6293
HKYI	40	-1072.4064	2224.8128	0.1470	0.7763
GTRG	44	-1068.5180	2225.0360	0.1315	0.9078
GTRIG	45	-1068.4600	2226.9200	0.0513	0.9591
GTRI	44	-1069.6865	2227.3731	0.0409	1.0000
K2PG	37	-1085.8491	2245.6981	0.0000	1.0000
HKY	39	-1085.0901	2248.1803	0.0000	1.0000
GTR	43	-1082.4304	2250.8609	0.0000	1.0000
K2PIG	38	-1088.0228	2252.0455	0.0000	1.0000
SYMG	41	-1085.2315	2252.4630	0.0000	1.0000
K2PI	37	-1089.2696	2252.5391	0.0000	1.0000
SYMIG	42	-1087.4196	2258.8393	0.0000	1.0000
SYMI	41	-1088.5707	2259.1415	0.0000	1.0000
K2P	36	-1099.9355	2271.8710	0.0000	1.0000
SYM	40	-1098.9950	2277.9901	0.0000	1.0000
F81G	39	-1100.7968	2279.5935	0.0000	1.0000
F81IG	40	-1100.6994	2281.3988	0.0000	1.0000
F81I	39	-1102.0413	2282.0827	0.0000	1.0000
JC69G	36	-1112.8591	2297.7182	0.0000	1.0000
F81	38	-1114.2890	2304.5780	0.0000	1.0000
JC69IG	37	-1115.6478	2305.2956	0.0000	1.0000
JC69I	36	-1116.9278	2305.8555	0.0000	1.0000
JC69	35	-1126.5055	2323.0110	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1071.2885	2375.3804	0.7197	0.7197
HKYI	40	-1072.4064	2377.6162	0.2353	0.9551
HKYIG	41	-1071.2062	2381.0359	0.0426	0.9977
K2PG	37	-1085.8491	2387.0412	0.0021	0.9998
GTRG	44	-1068.5180	2393.1197	0.0001	0.9999
K2PI	37	-1089.2696	2393.8822	0.0001	0.9999
GTRI	44	-1069.6865	2395.4567	0.0000	1.0000
HKY	39	-1085.0901	2397.1635	0.0000	1.0000
K2PIG	38	-1088.0228	2397.2087	0.0000	1.0000
GTRIG	45	-1068.4600	2398.8237	0.0000	1.0000
SYMG	41	-1085.2315	2409.0864	0.0000	1.0000
K2P	36	-1099.9355	2409.3940	0.0000	1.0000
GTR	43	-1082.4304	2415.1244	0.0000	1.0000
SYMI	41	-1088.5707	2415.7649	0.0000	1.0000
SYMIG	42	-1087.4196	2419.2827	0.0000	1.0000
F81G	39	-1100.7968	2428.5768	0.0000	1.0000
SYM	40	-1098.9950	2430.7934	0.0000	1.0000
F81I	39	-1102.0413	2431.0659	0.0000	1.0000
F81IG	40	-1100.6994	2434.2021	0.0000	1.0000
JC69G	36	-1112.8591	2435.2411	0.0000	1.0000
JC69I	36	-1116.9278	2443.3785	0.0000	1.0000
JC69IG	37	-1115.6478	2446.6386	0.0000	1.0000
F81	38	-1114.2890	2449.7412	0.0000	1.0000
JC69	35	-1126.5055	2456.7139	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/chrZ_4746.phylip = 1 - 337;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_4746.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/chrZ_4746.phylip = 1 - 337;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_4746.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 (HKYG)]

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

