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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	44	-1076.9733	2252.4229	0.2307	0.2307
GTRI	44	-1077.0987	2252.6737	0.2036	0.4343
HKYG	40	-1082.4383	2253.4630	0.1372	0.5715
HKYI	40	-1082.6126	2253.8116	0.1152	0.6867
GTRIG	45	-1076.8417	2254.6649	0.0752	0.7619
K2PG	37	-1086.8981	2255.1002	0.0605	0.8224
HKYIG	41	-1082.1154	2255.2701	0.0556	0.8780
K2PI	37	-1087.0733	2255.4506	0.0508	0.9288
K2PIG	38	-1086.4922	2256.7032	0.0271	0.9559
SYMG	41	-1083.1508	2257.3410	0.0197	0.9756
SYMI	41	-1083.3125	2257.6644	0.0168	0.9924
SYMIG	42	-1082.8761	2259.2574	0.0076	1.0000
GTR	43	-1095.5937	2287.1715	0.0000	1.0000
HKY	39	-1101.1003	2288.3469	0.0000	1.0000
K2P	36	-1106.0496	2291.0007	0.0000	1.0000
SYM	40	-1101.7531	2292.0927	0.0000	1.0000
F81G	39	-1110.1643	2306.4748	0.0000	1.0000
F81I	39	-1110.3326	2306.8114	0.0000	1.0000
JC69G	36	-1114.6950	2308.2915	0.0000	1.0000
F81IG	40	-1109.8692	2308.3247	0.0000	1.0000
JC69I	36	-1114.8777	2308.6569	0.0000	1.0000
JC69IG	37	-1114.3115	2309.9269	0.0000	1.0000
F81	38	-1128.6647	2341.0481	0.0000	1.0000
JC69	35	-1133.4440	2343.3996	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1076.9733	2241.9467	0.3257	0.3257
GTRI	44	-1077.0987	2242.1975	0.2873	0.6131
GTRIG	45	-1076.8417	2243.6835	0.1367	0.7497
HKYG	40	-1082.4383	2244.8766	0.0753	0.8250
HKYI	40	-1082.6126	2245.2252	0.0632	0.8882
HKYIG	41	-1082.1154	2246.2307	0.0382	0.9265
K2PG	37	-1086.8981	2247.7963	0.0175	0.9440
K2PI	37	-1087.0733	2248.1467	0.0147	0.9586
SYMG	41	-1083.1508	2248.3016	0.0136	0.9722
SYMI	41	-1083.3125	2248.6250	0.0116	0.9838
K2PIG	38	-1086.4922	2248.9844	0.0097	0.9934
SYMIG	42	-1082.8761	2249.7522	0.0066	1.0000
GTR	43	-1095.5937	2277.1873	0.0000	1.0000
HKY	39	-1101.1003	2280.2007	0.0000	1.0000
SYM	40	-1101.7531	2283.5063	0.0000	1.0000
K2P	36	-1106.0496	2284.0992	0.0000	1.0000
F81G	39	-1110.1643	2298.3285	0.0000	1.0000
F81I	39	-1110.3326	2298.6651	0.0000	1.0000
F81IG	40	-1109.8692	2299.7384	0.0000	1.0000
JC69G	36	-1114.6950	2301.3899	0.0000	1.0000
JC69I	36	-1114.8777	2301.7553	0.0000	1.0000
JC69IG	37	-1114.3115	2302.6230	0.0000	1.0000
F81	38	-1128.6647	2333.3293	0.0000	1.0000
JC69	35	-1133.4440	2336.8880	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1086.8981	2397.5491	0.5177	0.5177
K2PI	37	-1087.0733	2397.8995	0.4345	0.9522
K2PIG	38	-1086.4922	2402.7846	0.0378	0.9899
HKYG	40	-1082.4383	2406.7715	0.0051	0.9951
HKYI	40	-1082.6126	2407.1201	0.0043	0.9994
HKYIG	41	-1082.1154	2412.1730	0.0003	0.9998
SYMG	41	-1083.1508	2414.2439	0.0001	0.9999
SYMI	41	-1083.3125	2414.5673	0.0001	1.0000
SYMIG	42	-1082.8761	2419.7418	0.0000	1.0000
GTRG	44	-1076.9733	2420.0311	0.0000	1.0000
GTRI	44	-1077.0987	2420.2819	0.0000	1.0000
GTRIG	45	-1076.8417	2425.8152	0.0000	1.0000
K2P	36	-1106.0496	2429.8046	0.0000	1.0000
HKY	39	-1101.1003	2438.0482	0.0000	1.0000
SYM	40	-1101.7531	2445.4012	0.0000	1.0000
JC69G	36	-1114.6950	2447.0953	0.0000	1.0000
JC69I	36	-1114.8777	2447.4607	0.0000	1.0000
GTR	43	-1095.5937	2451.2243	0.0000	1.0000
JC69IG	37	-1114.3115	2452.3758	0.0000	1.0000
F81G	39	-1110.1643	2456.1761	0.0000	1.0000
F81I	39	-1110.3326	2456.5127	0.0000	1.0000
F81IG	40	-1109.8692	2461.6332	0.0000	1.0000
JC69	35	-1133.4440	2478.5460	0.0000	1.0000
F81	38	-1128.6647	2487.1295	0.0000	1.0000



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

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

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


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr14_715.phylip = 1 - 423;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr14_715.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,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/chr14_715.phylip = 1 - 423;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr14_715.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

