
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:45 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_10269.phylip" (19 taxa, 832 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	-3312.2833	6719.8338	0.6853	0.6853
GTRG	44	-3314.1894	6721.4105	0.3115	0.9968
HKYIG	41	-3322.7998	6731.9591	0.0016	0.9984
HKYG	40	-3324.5804	6733.3074	0.0008	0.9992
GTRI	44	-3320.1752	6733.3822	0.0008	1.0000
HKYI	40	-3329.0837	6742.3139	0.0000	1.0000
K2PIG	38	-3347.4462	6774.6300	0.0000	1.0000
K2PG	37	-3348.6489	6774.8393	0.0000	1.0000
SYMIG	42	-3346.2011	6780.9801	0.0000	1.0000
SYMG	41	-3348.3877	6783.1350	0.0000	1.0000
K2PI	37	-3354.0855	6785.7127	0.0000	1.0000
SYMI	41	-3352.6313	6791.6222	0.0000	1.0000
F81IG	40	-3381.7331	6847.6128	0.0000	1.0000
F81G	39	-3383.5687	6849.0768	0.0000	1.0000
F81I	39	-3389.0034	6859.9462	0.0000	1.0000
JC69IG	37	-3400.3942	6878.3299	0.0000	1.0000
JC69G	36	-3402.3811	6880.1131	0.0000	1.0000
JC69I	36	-3407.3275	6890.0058	0.0000	1.0000
GTR	43	-3458.2865	7007.3751	0.0000	1.0000
HKY	39	-3476.6665	7035.2724	0.0000	1.0000
K2P	36	-3505.9083	7087.1676	0.0000	1.0000
SYM	40	-3502.9402	7090.0271	0.0000	1.0000
F81	38	-3526.5189	7132.7755	0.0000	1.0000
JC69	35	-3547.1062	7167.3783	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-3312.2833	6714.5667	0.7106	0.7106
GTRG	44	-3314.1894	6716.3787	0.2872	0.9977
HKYIG	41	-3322.7998	6727.5996	0.0011	0.9988
GTRI	44	-3320.1752	6728.3504	0.0007	0.9995
HKYG	40	-3324.5804	6729.1608	0.0005	1.0000
HKYI	40	-3329.0837	6738.1673	0.0000	1.0000
K2PIG	38	-3347.4462	6770.8923	0.0000	1.0000
K2PG	37	-3348.6489	6771.2977	0.0000	1.0000
SYMIG	42	-3346.2011	6776.4022	0.0000	1.0000
SYMG	41	-3348.3877	6778.7755	0.0000	1.0000
K2PI	37	-3354.0855	6782.1711	0.0000	1.0000
SYMI	41	-3352.6313	6787.2627	0.0000	1.0000
F81IG	40	-3381.7331	6843.4662	0.0000	1.0000
F81G	39	-3383.5687	6845.1374	0.0000	1.0000
F81I	39	-3389.0034	6856.0068	0.0000	1.0000
JC69IG	37	-3400.3942	6874.7884	0.0000	1.0000
JC69G	36	-3402.3811	6876.7622	0.0000	1.0000
JC69I	36	-3407.3275	6886.6549	0.0000	1.0000
GTR	43	-3458.2865	7002.5731	0.0000	1.0000
HKY	39	-3476.6665	7031.3330	0.0000	1.0000
K2P	36	-3505.9083	7083.8166	0.0000	1.0000
SYM	40	-3502.9402	7085.8805	0.0000	1.0000
F81	38	-3526.5189	7129.0378	0.0000	1.0000
JC69	35	-3547.1062	7164.2125	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-3324.5804	6918.1141	0.7844	0.7844
HKYIG	41	-3322.7998	6921.2767	0.1613	0.9457
GTRG	44	-3314.1894	6924.2273	0.0369	0.9826
HKYI	40	-3329.0837	6927.1206	0.0087	0.9913
GTRIG	45	-3312.2833	6927.1391	0.0086	0.9999
GTRI	44	-3320.1752	6936.1990	0.0001	1.0000
K2PG	37	-3348.6489	6946.0795	0.0000	1.0000
K2PIG	38	-3347.4462	6950.3980	0.0000	1.0000
K2PI	37	-3354.0855	6956.9529	0.0000	1.0000
SYMG	41	-3348.3877	6972.4526	0.0000	1.0000
SYMIG	42	-3346.2011	6974.8032	0.0000	1.0000
SYMI	41	-3352.6313	6980.9398	0.0000	1.0000
F81G	39	-3383.5687	7029.3669	0.0000	1.0000
F81IG	40	-3381.7331	7032.4195	0.0000	1.0000
F81I	39	-3389.0034	7040.2363	0.0000	1.0000
JC69G	36	-3402.3811	7046.8201	0.0000	1.0000
JC69IG	37	-3400.3942	7049.5702	0.0000	1.0000
JC69I	36	-3407.3275	7056.7129	0.0000	1.0000
GTR	43	-3458.2865	7205.6979	0.0000	1.0000
HKY	39	-3476.6665	7215.5625	0.0000	1.0000
K2P	36	-3505.9083	7253.8746	0.0000	1.0000
SYM	40	-3502.9402	7274.8338	0.0000	1.0000
F81	38	-3526.5189	7308.5434	0.0000	1.0000
JC69	35	-3547.1062	7329.5466	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/chr7_10269.phylip = 1 - 832;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10269.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/chr7_10269.phylip = 1 - 832;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10269.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/chr7_10269.phylip = 1 - 832;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10269.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

