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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	45	-4384.8993	8863.4430	0.7289	0.7289
GTRG	44	-4386.9694	8865.4216	0.2711	1.0000
HKYIG	41	-4405.2615	8895.5441	0.0000	1.0000
HKYG	40	-4407.8554	8898.5854	0.0000	1.0000
GTRI	44	-4406.6991	8904.8811	0.0000	1.0000
HKYI	40	-4423.6764	8930.2275	0.0000	1.0000
SYMIG	42	-4437.9100	8962.9911	0.0000	1.0000
SYMG	41	-4440.4035	8965.8280	0.0000	1.0000
K2PIG	38	-4451.4392	8981.4716	0.0000	1.0000
K2PG	37	-4454.2712	8985.0005	0.0000	1.0000
SYMI	41	-4456.5283	8998.0776	0.0000	1.0000
K2PI	37	-4468.9784	9014.4149	0.0000	1.0000
GTR	43	-4547.7901	9184.9053	0.0000	1.0000
F81G	39	-4558.7055	9198.1430	0.0000	1.0000
F81IG	40	-4563.7534	9210.3815	0.0000	1.0000
HKY	39	-4569.9362	9220.6044	0.0000	1.0000
F81I	39	-4571.1050	9222.9420	0.0000	1.0000
JC69G	36	-4594.6647	9263.6560	0.0000	1.0000
SYM	40	-4594.6443	9272.1633	0.0000	1.0000
JC69IG	37	-4599.6435	9275.7450	0.0000	1.0000
JC69I	36	-4610.1094	9294.5455	0.0000	1.0000
K2P	36	-4613.4346	9301.1959	0.0000	1.0000
F81	38	-4710.5421	9499.6774	0.0000	1.0000
JC69	35	-4744.1613	9560.5216	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-4384.8993	8859.7987	0.7446	0.7446
GTRG	44	-4386.9694	8861.9388	0.2554	1.0000
HKYIG	41	-4405.2615	8892.5230	0.0000	1.0000
HKYG	40	-4407.8554	8895.7107	0.0000	1.0000
GTRI	44	-4406.6991	8901.3982	0.0000	1.0000
HKYI	40	-4423.6764	8927.3528	0.0000	1.0000
SYMIG	42	-4437.9100	8959.8199	0.0000	1.0000
SYMG	41	-4440.4035	8962.8070	0.0000	1.0000
K2PIG	38	-4451.4392	8978.8784	0.0000	1.0000
K2PG	37	-4454.2712	8982.5424	0.0000	1.0000
SYMI	41	-4456.5283	8995.0565	0.0000	1.0000
K2PI	37	-4468.9784	9011.9569	0.0000	1.0000
GTR	43	-4547.7901	9181.5802	0.0000	1.0000
F81G	39	-4558.7055	9195.4110	0.0000	1.0000
F81IG	40	-4563.7534	9207.5069	0.0000	1.0000
HKY	39	-4569.9362	9217.8723	0.0000	1.0000
F81I	39	-4571.1050	9220.2100	0.0000	1.0000
JC69G	36	-4594.6647	9261.3294	0.0000	1.0000
SYM	40	-4594.6443	9269.2886	0.0000	1.0000
JC69IG	37	-4599.6435	9273.2870	0.0000	1.0000
JC69I	36	-4610.1094	9292.2188	0.0000	1.0000
K2P	36	-4613.4346	9298.8692	0.0000	1.0000
F81	38	-4710.5421	9497.0842	0.0000	1.0000
JC69	35	-4744.1613	9558.3227	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-4386.9694	9085.2372	0.8116	0.8116
GTRIG	45	-4384.8993	9088.1720	0.1871	0.9987
HKYG	40	-4407.8554	9098.7092	0.0010	0.9996
HKYIG	41	-4405.2615	9100.5965	0.0004	1.0000
GTRI	44	-4406.6991	9124.6966	0.0000	1.0000
HKYI	40	-4423.6764	9130.3513	0.0000	1.0000
K2PG	37	-4454.2712	9170.3161	0.0000	1.0000
SYMG	41	-4440.4035	9170.8805	0.0000	1.0000
K2PIG	38	-4451.4392	9171.7270	0.0000	1.0000
SYMIG	42	-4437.9100	9172.9684	0.0000	1.0000
K2PI	37	-4468.9784	9199.7305	0.0000	1.0000
SYMI	41	-4456.5283	9203.1300	0.0000	1.0000
F81G	39	-4558.7055	9393.3345	0.0000	1.0000
GTR	43	-4547.7901	9399.8036	0.0000	1.0000
F81IG	40	-4563.7534	9410.5054	0.0000	1.0000
HKY	39	-4569.9362	9415.7959	0.0000	1.0000
F81I	39	-4571.1050	9418.1336	0.0000	1.0000
JC69G	36	-4594.6647	9444.0281	0.0000	1.0000
JC69IG	37	-4599.6435	9461.0606	0.0000	1.0000
SYM	40	-4594.6443	9472.2871	0.0000	1.0000
JC69I	36	-4610.1094	9474.9175	0.0000	1.0000
K2P	36	-4613.4346	9481.5679	0.0000	1.0000
F81	38	-4710.5421	9689.9328	0.0000	1.0000
JC69	35	-4744.1613	9735.9464	0.0000	1.0000



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

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

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


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

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

