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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	40	-3437.5237	6959.7874	0.5446	0.5446
GTRG	44	-3433.9522	6961.6601	0.2135	0.7581
HKYIG	41	-3437.5356	6962.0552	0.1752	0.9333
GTRIG	45	-3433.9832	6963.9925	0.0665	0.9998
HKYI	40	-3445.7335	6976.2069	0.0001	0.9999
GTRI	44	-3442.3111	6978.3779	0.0001	1.0000
SYMG	41	-3488.9892	7064.9625	0.0000	1.0000
SYMIG	42	-3488.9687	7067.1721	0.0000	1.0000
SYMI	41	-3495.7734	7078.5309	0.0000	1.0000
K2PG	37	-3504.5555	7087.1570	0.0000	1.0000
K2PIG	38	-3504.4566	7089.1841	0.0000	1.0000
K2PI	37	-3510.0947	7098.2354	0.0000	1.0000
HKY	39	-3516.3510	7115.2042	0.0000	1.0000
GTR	43	-3512.4060	7116.3040	0.0000	1.0000
F81G	39	-3518.1459	7118.7940	0.0000	1.0000
F81IG	40	-3518.0132	7120.7663	0.0000	1.0000
F81I	39	-3523.4292	7129.3605	0.0000	1.0000
SYM	40	-3561.1230	7206.9859	0.0000	1.0000
JC69G	36	-3566.4947	7208.8170	0.0000	1.0000
JC69IG	37	-3566.2941	7210.6342	0.0000	1.0000
JC69I	36	-3570.6713	7217.1701	0.0000	1.0000
K2P	36	-3575.0743	7225.9762	0.0000	1.0000
F81	38	-3586.4574	7253.1857	0.0000	1.0000
JC69	35	-3631.3646	7336.3447	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-3437.5237	6955.0475	0.4449	0.4449
GTRG	44	-3433.9522	6955.9043	0.2898	0.7347
HKYIG	41	-3437.5356	6957.0711	0.1617	0.8964
GTRIG	45	-3433.9832	6957.9663	0.1034	0.9998
HKYI	40	-3445.7335	6971.4671	0.0001	0.9999
GTRI	44	-3442.3111	6972.6221	0.0001	1.0000
SYMG	41	-3488.9892	7059.9785	0.0000	1.0000
SYMIG	42	-3488.9687	7061.9374	0.0000	1.0000
SYMI	41	-3495.7734	7073.5468	0.0000	1.0000
K2PG	37	-3504.5555	7083.1109	0.0000	1.0000
K2PIG	38	-3504.4566	7084.9132	0.0000	1.0000
K2PI	37	-3510.0947	7094.1894	0.0000	1.0000
HKY	39	-3516.3510	7110.7021	0.0000	1.0000
GTR	43	-3512.4060	7110.8120	0.0000	1.0000
F81G	39	-3518.1459	7114.2919	0.0000	1.0000
F81IG	40	-3518.0132	7116.0264	0.0000	1.0000
F81I	39	-3523.4292	7124.8584	0.0000	1.0000
SYM	40	-3561.1230	7202.2460	0.0000	1.0000
JC69G	36	-3566.4947	7204.9895	0.0000	1.0000
JC69IG	37	-3566.2941	7206.5881	0.0000	1.0000
JC69I	36	-3570.6713	7213.3425	0.0000	1.0000
K2P	36	-3575.0743	7222.1486	0.0000	1.0000
F81	38	-3586.4574	7248.9148	0.0000	1.0000
JC69	35	-3631.3646	7332.7292	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-3437.5237	7138.9333	0.9645	0.9645
HKYIG	41	-3437.5356	7145.5541	0.0352	0.9997
HKYI	40	-3445.7335	7155.3529	0.0003	0.9999
GTRG	44	-3433.9522	7158.1787	0.0001	1.0000
GTRIG	45	-3433.9832	7164.8379	0.0000	1.0000
GTRI	44	-3442.3111	7174.8965	0.0000	1.0000
SYMG	41	-3488.9892	7248.4614	0.0000	1.0000
K2PG	37	-3504.5555	7253.2053	0.0000	1.0000
SYMIG	42	-3488.9687	7255.0175	0.0000	1.0000
K2PIG	38	-3504.4566	7259.6047	0.0000	1.0000
SYMI	41	-3495.7734	7262.0298	0.0000	1.0000
K2PI	37	-3510.0947	7264.2838	0.0000	1.0000
HKY	39	-3516.3510	7289.9907	0.0000	1.0000
F81G	39	-3518.1459	7293.5805	0.0000	1.0000
F81IG	40	-3518.0132	7299.9122	0.0000	1.0000
F81I	39	-3523.4292	7304.1470	0.0000	1.0000
GTR	43	-3512.4060	7308.4892	0.0000	1.0000
JC69G	36	-3566.4947	7370.4867	0.0000	1.0000
JC69IG	37	-3566.2941	7376.6825	0.0000	1.0000
JC69I	36	-3570.6713	7378.8398	0.0000	1.0000
SYM	40	-3561.1230	7386.1318	0.0000	1.0000
K2P	36	-3575.0743	7387.6458	0.0000	1.0000
F81	38	-3586.4574	7423.6064	0.0000	1.0000
JC69	35	-3631.3646	7493.6293	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr23_707.phylip = 1 - 733;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr23_707.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;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr23_707.phylip = 1 - 733;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr23_707.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;


[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/chr23_707.phylip = 1 - 733;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr23_707.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

