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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: GTRI


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
GTRI	44	-2363.4601	4820.4978	0.4272	0.4272
GTRG	44	-2363.6235	4820.8244	0.3628	0.7900
GTRIG	45	-2363.0394	4821.9180	0.2100	1.0000
GTR	43	-2381.7380	4854.7980	0.0000	1.0000
HKYG	40	-2412.5618	4909.7174	0.0000	1.0000
HKYI	40	-2412.8699	4910.3337	0.0000	1.0000
HKYIG	41	-2412.6006	4912.0315	0.0000	1.0000
SYMG	41	-2425.9633	4938.7570	0.0000	1.0000
SYMI	41	-2426.1440	4939.1182	0.0000	1.0000
SYMIG	42	-2425.5346	4940.1423	0.0000	1.0000
HKY	39	-2434.8964	4952.1565	0.0000	1.0000
SYM	40	-2440.6428	4965.8794	0.0000	1.0000
F81G	39	-2462.7928	5007.9491	0.0000	1.0000
K2PG	37	-2465.5069	5008.9356	0.0000	1.0000
K2PI	37	-2466.2557	5010.4332	0.0000	1.0000
F81I	39	-2464.0590	5010.4817	0.0000	1.0000
K2PIG	38	-2465.8773	5011.8943	0.0000	1.0000
F81IG	40	-2463.9045	5012.4028	0.0000	1.0000
K2P	36	-2485.3430	5046.3963	0.0000	1.0000
F81	38	-2484.5380	5049.2156	0.0000	1.0000
JC69G	36	-2506.8686	5089.4474	0.0000	1.0000
JC69I	36	-2507.5559	5090.8221	0.0000	1.0000
JC69IG	37	-2507.2475	5092.4169	0.0000	1.0000
JC69	35	-2525.5123	5124.5294	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	44	-2363.4601	4814.9203	0.4150	0.4150
GTRG	44	-2363.6235	4815.2470	0.3525	0.7675
GTRIG	45	-2363.0394	4816.0788	0.2325	1.0000
GTR	43	-2381.7380	4849.4759	0.0000	1.0000
HKYG	40	-2412.5618	4905.1235	0.0000	1.0000
HKYI	40	-2412.8699	4905.7398	0.0000	1.0000
HKYIG	41	-2412.6006	4907.2012	0.0000	1.0000
SYMG	41	-2425.9633	4933.9267	0.0000	1.0000
SYMI	41	-2426.1440	4934.2879	0.0000	1.0000
SYMIG	42	-2425.5346	4935.0693	0.0000	1.0000
HKY	39	-2434.8964	4947.7928	0.0000	1.0000
SYM	40	-2440.6428	4961.2855	0.0000	1.0000
F81G	39	-2462.7928	5003.5855	0.0000	1.0000
K2PG	37	-2465.5069	5005.0137	0.0000	1.0000
F81I	39	-2464.0590	5006.1180	0.0000	1.0000
K2PI	37	-2466.2557	5006.5113	0.0000	1.0000
K2PIG	38	-2465.8773	5007.7546	0.0000	1.0000
F81IG	40	-2463.9045	5007.8090	0.0000	1.0000
K2P	36	-2485.3430	5042.6860	0.0000	1.0000
F81	38	-2484.5380	5045.0760	0.0000	1.0000
JC69G	36	-2506.8686	5085.7371	0.0000	1.0000
JC69I	36	-2507.5559	5087.1118	0.0000	1.0000
JC69IG	37	-2507.2475	5088.4950	0.0000	1.0000
JC69	35	-2525.5123	5121.0246	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	44	-2363.4601	5018.4959	0.5250	0.5250
GTRG	44	-2363.6235	5018.8226	0.4459	0.9709
GTRIG	45	-2363.0394	5024.2811	0.0291	1.0000
GTR	43	-2381.7380	5048.4248	0.0000	1.0000
HKYG	40	-2412.5618	5090.1922	0.0000	1.0000
HKYI	40	-2412.8699	5090.8085	0.0000	1.0000
HKYIG	41	-2412.6006	5096.8966	0.0000	1.0000
SYMG	41	-2425.9633	5123.6221	0.0000	1.0000
SYMI	41	-2426.1440	5123.9833	0.0000	1.0000
HKY	39	-2434.8964	5128.2348	0.0000	1.0000
SYMIG	42	-2425.5346	5129.3914	0.0000	1.0000
SYM	40	-2440.6428	5146.3542	0.0000	1.0000
K2PG	37	-2465.5069	5176.2023	0.0000	1.0000
K2PI	37	-2466.2557	5177.6999	0.0000	1.0000
K2PIG	38	-2465.8773	5183.5699	0.0000	1.0000
F81G	39	-2462.7928	5184.0275	0.0000	1.0000
F81I	39	-2464.0590	5186.5600	0.0000	1.0000
F81IG	40	-2463.9045	5192.8777	0.0000	1.0000
K2P	36	-2485.3430	5209.2479	0.0000	1.0000
F81	38	-2484.5380	5220.8912	0.0000	1.0000
JC69G	36	-2506.8686	5252.2989	0.0000	1.0000
JC69I	36	-2507.5559	5253.6737	0.0000	1.0000
JC69IG	37	-2507.2475	5259.6836	0.0000	1.0000
JC69	35	-2525.5123	5282.9597	0.0000	1.0000



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

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

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


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

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


[Mrbayes block for the best BIC model (GTRI)]

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


-------------------------------------------------------------
End of Output

