
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:39 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_3291.phylip" (19 taxa, 973 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	-3434.9961	6962.2594	0.9989	0.9989
GTRIG	45	-3440.6856	6975.8373	0.0011	1.0000
GTRI	44	-3446.4654	6985.1980	0.0000	1.0000
HKYG	40	-3455.3260	6994.1713	0.0000	1.0000
HKYIG	41	-3456.4856	6998.6705	0.0000	1.0000
HKYI	40	-3461.6309	7006.7811	0.0000	1.0000
SYMG	41	-3474.7949	7035.2891	0.0000	1.0000
K2PG	37	-3481.5098	7040.0271	0.0000	1.0000
SYMIG	42	-3476.2217	7040.3273	0.0000	1.0000
K2PIG	38	-3482.7621	7044.6977	0.0000	1.0000
SYMI	41	-3481.2141	7048.1273	0.0000	1.0000
K2PI	37	-3487.7718	7052.5511	0.0000	1.0000
GTR	43	-3494.1846	7078.4423	0.0000	1.0000
HKY	39	-3508.8689	7099.0818	0.0000	1.0000
SYM	40	-3527.0058	7137.5310	0.0000	1.0000
K2P	36	-3534.0249	7142.8959	0.0000	1.0000
F81G	39	-3546.0755	7173.4951	0.0000	1.0000
F81IG	40	-3552.3345	7188.1883	0.0000	1.0000
F81I	39	-3556.5673	7194.4787	0.0000	1.0000
JC69G	36	-3570.0158	7214.8778	0.0000	1.0000
JC69IG	37	-3572.1672	7221.3418	0.0000	1.0000
JC69I	36	-3576.4459	7227.7380	0.0000	1.0000
F81	38	-3600.4881	7280.1496	0.0000	1.0000
JC69	35	-3618.5129	7309.7153	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-3434.9961	6957.9922	0.9987	0.9987
GTRIG	45	-3440.6856	6971.3713	0.0012	1.0000
GTRI	44	-3446.4654	6980.9307	0.0000	1.0000
HKYG	40	-3455.3260	6990.6520	0.0000	1.0000
HKYIG	41	-3456.4856	6994.9713	0.0000	1.0000
HKYI	40	-3461.6309	7003.2618	0.0000	1.0000
SYMG	41	-3474.7949	7031.5898	0.0000	1.0000
SYMIG	42	-3476.2217	7036.4434	0.0000	1.0000
K2PG	37	-3481.5098	7037.0196	0.0000	1.0000
K2PIG	38	-3482.7621	7041.5243	0.0000	1.0000
SYMI	41	-3481.2141	7044.4281	0.0000	1.0000
K2PI	37	-3487.7718	7049.5437	0.0000	1.0000
GTR	43	-3494.1846	7074.3691	0.0000	1.0000
HKY	39	-3508.8689	7095.7377	0.0000	1.0000
SYM	40	-3527.0058	7134.0116	0.0000	1.0000
K2P	36	-3534.0249	7140.0498	0.0000	1.0000
F81G	39	-3546.0755	7170.1511	0.0000	1.0000
F81IG	40	-3552.3345	7184.6690	0.0000	1.0000
F81I	39	-3556.5673	7191.1346	0.0000	1.0000
JC69G	36	-3570.0158	7212.0316	0.0000	1.0000
JC69IG	37	-3572.1672	7218.3344	0.0000	1.0000
JC69I	36	-3576.4459	7224.8918	0.0000	1.0000
F81	38	-3600.4881	7276.9761	0.0000	1.0000
JC69	35	-3618.5129	7307.0258	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-3434.9961	7172.7291	0.9985	0.9985
HKYG	40	-3455.3260	7185.8674	0.0014	0.9999
GTRIG	45	-3440.6856	7190.9886	0.0001	1.0000
HKYIG	41	-3456.4856	7195.0670	0.0000	1.0000
GTRI	44	-3446.4654	7195.6676	0.0000	1.0000
HKYI	40	-3461.6309	7198.4771	0.0000	1.0000
K2PG	37	-3481.5098	7217.5938	0.0000	1.0000
K2PIG	38	-3482.7621	7226.9789	0.0000	1.0000
K2PI	37	-3487.7718	7230.1179	0.0000	1.0000
SYMG	41	-3474.7949	7231.6856	0.0000	1.0000
SYMIG	42	-3476.2217	7241.4195	0.0000	1.0000
SYMI	41	-3481.2141	7244.5238	0.0000	1.0000
GTR	43	-3494.1846	7284.2256	0.0000	1.0000
HKY	39	-3508.8689	7286.0727	0.0000	1.0000
K2P	36	-3534.0249	7315.7436	0.0000	1.0000
SYM	40	-3527.0058	7329.2270	0.0000	1.0000
F81G	39	-3546.0755	7360.4861	0.0000	1.0000
F81IG	40	-3552.3345	7379.8844	0.0000	1.0000
F81I	39	-3556.5673	7381.4696	0.0000	1.0000
JC69G	36	-3570.0158	7387.7254	0.0000	1.0000
JC69IG	37	-3572.1672	7398.9086	0.0000	1.0000
JC69I	36	-3576.4459	7400.5856	0.0000	1.0000
F81	38	-3600.4881	7462.4307	0.0000	1.0000
JC69	35	-3618.5129	7477.8393	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/chr6_3291.phylip = 1 - 973;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3291.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/chr6_3291.phylip = 1 - 973;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3291.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/chr6_3291.phylip = 1 - 973;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3291.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

