
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:47 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_10316.phylip" (19 taxa, 670 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	-2146.5319	4387.3998	0.7384	0.7384
GTRIG	45	-2147.0733	4390.7812	0.1362	0.8745
GTRI	44	-2148.3053	4390.9466	0.1253	0.9999
HKYG	40	-2160.6195	4406.4536	0.0001	0.9999
HKYIG	41	-2160.1502	4407.7845	0.0000	1.0000
HKYI	40	-2161.6670	4408.5485	0.0000	1.0000
SYMG	41	-2161.3187	4410.1215	0.0000	1.0000
SYMI	41	-2162.9873	4413.4588	0.0000	1.0000
SYMIG	42	-2161.9067	4413.5741	0.0000	1.0000
K2PG	37	-2177.9295	4434.3083	0.0000	1.0000
K2PIG	38	-2177.3586	4435.4144	0.0000	1.0000
K2PI	37	-2179.0134	4436.4761	0.0000	1.0000
GTR	43	-2202.8036	4497.6520	0.0000	1.0000
HKY	39	-2215.9923	4514.9371	0.0000	1.0000
SYM	40	-2215.2350	4515.6847	0.0000	1.0000
F81G	39	-2219.0873	4521.1270	0.0000	1.0000
F81IG	40	-2218.7509	4522.7165	0.0000	1.0000
F81I	39	-2219.9090	4522.7704	0.0000	1.0000
K2P	36	-2232.4430	4541.0945	0.0000	1.0000
JC69G	36	-2235.7378	4547.6841	0.0000	1.0000
JC69I	36	-2236.4191	4549.0467	0.0000	1.0000
JC69IG	37	-2235.5939	4549.6371	0.0000	1.0000
F81	38	-2272.5216	4625.7405	0.0000	1.0000
JC69	35	-2287.0916	4648.1580	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2146.5319	4381.0638	0.7226	0.7226
GTRIG	45	-2147.0733	4384.1466	0.1547	0.8773
GTRI	44	-2148.3053	4384.6106	0.1227	0.9999
HKYG	40	-2160.6195	4401.2389	0.0000	1.0000
HKYIG	41	-2160.1502	4402.3004	0.0000	1.0000
HKYI	40	-2161.6670	4403.3339	0.0000	1.0000
SYMG	41	-2161.3187	4404.6374	0.0000	1.0000
SYMIG	42	-2161.9067	4407.8133	0.0000	1.0000
SYMI	41	-2162.9873	4407.9747	0.0000	1.0000
K2PG	37	-2177.9295	4429.8589	0.0000	1.0000
K2PIG	38	-2177.3586	4430.7171	0.0000	1.0000
K2PI	37	-2179.0134	4432.0267	0.0000	1.0000
GTR	43	-2202.8036	4491.6073	0.0000	1.0000
HKY	39	-2215.9923	4509.9847	0.0000	1.0000
SYM	40	-2215.2350	4510.4701	0.0000	1.0000
F81G	39	-2219.0873	4516.1746	0.0000	1.0000
F81IG	40	-2218.7509	4517.5019	0.0000	1.0000
F81I	39	-2219.9090	4517.8180	0.0000	1.0000
K2P	36	-2232.4430	4536.8859	0.0000	1.0000
JC69G	36	-2235.7378	4543.4755	0.0000	1.0000
JC69I	36	-2236.4191	4544.8381	0.0000	1.0000
JC69IG	37	-2235.5939	4545.1877	0.0000	1.0000
F81	38	-2272.5216	4621.0432	0.0000	1.0000
JC69	35	-2287.0916	4644.1832	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-2146.5319	4579.3841	0.5940	0.5940
HKYG	40	-2160.6195	4581.5300	0.2031	0.7971
GTRI	44	-2148.3053	4582.9308	0.1008	0.8980
HKYI	40	-2161.6670	4583.6250	0.0713	0.9692
GTRIG	45	-2147.0733	4586.9741	0.0134	0.9826
HKYIG	41	-2160.1502	4587.0988	0.0125	0.9951
SYMG	41	-2161.3187	4589.4358	0.0039	0.9990
SYMI	41	-2162.9873	4592.7731	0.0007	0.9998
K2PG	37	-2177.9295	4596.6282	0.0001	0.9999
SYMIG	42	-2161.9067	4597.1190	0.0001	1.0000
K2PI	37	-2179.0134	4598.7960	0.0000	1.0000
K2PIG	38	-2177.3586	4601.9937	0.0000	1.0000
GTR	43	-2202.8036	4685.4202	0.0000	1.0000
HKY	39	-2215.9923	4685.7685	0.0000	1.0000
SYM	40	-2215.2350	4690.7612	0.0000	1.0000
F81G	39	-2219.0873	4691.9584	0.0000	1.0000
F81I	39	-2219.9090	4693.6019	0.0000	1.0000
F81IG	40	-2218.7509	4697.7930	0.0000	1.0000
K2P	36	-2232.4430	4699.1479	0.0000	1.0000
JC69G	36	-2235.7378	4705.7375	0.0000	1.0000
JC69I	36	-2236.4191	4707.1001	0.0000	1.0000
JC69IG	37	-2235.5939	4711.9570	0.0000	1.0000
F81	38	-2272.5216	4792.3198	0.0000	1.0000
JC69	35	-2287.0916	4801.9379	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/chr7_10316.phylip = 1 - 670;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10316.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/chr7_10316.phylip = 1 - 670;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10316.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 (GTRG)]

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


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

