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

Minimum AIC  model: SYMIG
Minimum AICc model: K2PG
Minimum BIC  model: K2PG


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
K2PG	37	-1104.4290	2293.0464	0.3418	0.3418
K2PIG	38	-1103.2419	2293.2620	0.3068	0.6486
SYMG	41	-1099.9476	2294.5570	0.1606	0.8092
SYMIG	42	-1098.8384	2295.0052	0.1283	0.9375
K2PI	37	-1106.9996	2298.1876	0.0261	0.9637
HKYG	40	-1104.0821	2300.1789	0.0097	0.9733
HKYIG	41	-1102.8760	2300.4138	0.0086	0.9819
GTRG	44	-1099.0795	2300.8803	0.0068	0.9887
SYMI	41	-1103.1507	2300.9632	0.0065	0.9952
GTRIG	45	-1098.3203	2302.0883	0.0037	0.9990
HKYI	40	-1106.5965	2305.2077	0.0008	0.9997
GTRI	44	-1102.3126	2307.3463	0.0003	1.0000
JC69G	36	-1130.0838	2341.7850	0.0000	1.0000
JC69IG	37	-1128.8160	2341.8205	0.0000	1.0000
K2P	36	-1133.1013	2347.8199	0.0000	1.0000
JC69I	36	-1133.2609	2348.1391	0.0000	1.0000
SYM	40	-1128.3077	2348.6301	0.0000	1.0000
F81G	39	-1129.6522	2348.6912	0.0000	1.0000
F81IG	40	-1128.4012	2348.8170	0.0000	1.0000
GTR	43	-1127.4493	2354.9134	0.0000	1.0000
HKY	39	-1132.7879	2354.9626	0.0000	1.0000
F81I	39	-1132.8134	2355.0137	0.0000	1.0000
JC69	35	-1158.5282	2396.1211	0.0000	1.0000
F81	38	-1158.0049	2402.7880	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMIG	42	-1098.8384	2281.6768	0.2872	0.2872
SYMG	41	-1099.9476	2281.8952	0.2575	0.5447
K2PIG	38	-1103.2419	2282.4838	0.1919	0.7366
K2PG	37	-1104.4290	2282.8580	0.1591	0.8957
GTRG	44	-1099.0795	2286.1591	0.0305	0.9263
GTRIG	45	-1098.3203	2286.6406	0.0240	0.9503
HKYIG	41	-1102.8760	2287.7520	0.0138	0.9640
K2PI	37	-1106.9996	2287.9992	0.0122	0.9762
HKYG	40	-1104.0821	2288.1642	0.0112	0.9874
SYMI	41	-1103.1507	2288.3015	0.0105	0.9979
GTRI	44	-1102.3126	2292.6251	0.0012	0.9991
HKYI	40	-1106.5965	2293.1931	0.0009	1.0000
JC69IG	37	-1128.8160	2331.6321	0.0000	1.0000
JC69G	36	-1130.0838	2332.1677	0.0000	1.0000
SYM	40	-1128.3077	2336.6154	0.0000	1.0000
F81IG	40	-1128.4012	2336.8023	0.0000	1.0000
F81G	39	-1129.6522	2337.3044	0.0000	1.0000
K2P	36	-1133.1013	2338.2026	0.0000	1.0000
JC69I	36	-1133.2609	2338.5218	0.0000	1.0000
GTR	43	-1127.4493	2340.8985	0.0000	1.0000
HKY	39	-1132.7879	2343.5757	0.0000	1.0000
F81I	39	-1132.8134	2343.6268	0.0000	1.0000
JC69	35	-1158.5282	2387.0563	0.0000	1.0000
F81	38	-1158.0049	2392.0098	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1104.4290	2421.5856	0.7919	0.7919
K2PIG	38	-1103.2419	2424.9608	0.1465	0.9383
K2PI	37	-1106.9996	2426.7268	0.0606	0.9989
SYMG	41	-1099.9476	2435.6204	0.0007	0.9996
HKYG	40	-1104.0821	2438.1399	0.0002	0.9998
SYMIG	42	-1098.8384	2439.1513	0.0001	0.9999
HKYIG	41	-1102.8760	2441.4772	0.0000	1.0000
SYMI	41	-1103.1507	2442.0266	0.0000	1.0000
HKYI	40	-1106.5965	2443.1688	0.0000	1.0000
GTRG	44	-1099.0795	2451.1324	0.0000	1.0000
GTRIG	45	-1098.3203	2455.3632	0.0000	1.0000
GTRI	44	-1102.3126	2457.5984	0.0000	1.0000
JC69G	36	-1130.0838	2467.1458	0.0000	1.0000
JC69IG	37	-1128.8160	2470.3596	0.0000	1.0000
K2P	36	-1133.1013	2473.1807	0.0000	1.0000
JC69I	36	-1133.2609	2473.4999	0.0000	1.0000
F81G	39	-1129.6522	2483.5307	0.0000	1.0000
SYM	40	-1128.3077	2486.5912	0.0000	1.0000
F81IG	40	-1128.4012	2486.7781	0.0000	1.0000
HKY	39	-1132.7879	2489.8020	0.0000	1.0000
F81I	39	-1132.8134	2489.8531	0.0000	1.0000
GTR	43	-1127.4493	2502.1224	0.0000	1.0000
JC69	35	-1158.5282	2518.2851	0.0000	1.0000
F81	38	-1158.0049	2534.4867	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_18714.phylip = 1 - 314;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_18714.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Fixed(Equal) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

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


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

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


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

