
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:27 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/chrZ_6036.phylip" (19 taxa, 476 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	-1545.3464	3178.2330	0.5648	0.5648
HKYIG	41	-1544.5245	3178.9845	0.3879	0.9527
HKYI	40	-1548.6788	3184.8977	0.0202	0.9729
GTRG	44	-1544.1140	3185.4160	0.0156	0.9885
GTRIG	45	-1543.2935	3186.2150	0.0104	0.9989
GTRI	44	-1546.7772	3190.7424	0.0011	1.0000
K2PG	37	-1561.3520	3203.1242	0.0000	1.0000
K2PIG	38	-1560.6169	3204.0165	0.0000	1.0000
K2PI	37	-1563.3793	3207.1788	0.0000	1.0000
SYMG	41	-1558.8510	3207.6374	0.0000	1.0000
SYMIG	42	-1558.0736	3208.4889	0.0000	1.0000
SYMI	41	-1560.7939	3211.5234	0.0000	1.0000
HKY	39	-1584.3342	3253.8244	0.0000	1.0000
GTR	43	-1581.8547	3258.4687	0.0000	1.0000
F81G	39	-1595.0074	3275.1707	0.0000	1.0000
F81IG	40	-1594.5347	3276.6096	0.0000	1.0000
K2P	36	-1600.0239	3278.1162	0.0000	1.0000
F81I	39	-1598.2831	3281.7222	0.0000	1.0000
SYM	40	-1598.5316	3284.6035	0.0000	1.0000
JC69G	36	-1608.2946	3294.6575	0.0000	1.0000
JC69IG	37	-1607.8519	3296.1239	0.0000	1.0000
JC69I	36	-1609.9426	3297.9535	0.0000	1.0000
F81	38	-1633.2266	3349.2358	0.0000	1.0000
JC69	35	-1644.7676	3365.2624	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1545.3464	3170.6928	0.5019	0.5019
HKYIG	41	-1544.5245	3171.0490	0.4201	0.9220
GTRG	44	-1544.1140	3176.2281	0.0315	0.9535
GTRIG	45	-1543.2935	3176.5871	0.0263	0.9799
HKYI	40	-1548.6788	3177.3575	0.0179	0.9978
GTRI	44	-1546.7772	3181.5545	0.0022	1.0000
K2PG	37	-1561.3520	3196.7041	0.0000	1.0000
K2PIG	38	-1560.6169	3197.2339	0.0000	1.0000
SYMG	41	-1558.8510	3199.7020	0.0000	1.0000
SYMIG	42	-1558.0736	3200.1471	0.0000	1.0000
K2PI	37	-1563.3793	3200.7587	0.0000	1.0000
SYMI	41	-1560.7939	3203.5879	0.0000	1.0000
HKY	39	-1584.3342	3246.6685	0.0000	1.0000
GTR	43	-1581.8547	3249.7094	0.0000	1.0000
F81G	39	-1595.0074	3268.0147	0.0000	1.0000
F81IG	40	-1594.5347	3269.0694	0.0000	1.0000
K2P	36	-1600.0239	3272.0479	0.0000	1.0000
F81I	39	-1598.2831	3274.5662	0.0000	1.0000
SYM	40	-1598.5316	3277.0633	0.0000	1.0000
JC69G	36	-1608.2946	3288.5892	0.0000	1.0000
JC69IG	37	-1607.8519	3289.7038	0.0000	1.0000
JC69I	36	-1609.9426	3291.8852	0.0000	1.0000
F81	38	-1633.2266	3342.4532	0.0000	1.0000
JC69	35	-1644.7676	3359.5351	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1545.3464	3337.3095	0.8761	0.8761
HKYIG	41	-1544.5245	3341.8311	0.0913	0.9675
HKYI	40	-1548.6788	3343.9742	0.0313	0.9987
K2PG	37	-1561.3520	3350.8245	0.0010	0.9998
K2PI	37	-1563.3793	3354.8792	0.0001	0.9999
K2PIG	38	-1560.6169	3355.5198	0.0001	1.0000
GTRG	44	-1544.1140	3359.5065	0.0000	1.0000
GTRIG	45	-1543.2935	3364.0309	0.0000	1.0000
GTRI	44	-1546.7772	3364.8328	0.0000	1.0000
SYMG	41	-1558.8510	3370.4841	0.0000	1.0000
SYMI	41	-1560.7939	3374.3700	0.0000	1.0000
SYMIG	42	-1558.0736	3375.0947	0.0000	1.0000
HKY	39	-1584.3342	3409.1198	0.0000	1.0000
K2P	36	-1600.0239	3422.0029	0.0000	1.0000
GTR	43	-1581.8547	3428.8224	0.0000	1.0000
F81G	39	-1595.0074	3430.4660	0.0000	1.0000
F81IG	40	-1594.5347	3435.6861	0.0000	1.0000
F81I	39	-1598.2831	3437.0175	0.0000	1.0000
JC69G	36	-1608.2946	3438.5442	0.0000	1.0000
JC69I	36	-1609.9426	3441.8402	0.0000	1.0000
SYM	40	-1598.5316	3443.6800	0.0000	1.0000
JC69IG	37	-1607.8519	3443.8242	0.0000	1.0000
F81	38	-1633.2266	3500.7391	0.0000	1.0000
JC69	35	-1644.7676	3505.3247	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/chrZ_6036.phylip = 1 - 476;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6036.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/chrZ_6036.phylip = 1 - 476;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6036.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/chrZ_6036.phylip = 1 - 476;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6036.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

