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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	45	-4815.5583	9724.5868	0.6724	0.6724
GTRG	44	-4817.3540	9726.0247	0.3276	1.0000
GTRI	44	-4828.0600	9747.4367	0.0000	1.0000
HKYG	40	-4835.3841	9753.5062	0.0000	1.0000
HKYIG	41	-4835.1238	9755.1247	0.0000	1.0000
HKYI	40	-4849.8381	9782.4141	0.0000	1.0000
F81G	39	-4920.3262	9921.2545	0.0000	1.0000
F81IG	40	-4920.0925	9922.9230	0.0000	1.0000
SYMG	41	-4922.6128	9930.1027	0.0000	1.0000
GTR	43	-4923.5091	9936.1847	0.0000	1.0000
K2PG	37	-4930.1229	9936.5872	0.0000	1.0000
SYMIG	42	-4924.8583	9936.7366	0.0000	1.0000
K2PIG	38	-4932.2387	9942.9473	0.0000	1.0000
F81I	39	-4931.8224	9944.2471	0.0000	1.0000
SYMI	41	-4935.2249	9955.3270	0.0000	1.0000
K2PI	37	-4941.6359	9959.6132	0.0000	1.0000
HKY	39	-4943.1630	9966.9281	0.0000	1.0000
JC69G	36	-4992.0548	10058.3259	0.0000	1.0000
JC69IG	37	-4993.3646	10063.0706	0.0000	1.0000
JC69I	36	-5002.5877	10079.3917	0.0000	1.0000
F81	38	-5015.7443	10109.9585	0.0000	1.0000
SYM	40	-5021.3159	10125.3698	0.0000	1.0000
K2P	36	-5027.9543	10130.1249	0.0000	1.0000
JC69	35	-5085.9967	10244.0882	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-4815.5583	9721.1165	0.6891	0.6891
GTRG	44	-4817.3540	9722.7081	0.3109	1.0000
GTRI	44	-4828.0600	9744.1201	0.0000	1.0000
HKYG	40	-4835.3841	9750.7683	0.0000	1.0000
HKYIG	41	-4835.1238	9752.2475	0.0000	1.0000
HKYI	40	-4849.8381	9779.6762	0.0000	1.0000
F81G	39	-4920.3262	9918.6524	0.0000	1.0000
F81IG	40	-4920.0925	9920.1851	0.0000	1.0000
SYMG	41	-4922.6128	9927.2255	0.0000	1.0000
GTR	43	-4923.5091	9933.0182	0.0000	1.0000
SYMIG	42	-4924.8583	9933.7166	0.0000	1.0000
K2PG	37	-4930.1229	9934.2459	0.0000	1.0000
K2PIG	38	-4932.2387	9940.4773	0.0000	1.0000
F81I	39	-4931.8224	9941.6449	0.0000	1.0000
SYMI	41	-4935.2249	9952.4498	0.0000	1.0000
K2PI	37	-4941.6359	9957.2719	0.0000	1.0000
HKY	39	-4943.1630	9964.3260	0.0000	1.0000
JC69G	36	-4992.0548	10056.1096	0.0000	1.0000
JC69IG	37	-4993.3646	10060.7292	0.0000	1.0000
JC69I	36	-5002.5877	10077.1754	0.0000	1.0000
F81	38	-5015.7443	10107.4885	0.0000	1.0000
SYM	40	-5021.3159	10122.6319	0.0000	1.0000
K2P	36	-5027.9543	10127.9086	0.0000	1.0000
JC69	35	-5085.9967	10241.9935	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-4817.3540	9948.0787	0.8370	0.8370
GTRIG	45	-4815.5583	9951.6092	0.1433	0.9803
HKYG	40	-4835.3841	9955.6507	0.0190	0.9993
HKYIG	41	-4835.1238	9962.2520	0.0007	1.0000
GTRI	44	-4828.0600	9969.4907	0.0000	1.0000
HKYI	40	-4849.8381	9984.5586	0.0000	1.0000
F81G	39	-4920.3262	10118.4127	0.0000	1.0000
K2PG	37	-4930.1229	10123.7621	0.0000	1.0000
F81IG	40	-4920.0925	10125.0675	0.0000	1.0000
K2PIG	38	-4932.2387	10135.1156	0.0000	1.0000
SYMG	41	-4922.6128	10137.2300	0.0000	1.0000
F81I	39	-4931.8224	10141.4052	0.0000	1.0000
K2PI	37	-4941.6359	10146.7881	0.0000	1.0000
SYMIG	42	-4924.8583	10148.8431	0.0000	1.0000
GTR	43	-4923.5091	10153.2668	0.0000	1.0000
SYMI	41	-4935.2249	10162.4543	0.0000	1.0000
HKY	39	-4943.1630	10164.0863	0.0000	1.0000
JC69G	36	-4992.0548	10240.5038	0.0000	1.0000
JC69IG	37	-4993.3646	10250.2455	0.0000	1.0000
JC69I	36	-5002.5877	10261.5696	0.0000	1.0000
F81	38	-5015.7443	10302.1268	0.0000	1.0000
K2P	36	-5027.9543	10312.3027	0.0000	1.0000
SYM	40	-5021.3159	10327.5143	0.0000	1.0000
JC69	35	-5085.9967	10421.2656	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_675.phylip = 1 - 1239;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_675.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_675.phylip = 1 - 1239;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_675.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.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_675.phylip = 1 - 1239;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_675.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

