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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	-1537.3035	3172.8660	0.8520	0.8520
GTRIG	45	-1537.8116	3176.3764	0.1473	0.9993
HKYG	40	-1550.3086	3189.0274	0.0003	0.9996
GTRI	44	-1545.4580	3189.1751	0.0002	0.9998
HKYIG	41	-1549.5283	3189.9101	0.0002	1.0000
HKYI	40	-1554.9886	3198.3874	0.0000	1.0000
SYMG	41	-1559.0708	3208.9951	0.0000	1.0000
SYMIG	42	-1558.2745	3209.8583	0.0000	1.0000
K2PG	37	-1566.0770	3213.3093	0.0000	1.0000
K2PIG	38	-1565.1992	3213.9596	0.0000	1.0000
SYMI	41	-1563.6065	3218.0666	0.0000	1.0000
K2PI	37	-1569.2692	3219.6936	0.0000	1.0000
F81G	39	-1578.5303	3243.0402	0.0000	1.0000
F81IG	40	-1578.0966	3244.6035	0.0000	1.0000
F81I	39	-1581.3230	3248.6256	0.0000	1.0000
GTR	43	-1578.0694	3251.9165	0.0000	1.0000
HKY	39	-1589.7422	3265.4640	0.0000	1.0000
JC69G	36	-1597.8394	3274.4403	0.0000	1.0000
JC69IG	37	-1598.1762	3277.5077	0.0000	1.0000
JC69I	36	-1600.7333	3280.2281	0.0000	1.0000
SYM	40	-1599.5174	3287.4451	0.0000	1.0000
K2P	36	-1605.1659	3289.0931	0.0000	1.0000
F81	38	-1619.4660	3322.4933	0.0000	1.0000
JC69	35	-1637.0113	3350.4023	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1537.3035	3162.6070	0.8184	0.8184
GTRIG	45	-1537.8116	3165.6231	0.1811	0.9996
GTRI	44	-1545.4580	3178.9160	0.0002	0.9998
HKYG	40	-1550.3086	3180.6172	0.0001	0.9999
HKYIG	41	-1549.5283	3181.0566	0.0001	1.0000
HKYI	40	-1554.9886	3189.9772	0.0000	1.0000
SYMG	41	-1559.0708	3200.1416	0.0000	1.0000
SYMIG	42	-1558.2745	3200.5490	0.0000	1.0000
K2PG	37	-1566.0770	3206.1541	0.0000	1.0000
K2PIG	38	-1565.1992	3206.3983	0.0000	1.0000
SYMI	41	-1563.6065	3209.2131	0.0000	1.0000
K2PI	37	-1569.2692	3212.5383	0.0000	1.0000
F81G	39	-1578.5303	3235.0606	0.0000	1.0000
F81IG	40	-1578.0966	3236.1933	0.0000	1.0000
F81I	39	-1581.3230	3240.6461	0.0000	1.0000
GTR	43	-1578.0694	3242.1387	0.0000	1.0000
HKY	39	-1589.7422	3257.4845	0.0000	1.0000
JC69G	36	-1597.8394	3267.6789	0.0000	1.0000
JC69IG	37	-1598.1762	3270.3525	0.0000	1.0000
JC69I	36	-1600.7333	3273.4666	0.0000	1.0000
SYM	40	-1599.5174	3279.0349	0.0000	1.0000
K2P	36	-1605.1659	3282.3317	0.0000	1.0000
F81	38	-1619.4660	3314.9320	0.0000	1.0000
JC69	35	-1637.0113	3344.0225	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-1537.3035	3341.5157	0.6687	0.6687
HKYG	40	-1550.3086	3343.2615	0.2794	0.9481
HKYIG	41	-1549.5283	3347.7670	0.0294	0.9774
GTRIG	45	-1537.8116	3348.5980	0.0194	0.9968
HKYI	40	-1554.9886	3352.6215	0.0026	0.9994
K2PG	37	-1566.0770	3356.6001	0.0004	0.9997
GTRI	44	-1545.4580	3357.8248	0.0002	0.9999
K2PIG	38	-1565.1992	3360.9104	0.0000	1.0000
K2PI	37	-1569.2692	3362.9843	0.0000	1.0000
SYMG	41	-1559.0708	3366.8521	0.0000	1.0000
SYMIG	42	-1558.2745	3371.3255	0.0000	1.0000
SYMI	41	-1563.6065	3375.9235	0.0000	1.0000
F81G	39	-1578.5303	3393.6388	0.0000	1.0000
F81IG	40	-1578.0966	3398.8376	0.0000	1.0000
F81I	39	-1581.3230	3399.2243	0.0000	1.0000
JC69G	36	-1597.8394	3414.0588	0.0000	1.0000
HKY	39	-1589.7422	3416.0627	0.0000	1.0000
GTR	43	-1578.0694	3416.9814	0.0000	1.0000
JC69I	36	-1600.7333	3419.8465	0.0000	1.0000
JC69IG	37	-1598.1762	3420.7985	0.0000	1.0000
K2P	36	-1605.1659	3428.7116	0.0000	1.0000
SYM	40	-1599.5174	3441.6792	0.0000	1.0000
F81	38	-1619.4660	3469.4441	0.0000	1.0000
JC69	35	-1637.0113	3486.3363	0.0000	1.0000



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

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

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

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

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


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

