
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:11 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_4747.phylip" (19 taxa, 931 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	-3287.9140	6668.2975	0.6914	0.6914
GTRIG	45	-3287.6327	6669.9433	0.3036	0.9951
GTRI	44	-3292.8558	6678.1812	0.0049	1.0000
HKYG	40	-3308.7818	6701.2490	0.0000	1.0000
HKYIG	41	-3308.7207	6703.3154	0.0000	1.0000
HKYI	40	-3312.1190	6707.9235	0.0000	1.0000
GTR	43	-3326.2212	6742.7085	0.0000	1.0000
SYMG	41	-3333.9038	6753.6816	0.0000	1.0000
SYMIG	42	-3333.8659	6755.7993	0.0000	1.0000
K2PG	37	-3340.7838	6758.7165	0.0000	1.0000
SYMI	41	-3336.6782	6759.2305	0.0000	1.0000
K2PIG	38	-3340.7498	6760.8225	0.0000	1.0000
K2PI	37	-3343.4461	6764.0410	0.0000	1.0000
HKY	39	-3344.1135	6769.7287	0.0000	1.0000
SYM	40	-3370.2211	6824.1276	0.0000	1.0000
K2P	36	-3377.0419	6829.0638	0.0000	1.0000
F81G	39	-3395.1309	6871.7634	0.0000	1.0000
F81IG	40	-3395.1221	6873.9295	0.0000	1.0000
F81I	39	-3397.5957	6876.6931	0.0000	1.0000
JC69G	36	-3419.5528	6914.0854	0.0000	1.0000
JC69IG	37	-3419.5470	6916.2429	0.0000	1.0000
JC69I	36	-3421.7730	6918.5259	0.0000	1.0000
F81	38	-3427.8458	6935.0144	0.0000	1.0000
JC69	35	-3453.1582	6979.1320	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-3287.9140	6663.8280	0.6691	0.6691
GTRIG	45	-3287.6327	6665.2654	0.3261	0.9952
GTRI	44	-3292.8558	6673.7117	0.0048	1.0000
HKYG	40	-3308.7818	6697.5636	0.0000	1.0000
HKYIG	41	-3308.7207	6699.4414	0.0000	1.0000
HKYI	40	-3312.1190	6704.2381	0.0000	1.0000
GTR	43	-3326.2212	6738.4424	0.0000	1.0000
SYMG	41	-3333.9038	6749.8076	0.0000	1.0000
SYMIG	42	-3333.8659	6751.7318	0.0000	1.0000
SYMI	41	-3336.6782	6755.3564	0.0000	1.0000
K2PG	37	-3340.7838	6755.5676	0.0000	1.0000
K2PIG	38	-3340.7498	6757.4996	0.0000	1.0000
K2PI	37	-3343.4461	6760.8921	0.0000	1.0000
HKY	39	-3344.1135	6766.2270	0.0000	1.0000
SYM	40	-3370.2211	6820.4422	0.0000	1.0000
K2P	36	-3377.0419	6826.0839	0.0000	1.0000
F81G	39	-3395.1309	6868.2617	0.0000	1.0000
F81IG	40	-3395.1221	6870.2441	0.0000	1.0000
F81I	39	-3397.5957	6873.1915	0.0000	1.0000
JC69G	36	-3419.5528	6911.1056	0.0000	1.0000
JC69IG	37	-3419.5470	6913.0939	0.0000	1.0000
JC69I	36	-3421.7730	6915.5460	0.0000	1.0000
F81	38	-3427.8458	6931.6915	0.0000	1.0000
JC69	35	-3453.1582	6976.3163	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-3287.9140	6876.6234	0.9511	0.9511
GTRIG	45	-3287.6327	6882.8970	0.0413	0.9924
GTRI	44	-3292.8558	6886.5071	0.0068	0.9992
HKYG	40	-3308.7818	6891.0140	0.0007	0.9999
HKYI	40	-3312.1190	6897.6885	0.0000	1.0000
HKYIG	41	-3308.7207	6897.7280	0.0000	1.0000
K2PG	37	-3340.7838	6934.5092	0.0000	1.0000
K2PI	37	-3343.4461	6939.8337	0.0000	1.0000
K2PIG	38	-3340.7498	6941.2775	0.0000	1.0000
GTR	43	-3326.2212	6946.4015	0.0000	1.0000
SYMG	41	-3333.9038	6948.0942	0.0000	1.0000
SYMI	41	-3336.6782	6953.6431	0.0000	1.0000
HKY	39	-3344.1135	6954.8411	0.0000	1.0000
SYMIG	42	-3333.8659	6954.8547	0.0000	1.0000
K2P	36	-3377.0419	7000.1892	0.0000	1.0000
SYM	40	-3370.2211	7013.8926	0.0000	1.0000
F81G	39	-3395.1309	7056.8758	0.0000	1.0000
F81I	39	-3397.5957	7061.8056	0.0000	1.0000
F81IG	40	-3395.1221	7063.6945	0.0000	1.0000
JC69G	36	-3419.5528	7085.2109	0.0000	1.0000
JC69I	36	-3421.7730	7089.6514	0.0000	1.0000
JC69IG	37	-3419.5470	7092.0355	0.0000	1.0000
F81	38	-3427.8458	7115.4694	0.0000	1.0000
JC69	35	-3453.1582	7145.5854	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/chr7_4747.phylip = 1 - 931;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4747.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/chr7_4747.phylip = 1 - 931;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4747.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/chr7_4747.phylip = 1 - 931;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_4747.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

