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

Minimum AIC  model: HKYG
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	-1107.3172	2298.5707	0.3214	0.3214
HKYG	40	-1103.9041	2299.5225	0.1997	0.5211
K2PIG	38	-1106.8107	2300.1320	0.1472	0.6683
HKYIG	41	-1103.3680	2301.0800	0.0917	0.7600
K2PI	37	-1108.6599	2301.2562	0.0839	0.8439
SYMG	41	-1103.8228	2301.9897	0.0582	0.9021
HKYI	40	-1105.2829	2302.2800	0.0503	0.9524
SYMIG	42	-1103.4284	2303.8496	0.0229	0.9753
SYMI	41	-1105.1549	2304.6538	0.0154	0.9907
GTRG	44	-1102.1593	2306.6664	0.0056	0.9963
GTRIG	45	-1101.6812	2308.4169	0.0023	0.9986
GTRI	44	-1103.5539	2309.4557	0.0014	1.0000
JC69G	36	-1128.2985	2337.9772	0.0000	1.0000
F81G	39	-1124.9047	2338.9125	0.0000	1.0000
JC69IG	37	-1127.7149	2339.3662	0.0000	1.0000
F81IG	40	-1124.2519	2340.2181	0.0000	1.0000
JC69I	36	-1129.8374	2341.0550	0.0000	1.0000
F81I	39	-1126.5446	2342.1923	0.0000	1.0000
SYM	40	-1138.8456	2369.4054	0.0000	1.0000
K2P	36	-1144.1510	2369.6824	0.0000	1.0000
HKY	39	-1140.9952	2371.0937	0.0000	1.0000
GTR	43	-1137.5289	2374.7184	0.0000	1.0000
JC69	35	-1163.7029	2406.2479	0.0000	1.0000
F81	38	-1160.6013	2407.7132	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1103.9041	2287.8082	0.2479	0.2479
K2PG	37	-1107.3172	2288.6343	0.1640	0.4119
HKYIG	41	-1103.3680	2288.7359	0.1559	0.5678
K2PIG	38	-1106.8107	2289.6213	0.1001	0.6679
SYMG	41	-1103.8228	2289.6456	0.0989	0.7668
HKYI	40	-1105.2829	2290.5657	0.0624	0.8292
SYMIG	42	-1103.4284	2290.8568	0.0540	0.8832
K2PI	37	-1108.6599	2291.3198	0.0428	0.9260
SYMI	41	-1105.1549	2292.3097	0.0261	0.9521
GTRG	44	-1102.1593	2292.3186	0.0260	0.9781
GTRIG	45	-1101.6812	2293.3623	0.0154	0.9936
GTRI	44	-1103.5539	2295.1078	0.0064	1.0000
F81G	39	-1124.9047	2327.8093	0.0000	1.0000
F81IG	40	-1124.2519	2328.5038	0.0000	1.0000
JC69G	36	-1128.2985	2328.5970	0.0000	1.0000
JC69IG	37	-1127.7149	2329.4298	0.0000	1.0000
F81I	39	-1126.5446	2331.0891	0.0000	1.0000
JC69I	36	-1129.8374	2331.6747	0.0000	1.0000
SYM	40	-1138.8456	2357.6911	0.0000	1.0000
HKY	39	-1140.9952	2359.9905	0.0000	1.0000
K2P	36	-1144.1510	2360.3021	0.0000	1.0000
GTR	43	-1137.5289	2361.0577	0.0000	1.0000
F81	38	-1160.6013	2397.2026	0.0000	1.0000
JC69	35	-1163.7029	2397.4058	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1107.3172	2428.1777	0.7346	0.7346
K2PI	37	-1108.6599	2430.8631	0.1918	0.9264
K2PIG	38	-1106.8107	2432.9361	0.0680	0.9945
HKYG	40	-1103.9041	2438.6659	0.0039	0.9983
HKYI	40	-1105.2829	2441.4234	0.0010	0.9993
HKYIG	41	-1103.3680	2443.3650	0.0004	0.9997
SYMG	41	-1103.8228	2444.2747	0.0002	0.9999
SYMI	41	-1105.1549	2446.9388	0.0001	1.0000
SYMIG	42	-1103.4284	2449.2573	0.0000	1.0000
GTRG	44	-1102.1593	2458.2620	0.0000	1.0000
GTRI	44	-1103.5539	2461.0512	0.0000	1.0000
GTRIG	45	-1101.6812	2463.0772	0.0000	1.0000
JC69G	36	-1128.2985	2464.3688	0.0000	1.0000
JC69I	36	-1129.8374	2467.4466	0.0000	1.0000
JC69IG	37	-1127.7149	2468.9731	0.0000	1.0000
F81G	39	-1124.9047	2474.8955	0.0000	1.0000
F81I	39	-1126.5446	2478.1753	0.0000	1.0000
F81IG	40	-1124.2519	2479.3615	0.0000	1.0000
K2P	36	-1144.1510	2496.0740	0.0000	1.0000
HKY	39	-1140.9952	2507.0767	0.0000	1.0000
SYM	40	-1138.8456	2508.5488	0.0000	1.0000
GTR	43	-1137.5289	2523.2297	0.0000	1.0000
JC69	35	-1163.7029	2529.4063	0.0000	1.0000
F81	38	-1160.6013	2540.5173	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/chr2_26954.phylip = 1 - 321;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_26954.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 (K2PG)]

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

