
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:11 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_32302.phylip" (19 taxa, 596 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	-1646.4192	3378.7482	0.4426	0.4426
HKYG	40	-1646.5468	3379.0035	0.3896	0.8322
HKYIG	41	-1646.5467	3381.3099	0.1230	0.9551
GTRI	44	-1644.7774	3384.7418	0.0221	0.9772
GTRG	44	-1645.0182	3385.2233	0.0174	0.9946
GTRIG	45	-1645.0181	3387.5635	0.0054	1.0000
HKY	39	-1672.8689	3429.3493	0.0000	1.0000
GTR	43	-1672.2578	3437.3707	0.0000	1.0000
K2PG	37	-1688.1566	3455.3527	0.0000	1.0000
K2PI	37	-1688.2629	3455.5652	0.0000	1.0000
K2PIG	38	-1688.1620	3457.6453	0.0000	1.0000
SYMG	41	-1687.1510	3462.5186	0.0000	1.0000
SYMI	41	-1687.1674	3462.5514	0.0000	1.0000
F81I	39	-1690.4446	3464.5008	0.0000	1.0000
F81G	39	-1690.4825	3464.5765	0.0000	1.0000
SYMIG	42	-1687.1509	3464.8335	0.0000	1.0000
F81IG	40	-1690.4824	3466.8747	0.0000	1.0000
K2P	36	-1711.4725	3499.7106	0.0000	1.0000
SYM	40	-1710.5166	3506.9432	0.0000	1.0000
F81	38	-1714.4294	3510.1802	0.0000	1.0000
JC69G	36	-1724.6366	3526.0389	0.0000	1.0000
JC69I	36	-1724.6554	3526.0766	0.0000	1.0000
JC69IG	37	-1724.6365	3528.3124	0.0000	1.0000
JC69	35	-1747.6316	3569.7631	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1646.4192	3372.8383	0.4166	0.4166
HKYG	40	-1646.5468	3373.0936	0.3667	0.7833
HKYIG	41	-1646.5467	3375.0933	0.1349	0.9182
GTRI	44	-1644.7774	3377.5549	0.0394	0.9576
GTRG	44	-1645.0182	3378.0363	0.0310	0.9886
GTRIG	45	-1645.0181	3380.0362	0.0114	1.0000
HKY	39	-1672.8689	3423.7378	0.0000	1.0000
GTR	43	-1672.2578	3430.5156	0.0000	1.0000
K2PG	37	-1688.1566	3450.3132	0.0000	1.0000
K2PI	37	-1688.2629	3450.5257	0.0000	1.0000
K2PIG	38	-1688.1620	3452.3240	0.0000	1.0000
SYMG	41	-1687.1510	3456.3020	0.0000	1.0000
SYMI	41	-1687.1674	3456.3348	0.0000	1.0000
SYMIG	42	-1687.1509	3458.3018	0.0000	1.0000
F81I	39	-1690.4446	3458.8893	0.0000	1.0000
F81G	39	-1690.4825	3458.9650	0.0000	1.0000
F81IG	40	-1690.4824	3460.9648	0.0000	1.0000
K2P	36	-1711.4725	3494.9449	0.0000	1.0000
SYM	40	-1710.5166	3501.0333	0.0000	1.0000
F81	38	-1714.4294	3504.8588	0.0000	1.0000
JC69G	36	-1724.6366	3521.2733	0.0000	1.0000
JC69I	36	-1724.6554	3521.3109	0.0000	1.0000
JC69IG	37	-1724.6365	3523.2730	0.0000	1.0000
JC69	35	-1747.6316	3565.2631	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1646.4192	3548.4479	0.5219	0.5219
HKYG	40	-1646.5468	3548.7032	0.4593	0.9812
HKYIG	41	-1646.5467	3555.0932	0.0188	1.0000
GTRI	44	-1644.7774	3570.7254	0.0000	1.0000
GTRG	44	-1645.0182	3571.2069	0.0000	1.0000
GTRIG	45	-1645.0181	3577.5970	0.0000	1.0000
HKY	39	-1672.8689	3594.9572	0.0000	1.0000
K2PG	37	-1688.1566	3612.7521	0.0000	1.0000
K2PI	37	-1688.2629	3612.9646	0.0000	1.0000
K2PIG	38	-1688.1620	3619.1531	0.0000	1.0000
GTR	43	-1672.2578	3619.2960	0.0000	1.0000
F81I	39	-1690.4446	3630.1086	0.0000	1.0000
F81G	39	-1690.4825	3630.1844	0.0000	1.0000
SYMG	41	-1687.1510	3636.3019	0.0000	1.0000
SYMI	41	-1687.1674	3636.3346	0.0000	1.0000
F81IG	40	-1690.4824	3636.5744	0.0000	1.0000
SYMIG	42	-1687.1509	3642.6919	0.0000	1.0000
K2P	36	-1711.4725	3652.9936	0.0000	1.0000
F81	38	-1714.4294	3671.6880	0.0000	1.0000
SYM	40	-1710.5166	3676.6429	0.0000	1.0000
JC69G	36	-1724.6366	3679.3219	0.0000	1.0000
JC69I	36	-1724.6554	3679.3596	0.0000	1.0000
JC69IG	37	-1724.6365	3685.7119	0.0000	1.0000
JC69	35	-1747.6316	3718.9216	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_32302.phylip = 1 - 596;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_32302.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_32302.phylip = 1 - 596;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_32302.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_32302.phylip = 1 - 596;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr1_32302.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

