
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:00 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/chr8_4260.phylip" (19 taxa, 1190 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	-5775.6871	11642.8327	0.7466	0.7466
GTRIG	45	-5775.6875	11644.9938	0.2534	1.0000
GTRI	44	-5795.9302	11683.3189	0.0000	1.0000
SYMG	41	-5805.4787	11695.9574	0.0000	1.0000
SYMIG	42	-5805.4789	11698.1069	0.0000	1.0000
HKYG	40	-5811.8317	11706.5181	0.0000	1.0000
HKYIG	41	-5811.8321	11708.6641	0.0000	1.0000
K2PG	37	-5824.8190	11726.0790	0.0000	1.0000
K2PIG	38	-5824.8193	11728.2137	0.0000	1.0000
SYMI	41	-5823.1374	11731.2747	0.0000	1.0000
HKYI	40	-5830.3487	11743.5520	0.0000	1.0000
K2PI	37	-5843.5292	11763.4994	0.0000	1.0000
GTR	43	-5849.1125	11787.5269	0.0000	1.0000
SYM	40	-5876.3433	11835.5413	0.0000	1.0000
HKY	39	-5883.2209	11847.1548	0.0000	1.0000
K2P	36	-5898.2345	11870.7795	0.0000	1.0000
F81G	39	-5936.2722	11953.2574	0.0000	1.0000
F81IG	40	-5936.2726	11955.3998	0.0000	1.0000
JC69G	36	-5944.9371	11964.1848	0.0000	1.0000
JC69IG	37	-5944.9375	11966.3160	0.0000	1.0000
F81I	39	-5953.1175	11986.9481	0.0000	1.0000
JC69I	36	-5961.9772	11998.2649	0.0000	1.0000
F81	38	-6002.7401	12084.0554	0.0000	1.0000
JC69	35	-6012.5334	12097.2504	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-5775.6871	11639.3741	0.7311	0.7311
GTRIG	45	-5775.6875	11641.3749	0.2689	1.0000
GTRI	44	-5795.9302	11679.8604	0.0000	1.0000
SYMG	41	-5805.4787	11692.9574	0.0000	1.0000
SYMIG	42	-5805.4789	11694.9579	0.0000	1.0000
HKYG	40	-5811.8317	11703.6635	0.0000	1.0000
HKYIG	41	-5811.8321	11705.6641	0.0000	1.0000
K2PG	37	-5824.8190	11723.6380	0.0000	1.0000
K2PIG	38	-5824.8193	11725.6385	0.0000	1.0000
SYMI	41	-5823.1374	11728.2747	0.0000	1.0000
HKYI	40	-5830.3487	11740.6974	0.0000	1.0000
K2PI	37	-5843.5292	11761.0584	0.0000	1.0000
GTR	43	-5849.1125	11784.2250	0.0000	1.0000
SYM	40	-5876.3433	11832.6867	0.0000	1.0000
HKY	39	-5883.2209	11844.4418	0.0000	1.0000
K2P	36	-5898.2345	11868.4690	0.0000	1.0000
F81G	39	-5936.2722	11950.5443	0.0000	1.0000
F81IG	40	-5936.2726	11952.5451	0.0000	1.0000
JC69G	36	-5944.9371	11961.8743	0.0000	1.0000
JC69IG	37	-5944.9375	11963.8750	0.0000	1.0000
F81I	39	-5953.1175	11984.2351	0.0000	1.0000
JC69I	36	-5961.9772	11995.9544	0.0000	1.0000
F81	38	-6002.7401	12081.4802	0.0000	1.0000
JC69	35	-6012.5334	12095.0667	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-5775.6871	11862.9693	0.9718	0.9718
GTRIG	45	-5775.6875	11870.0518	0.0282	1.0000
SYMG	41	-5805.4787	11901.3074	0.0000	1.0000
GTRI	44	-5795.9302	11903.4556	0.0000	1.0000
HKYG	40	-5811.8317	11906.9318	0.0000	1.0000
SYMIG	42	-5805.4789	11908.3896	0.0000	1.0000
K2PG	37	-5824.8190	11911.6612	0.0000	1.0000
HKYIG	41	-5811.8321	11914.0142	0.0000	1.0000
K2PIG	38	-5824.8193	11918.7434	0.0000	1.0000
SYMI	41	-5823.1374	11936.6248	0.0000	1.0000
HKYI	40	-5830.3487	11943.9657	0.0000	1.0000
K2PI	37	-5843.5292	11949.0816	0.0000	1.0000
GTR	43	-5849.1125	12002.7385	0.0000	1.0000
SYM	40	-5876.3433	12035.9550	0.0000	1.0000
HKY	39	-5883.2209	12042.6284	0.0000	1.0000
K2P	36	-5898.2345	12051.4105	0.0000	1.0000
JC69G	36	-5944.9371	12144.8158	0.0000	1.0000
F81G	39	-5936.2722	12148.7310	0.0000	1.0000
JC69IG	37	-5944.9375	12151.8982	0.0000	1.0000
F81IG	40	-5936.2726	12155.8135	0.0000	1.0000
JC69I	36	-5961.9772	12178.8959	0.0000	1.0000
F81I	39	-5953.1175	12182.4217	0.0000	1.0000
JC69	35	-6012.5334	12272.9265	0.0000	1.0000
F81	38	-6002.7401	12274.5852	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/chr8_4260.phylip = 1 - 1190;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_4260.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/chr8_4260.phylip = 1 - 1190;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_4260.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/chr8_4260.phylip = 1 - 1190;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_4260.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

