
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:14 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_496.phylip" (19 taxa, 669 characters)

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


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	-2991.5479	6077.4419	0.7496	0.7496
GTRIG	45	-2991.5072	6079.6597	0.2473	0.9969
SYMG	41	-3000.7730	6089.0389	0.0023	0.9992
SYMIG	42	-3000.6953	6091.1606	0.0008	1.0000
K2PG	37	-3014.9259	6108.3082	0.0000	1.0000
HKYG	40	-3012.3780	6109.9790	0.0000	1.0000
GTRI	44	-3009.1477	6112.6415	0.0000	1.0000
K2PIG	38	-3018.5728	6117.8504	0.0000	1.0000
HKYIG	41	-3016.0165	6119.5258	0.0000	1.0000
SYMI	41	-3016.9865	6121.4658	0.0000	1.0000
K2PI	37	-3030.6704	6139.7972	0.0000	1.0000
HKYI	40	-3028.0915	6141.4059	0.0000	1.0000
GTR	43	-3054.2154	6200.4852	0.0000	1.0000
SYM	40	-3062.7504	6210.7236	0.0000	1.0000
HKY	39	-3074.9649	6232.8901	0.0000	1.0000
K2P	36	-3078.3620	6232.9392	0.0000	1.0000
F81IG	40	-3111.0490	6307.3209	0.0000	1.0000
JC69G	36	-3117.7145	6311.6441	0.0000	1.0000
JC69IG	37	-3117.6651	6313.7866	0.0000	1.0000
F81G	39	-3116.1977	6315.3557	0.0000	1.0000
JC69I	36	-3131.2671	6338.7493	0.0000	1.0000
F81I	39	-3129.7963	6342.5529	0.0000	1.0000
JC69	35	-3171.4997	6416.9805	0.0000	1.0000
F81	38	-3169.3237	6419.3521	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-2991.5479	6071.0957	0.7215	0.7215
GTRIG	45	-2991.5072	6073.0144	0.2765	0.9980
SYMG	41	-3000.7730	6083.5461	0.0014	0.9994
SYMIG	42	-3000.6953	6085.3906	0.0006	1.0000
K2PG	37	-3014.9259	6103.8518	0.0000	1.0000
HKYG	40	-3012.3780	6104.7560	0.0000	1.0000
GTRI	44	-3009.1477	6106.2953	0.0000	1.0000
K2PIG	38	-3018.5728	6113.1457	0.0000	1.0000
HKYIG	41	-3016.0165	6114.0330	0.0000	1.0000
SYMI	41	-3016.9865	6115.9730	0.0000	1.0000
K2PI	37	-3030.6704	6135.3408	0.0000	1.0000
HKYI	40	-3028.0915	6136.1829	0.0000	1.0000
GTR	43	-3054.2154	6194.4308	0.0000	1.0000
SYM	40	-3062.7504	6205.5007	0.0000	1.0000
HKY	39	-3074.9649	6227.9298	0.0000	1.0000
K2P	36	-3078.3620	6228.7240	0.0000	1.0000
F81IG	40	-3111.0490	6302.0980	0.0000	1.0000
JC69G	36	-3117.7145	6307.4290	0.0000	1.0000
JC69IG	37	-3117.6651	6309.3302	0.0000	1.0000
F81G	39	-3116.1977	6310.3955	0.0000	1.0000
JC69I	36	-3131.2671	6334.5341	0.0000	1.0000
F81I	39	-3129.7963	6337.5927	0.0000	1.0000
JC69	35	-3171.4997	6412.9995	0.0000	1.0000
F81	38	-3169.3237	6414.6473	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMG	41	-3000.7730	6268.2832	0.5071	0.5071
GTRG	44	-2991.5479	6269.3502	0.2974	0.8046
K2PG	37	-3014.9259	6270.5658	0.1620	0.9665
SYMIG	42	-3000.6953	6274.6335	0.0212	0.9877
GTRIG	45	-2991.5072	6275.7747	0.0120	0.9997
K2PIG	38	-3018.5728	6284.3655	0.0002	0.9999
HKYG	40	-3012.3780	6284.9874	0.0001	1.0000
HKYIG	41	-3016.0165	6298.7701	0.0000	1.0000
SYMI	41	-3016.9865	6300.7101	0.0000	1.0000
K2PI	37	-3030.6704	6302.0548	0.0000	1.0000
GTRI	44	-3009.1477	6304.5498	0.0000	1.0000
HKYI	40	-3028.0915	6316.4143	0.0000	1.0000
SYM	40	-3062.7504	6385.7321	0.0000	1.0000
GTR	43	-3054.2154	6388.1795	0.0000	1.0000
K2P	36	-3078.3620	6390.9322	0.0000	1.0000
HKY	39	-3074.9649	6403.6554	0.0000	1.0000
JC69G	36	-3117.7145	6469.6372	0.0000	1.0000
JC69IG	37	-3117.6651	6476.0442	0.0000	1.0000
F81IG	40	-3111.0490	6482.3293	0.0000	1.0000
F81G	39	-3116.1977	6486.1211	0.0000	1.0000
JC69I	36	-3131.2671	6496.7424	0.0000	1.0000
F81I	39	-3129.7963	6513.3183	0.0000	1.0000
JC69	35	-3171.4997	6570.7019	0.0000	1.0000
F81	38	-3169.3237	6585.8671	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_496.phylip = 1 - 669;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_496.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_496.phylip = 1 - 669;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_496.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 (SYMG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_496.phylip = 1 - 669;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_496.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=Fixed(Equal);
END;


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

