
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:19 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/chr3_24921.phylip" (19 taxa, 306 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	-998.2487	2099.6698	0.4202	0.4202
GTRI	44	-998.4151	2100.0026	0.3558	0.7760
GTRIG	45	-998.2262	2102.3754	0.1086	0.8846
HKYG	40	-1006.5977	2105.5728	0.0220	0.9066
HKYI	40	-1006.8496	2106.0766	0.0171	0.9236
SYMG	41	-1005.5560	2106.1575	0.0164	0.9400
K2PG	37	-1010.8485	2106.1896	0.0161	0.9562
SYMI	41	-1005.5989	2106.2433	0.0157	0.9719
K2PI	37	-1011.0161	2106.5247	0.0136	0.9855
HKYIG	41	-1006.5774	2108.2003	0.0059	0.9914
K2PIG	38	-1010.8483	2108.7976	0.0044	0.9958
SYMIG	42	-1005.5720	2108.8779	0.0042	1.0000
GTR	43	-1010.9646	2122.3720	0.0000	1.0000
SYM	40	-1019.1901	2130.7575	0.0000	1.0000
HKY	39	-1020.6786	2131.0865	0.0000	1.0000
K2P	36	-1026.1393	2134.1820	0.0000	1.0000
JC69G	36	-1042.6656	2167.2346	0.0000	1.0000
F81G	39	-1038.7585	2167.2463	0.0000	1.0000
F81I	39	-1038.8879	2167.5050	0.0000	1.0000
JC69I	36	-1042.8069	2167.5172	0.0000	1.0000
JC69IG	37	-1042.6191	2169.7308	0.0000	1.0000
F81IG	40	-1038.7099	2169.7971	0.0000	1.0000
F81	38	-1052.5696	2192.2404	0.0000	1.0000
JC69	35	-1057.1547	2193.6428	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-998.2487	2084.4973	0.4366	0.4366
GTRI	44	-998.4151	2084.8302	0.3697	0.8063
GTRIG	45	-998.2262	2086.4524	0.1643	0.9706
SYMG	41	-1005.5560	2093.1121	0.0059	0.9765
HKYG	40	-1006.5977	2093.1955	0.0056	0.9822
SYMI	41	-1005.5989	2093.1978	0.0056	0.9878
HKYI	40	-1006.8496	2093.6992	0.0044	0.9922
SYMIG	42	-1005.5720	2095.1441	0.0021	0.9943
HKYIG	41	-1006.5774	2095.1548	0.0021	0.9964
K2PG	37	-1010.8485	2095.6971	0.0016	0.9980
K2PI	37	-1011.0161	2096.0322	0.0014	0.9994
K2PIG	38	-1010.8483	2097.6965	0.0006	1.0000
GTR	43	-1010.9646	2107.9293	0.0000	1.0000
SYM	40	-1019.1901	2118.3801	0.0000	1.0000
HKY	39	-1020.6786	2119.3572	0.0000	1.0000
K2P	36	-1026.1393	2124.2786	0.0000	1.0000
F81G	39	-1038.7585	2155.5170	0.0000	1.0000
F81I	39	-1038.8879	2155.7757	0.0000	1.0000
JC69G	36	-1042.6656	2157.3312	0.0000	1.0000
F81IG	40	-1038.7099	2157.4198	0.0000	1.0000
JC69I	36	-1042.8069	2157.6138	0.0000	1.0000
JC69IG	37	-1042.6191	2159.2382	0.0000	1.0000
F81	38	-1052.5696	2181.1393	0.0000	1.0000
JC69	35	-1057.1547	2184.3094	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1010.8485	2233.4697	0.5175	0.5175
K2PI	37	-1011.0161	2233.8048	0.4377	0.9552
K2PIG	38	-1010.8483	2239.1928	0.0296	0.9847
HKYG	40	-1006.5977	2242.1389	0.0068	0.9915
HKYI	40	-1006.8496	2242.6426	0.0053	0.9968
SYMG	41	-1005.5560	2245.7790	0.0011	0.9979
SYMI	41	-1005.5989	2245.8648	0.0011	0.9990
HKYIG	41	-1006.5774	2247.8218	0.0004	0.9994
GTRG	44	-998.2487	2248.3351	0.0003	0.9997
GTRI	44	-998.4151	2248.6679	0.0003	0.9999
SYMIG	42	-1005.5720	2251.5347	0.0001	1.0000
GTRIG	45	-998.2262	2254.0137	0.0000	1.0000
K2P	36	-1026.1393	2258.3277	0.0000	1.0000
HKY	39	-1020.6786	2264.5770	0.0000	1.0000
SYM	40	-1019.1901	2267.3235	0.0000	1.0000
GTR	43	-1010.9646	2268.0434	0.0000	1.0000
JC69G	36	-1042.6656	2291.3803	0.0000	1.0000
JC69I	36	-1042.8069	2291.6629	0.0000	1.0000
JC69IG	37	-1042.6191	2297.0109	0.0000	1.0000
F81G	39	-1038.7585	2300.7368	0.0000	1.0000
F81I	39	-1038.8879	2300.9955	0.0000	1.0000
F81IG	40	-1038.7099	2306.3632	0.0000	1.0000
JC69	35	-1057.1547	2314.6349	0.0000	1.0000
F81	38	-1052.5696	2322.6355	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/chr3_24921.phylip = 1 - 306;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_24921.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/chr3_24921.phylip = 1 - 306;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_24921.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/chr3_24921.phylip = 1 - 306;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_24921.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

