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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	40	-2371.0930	4827.5368	0.4553	0.4553
GTRG	44	-2367.0987	4828.6999	0.2545	0.7098
HKYIG	41	-2370.8950	4829.4175	0.1778	0.8876
GTRIG	45	-2366.8974	4830.6041	0.0982	0.9858
K2PG	37	-2378.8498	4836.2645	0.0058	0.9916
HKYI	40	-2375.8220	4836.9948	0.0040	0.9956
K2PIG	38	-2378.6548	4838.1292	0.0023	0.9979
GTRI	44	-2372.4404	4839.3832	0.0012	0.9991
SYMG	41	-2376.6457	4840.9188	0.0006	0.9997
SYMIG	42	-2376.4896	4842.8907	0.0002	0.9999
K2PI	37	-2383.0737	4844.7123	0.0001	1.0000
SYMI	41	-2381.1817	4849.9909	0.0000	1.0000
HKY	39	-2402.6868	4888.4551	0.0000	1.0000
GTR	43	-2399.2871	4890.7774	0.0000	1.0000
K2P	36	-2410.4559	4897.2295	0.0000	1.0000
SYM	40	-2407.5913	4900.5333	0.0000	1.0000
F81G	39	-2448.3540	4979.7894	0.0000	1.0000
F81IG	40	-2448.1623	4981.6753	0.0000	1.0000
JC69G	36	-2455.3024	4986.9224	0.0000	1.0000
JC69IG	37	-2455.1331	4988.8311	0.0000	1.0000
F81I	39	-2455.3640	4993.8095	0.0000	1.0000
JC69I	36	-2462.0775	5000.4727	0.0000	1.0000
F81	38	-2481.8187	5044.4570	0.0000	1.0000
JC69	35	-2488.1329	5050.3435	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-2371.0930	4822.1860	0.3427	0.3427
GTRG	44	-2367.0987	4822.1975	0.3407	0.6834
HKYIG	41	-2370.8950	4823.7900	0.1537	0.8370
GTRIG	45	-2366.8974	4823.7948	0.1533	0.9903
HKYI	40	-2375.8220	4831.6441	0.0030	0.9934
K2PG	37	-2378.8498	4831.6995	0.0029	0.9963
GTRI	44	-2372.4404	4832.8808	0.0016	0.9979
K2PIG	38	-2378.6548	4833.3096	0.0013	0.9993
SYMG	41	-2376.6457	4835.2913	0.0005	0.9997
SYMIG	42	-2376.4896	4836.9791	0.0002	1.0000
K2PI	37	-2383.0737	4840.1473	0.0000	1.0000
SYMI	41	-2381.1817	4844.3634	0.0000	1.0000
HKY	39	-2402.6868	4883.3737	0.0000	1.0000
GTR	43	-2399.2871	4884.5741	0.0000	1.0000
K2P	36	-2410.4559	4892.9119	0.0000	1.0000
SYM	40	-2407.5913	4895.1826	0.0000	1.0000
F81G	39	-2448.3540	4974.7080	0.0000	1.0000
F81IG	40	-2448.1623	4976.3245	0.0000	1.0000
JC69G	36	-2455.3024	4982.6047	0.0000	1.0000
JC69IG	37	-2455.1331	4984.2661	0.0000	1.0000
F81I	39	-2455.3640	4988.7280	0.0000	1.0000
JC69I	36	-2462.0775	4996.1550	0.0000	1.0000
F81	38	-2481.8187	5039.6375	0.0000	1.0000
JC69	35	-2488.1329	5046.2658	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-2378.8498	4997.5745	0.8266	0.8266
HKYG	40	-2371.0930	5001.5103	0.1155	0.9421
K2PIG	38	-2378.6548	5003.6677	0.0393	0.9813
K2PI	37	-2383.0737	5006.0223	0.0121	0.9934
HKYIG	41	-2370.8950	5007.5974	0.0055	0.9989
HKYI	40	-2375.8220	5010.9684	0.0010	1.0000
SYMG	41	-2376.6457	5019.0987	0.0000	1.0000
GTRG	44	-2367.0987	5019.4542	0.0000	1.0000
SYMIG	42	-2376.4896	5025.2696	0.0000	1.0000
GTRIG	45	-2366.8974	5025.5347	0.0000	1.0000
SYMI	41	-2381.1817	5028.1708	0.0000	1.0000
GTRI	44	-2372.4404	5030.1375	0.0000	1.0000
K2P	36	-2410.4559	5054.3037	0.0000	1.0000
HKY	39	-2402.6868	5058.2149	0.0000	1.0000
SYM	40	-2407.5913	5074.5069	0.0000	1.0000
GTR	43	-2399.2871	5077.3478	0.0000	1.0000
JC69G	36	-2455.3024	5143.9966	0.0000	1.0000
F81G	39	-2448.3540	5149.5492	0.0000	1.0000
JC69IG	37	-2455.1331	5150.1411	0.0000	1.0000
F81IG	40	-2448.1623	5155.6488	0.0000	1.0000
JC69I	36	-2462.0775	5157.5469	0.0000	1.0000
F81I	39	-2455.3640	5163.5692	0.0000	1.0000
JC69	35	-2488.1329	5203.1746	0.0000	1.0000
F81	38	-2481.8187	5209.9955	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4746.phylip = 1 - 654;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4746.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4746.phylip = 1 - 654;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4746.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(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/chr7_4746.phylip = 1 - 654;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4746.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

