
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:50 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/chr7_6308.phylip" (19 taxa, 338 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	-1196.6702	2484.3842	0.4607	0.4607
HKYI	40	-1196.8865	2484.8169	0.3710	0.8317
HKYIG	41	-1196.5977	2486.8305	0.1356	0.9673
GTRG	44	-1194.9462	2491.4077	0.0137	0.9810
GTRI	44	-1194.9942	2491.5037	0.0131	0.9941
GTRIG	45	-1194.8121	2493.8024	0.0042	0.9983
K2PG	37	-1206.8253	2497.0239	0.0008	0.9991
K2PI	37	-1207.2362	2497.8457	0.0005	0.9996
K2PIG	38	-1206.7913	2499.4956	0.0002	0.9999
HKY	39	-1206.4434	2501.3566	0.0001	1.0000
SYMG	41	-1206.0964	2505.8280	0.0000	1.0000
SYMI	41	-1206.5236	2506.6824	0.0000	1.0000
SYMIG	42	-1206.0661	2508.3763	0.0000	1.0000
GTR	43	-1204.8235	2508.5176	0.0000	1.0000
K2P	36	-1215.7407	2512.3319	0.0000	1.0000
SYM	40	-1215.1390	2521.3217	0.0000	1.0000
F81G	39	-1217.8266	2524.1231	0.0000	1.0000
F81I	39	-1218.1578	2524.7853	0.0000	1.0000
F81IG	40	-1217.7887	2526.6211	0.0000	1.0000
JC69G	36	-1227.1174	2535.0852	0.0000	1.0000
JC69I	36	-1227.5635	2535.9776	0.0000	1.0000
JC69IG	37	-1227.0883	2537.5500	0.0000	1.0000
F81	38	-1226.6080	2539.1290	0.0000	1.0000
JC69	35	-1235.8210	2549.9863	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1196.6702	2473.3404	0.4087	0.4087
HKYI	40	-1196.8865	2473.7731	0.3292	0.7379
HKYIG	41	-1196.5977	2475.1954	0.1617	0.8996
GTRG	44	-1194.9462	2477.8923	0.0420	0.9416
GTRI	44	-1194.9942	2477.9883	0.0400	0.9816
GTRIG	45	-1194.8121	2479.6243	0.0177	0.9993
K2PG	37	-1206.8253	2487.6506	0.0003	0.9996
K2PI	37	-1207.2362	2488.4724	0.0002	0.9998
K2PIG	38	-1206.7913	2489.5825	0.0001	0.9999
HKY	39	-1206.4434	2490.8868	0.0001	1.0000
SYMG	41	-1206.0964	2494.1929	0.0000	1.0000
SYMI	41	-1206.5236	2495.0472	0.0000	1.0000
GTR	43	-1204.8235	2495.6469	0.0000	1.0000
SYMIG	42	-1206.0661	2496.1322	0.0000	1.0000
K2P	36	-1215.7407	2503.4814	0.0000	1.0000
SYM	40	-1215.1390	2510.2780	0.0000	1.0000
F81G	39	-1217.8266	2513.6533	0.0000	1.0000
F81I	39	-1218.1578	2514.3155	0.0000	1.0000
F81IG	40	-1217.7887	2515.5774	0.0000	1.0000
JC69G	36	-1227.1174	2526.2347	0.0000	1.0000
JC69I	36	-1227.5635	2527.1271	0.0000	1.0000
JC69IG	37	-1227.0883	2528.1767	0.0000	1.0000
F81	38	-1226.6080	2529.2160	0.0000	1.0000
JC69	35	-1235.8210	2541.6419	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1196.6702	2626.2622	0.4384	0.4384
HKYI	40	-1196.8865	2626.6949	0.3531	0.7915
K2PG	37	-1206.8253	2629.1033	0.1059	0.8974
K2PI	37	-1207.2362	2629.9251	0.0702	0.9676
HKYIG	41	-1196.5977	2631.9402	0.0256	0.9933
K2PIG	38	-1206.7913	2634.8583	0.0060	0.9992
HKY	39	-1206.4434	2639.9856	0.0005	0.9997
K2P	36	-1215.7407	2641.1111	0.0003	1.0000
GTRG	44	-1194.9462	2646.1063	0.0000	1.0000
GTRI	44	-1194.9942	2646.2024	0.0000	1.0000
SYMG	41	-1206.0964	2650.9378	0.0000	1.0000
GTRIG	45	-1194.8121	2651.6613	0.0000	1.0000
SYMI	41	-1206.5236	2651.7921	0.0000	1.0000
SYMIG	42	-1206.0661	2656.7001	0.0000	1.0000
GTR	43	-1204.8235	2660.0379	0.0000	1.0000
F81G	39	-1217.8266	2662.7521	0.0000	1.0000
SYM	40	-1215.1390	2663.1998	0.0000	1.0000
F81I	39	-1218.1578	2663.4143	0.0000	1.0000
JC69G	36	-1227.1174	2663.8644	0.0000	1.0000
JC69I	36	-1227.5635	2664.7567	0.0000	1.0000
F81IG	40	-1217.7887	2668.4992	0.0000	1.0000
JC69IG	37	-1227.0883	2669.6294	0.0000	1.0000
F81	38	-1226.6080	2674.4917	0.0000	1.0000
JC69	35	-1235.8210	2675.4485	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/chr7_6308.phylip = 1 - 338;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_6308.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/chr7_6308.phylip = 1 - 338;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_6308.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/chr7_6308.phylip = 1 - 338;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_6308.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

