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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	44	-4307.6605	8708.1620	0.2382	0.2382
SYMG	41	-4311.0028	8708.2004	0.2337	0.4720
K2PG	37	-4315.7514	8708.9112	0.1638	0.6358
SYMIG	42	-4310.8561	8710.1170	0.0896	0.7254
HKYG	40	-4313.0877	8710.1657	0.0875	0.8129
GTRIG	45	-4307.5661	8710.1995	0.0860	0.8989
K2PIG	38	-4315.5572	8710.7116	0.0666	0.9655
HKYIG	41	-4312.9152	8712.0252	0.0345	1.0000
GTRI	44	-4338.4861	8769.8133	0.0000	1.0000
SYMI	41	-4343.5629	8773.3208	0.0000	1.0000
K2PI	37	-4350.3013	8778.0110	0.0000	1.0000
HKYI	40	-4347.0501	8778.0905	0.0000	1.0000
JC69G	36	-4395.6257	8866.4766	0.0000	1.0000
F81G	39	-4392.5712	8866.9334	0.0000	1.0000
JC69IG	37	-4395.4793	8868.3672	0.0000	1.0000
F81IG	40	-4392.4337	8868.8576	0.0000	1.0000
F81I	39	-4422.0780	8925.9469	0.0000	1.0000
JC69I	36	-4425.7143	8926.6538	0.0000	1.0000
GTR	43	-4445.9898	8982.5999	0.0000	1.0000
SYM	40	-4454.8524	8993.6950	0.0000	1.0000
HKY	39	-4459.6592	9001.1094	0.0000	1.0000
K2P	36	-4465.2392	9005.7036	0.0000	1.0000
F81	38	-4530.9120	9141.4211	0.0000	1.0000
JC69	35	-4535.5606	9144.1684	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-4307.6605	8703.3209	0.3079	0.3079
SYMG	41	-4311.0028	8704.0055	0.2187	0.5266
GTRIG	45	-4307.5661	8705.1322	0.1245	0.6510
K2PG	37	-4315.7514	8705.5027	0.1034	0.7545
SYMIG	42	-4310.8561	8705.7121	0.0931	0.8476
HKYG	40	-4313.0877	8706.1754	0.0739	0.9215
K2PIG	38	-4315.5572	8707.1145	0.0462	0.9677
HKYIG	41	-4312.9152	8707.8303	0.0323	1.0000
GTRI	44	-4338.4861	8764.9722	0.0000	1.0000
SYMI	41	-4343.5629	8769.1259	0.0000	1.0000
HKYI	40	-4347.0501	8774.1003	0.0000	1.0000
K2PI	37	-4350.3013	8774.6026	0.0000	1.0000
F81G	39	-4392.5712	8863.1424	0.0000	1.0000
JC69G	36	-4395.6257	8863.2514	0.0000	1.0000
F81IG	40	-4392.4337	8864.8673	0.0000	1.0000
JC69IG	37	-4395.4793	8864.9587	0.0000	1.0000
F81I	39	-4422.0780	8922.1559	0.0000	1.0000
JC69I	36	-4425.7143	8923.4286	0.0000	1.0000
GTR	43	-4445.9898	8977.9796	0.0000	1.0000
SYM	40	-4454.8524	8989.7047	0.0000	1.0000
HKY	39	-4459.6592	8997.3184	0.0000	1.0000
K2P	36	-4465.2392	9002.4784	0.0000	1.0000
F81	38	-4530.9120	9137.8240	0.0000	1.0000
JC69	35	-4535.5606	9141.1212	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-4315.7514	8881.6381	0.9596	0.9596
K2PIG	38	-4315.5572	8888.0102	0.0397	0.9993
HKYG	40	-4313.0877	8896.5920	0.0005	0.9998
SYMG	41	-4311.0028	8899.1825	0.0001	1.0000
HKYIG	41	-4312.9152	8903.0073	0.0000	1.0000
SYMIG	42	-4310.8561	8905.6495	0.0000	1.0000
GTRG	44	-4307.6605	8912.7791	0.0000	1.0000
GTRIG	45	-4307.5661	8919.3508	0.0000	1.0000
K2PI	37	-4350.3013	8950.7379	0.0000	1.0000
SYMI	41	-4343.5629	8964.3029	0.0000	1.0000
HKYI	40	-4347.0501	8964.5168	0.0000	1.0000
GTRI	44	-4338.4861	8974.4305	0.0000	1.0000
JC69G	36	-4395.6257	9034.6263	0.0000	1.0000
JC69IG	37	-4395.4793	9041.0940	0.0000	1.0000
F81G	39	-4392.5712	9048.7986	0.0000	1.0000
F81IG	40	-4392.4337	9055.2839	0.0000	1.0000
JC69I	36	-4425.7143	9094.8036	0.0000	1.0000
F81I	39	-4422.0780	9107.8121	0.0000	1.0000
K2P	36	-4465.2392	9173.8534	0.0000	1.0000
SYM	40	-4454.8524	9180.1213	0.0000	1.0000
GTR	43	-4445.9898	9182.6775	0.0000	1.0000
HKY	39	-4459.6592	9182.9745	0.0000	1.0000
JC69	35	-4535.5606	9307.7357	0.0000	1.0000
F81	38	-4530.9120	9318.7198	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/chr6_3036.phylip = 1 - 863;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3036.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/chr6_3036.phylip = 1 - 863;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3036.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 (K2PG)]

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

