
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_11572.phylip" (19 taxa, 859 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	44	-3831.1623	7755.1894	0.4343	0.4343
GTRIG	45	-3830.4403	7755.9729	0.2935	0.7279
HKYG	40	-3836.5748	7757.1594	0.1622	0.8901
HKYIG	41	-3835.8667	7757.9489	0.1093	0.9994
GTRI	44	-3837.6914	7768.2476	0.0006	1.0000
HKYI	40	-3849.2835	7782.5768	0.0000	1.0000
SYMG	41	-3873.8196	7833.8546	0.0000	1.0000
K2PG	37	-3881.3693	7840.1636	0.0000	1.0000
K2PIG	38	-3880.5702	7840.7551	0.0000	1.0000
SYMIG	42	-3879.0866	7846.5997	0.0000	1.0000
K2PI	37	-3885.5350	7848.4952	0.0000	1.0000
SYMI	41	-3884.0140	7854.2435	0.0000	1.0000
HKY	39	-3912.6040	7907.0174	0.0000	1.0000
F81G	39	-3913.8121	7909.4338	0.0000	1.0000
F81IG	40	-3912.8590	7909.7279	0.0000	1.0000
F81I	39	-3915.3594	7912.5283	0.0000	1.0000
GTR	43	-3916.9056	7924.4541	0.0000	1.0000
JC69G	36	-3942.6407	7960.5222	0.0000	1.0000
JC69IG	37	-3941.6852	7960.7955	0.0000	1.0000
JC69I	36	-3945.6573	7966.5555	0.0000	1.0000
SYM	40	-3950.5266	7985.0630	0.0000	1.0000
K2P	36	-3958.6294	7992.4997	0.0000	1.0000
F81	38	-3987.9095	8055.4337	0.0000	1.0000
JC69	35	-4016.6338	8106.3295	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-3831.1623	7750.3246	0.4578	0.4578
GTRIG	45	-3830.4403	7750.8807	0.3467	0.8046
HKYG	40	-3836.5748	7753.1496	0.1115	0.9161
HKYIG	41	-3835.8667	7753.7334	0.0833	0.9993
GTRI	44	-3837.6914	7763.3828	0.0007	1.0000
HKYI	40	-3849.2835	7778.5670	0.0000	1.0000
SYMG	41	-3873.8196	7829.6392	0.0000	1.0000
K2PG	37	-3881.3693	7836.7385	0.0000	1.0000
K2PIG	38	-3880.5702	7837.1404	0.0000	1.0000
SYMIG	42	-3879.0866	7842.1732	0.0000	1.0000
K2PI	37	-3885.5350	7845.0701	0.0000	1.0000
SYMI	41	-3884.0140	7850.0280	0.0000	1.0000
HKY	39	-3912.6040	7903.2079	0.0000	1.0000
F81G	39	-3913.8121	7905.6243	0.0000	1.0000
F81IG	40	-3912.8590	7905.7181	0.0000	1.0000
F81I	39	-3915.3594	7908.7188	0.0000	1.0000
GTR	43	-3916.9056	7919.8111	0.0000	1.0000
JC69G	36	-3942.6407	7957.2813	0.0000	1.0000
JC69IG	37	-3941.6852	7957.3704	0.0000	1.0000
JC69I	36	-3945.6573	7963.3147	0.0000	1.0000
SYM	40	-3950.5266	7981.0532	0.0000	1.0000
K2P	36	-3958.6294	7989.2588	0.0000	1.0000
F81	38	-3987.9095	8051.8191	0.0000	1.0000
JC69	35	-4016.6338	8103.2675	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-3836.5748	7943.3803	0.9349	0.9349
HKYIG	41	-3835.8667	7948.7200	0.0648	0.9997
GTRG	44	-3831.1623	7959.5784	0.0003	1.0000
GTRIG	45	-3830.4403	7964.8903	0.0000	1.0000
HKYI	40	-3849.2835	7968.7978	0.0000	1.0000
GTRI	44	-3837.6914	7972.6366	0.0000	1.0000
K2PG	37	-3881.3693	8012.7020	0.0000	1.0000
K2PIG	38	-3880.5702	8017.8597	0.0000	1.0000
K2PI	37	-3885.5350	8021.0335	0.0000	1.0000
SYMG	41	-3873.8196	8024.6257	0.0000	1.0000
SYMIG	42	-3879.0866	8041.9155	0.0000	1.0000
SYMI	41	-3884.0140	8045.0146	0.0000	1.0000
HKY	39	-3912.6040	8088.6829	0.0000	1.0000
F81G	39	-3913.8121	8091.0993	0.0000	1.0000
F81I	39	-3915.3594	8094.1938	0.0000	1.0000
F81IG	40	-3912.8590	8095.9488	0.0000	1.0000
GTR	43	-3916.9056	8124.3092	0.0000	1.0000
JC69G	36	-3942.6407	8128.4890	0.0000	1.0000
JC69IG	37	-3941.6852	8133.3338	0.0000	1.0000
JC69I	36	-3945.6573	8134.5223	0.0000	1.0000
K2P	36	-3958.6294	8160.4665	0.0000	1.0000
SYM	40	-3950.5266	8171.2840	0.0000	1.0000
F81	38	-3987.9095	8232.5383	0.0000	1.0000
JC69	35	-4016.6338	8269.7195	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/chrZ_11572.phylip = 1 - 859;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_11572.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/chrZ_11572.phylip = 1 - 859;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_11572.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 (HKYG)]

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

