
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:04 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/chr2_21358.phylip" (19 taxa, 1067 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	-4852.4553	9796.7854	0.7146	0.7146
GTRIG	45	-4852.2833	9798.6215	0.2853	1.0000
HKYG	40	-4867.1570	9817.5109	0.0000	1.0000
HKYIG	41	-4866.9750	9819.3100	0.0000	1.0000
GTRI	44	-4865.2902	9822.4551	0.0000	1.0000
HKYI	40	-4877.7700	9838.7369	0.0000	1.0000
SYMIG	42	-4899.1126	9885.7526	0.0000	1.0000
SYMG	41	-4900.2242	9885.8084	0.0000	1.0000
SYMI	41	-4905.5360	9896.4320	0.0000	1.0000
K2PG	37	-4911.5777	9899.8881	0.0000	1.0000
K2PIG	38	-4911.2566	9901.3964	0.0000	1.0000
K2PI	37	-4919.3902	9915.5131	0.0000	1.0000
GTR	43	-4973.9747	10037.6483	0.0000	1.0000
HKY	39	-4987.2015	10055.4411	0.0000	1.0000
F81IG	40	-4995.0001	10073.1972	0.0000	1.0000
F81G	39	-4997.0738	10075.1856	0.0000	1.0000
F81I	39	-5003.2229	10087.4839	0.0000	1.0000
SYM	40	-5016.4622	10116.1212	0.0000	1.0000
K2P	36	-5033.1744	10140.9352	0.0000	1.0000
JC69G	36	-5034.0721	10142.7306	0.0000	1.0000
JC69IG	37	-5033.4463	10143.6253	0.0000	1.0000
JC69I	36	-5039.3540	10153.2944	0.0000	1.0000
F81	38	-5113.8291	10306.5415	0.0000	1.0000
JC69	35	-5148.0597	10368.5636	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-4852.4553	9792.9106	0.6959	0.6959
GTRIG	45	-4852.2833	9794.5667	0.3041	1.0000
HKYG	40	-4867.1570	9814.3140	0.0000	1.0000
HKYIG	41	-4866.9750	9815.9500	0.0000	1.0000
GTRI	44	-4865.2902	9818.5803	0.0000	1.0000
HKYI	40	-4877.7700	9835.5400	0.0000	1.0000
SYMIG	42	-4899.1126	9882.2253	0.0000	1.0000
SYMG	41	-4900.2242	9882.4484	0.0000	1.0000
SYMI	41	-4905.5360	9893.0720	0.0000	1.0000
K2PG	37	-4911.5777	9897.1553	0.0000	1.0000
K2PIG	38	-4911.2566	9898.5131	0.0000	1.0000
K2PI	37	-4919.3902	9912.7804	0.0000	1.0000
GTR	43	-4973.9747	10033.9494	0.0000	1.0000
HKY	39	-4987.2015	10052.4031	0.0000	1.0000
F81IG	40	-4995.0001	10070.0003	0.0000	1.0000
F81G	39	-4997.0738	10072.1476	0.0000	1.0000
F81I	39	-5003.2229	10084.4459	0.0000	1.0000
SYM	40	-5016.4622	10112.9243	0.0000	1.0000
K2P	36	-5033.1744	10138.3488	0.0000	1.0000
JC69G	36	-5034.0721	10140.1442	0.0000	1.0000
JC69IG	37	-5033.4463	10140.8926	0.0000	1.0000
JC69I	36	-5039.3540	10150.7080	0.0000	1.0000
F81	38	-5113.8291	10303.6583	0.0000	1.0000
JC69	35	-5148.0597	10366.1194	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-4852.4553	10011.7053	0.6566	0.6566
HKYG	40	-4867.1570	10013.2183	0.3082	0.9648
GTRIG	45	-4852.2833	10018.3340	0.0239	0.9887
HKYIG	41	-4866.9750	10019.8269	0.0113	1.0000
HKYI	40	-4877.7700	10034.4443	0.0000	1.0000
GTRI	44	-4865.2902	10037.3750	0.0000	1.0000
K2PG	37	-4911.5777	10081.1418	0.0000	1.0000
SYMG	41	-4900.2242	10086.3252	0.0000	1.0000
K2PIG	38	-4911.2566	10087.4722	0.0000	1.0000
SYMIG	42	-4899.1126	10091.0747	0.0000	1.0000
K2PI	37	-4919.3902	10096.7668	0.0000	1.0000
SYMI	41	-4905.5360	10096.9489	0.0000	1.0000
HKY	39	-4987.2015	10246.3347	0.0000	1.0000
GTR	43	-4973.9747	10247.7714	0.0000	1.0000
F81G	39	-4997.0738	10266.0792	0.0000	1.0000
F81IG	40	-4995.0001	10268.9046	0.0000	1.0000
F81I	39	-5003.2229	10278.3775	0.0000	1.0000
SYM	40	-5016.4622	10311.8286	0.0000	1.0000
K2P	36	-5033.1744	10317.3627	0.0000	1.0000
JC69G	36	-5034.0721	10319.1580	0.0000	1.0000
JC69IG	37	-5033.4463	10324.8790	0.0000	1.0000
JC69I	36	-5039.3540	10329.7218	0.0000	1.0000
F81	38	-5113.8291	10492.6173	0.0000	1.0000
JC69	35	-5148.0597	10540.1606	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/chr2_21358.phylip = 1 - 1067;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21358.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/chr2_21358.phylip = 1 - 1067;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21358.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/chr2_21358.phylip = 1 - 1067;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_21358.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

