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

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


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	-1049.5822	2189.2567	0.3236	0.3236
HKYG	40	-1049.9823	2190.0569	0.2169	0.5405
HKY	39	-1051.3960	2190.3625	0.1861	0.7266
HKYIG	41	-1049.5949	2191.8193	0.0898	0.8165
GTRI	44	-1046.0747	2192.4858	0.0644	0.8808
GTR	43	-1047.6792	2193.1099	0.0471	0.9280
GTRG	44	-1046.4066	2193.1496	0.0462	0.9742
GTRIG	45	-1046.0852	2195.1080	0.0174	0.9915
K2PI	37	-1057.8450	2198.2631	0.0036	0.9951
K2PG	37	-1058.3095	2199.1921	0.0023	0.9974
K2P	36	-1060.0468	2200.1909	0.0014	0.9987
K2PIG	38	-1057.8596	2200.7835	0.0010	0.9997
SYMI	41	-1056.3576	2205.3449	0.0001	0.9998
SYMG	41	-1056.7876	2206.2048	0.0001	0.9999
SYM	40	-1058.2803	2206.6529	0.0001	1.0000
SYMIG	42	-1056.3714	2207.9254	0.0000	1.0000
F81I	39	-1065.6870	2218.9447	0.0000	1.0000
F81G	39	-1066.1614	2219.8934	0.0000	1.0000
F81	38	-1067.5629	2220.1900	0.0000	1.0000
F81IG	40	-1065.7023	2221.4969	0.0000	1.0000
JC69I	36	-1073.0324	2226.1621	0.0000	1.0000
JC69G	36	-1073.5375	2227.1723	0.0000	1.0000
JC69	35	-1075.1275	2227.8913	0.0000	1.0000
JC69IG	37	-1073.0487	2228.6706	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1049.5822	2179.1644	0.2442	0.2442
HKYG	40	-1049.9823	2179.9646	0.1637	0.4078
GTRI	44	-1046.0747	2180.1494	0.1492	0.5570
HKY	39	-1051.3960	2180.7920	0.1082	0.6652
GTRG	44	-1046.4066	2180.8132	0.1071	0.7723
HKYIG	41	-1049.5949	2181.1897	0.0887	0.8610
GTR	43	-1047.6792	2181.3583	0.0815	0.9425
GTRIG	45	-1046.0852	2182.1705	0.0543	0.9969
K2PI	37	-1057.8450	2189.6899	0.0013	0.9981
K2PG	37	-1058.3095	2190.6190	0.0008	0.9989
K2PIG	38	-1057.8596	2191.7193	0.0005	0.9994
K2P	36	-1060.0468	2192.0937	0.0004	0.9998
SYMI	41	-1056.3576	2194.7153	0.0001	0.9999
SYMG	41	-1056.7876	2195.5751	0.0001	0.9999
SYM	40	-1058.2803	2196.5606	0.0000	1.0000
SYMIG	42	-1056.3714	2196.7428	0.0000	1.0000
F81I	39	-1065.6870	2209.3741	0.0000	1.0000
F81G	39	-1066.1614	2210.3229	0.0000	1.0000
F81	38	-1067.5629	2211.1258	0.0000	1.0000
F81IG	40	-1065.7023	2211.4046	0.0000	1.0000
JC69I	36	-1073.0324	2218.0648	0.0000	1.0000
JC69G	36	-1073.5375	2219.0751	0.0000	1.0000
JC69IG	37	-1073.0487	2220.0975	0.0000	1.0000
JC69	35	-1075.1275	2220.2550	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	36	-1060.0468	2332.5885	0.3267	0.3267
HKY	39	-1051.3960	2332.9947	0.2666	0.5933
K2PI	37	-1057.8450	2334.0874	0.1544	0.7477
K2PG	37	-1058.3095	2335.0164	0.0970	0.8447
HKYI	40	-1049.5822	2335.2697	0.0855	0.9302
HKYG	40	-1049.9823	2336.0700	0.0573	0.9875
K2PIG	38	-1057.8596	2340.0193	0.0080	0.9954
HKYIG	41	-1049.5949	2341.1977	0.0044	0.9999
GTR	43	-1047.6792	2349.1716	0.0001	0.9999
GTRI	44	-1046.0747	2351.8652	0.0000	1.0000
GTRG	44	-1046.4066	2352.5290	0.0000	1.0000
SYM	40	-1058.2803	2352.6659	0.0000	1.0000
SYMI	41	-1056.3576	2354.7233	0.0000	1.0000
SYMG	41	-1056.7876	2355.5831	0.0000	1.0000
JC69	35	-1075.1275	2356.8471	0.0000	1.0000
GTRIG	45	-1046.0852	2357.7890	0.0000	1.0000
JC69I	36	-1073.0324	2358.5596	0.0000	1.0000
F81	38	-1067.5629	2359.4258	0.0000	1.0000
JC69G	36	-1073.5375	2359.5699	0.0000	1.0000
SYMIG	42	-1056.3714	2360.6534	0.0000	1.0000
F81I	39	-1065.6870	2361.5768	0.0000	1.0000
F81G	39	-1066.1614	2362.5256	0.0000	1.0000
JC69IG	37	-1073.0487	2364.4949	0.0000	1.0000
F81IG	40	-1065.7023	2367.5100	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/chr2_24763.phylip = 1 - 366;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_24763.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/chr2_24763.phylip = 1 - 366;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_24763.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 (K2P)]

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


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

