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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	44	-2193.3433	4480.9924	0.5184	0.5184
GTRI	44	-2194.0840	4482.4738	0.2472	0.7656
GTRIG	45	-2193.3504	4483.3037	0.1632	0.9288
HKYG	40	-2200.3654	4485.9207	0.0441	0.9729
HKYIG	41	-2200.3547	4488.1674	0.0143	0.9873
HKYI	40	-2201.6081	4488.4061	0.0127	1.0000
K2PG	37	-2226.3080	4531.0443	0.0000	1.0000
K2PIG	38	-2226.2020	4533.0791	0.0000	1.0000
K2PI	37	-2228.2521	4534.9326	0.0000	1.0000
SYMG	41	-2224.1827	4535.8235	0.0000	1.0000
SYMI	41	-2224.8928	4537.2435	0.0000	1.0000
SYMIG	42	-2224.1211	4537.9754	0.0000	1.0000
GTR	43	-2225.7455	4543.5068	0.0000	1.0000
HKY	39	-2234.0166	4550.9621	0.0000	1.0000
F81G	39	-2248.2160	4579.3609	0.0000	1.0000
F81IG	40	-2248.1788	4581.5474	0.0000	1.0000
F81I	39	-2249.6576	4582.2440	0.0000	1.0000
SYM	40	-2254.2739	4593.7377	0.0000	1.0000
K2P	36	-2260.0926	4596.3740	0.0000	1.0000
JC69G	36	-2269.7406	4615.6699	0.0000	1.0000
JC69IG	37	-2269.6917	4617.8118	0.0000	1.0000
JC69I	36	-2271.6030	4619.3947	0.0000	1.0000
F81	38	-2279.8685	4640.4121	0.0000	1.0000
JC69	35	-2300.1532	4674.2624	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2193.3433	4474.6867	0.5200	0.5200
GTRI	44	-2194.0840	4476.1680	0.2479	0.7680
GTRIG	45	-2193.3504	4476.7009	0.1900	0.9579
HKYG	40	-2200.3654	4480.7308	0.0253	0.9833
HKYIG	41	-2200.3547	4482.7094	0.0094	0.9927
HKYI	40	-2201.6081	4483.2162	0.0073	1.0000
K2PG	37	-2226.3080	4526.6159	0.0000	1.0000
K2PIG	38	-2226.2020	4528.4040	0.0000	1.0000
SYMG	41	-2224.1827	4530.3655	0.0000	1.0000
K2PI	37	-2228.2521	4530.5043	0.0000	1.0000
SYMI	41	-2224.8928	4531.7855	0.0000	1.0000
SYMIG	42	-2224.1211	4532.2421	0.0000	1.0000
GTR	43	-2225.7455	4537.4909	0.0000	1.0000
HKY	39	-2234.0166	4546.0332	0.0000	1.0000
F81G	39	-2248.2160	4574.4320	0.0000	1.0000
F81IG	40	-2248.1788	4576.3576	0.0000	1.0000
F81I	39	-2249.6576	4577.3151	0.0000	1.0000
SYM	40	-2254.2739	4588.5478	0.0000	1.0000
K2P	36	-2260.0926	4592.1853	0.0000	1.0000
JC69G	36	-2269.7406	4611.4813	0.0000	1.0000
JC69IG	37	-2269.6917	4613.3834	0.0000	1.0000
JC69I	36	-2271.6030	4615.2060	0.0000	1.0000
F81	38	-2279.8685	4635.7371	0.0000	1.0000
JC69	35	-2300.1532	4670.3063	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2200.3654	4661.2006	0.7511	0.7511
HKYI	40	-2201.6081	4663.6860	0.2168	0.9679
HKYIG	41	-2200.3547	4667.6910	0.0293	0.9972
GTRG	44	-2193.3433	4673.2035	0.0019	0.9990
GTRI	44	-2194.0840	4674.6848	0.0009	0.9999
GTRIG	45	-2193.3504	4679.7294	0.0001	1.0000
K2PG	37	-2226.3080	4693.5505	0.0000	1.0000
K2PI	37	-2228.2521	4697.4388	0.0000	1.0000
K2PIG	38	-2226.2020	4699.8504	0.0000	1.0000
SYMG	41	-2224.1827	4715.3471	0.0000	1.0000
SYMI	41	-2224.8928	4716.7671	0.0000	1.0000
SYMIG	42	-2224.1211	4721.7354	0.0000	1.0000
HKY	39	-2234.0166	4721.9912	0.0000	1.0000
GTR	43	-2225.7455	4731.4960	0.0000	1.0000
F81G	39	-2248.2160	4750.3901	0.0000	1.0000
F81I	39	-2249.6576	4753.2732	0.0000	1.0000
K2P	36	-2260.0926	4754.6081	0.0000	1.0000
F81IG	40	-2248.1788	4756.8274	0.0000	1.0000
SYM	40	-2254.2739	4769.0176	0.0000	1.0000
JC69G	36	-2269.7406	4773.9041	0.0000	1.0000
JC69I	36	-2271.6030	4777.6289	0.0000	1.0000
JC69IG	37	-2269.6917	4780.3180	0.0000	1.0000
F81	38	-2279.8685	4807.1834	0.0000	1.0000
JC69	35	-2300.1532	4828.2174	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/chr1_5292.phylip = 1 - 673;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_5292.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/chr1_5292.phylip = 1 - 673;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_5292.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 (HKYG)]

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

