
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:02 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_23690.phylip" (19 taxa, 947 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	45	-4379.4071	8853.4090	0.9727	0.9727
SYMIG	42	-4386.2845	8860.5645	0.0272	0.9999
GTRG	44	-4390.1721	8872.7344	0.0001	1.0000
HKYG	40	-4396.2845	8876.1893	0.0000	1.0000
K2PG	37	-4399.7348	8876.5631	0.0000	1.0000
HKYIG	41	-4395.9804	8877.7663	0.0000	1.0000
K2PIG	38	-4399.4369	8878.1381	0.0000	1.0000
SYMG	41	-4396.9723	8879.7501	0.0000	1.0000
GTRI	44	-4401.6403	8895.6709	0.0000	1.0000
K2PI	37	-4409.8815	8896.8565	0.0000	1.0000
SYMI	41	-4407.6510	8901.1076	0.0000	1.0000
HKYI	40	-4416.8242	8917.2686	0.0000	1.0000
K2P	36	-4500.5644	9076.0562	0.0000	1.0000
GTR	43	-4494.7314	9079.6533	0.0000	1.0000
SYM	40	-4503.7906	9091.2015	0.0000	1.0000
HKY	39	-4512.1914	9105.8227	0.0000	1.0000
JC69IG	37	-4575.1512	9227.3960	0.0000	1.0000
F81IG	40	-4571.9735	9227.5672	0.0000	1.0000
JC69G	36	-4586.0452	9247.0179	0.0000	1.0000
F81G	39	-4583.1347	9247.7093	0.0000	1.0000
JC69I	36	-4594.5424	9264.0123	0.0000	1.0000
F81I	39	-4591.9409	9265.3217	0.0000	1.0000
F81	38	-4680.6424	9440.5492	0.0000	1.0000
JC69	35	-4684.4560	9441.6782	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-4379.4071	8848.8141	0.9796	0.9796
SYMIG	42	-4386.2845	8856.5690	0.0203	0.9999
GTRG	44	-4390.1721	8868.3442	0.0001	1.0000
HKYG	40	-4396.2845	8872.5690	0.0000	1.0000
K2PG	37	-4399.7348	8873.4696	0.0000	1.0000
HKYIG	41	-4395.9804	8873.9608	0.0000	1.0000
K2PIG	38	-4399.4369	8874.8738	0.0000	1.0000
SYMG	41	-4396.9723	8875.9446	0.0000	1.0000
GTRI	44	-4401.6403	8891.2807	0.0000	1.0000
K2PI	37	-4409.8815	8893.7630	0.0000	1.0000
SYMI	41	-4407.6510	8897.3021	0.0000	1.0000
HKYI	40	-4416.8242	8913.6483	0.0000	1.0000
K2P	36	-4500.5644	9073.1287	0.0000	1.0000
GTR	43	-4494.7314	9075.4628	0.0000	1.0000
SYM	40	-4503.7906	9087.5812	0.0000	1.0000
HKY	39	-4512.1914	9102.3828	0.0000	1.0000
F81IG	40	-4571.9735	9223.9469	0.0000	1.0000
JC69IG	37	-4575.1512	9224.3025	0.0000	1.0000
JC69G	36	-4586.0452	9244.0905	0.0000	1.0000
F81G	39	-4583.1347	9244.2694	0.0000	1.0000
JC69I	36	-4594.5424	9261.0848	0.0000	1.0000
F81I	39	-4591.9409	9261.8818	0.0000	1.0000
F81	38	-4680.6424	9437.2849	0.0000	1.0000
JC69	35	-4684.4560	9438.9120	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-4399.7348	9053.0417	0.9338	0.9338
K2PIG	38	-4399.4369	9059.2992	0.0409	0.9746
SYMIG	42	-4386.2845	9060.4075	0.0235	0.9981
HKYG	40	-4396.2845	9066.7010	0.0010	0.9991
GTRIG	45	-4379.4071	9067.2126	0.0008	0.9999
HKYIG	41	-4395.9804	9072.9460	0.0000	0.9999
K2PI	37	-4409.8815	9073.3351	0.0000	1.0000
SYMG	41	-4396.9723	9074.9299	0.0000	1.0000
GTRG	44	-4390.1721	9081.8893	0.0000	1.0000
SYMI	41	-4407.6510	9096.2873	0.0000	1.0000
GTRI	44	-4401.6403	9104.8258	0.0000	1.0000
HKYI	40	-4416.8242	9107.7803	0.0000	1.0000
K2P	36	-4500.5644	9247.8475	0.0000	1.0000
SYM	40	-4503.7906	9281.7132	0.0000	1.0000
GTR	43	-4494.7314	9284.1547	0.0000	1.0000
HKY	39	-4512.1914	9291.6614	0.0000	1.0000
JC69IG	37	-4575.1512	9403.8746	0.0000	1.0000
F81IG	40	-4571.9735	9418.0789	0.0000	1.0000
JC69G	36	-4586.0452	9418.8092	0.0000	1.0000
F81G	39	-4583.1347	9433.5480	0.0000	1.0000
JC69I	36	-4594.5424	9435.8035	0.0000	1.0000
F81I	39	-4591.9409	9451.1605	0.0000	1.0000
JC69	35	-4684.4560	9608.7775	0.0000	1.0000
F81	38	-4680.6424	9621.7102	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/chr2_23690.phylip = 1 - 947;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_23690.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/chr2_23690.phylip = 1 - 947;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_23690.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 (K2PG)]

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

