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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	40	-1701.4121	3487.9174	0.3486	0.3486
HKYG	40	-1701.4357	3487.9645	0.3405	0.6892
GTRG	44	-1698.1670	3490.5216	0.0948	0.7840
GTRI	44	-1698.2937	3490.7748	0.0835	0.8675
HKYIG	41	-1701.7914	3490.9390	0.0770	0.9445
GTRIG	45	-1697.5586	3491.5960	0.0554	0.9999
F81I	39	-1711.4348	3505.7068	0.0000	0.9999
F81G	39	-1711.4728	3505.7827	0.0000	1.0000
F81IG	40	-1712.1953	3509.4837	0.0000	1.0000
K2PI	37	-1719.8545	3518.0551	0.0000	1.0000
K2PG	37	-1719.8832	3518.1126	0.0000	1.0000
K2PIG	38	-1719.8832	3520.3546	0.0000	1.0000
SYMI	41	-1719.3013	3525.9588	0.0000	1.0000
SYMG	41	-1719.3302	3526.0165	0.0000	1.0000
SYMIG	42	-1719.5499	3528.7260	0.0000	1.0000
HKY	39	-1724.8413	3532.5198	0.0000	1.0000
JC69I	36	-1728.4566	3533.0244	0.0000	1.0000
JC69G	36	-1728.4964	3533.1039	0.0000	1.0000
GTR	43	-1722.2247	3536.3527	0.0000	1.0000
JC69IG	37	-1729.7223	3537.7909	0.0000	1.0000
F81	38	-1734.3284	3549.2450	0.0000	1.0000
K2P	36	-1743.0277	3562.1665	0.0000	1.0000
SYM	40	-1742.2322	3569.5576	0.0000	1.0000
JC69	35	-1751.2805	3576.4440	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1701.4121	3482.8242	0.2915	0.2915
HKYG	40	-1701.4357	3482.8714	0.2847	0.5762
GTRG	44	-1698.1670	3484.3341	0.1370	0.7132
GTRI	44	-1698.2937	3484.5873	0.1207	0.8339
GTRIG	45	-1697.5586	3485.1171	0.0926	0.9265
HKYIG	41	-1701.7914	3485.5828	0.0734	0.9999
F81I	39	-1711.4348	3500.8696	0.0000	1.0000
F81G	39	-1711.4728	3500.9455	0.0000	1.0000
F81IG	40	-1712.1953	3504.3905	0.0000	1.0000
K2PI	37	-1719.8545	3513.7089	0.0000	1.0000
K2PG	37	-1719.8832	3513.7663	0.0000	1.0000
K2PIG	38	-1719.8832	3515.7663	0.0000	1.0000
SYMI	41	-1719.3013	3520.6027	0.0000	1.0000
SYMG	41	-1719.3302	3520.6603	0.0000	1.0000
SYMIG	42	-1719.5499	3523.0999	0.0000	1.0000
HKY	39	-1724.8413	3527.6826	0.0000	1.0000
JC69I	36	-1728.4566	3528.9133	0.0000	1.0000
JC69G	36	-1728.4964	3528.9928	0.0000	1.0000
GTR	43	-1722.2247	3530.4494	0.0000	1.0000
JC69IG	37	-1729.7223	3533.4447	0.0000	1.0000
F81	38	-1734.3284	3544.6568	0.0000	1.0000
K2P	36	-1743.0277	3558.0553	0.0000	1.0000
SYM	40	-1742.2322	3564.4645	0.0000	1.0000
JC69	35	-1751.2805	3572.5611	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1701.4121	3664.0010	0.4986	0.4986
HKYG	40	-1701.4357	3664.0481	0.4870	0.9856
HKYIG	41	-1701.7914	3671.2890	0.0130	0.9986
F81I	39	-1711.4348	3677.5170	0.0006	0.9992
F81G	39	-1711.4728	3677.5929	0.0006	0.9998
K2PI	37	-1719.8545	3681.2974	0.0001	0.9998
K2PG	37	-1719.8832	3681.3548	0.0001	0.9999
GTRG	44	-1698.1670	3683.6285	0.0000	1.0000
GTRI	44	-1698.2937	3683.8818	0.0000	1.0000
F81IG	40	-1712.1953	3685.5673	0.0000	1.0000
K2PIG	38	-1719.8832	3687.8843	0.0000	1.0000
GTRIG	45	-1697.5586	3688.9410	0.0000	1.0000
JC69I	36	-1728.4566	3691.9724	0.0000	1.0000
JC69G	36	-1728.4964	3692.0518	0.0000	1.0000
JC69IG	37	-1729.7223	3701.0332	0.0000	1.0000
HKY	39	-1724.8413	3704.3299	0.0000	1.0000
SYMI	41	-1719.3013	3706.3089	0.0000	1.0000
SYMG	41	-1719.3302	3706.3665	0.0000	1.0000
SYMIG	42	-1719.5499	3713.3355	0.0000	1.0000
F81	38	-1734.3284	3716.7747	0.0000	1.0000
K2P	36	-1743.0277	3721.1144	0.0000	1.0000
GTR	43	-1722.2247	3725.2144	0.0000	1.0000
JC69	35	-1751.2805	3731.0907	0.0000	1.0000
SYM	40	-1742.2322	3745.6412	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_32365.phylip = 1 - 685;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_32365.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_32365.phylip = 1 - 685;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_32365.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_32365.phylip = 1 - 685;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_32365.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

