
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:37 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/chr5_11794.phylip" (19 taxa, 453 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	-1106.9760	2301.9132	0.3759	0.3759
HKYG	40	-1107.2043	2302.3697	0.2992	0.6752
GTRI	44	-1103.3076	2304.3211	0.1128	0.7880
GTRG	44	-1103.5134	2304.7327	0.0918	0.8798
HKYIG	41	-1107.2360	2304.8517	0.0865	0.9663
GTRIG	45	-1103.5391	2307.2501	0.0261	0.9923
HKY	39	-1112.4119	2310.3782	0.0055	0.9978
GTR	43	-1108.4740	2312.1997	0.0022	1.0000
F81I	39	-1121.5148	2328.5840	0.0000	1.0000
F81G	39	-1121.8598	2329.2741	0.0000	1.0000
F81IG	40	-1121.9208	2331.8028	0.0000	1.0000
F81	38	-1127.4612	2338.0818	0.0000	1.0000
K2PI	37	-1156.9664	2394.7087	0.0000	1.0000
K2PG	37	-1157.1588	2395.0936	0.0000	1.0000
K2PIG	38	-1157.2000	2397.5594	0.0000	1.0000
SYMI	41	-1155.3571	2401.0938	0.0000	1.0000
SYMG	41	-1155.4700	2401.3195	0.0000	1.0000
SYMIG	42	-1155.4976	2403.8050	0.0000	1.0000
K2P	36	-1162.9887	2404.3813	0.0000	1.0000
SYM	40	-1161.0492	2410.0595	0.0000	1.0000
JC69I	36	-1166.9949	2412.3936	0.0000	1.0000
JC69G	36	-1167.2828	2412.9695	0.0000	1.0000
JC69IG	37	-1167.3427	2415.4612	0.0000	1.0000
JC69	35	-1173.0940	2422.2312	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1106.9760	2293.9520	0.2767	0.2767
HKYG	40	-1107.2043	2294.4085	0.2202	0.4969
GTRI	44	-1103.3076	2294.6153	0.1986	0.6955
GTRG	44	-1103.5134	2295.0269	0.1617	0.8572
HKYIG	41	-1107.2360	2296.4721	0.0785	0.9357
GTRIG	45	-1103.5391	2297.0781	0.0580	0.9936
HKY	39	-1112.4119	2302.8237	0.0033	0.9969
GTR	43	-1108.4740	2302.9479	0.0031	1.0000
F81I	39	-1121.5148	2321.0295	0.0000	1.0000
F81G	39	-1121.8598	2321.7196	0.0000	1.0000
F81IG	40	-1121.9208	2323.8417	0.0000	1.0000
F81	38	-1127.4612	2330.9224	0.0000	1.0000
K2PI	37	-1156.9664	2387.9328	0.0000	1.0000
K2PG	37	-1157.1588	2388.3177	0.0000	1.0000
K2PIG	38	-1157.2000	2390.4000	0.0000	1.0000
SYMI	41	-1155.3571	2392.7142	0.0000	1.0000
SYMG	41	-1155.4700	2392.9399	0.0000	1.0000
SYMIG	42	-1155.4976	2394.9952	0.0000	1.0000
K2P	36	-1162.9887	2397.9774	0.0000	1.0000
SYM	40	-1161.0492	2402.0984	0.0000	1.0000
JC69I	36	-1166.9949	2405.9897	0.0000	1.0000
JC69G	36	-1167.2828	2406.5657	0.0000	1.0000
JC69IG	37	-1167.3427	2408.6853	0.0000	1.0000
JC69	35	-1173.0940	2416.1880	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1106.9760	2458.5877	0.5194	0.5194
HKYG	40	-1107.2043	2459.0442	0.4134	0.9328
HKY	39	-1112.4119	2463.3435	0.0482	0.9810
HKYIG	41	-1107.2360	2465.2237	0.0188	0.9998
GTRI	44	-1103.3076	2475.7145	0.0001	0.9999
GTRG	44	-1103.5134	2476.1261	0.0001	1.0000
GTR	43	-1108.4740	2479.9313	0.0000	1.0000
F81I	39	-1121.5148	2481.5493	0.0000	1.0000
F81G	39	-1121.8598	2482.2394	0.0000	1.0000
GTRIG	45	-1103.5391	2482.2933	0.0000	1.0000
F81	38	-1127.4612	2487.3263	0.0000	1.0000
F81IG	40	-1121.9208	2488.4774	0.0000	1.0000
K2PI	37	-1156.9664	2540.2208	0.0000	1.0000
K2PG	37	-1157.1588	2540.6057	0.0000	1.0000
K2P	36	-1162.9887	2546.1496	0.0000	1.0000
K2PIG	38	-1157.2000	2546.8039	0.0000	1.0000
JC69I	36	-1166.9949	2554.1618	0.0000	1.0000
JC69G	36	-1167.2828	2554.7378	0.0000	1.0000
JC69	35	-1173.0940	2560.2442	0.0000	1.0000
JC69IG	37	-1167.3427	2560.9733	0.0000	1.0000
SYMI	41	-1155.3571	2561.4658	0.0000	1.0000
SYMG	41	-1155.4700	2561.6915	0.0000	1.0000
SYM	40	-1161.0492	2566.7340	0.0000	1.0000
SYMIG	42	-1155.4976	2567.8627	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/chr5_11794.phylip = 1 - 453;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11794.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/chr5_11794.phylip = 1 - 453;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11794.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/chr5_11794.phylip = 1 - 453;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11794.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

