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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	-2372.6402	4841.6348	0.8537	0.8537
GTRIG	45	-2373.2934	4845.3394	0.1339	0.9876
GTRI	44	-2376.8735	4850.1014	0.0124	1.0000
HKYG	40	-2393.3828	4873.6274	0.0000	1.0000
HKYIG	41	-2393.5055	4876.2311	0.0000	1.0000
HKYI	40	-2394.8948	4876.6515	0.0000	1.0000
SYMG	41	-2395.3588	4879.9377	0.0000	1.0000
SYMIG	42	-2396.0960	4883.7802	0.0000	1.0000
SYMI	41	-2397.4980	4884.2161	0.0000	1.0000
K2PI	37	-2409.2723	4898.3908	0.0000	1.0000
K2PIG	38	-2408.2512	4898.6774	0.0000	1.0000
K2PG	37	-2409.4745	4898.7953	0.0000	1.0000
F81I	39	-2436.3622	4957.2379	0.0000	1.0000
F81IG	40	-2435.7746	4958.4111	0.0000	1.0000
F81G	39	-2437.3836	4959.2808	0.0000	1.0000
GTR	43	-2433.1929	4960.3521	0.0000	1.0000
JC69G	36	-2443.3172	4964.1614	0.0000	1.0000
JC69I	36	-2444.1939	4965.9148	0.0000	1.0000
JC69IG	37	-2443.5842	4967.0145	0.0000	1.0000
HKY	39	-2451.2259	4986.9654	0.0000	1.0000
SYM	40	-2453.1515	4993.1649	0.0000	1.0000
K2P	36	-2471.8030	5021.1329	0.0000	1.0000
F81	38	-2496.1515	5074.4779	0.0000	1.0000
JC69	35	-2499.7169	5074.6512	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2372.6402	4833.2803	0.8292	0.8292
GTRIG	45	-2373.2934	4836.5867	0.1587	0.9880
GTRI	44	-2376.8735	4841.7470	0.0120	1.0000
HKYG	40	-2393.3828	4866.7655	0.0000	1.0000
HKYIG	41	-2393.5055	4869.0110	0.0000	1.0000
HKYI	40	-2394.8948	4869.7895	0.0000	1.0000
SYMG	41	-2395.3588	4872.7176	0.0000	1.0000
SYMIG	42	-2396.0960	4876.1920	0.0000	1.0000
SYMI	41	-2397.4980	4876.9960	0.0000	1.0000
K2PIG	38	-2408.2512	4892.5024	0.0000	1.0000
K2PI	37	-2409.2723	4892.5446	0.0000	1.0000
K2PG	37	-2409.4745	4892.9491	0.0000	1.0000
F81I	39	-2436.3622	4950.7244	0.0000	1.0000
F81IG	40	-2435.7746	4951.5492	0.0000	1.0000
GTR	43	-2433.1929	4952.3858	0.0000	1.0000
F81G	39	-2437.3836	4952.7672	0.0000	1.0000
JC69G	36	-2443.3172	4958.6344	0.0000	1.0000
JC69I	36	-2444.1939	4960.3879	0.0000	1.0000
JC69IG	37	-2443.5842	4961.1683	0.0000	1.0000
HKY	39	-2451.2259	4980.4518	0.0000	1.0000
SYM	40	-2453.1515	4986.3030	0.0000	1.0000
K2P	36	-2471.8030	5015.6059	0.0000	1.0000
F81	38	-2496.1515	5068.3029	0.0000	1.0000
JC69	35	-2499.7169	5069.4338	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-2372.6402	5020.3641	0.9637	0.9637
GTRIG	45	-2373.2934	5027.9224	0.0220	0.9857
GTRI	44	-2376.8735	5028.8307	0.0140	0.9997
HKYG	40	-2393.3828	5036.8417	0.0003	0.9999
HKYI	40	-2394.8948	5039.8657	0.0001	1.0000
HKYIG	41	-2393.5055	5043.3390	0.0000	1.0000
SYMG	41	-2395.3588	5047.0456	0.0000	1.0000
K2PI	37	-2409.2723	5049.8651	0.0000	1.0000
K2PG	37	-2409.4745	5050.2695	0.0000	1.0000
SYMI	41	-2397.4980	5051.3241	0.0000	1.0000
K2PIG	38	-2408.2512	5054.0748	0.0000	1.0000
SYMIG	42	-2396.0960	5054.7719	0.0000	1.0000
JC69G	36	-2443.3172	5111.7029	0.0000	1.0000
JC69I	36	-2444.1939	5113.4564	0.0000	1.0000
F81I	39	-2436.3622	5116.5486	0.0000	1.0000
JC69IG	37	-2443.5842	5118.4888	0.0000	1.0000
F81G	39	-2437.3836	5118.5915	0.0000	1.0000
F81IG	40	-2435.7746	5121.6254	0.0000	1.0000
GTR	43	-2433.1929	5135.2176	0.0000	1.0000
HKY	39	-2451.2259	5146.2761	0.0000	1.0000
SYM	40	-2453.1515	5156.3791	0.0000	1.0000
K2P	36	-2471.8030	5168.6745	0.0000	1.0000
JC69	35	-2499.7169	5218.2504	0.0000	1.0000
F81	38	-2496.1515	5229.8753	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/chr2_13034.phylip = 1 - 519;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_13034.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/chr2_13034.phylip = 1 - 519;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_13034.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 (GTRG)]

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


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

