
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:22 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_9156.phylip" (19 taxa, 946 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	-3802.3551	7699.3101	0.6328	0.6328
GTRG	44	-3804.0024	7700.4000	0.3669	0.9997
SYMIG	42	-3813.7942	7715.5885	0.0002	0.9999
SYMG	41	-3815.3429	7716.4956	0.0001	1.0000
GTRI	44	-3817.2787	7726.9525	0.0000	1.0000
HKYIG	41	-3823.1485	7732.1067	0.0000	1.0000
HKYG	40	-3824.5536	7732.7314	0.0000	1.0000
K2PIG	38	-3827.2708	7733.8096	0.0000	1.0000
K2PG	37	-3828.7093	7734.5155	0.0000	1.0000
SYMI	41	-3827.2801	7740.3699	0.0000	1.0000
HKYI	40	-3835.7212	7755.0667	0.0000	1.0000
K2PI	37	-3839.7768	7756.6504	0.0000	1.0000
GTR	43	-3908.6268	7907.4486	0.0000	1.0000
SYM	40	-3917.0227	7917.6696	0.0000	1.0000
F81IG	40	-3921.5275	7926.6794	0.0000	1.0000
F81G	39	-3922.6371	7926.7180	0.0000	1.0000
JC69IG	37	-3925.2959	7927.6887	0.0000	1.0000
JC69G	36	-3926.3827	7927.6960	0.0000	1.0000
HKY	39	-3925.1294	7931.7025	0.0000	1.0000
K2P	36	-3929.5375	7934.0057	0.0000	1.0000
F81I	39	-3931.9503	7945.3442	0.0000	1.0000
JC69I	36	-3935.5063	7945.9433	0.0000	1.0000
F81	38	-4017.4334	8114.1347	0.0000	1.0000
JC69	35	-4021.3064	8115.3820	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-3802.3551	7694.7101	0.6563	0.6563
GTRG	44	-3804.0024	7696.0049	0.3435	0.9998
SYMIG	42	-3813.7942	7711.5885	0.0001	0.9999
SYMG	41	-3815.3429	7712.6858	0.0001	1.0000
GTRI	44	-3817.2787	7722.5574	0.0000	1.0000
HKYIG	41	-3823.1485	7728.2970	0.0000	1.0000
HKYG	40	-3824.5536	7729.1071	0.0000	1.0000
K2PIG	38	-3827.2708	7730.5417	0.0000	1.0000
K2PG	37	-3828.7093	7731.4186	0.0000	1.0000
SYMI	41	-3827.2801	7736.5601	0.0000	1.0000
HKYI	40	-3835.7212	7751.4424	0.0000	1.0000
K2PI	37	-3839.7768	7753.5535	0.0000	1.0000
GTR	43	-3908.6268	7903.2535	0.0000	1.0000
SYM	40	-3917.0227	7914.0453	0.0000	1.0000
F81IG	40	-3921.5275	7923.0551	0.0000	1.0000
F81G	39	-3922.6371	7923.2743	0.0000	1.0000
JC69IG	37	-3925.2959	7924.5918	0.0000	1.0000
JC69G	36	-3926.3827	7924.7653	0.0000	1.0000
HKY	39	-3925.1294	7928.2588	0.0000	1.0000
K2P	36	-3929.5375	7931.0750	0.0000	1.0000
F81I	39	-3931.9503	7941.9005	0.0000	1.0000
JC69I	36	-3935.5063	7943.0126	0.0000	1.0000
F81	38	-4017.4334	8110.8668	0.0000	1.0000
JC69	35	-4021.3064	8112.6128	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-3804.0024	7909.5036	0.4717	0.4717
K2PG	37	-3828.7093	7910.9516	0.2287	0.7004
SYMG	41	-3815.3429	7911.6278	0.1631	0.8635
GTRIG	45	-3802.3551	7913.0611	0.0796	0.9431
K2PIG	38	-3827.2708	7914.9269	0.0313	0.9745
SYMIG	42	-3813.7942	7915.3827	0.0249	0.9994
HKYG	40	-3824.5536	7923.1968	0.0005	0.9999
HKYIG	41	-3823.1485	7927.2389	0.0001	1.0000
K2PI	37	-3839.7768	7933.0865	0.0000	1.0000
SYMI	41	-3827.2801	7935.5021	0.0000	1.0000
GTRI	44	-3817.2787	7936.0561	0.0000	1.0000
HKYI	40	-3835.7212	7945.5321	0.0000	1.0000
JC69G	36	-3926.3827	8099.4461	0.0000	1.0000
JC69IG	37	-3925.2959	8104.1248	0.0000	1.0000
K2P	36	-3929.5375	8105.7557	0.0000	1.0000
SYM	40	-3917.0227	8108.1350	0.0000	1.0000
GTR	43	-3908.6268	8111.9000	0.0000	1.0000
F81G	39	-3922.6371	8112.5117	0.0000	1.0000
F81IG	40	-3921.5275	8117.1448	0.0000	1.0000
HKY	39	-3925.1294	8117.4963	0.0000	1.0000
JC69I	36	-3935.5063	8117.6933	0.0000	1.0000
F81I	39	-3931.9503	8131.1380	0.0000	1.0000
JC69	35	-4021.3064	8282.4413	0.0000	1.0000
F81	38	-4017.4334	8295.2520	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_9156.phylip = 1 - 946;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_9156.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_9156.phylip = 1 - 946;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_9156.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_9156.phylip = 1 - 946;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr8_9156.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

