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

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
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
HKYIG	41	-1334.5378	2760.9438	0.5026	0.5026
HKYG	40	-1335.9643	2761.3000	0.4206	0.9232
HKYI	40	-1338.1240	2765.6195	0.0485	0.9718
GTRIG	45	-1333.5947	2769.1895	0.0081	0.9799
K2PG	37	-1343.6813	2769.3285	0.0076	0.9875
K2PIG	38	-1342.5879	2769.5963	0.0066	0.9941
GTRG	44	-1335.6384	2770.7220	0.0038	0.9979
GTRI	44	-1336.8973	2773.2398	0.0011	0.9990
K2PI	37	-1345.9273	2773.8206	0.0008	0.9998
SYMIG	42	-1341.4062	2777.1917	0.0001	0.9999
SYMG	41	-1344.3163	2780.5007	0.0000	1.0000
SYMI	41	-1344.5667	2781.0016	0.0000	1.0000
F81G	39	-1383.3566	2853.6021	0.0000	1.0000
F81IG	40	-1384.0104	2857.3922	0.0000	1.0000
F81I	39	-1385.5498	2857.9884	0.0000	1.0000
HKY	39	-1385.5753	2858.0395	0.0000	1.0000
JC69G	36	-1390.3287	2860.1829	0.0000	1.0000
K2P	36	-1391.0837	2861.6927	0.0000	1.0000
JC69IG	37	-1390.7125	2863.3911	0.0000	1.0000
GTR	43	-1383.2528	2863.4106	0.0000	1.0000
JC69I	36	-1392.3863	2864.2981	0.0000	1.0000
SYM	40	-1393.3989	2876.1693	0.0000	1.0000
F81	38	-1429.4456	2943.3116	0.0000	1.0000
JC69	35	-1437.1266	2951.3519	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	41	-1334.5378	2751.0756	0.5526	0.5526
HKYG	40	-1335.9643	2751.9286	0.3607	0.9133
HKYI	40	-1338.1240	2756.2481	0.0416	0.9549
GTRIG	45	-1333.5947	2757.1895	0.0260	0.9809
GTRG	44	-1335.6384	2759.2769	0.0092	0.9900
K2PIG	38	-1342.5879	2761.1758	0.0035	0.9936
K2PG	37	-1343.6813	2761.3625	0.0032	0.9968
GTRI	44	-1336.8973	2761.7947	0.0026	0.9994
K2PI	37	-1345.9273	2765.8546	0.0003	0.9997
SYMIG	42	-1341.4062	2766.8124	0.0002	0.9999
SYMG	41	-1344.3163	2770.6325	0.0000	1.0000
SYMI	41	-1344.5667	2771.1334	0.0000	1.0000
F81G	39	-1383.3566	2844.7132	0.0000	1.0000
F81IG	40	-1384.0104	2848.0208	0.0000	1.0000
F81I	39	-1385.5498	2849.0996	0.0000	1.0000
HKY	39	-1385.5753	2849.1506	0.0000	1.0000
GTR	43	-1383.2528	2852.5057	0.0000	1.0000
JC69G	36	-1390.3287	2852.6574	0.0000	1.0000
K2P	36	-1391.0837	2854.1673	0.0000	1.0000
JC69IG	37	-1390.7125	2855.4251	0.0000	1.0000
JC69I	36	-1392.3863	2856.7726	0.0000	1.0000
SYM	40	-1393.3989	2866.7979	0.0000	1.0000
F81	38	-1429.4456	2934.8911	0.0000	1.0000
JC69	35	-1437.1266	2944.2533	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1343.6813	2908.2047	0.6090	0.6090
HKYG	40	-1335.9643	2910.6769	0.1769	0.7860
K2PIG	38	-1342.5879	2911.9867	0.0919	0.8779
K2PI	37	-1345.9273	2912.6968	0.0644	0.9423
HKYIG	41	-1334.5378	2913.7926	0.0373	0.9796
HKYI	40	-1338.1240	2914.9964	0.0204	1.0000
SYMG	41	-1344.3163	2933.3495	0.0000	1.0000
SYMIG	42	-1341.4062	2933.4981	0.0000	1.0000
SYMI	41	-1344.5667	2933.8504	0.0000	1.0000
GTRG	44	-1335.6384	2933.9000	0.0000	1.0000
GTRIG	45	-1333.5947	2935.7813	0.0000	1.0000
GTRI	44	-1336.8973	2936.4178	0.0000	1.0000
JC69G	36	-1390.3287	2995.5309	0.0000	1.0000
K2P	36	-1391.0837	2997.0408	0.0000	1.0000
F81G	39	-1383.3566	2999.4928	0.0000	1.0000
JC69I	36	-1392.3863	2999.6461	0.0000	1.0000
JC69IG	37	-1390.7125	3002.2673	0.0000	1.0000
F81I	39	-1385.5498	3003.8792	0.0000	1.0000
HKY	39	-1385.5753	3003.9302	0.0000	1.0000
F81IG	40	-1384.0104	3006.7691	0.0000	1.0000
GTR	43	-1383.2528	3023.1601	0.0000	1.0000
SYM	40	-1393.3989	3025.5462	0.0000	1.0000
JC69	35	-1437.1266	3083.1581	0.0000	1.0000
F81	38	-1429.4456	3085.7020	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6695.phylip = 1 - 391;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6695.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYIG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6695.phylip = 1 - 391;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6695.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.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/chrZ_6695.phylip = 1 - 391;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6695.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

