
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:19 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_483.phylip" (19 taxa, 381 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	-1513.0434	3115.7340	0.5284	0.5284
HKYIG	41	-1512.4533	3117.0659	0.2715	0.7999
GTRG	44	-1509.3132	3118.4121	0.1385	0.9383
GTRIG	45	-1508.9038	3120.1657	0.0576	0.9960
HKYI	40	-1518.1614	3125.9699	0.0032	0.9991
GTRI	44	-1514.3741	3128.5340	0.0009	1.0000
K2PG	37	-1532.9352	3148.0687	0.0000	1.0000
K2PIG	38	-1532.0778	3148.8224	0.0000	1.0000
SYMG	41	-1531.5943	3155.3479	0.0000	1.0000
SYMIG	42	-1530.7593	3156.2049	0.0000	1.0000
K2PI	37	-1538.6535	3159.5054	0.0000	1.0000
SYMI	41	-1537.4944	3167.1481	0.0000	1.0000
F81G	39	-1542.5622	3172.2740	0.0000	1.0000
F81IG	40	-1541.9690	3173.5850	0.0000	1.0000
F81I	39	-1548.0026	3183.1548	0.0000	1.0000
HKY	39	-1551.2867	3189.7229	0.0000	1.0000
GTR	43	-1547.1347	3191.4980	0.0000	1.0000
JC69G	36	-1560.1117	3199.9675	0.0000	1.0000
JC69IG	37	-1559.5672	3201.3327	0.0000	1.0000
JC69I	36	-1566.4204	3212.5850	0.0000	1.0000
K2P	36	-1572.9015	3225.5472	0.0000	1.0000
SYM	40	-1571.4876	3232.6222	0.0000	1.0000
F81	38	-1581.1412	3246.9491	0.0000	1.0000
JC69	35	-1600.1129	3277.5301	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1513.0434	3106.0869	0.3495	0.3495
GTRG	44	-1509.3132	3106.6264	0.2669	0.6164
HKYIG	41	-1512.4533	3106.9066	0.2320	0.8484
GTRIG	45	-1508.9038	3107.8075	0.1479	0.9962
HKYI	40	-1518.1614	3116.3229	0.0021	0.9983
GTRI	44	-1514.3741	3116.7482	0.0017	1.0000
K2PG	37	-1532.9352	3139.8705	0.0000	1.0000
K2PIG	38	-1532.0778	3140.1557	0.0000	1.0000
SYMG	41	-1531.5943	3145.1886	0.0000	1.0000
SYMIG	42	-1530.7593	3145.5185	0.0000	1.0000
K2PI	37	-1538.6535	3151.3071	0.0000	1.0000
SYMI	41	-1537.4944	3156.9888	0.0000	1.0000
F81G	39	-1542.5622	3163.1244	0.0000	1.0000
F81IG	40	-1541.9690	3163.9379	0.0000	1.0000
F81I	39	-1548.0026	3174.0052	0.0000	1.0000
GTR	43	-1547.1347	3180.2695	0.0000	1.0000
HKY	39	-1551.2867	3180.5734	0.0000	1.0000
JC69G	36	-1560.1117	3192.2234	0.0000	1.0000
JC69IG	37	-1559.5672	3193.1344	0.0000	1.0000
JC69I	36	-1566.4204	3204.8409	0.0000	1.0000
K2P	36	-1572.9015	3217.8030	0.0000	1.0000
SYM	40	-1571.4876	3222.9752	0.0000	1.0000
F81	38	-1581.1412	3238.2824	0.0000	1.0000
JC69	35	-1600.1129	3270.2258	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1513.0434	3263.7989	0.9101	0.9101
HKYIG	41	-1512.4533	3268.5614	0.0841	0.9942
HKYI	40	-1518.1614	3274.0348	0.0054	0.9997
GTRG	44	-1509.3132	3280.1096	0.0003	1.0000
GTRIG	45	-1508.9038	3285.2335	0.0000	1.0000
K2PG	37	-1532.9352	3285.7541	0.0000	1.0000
K2PIG	38	-1532.0778	3289.9821	0.0000	1.0000
GTRI	44	-1514.3741	3290.2314	0.0000	1.0000
K2PI	37	-1538.6535	3297.1907	0.0000	1.0000
SYMG	41	-1531.5943	3306.8434	0.0000	1.0000
SYMIG	42	-1530.7593	3311.1161	0.0000	1.0000
F81G	39	-1542.5622	3316.8936	0.0000	1.0000
SYMI	41	-1537.4944	3318.6435	0.0000	1.0000
F81IG	40	-1541.9690	3321.6499	0.0000	1.0000
F81I	39	-1548.0026	3327.7744	0.0000	1.0000
JC69G	36	-1560.1117	3334.1641	0.0000	1.0000
HKY	39	-1551.2867	3334.3426	0.0000	1.0000
JC69IG	37	-1559.5672	3339.0180	0.0000	1.0000
JC69I	36	-1566.4204	3346.7816	0.0000	1.0000
GTR	43	-1547.1347	3349.8099	0.0000	1.0000
K2P	36	-1572.9015	3359.7438	0.0000	1.0000
SYM	40	-1571.4876	3380.6872	0.0000	1.0000
F81	38	-1581.1412	3388.1088	0.0000	1.0000
JC69	35	-1600.1129	3408.2238	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_483.phylip = 1 - 381;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_483.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_483.phylip = 1 - 381;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_483.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_483.phylip = 1 - 381;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_483.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

