
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:57 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/chr7_2589.phylip" (19 taxa, 978 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	-4470.9632	9025.4269	0.5774	0.5774
HKYIG	41	-4470.7726	9027.2247	0.2350	0.8124
GTRG	44	-4468.0461	9028.3365	0.1348	0.9472
GTRIG	45	-4467.8846	9030.2113	0.0528	1.0000
HKYI	40	-4481.0071	9045.5148	0.0000	1.0000
GTRI	44	-4478.2195	9048.6834	0.0000	1.0000
F81G	39	-4556.6179	9194.5619	0.0000	1.0000
F81IG	40	-4556.4601	9196.4208	0.0000	1.0000
F81I	39	-4566.4145	9214.1552	0.0000	1.0000
SYMG	41	-4566.6114	9218.9023	0.0000	1.0000
SYMIG	42	-4566.4559	9220.7749	0.0000	1.0000
HKY	39	-4571.3028	9223.9318	0.0000	1.0000
GTR	43	-4568.7481	9227.5476	0.0000	1.0000
K2PG	37	-4577.9684	9232.9283	0.0000	1.0000
K2PIG	38	-4577.7981	9234.7526	0.0000	1.0000
SYMI	41	-4579.6483	9244.9762	0.0000	1.0000
K2PI	37	-4589.4554	9255.9023	0.0000	1.0000
JC69G	36	-4639.4952	9353.8215	0.0000	1.0000
JC69IG	37	-4639.3358	9355.6630	0.0000	1.0000
JC69I	36	-4650.2801	9375.3912	0.0000	1.0000
F81	38	-4648.1574	9375.4714	0.0000	1.0000
SYM	40	-4662.6653	9408.8311	0.0000	1.0000
K2P	36	-4676.0608	9426.9527	0.0000	1.0000
JC69	35	-4729.1465	9530.9682	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-4470.9632	9021.9263	0.5181	0.5181
HKYIG	41	-4470.7726	9023.5452	0.2306	0.7487
GTRG	44	-4468.0461	9024.0921	0.1754	0.9241
GTRIG	45	-4467.8846	9025.7692	0.0758	1.0000
HKYI	40	-4481.0071	9042.0142	0.0000	1.0000
GTRI	44	-4478.2195	9044.4390	0.0000	1.0000
F81G	39	-4556.6179	9191.2357	0.0000	1.0000
F81IG	40	-4556.4601	9192.9203	0.0000	1.0000
F81I	39	-4566.4145	9210.8290	0.0000	1.0000
SYMG	41	-4566.6114	9215.2228	0.0000	1.0000
SYMIG	42	-4566.4559	9216.9118	0.0000	1.0000
HKY	39	-4571.3028	9220.6056	0.0000	1.0000
GTR	43	-4568.7481	9223.4962	0.0000	1.0000
K2PG	37	-4577.9684	9229.9368	0.0000	1.0000
K2PIG	38	-4577.7981	9231.5961	0.0000	1.0000
SYMI	41	-4579.6483	9241.2967	0.0000	1.0000
K2PI	37	-4589.4554	9252.9108	0.0000	1.0000
JC69G	36	-4639.4952	9350.9905	0.0000	1.0000
JC69IG	37	-4639.3358	9352.6715	0.0000	1.0000
F81	38	-4648.1574	9372.3149	0.0000	1.0000
JC69I	36	-4650.2801	9372.5602	0.0000	1.0000
SYM	40	-4662.6653	9405.3305	0.0000	1.0000
K2P	36	-4676.0608	9424.1217	0.0000	1.0000
JC69	35	-4729.1465	9528.2930	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-4470.9632	9217.3467	0.9627	0.9627
HKYIG	41	-4470.7726	9223.8511	0.0372	0.9999
HKYI	40	-4481.0071	9237.4346	0.0000	1.0000
GTRG	44	-4468.0461	9239.0545	0.0000	1.0000
GTRIG	45	-4467.8846	9245.6172	0.0000	1.0000
GTRI	44	-4478.2195	9259.4014	0.0000	1.0000
F81G	39	-4556.6179	9381.7706	0.0000	1.0000
F81IG	40	-4556.4601	9388.3407	0.0000	1.0000
F81I	39	-4566.4145	9401.3639	0.0000	1.0000
K2PG	37	-4577.9684	9410.7007	0.0000	1.0000
HKY	39	-4571.3028	9411.1405	0.0000	1.0000
SYMG	41	-4566.6114	9415.5287	0.0000	1.0000
K2PIG	38	-4577.7981	9417.2455	0.0000	1.0000
SYMIG	42	-4566.4559	9422.1032	0.0000	1.0000
GTR	43	-4568.7481	9433.5731	0.0000	1.0000
K2PI	37	-4589.4554	9433.6747	0.0000	1.0000
SYMI	41	-4579.6483	9441.6026	0.0000	1.0000
JC69G	36	-4639.4952	9526.8688	0.0000	1.0000
JC69IG	37	-4639.3358	9533.4354	0.0000	1.0000
JC69I	36	-4650.2801	9548.4385	0.0000	1.0000
F81	38	-4648.1574	9557.9642	0.0000	1.0000
K2P	36	-4676.0608	9600.0000	0.0000	1.0000
SYM	40	-4662.6653	9600.7509	0.0000	1.0000
JC69	35	-4729.1465	9699.2859	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/chr7_2589.phylip = 1 - 978;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_2589.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/chr7_2589.phylip = 1 - 978;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_2589.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/chr7_2589.phylip = 1 - 978;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_2589.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

