
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:07 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/chr7_11484.phylip" (19 taxa, 611 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	-1509.8941	3105.5425	0.6131	0.6131
HKYIG	41	-1509.4956	3107.0438	0.2894	0.9026
HKYI	40	-1512.4900	3110.7343	0.0457	0.9483
GTRG	44	-1508.1626	3111.3217	0.0341	0.9824
GTRIG	45	-1507.8097	3112.9468	0.0151	0.9975
GTRI	44	-1510.7753	3116.5471	0.0025	1.0000
HKY	39	-1535.6645	3154.7931	0.0000	1.0000
GTR	43	-1532.9653	3158.6043	0.0000	1.0000
F81G	39	-1540.0318	3163.5277	0.0000	1.0000
F81I	39	-1540.4646	3164.3934	0.0000	1.0000
F81IG	40	-1539.7501	3165.2546	0.0000	1.0000
K2PG	37	-1545.8765	3170.6606	0.0000	1.0000
SYMG	41	-1542.2138	3172.4804	0.0000	1.0000
SYMIG	42	-1541.7163	3173.7918	0.0000	1.0000
K2PI	37	-1548.9658	3176.8392	0.0000	1.0000
K2PIG	38	-1547.9106	3177.0031	0.0000	1.0000
SYMI	41	-1545.3843	3178.8214	0.0000	1.0000
F81	38	-1561.9149	3205.0117	0.0000	1.0000
SYM	40	-1562.6859	3211.1262	0.0000	1.0000
K2P	36	-1567.7878	3212.2168	0.0000	1.0000
JC69G	36	-1570.1676	3216.9763	0.0000	1.0000
JC69I	36	-1570.7474	3218.1360	0.0000	1.0000
JC69IG	37	-1569.8351	3218.5777	0.0000	1.0000
JC69	35	-1590.6260	3255.6345	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1509.8941	3099.7882	0.5593	0.5593
HKYIG	41	-1509.4956	3100.9911	0.3065	0.8659
GTRG	44	-1508.1626	3104.3252	0.0579	0.9237
HKYI	40	-1512.4900	3104.9800	0.0417	0.9655
GTRIG	45	-1507.8097	3105.6193	0.0303	0.9958
GTRI	44	-1510.7753	3109.5506	0.0042	1.0000
HKY	39	-1535.6645	3149.3290	0.0000	1.0000
GTR	43	-1532.9653	3151.9306	0.0000	1.0000
F81G	39	-1540.0318	3158.0636	0.0000	1.0000
F81I	39	-1540.4646	3158.9293	0.0000	1.0000
F81IG	40	-1539.7501	3159.5002	0.0000	1.0000
K2PG	37	-1545.8765	3165.7531	0.0000	1.0000
SYMG	41	-1542.2138	3166.4276	0.0000	1.0000
SYMIG	42	-1541.7163	3167.4326	0.0000	1.0000
K2PIG	38	-1547.9106	3171.8213	0.0000	1.0000
K2PI	37	-1548.9658	3171.9317	0.0000	1.0000
SYMI	41	-1545.3843	3172.7686	0.0000	1.0000
F81	38	-1561.9149	3199.8299	0.0000	1.0000
SYM	40	-1562.6859	3205.3718	0.0000	1.0000
K2P	36	-1567.7878	3207.5756	0.0000	1.0000
JC69G	36	-1570.1676	3212.3351	0.0000	1.0000
JC69I	36	-1570.7474	3213.4948	0.0000	1.0000
JC69IG	37	-1569.8351	3213.6702	0.0000	1.0000
JC69	35	-1590.6260	3251.2519	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1509.8941	3276.3920	0.8812	0.8812
HKYI	40	-1512.4900	3281.5838	0.0657	0.9469
HKYIG	41	-1509.4956	3282.0101	0.0531	1.0000
GTRG	44	-1508.1626	3298.5895	0.0000	1.0000
GTRI	44	-1510.7753	3303.8149	0.0000	1.0000
GTRIG	45	-1507.8097	3304.2987	0.0000	1.0000
HKY	39	-1535.6645	3321.5178	0.0000	1.0000
K2PG	37	-1545.8765	3329.1116	0.0000	1.0000
F81G	39	-1540.0318	3330.2524	0.0000	1.0000
F81I	39	-1540.4646	3331.1181	0.0000	1.0000
K2PI	37	-1548.9658	3335.2903	0.0000	1.0000
F81IG	40	-1539.7501	3336.1041	0.0000	1.0000
K2PIG	38	-1547.9106	3339.5949	0.0000	1.0000
GTR	43	-1532.9653	3341.7797	0.0000	1.0000
SYMG	41	-1542.2138	3347.4466	0.0000	1.0000
SYMIG	42	-1541.7163	3352.8667	0.0000	1.0000
SYMI	41	-1545.3843	3353.7876	0.0000	1.0000
K2P	36	-1567.7878	3366.5191	0.0000	1.0000
F81	38	-1561.9149	3367.6035	0.0000	1.0000
JC69G	36	-1570.1676	3371.2786	0.0000	1.0000
JC69I	36	-1570.7474	3372.4383	0.0000	1.0000
JC69IG	37	-1569.8351	3377.0288	0.0000	1.0000
SYM	40	-1562.6859	3381.9757	0.0000	1.0000
JC69	35	-1590.6260	3405.7803	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/chr7_11484.phylip = 1 - 611;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_11484.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/chr7_11484.phylip = 1 - 611;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_11484.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/chr7_11484.phylip = 1 - 611;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_11484.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

