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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-3536.9932	7166.7864	0.5100	0.5100
GTRIG	45	-3535.9214	7166.8671	0.4898	0.9999
SYMG	41	-3549.9039	7185.9672	0.0000	0.9999
HKYG	40	-3551.0715	7186.0995	0.0000	0.9999
GTRI	44	-3546.7406	7186.2812	0.0000	0.9999
SYMIG	42	-3549.0406	7186.4488	0.0000	1.0000
HKYIG	41	-3550.2929	7186.7453	0.0000	1.0000
K2PIG	38	-3557.7051	7194.9769	0.0000	1.0000
K2PG	37	-3559.8992	7197.1783	0.0000	1.0000
SYMI	41	-3555.8420	7197.8433	0.0000	1.0000
HKYI	40	-3558.1305	7200.2176	0.0000	1.0000
K2PI	37	-3565.4258	7208.2314	0.0000	1.0000
GTR	43	-3630.1164	7350.8140	0.0000	1.0000
HKY	39	-3639.5672	7360.8934	0.0000	1.0000
SYM	40	-3638.5164	7360.9894	0.0000	1.0000
F81G	39	-3639.7758	7361.3106	0.0000	1.0000
F81IG	40	-3638.8745	7361.7056	0.0000	1.0000
F81I	39	-3642.4603	7366.6797	0.0000	1.0000
JC69G	36	-3645.7731	7366.7443	0.0000	1.0000
JC69IG	37	-3644.7336	7366.8470	0.0000	1.0000
K2P	36	-3647.1586	7369.5152	0.0000	1.0000
JC69I	36	-3647.8550	7370.9081	0.0000	1.0000
F81	38	-3725.6337	7530.8342	0.0000	1.0000
JC69	35	-3730.7326	7534.4868	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-3535.9214	7161.8428	0.5179	0.5179
GTRG	44	-3536.9932	7161.9864	0.4820	0.9999
GTRI	44	-3546.7406	7181.4812	0.0000	0.9999
SYMG	41	-3549.9039	7181.8077	0.0000	0.9999
SYMIG	42	-3549.0406	7182.0812	0.0000	1.0000
HKYG	40	-3551.0715	7182.1429	0.0000	1.0000
HKYIG	41	-3550.2929	7182.5859	0.0000	1.0000
K2PIG	38	-3557.7051	7191.4102	0.0000	1.0000
SYMI	41	-3555.8420	7193.6839	0.0000	1.0000
K2PG	37	-3559.8992	7193.7985	0.0000	1.0000
HKYI	40	-3558.1305	7196.2611	0.0000	1.0000
K2PI	37	-3565.4258	7204.8516	0.0000	1.0000
GTR	43	-3630.1164	7346.2329	0.0000	1.0000
SYM	40	-3638.5164	7357.0328	0.0000	1.0000
HKY	39	-3639.5672	7357.1343	0.0000	1.0000
F81G	39	-3639.7758	7357.5516	0.0000	1.0000
F81IG	40	-3638.8745	7357.7490	0.0000	1.0000
F81I	39	-3642.4603	7362.9207	0.0000	1.0000
JC69IG	37	-3644.7336	7363.4672	0.0000	1.0000
JC69G	36	-3645.7731	7363.5462	0.0000	1.0000
K2P	36	-3647.1586	7366.3171	0.0000	1.0000
JC69I	36	-3647.8550	7367.7100	0.0000	1.0000
F81	38	-3725.6337	7527.2674	0.0000	1.0000
JC69	35	-3730.7326	7531.4652	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-3559.8992	7370.2327	0.4701	0.4701
GTRG	44	-3536.9932	7371.8001	0.2147	0.6848
K2PIG	38	-3557.7051	7372.6129	0.1430	0.8278
HKYG	40	-3551.0715	7372.8827	0.1250	0.9528
GTRIG	45	-3535.9214	7376.4250	0.0213	0.9740
SYMG	41	-3549.9039	7377.3160	0.0136	0.9877
HKYIG	41	-3550.2929	7378.0941	0.0092	0.9969
K2PI	37	-3565.4258	7381.2858	0.0019	0.9987
SYMIG	42	-3549.0406	7382.3580	0.0011	0.9998
HKYI	40	-3558.1305	7387.0008	0.0001	1.0000
SYMI	41	-3555.8420	7389.1921	0.0000	1.0000
GTRI	44	-3546.7406	7391.2949	0.0000	1.0000
JC69G	36	-3645.7731	7535.2120	0.0000	1.0000
K2P	36	-3647.1586	7537.9829	0.0000	1.0000
JC69I	36	-3647.8550	7539.3757	0.0000	1.0000
JC69IG	37	-3644.7336	7539.9014	0.0000	1.0000
HKY	39	-3639.5672	7543.1056	0.0000	1.0000
F81G	39	-3639.7758	7543.5228	0.0000	1.0000
SYM	40	-3638.5164	7547.7725	0.0000	1.0000
F81IG	40	-3638.8745	7548.4887	0.0000	1.0000
F81I	39	-3642.4603	7548.8919	0.0000	1.0000
GTR	43	-3630.1164	7551.2781	0.0000	1.0000
JC69	35	-3730.7326	7698.3625	0.0000	1.0000
F81	38	-3725.6337	7708.4701	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/chr11_3201.phylip = 1 - 870;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_3201.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 (GTRG)]

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_3201.phylip = 1 - 870;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr11_3201.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

