
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:00 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/chr6_7142.phylip" (19 taxa, 874 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	45	-3221.1366	6537.2731	0.9237	0.9237
GTRG	44	-3225.3629	6543.5027	0.0410	0.9647
HKYG	40	-3229.9480	6543.8337	0.0347	0.9994
GTRI	44	-3230.3022	6553.3812	0.0003	0.9997
HKYI	40	-3234.7570	6553.4516	0.0003	1.0000
HKYIG	41	-3236.0444	6558.2282	0.0000	1.0000
K2PG	37	-3281.3528	6640.0692	0.0000	1.0000
K2PIG	38	-3280.5621	6640.6740	0.0000	1.0000
SYMG	41	-3279.8749	6645.8892	0.0000	1.0000
SYMIG	42	-3278.9908	6646.3282	0.0000	1.0000
K2PI	37	-3289.7815	6656.9267	0.0000	1.0000
SYMI	41	-3287.7598	6661.6590	0.0000	1.0000
F81IG	40	-3294.8716	6673.6807	0.0000	1.0000
F81G	39	-3297.6992	6677.1394	0.0000	1.0000
F81I	39	-3300.4280	6682.5970	0.0000	1.0000
GTR	43	-3331.3526	6753.2643	0.0000	1.0000
JC69IG	37	-3338.4740	6754.3117	0.0000	1.0000
JC69G	36	-3341.9469	6759.0766	0.0000	1.0000
JC69I	36	-3344.8577	6764.8982	0.0000	1.0000
HKY	39	-3345.0037	6771.7483	0.0000	1.0000
K2P	36	-3403.9250	6883.0327	0.0000	1.0000
SYM	40	-3402.1945	6888.3266	0.0000	1.0000
F81	38	-3405.0388	6889.6272	0.0000	1.0000
JC69	35	-3453.5969	6980.2009	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-3221.1366	6532.2731	0.9414	0.9414
GTRG	44	-3225.3629	6538.7259	0.0374	0.9787
HKYG	40	-3229.9480	6539.8961	0.0208	0.9995
GTRI	44	-3230.3022	6548.6044	0.0003	0.9998
HKYI	40	-3234.7570	6549.5140	0.0002	1.0000
HKYIG	41	-3236.0444	6554.0888	0.0000	1.0000
K2PG	37	-3281.3528	6636.7056	0.0000	1.0000
K2PIG	38	-3280.5621	6637.1243	0.0000	1.0000
SYMG	41	-3279.8749	6641.7498	0.0000	1.0000
SYMIG	42	-3278.9908	6641.9816	0.0000	1.0000
K2PI	37	-3289.7815	6653.5630	0.0000	1.0000
SYMI	41	-3287.7598	6657.5196	0.0000	1.0000
F81IG	40	-3294.8716	6669.7432	0.0000	1.0000
F81G	39	-3297.6992	6673.3984	0.0000	1.0000
F81I	39	-3300.4280	6678.8560	0.0000	1.0000
GTR	43	-3331.3526	6748.7052	0.0000	1.0000
JC69IG	37	-3338.4740	6750.9480	0.0000	1.0000
JC69G	36	-3341.9469	6755.8938	0.0000	1.0000
JC69I	36	-3344.8577	6761.7154	0.0000	1.0000
HKY	39	-3345.0037	6768.0073	0.0000	1.0000
K2P	36	-3403.9250	6879.8499	0.0000	1.0000
SYM	40	-3402.1945	6884.3890	0.0000	1.0000
F81	38	-3405.0388	6886.0775	0.0000	1.0000
JC69	35	-3453.5969	6977.1937	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-3229.9480	6730.8193	0.9914	0.9914
HKYI	40	-3234.7570	6740.4373	0.0081	0.9995
GTRIG	45	-3221.1366	6747.0617	0.0003	0.9998
GTRG	44	-3225.3629	6748.7414	0.0001	0.9999
HKYIG	41	-3236.0444	6749.7851	0.0001	1.0000
GTRI	44	-3230.3022	6758.6199	0.0000	1.0000
K2PG	37	-3281.3528	6813.3096	0.0000	1.0000
K2PIG	38	-3280.5621	6818.5013	0.0000	1.0000
K2PI	37	-3289.7815	6830.1670	0.0000	1.0000
SYMG	41	-3279.8749	6837.4461	0.0000	1.0000
SYMIG	42	-3278.9908	6842.4510	0.0000	1.0000
SYMI	41	-3287.7598	6853.2159	0.0000	1.0000
F81G	39	-3297.6992	6859.5485	0.0000	1.0000
F81IG	40	-3294.8716	6860.6664	0.0000	1.0000
F81I	39	-3300.4280	6865.0062	0.0000	1.0000
JC69IG	37	-3338.4740	6927.5520	0.0000	1.0000
JC69G	36	-3341.9469	6927.7247	0.0000	1.0000
JC69I	36	-3344.8577	6933.5463	0.0000	1.0000
GTR	43	-3331.3526	6953.9477	0.0000	1.0000
HKY	39	-3345.0037	6954.1574	0.0000	1.0000
K2P	36	-3403.9250	7051.6808	0.0000	1.0000
F81	38	-3405.0388	7067.4546	0.0000	1.0000
SYM	40	-3402.1945	7075.3122	0.0000	1.0000
JC69	35	-3453.5969	7144.2516	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7142.phylip = 1 - 874;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7142.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7142.phylip = 1 - 874;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7142.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.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/chr6_7142.phylip = 1 - 874;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_7142.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

