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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	40	-3750.8344	7586.2244	0.6708	0.6708
HKYIG	41	-3750.4397	7587.6694	0.3257	0.9965
GTRG	44	-3752.1415	7597.8137	0.0020	0.9986
GTRIG	45	-3751.6271	7599.0444	0.0011	0.9997
HKYI	40	-3758.5341	7601.6238	0.0003	1.0000
GTRI	44	-3756.9635	7607.4576	0.0000	1.0000
K2PG	37	-3772.6405	7623.1703	0.0000	1.0000
K2PIG	38	-3772.2941	7624.6934	0.0000	1.0000
SYMG	41	-3771.6782	7630.1464	0.0000	1.0000
SYMIG	42	-3771.3532	7631.7371	0.0000	1.0000
K2PI	37	-3780.3710	7638.6314	0.0000	1.0000
SYMI	41	-3779.6313	7646.0526	0.0000	1.0000
HKY	39	-3828.0876	7738.5026	0.0000	1.0000
GTR	43	-3824.4203	7740.1182	0.0000	1.0000
K2P	36	-3846.6975	7769.0746	0.0000	1.0000
SYM	40	-3843.3171	7771.1898	0.0000	1.0000
F81G	39	-3881.0805	7844.4884	0.0000	1.0000
F81IG	40	-3880.5276	7845.6107	0.0000	1.0000
F81I	39	-3886.6448	7855.6170	0.0000	1.0000
JC69G	36	-3892.7332	7861.1460	0.0000	1.0000
JC69IG	37	-3892.1397	7862.1688	0.0000	1.0000
JC69I	36	-3898.3996	7872.4787	0.0000	1.0000
F81	38	-3944.9271	7969.9594	0.0000	1.0000
JC69	35	-3955.8910	7985.2579	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-3750.8344	7581.6689	0.6433	0.6433
HKYIG	41	-3750.4397	7582.8794	0.3512	0.9945
GTRG	44	-3752.1415	7592.2830	0.0032	0.9977
GTRIG	45	-3751.6271	7593.2542	0.0020	0.9997
HKYI	40	-3758.5341	7597.0682	0.0003	1.0000
GTRI	44	-3756.9635	7601.9269	0.0000	1.0000
K2PG	37	-3772.6405	7619.2809	0.0000	1.0000
K2PIG	38	-3772.2941	7620.5881	0.0000	1.0000
SYMG	41	-3771.6782	7625.3565	0.0000	1.0000
SYMIG	42	-3771.3532	7626.7064	0.0000	1.0000
K2PI	37	-3780.3710	7634.7420	0.0000	1.0000
SYMI	41	-3779.6313	7641.2626	0.0000	1.0000
HKY	39	-3828.0876	7734.1753	0.0000	1.0000
GTR	43	-3824.4203	7734.8407	0.0000	1.0000
K2P	36	-3846.6975	7765.3951	0.0000	1.0000
SYM	40	-3843.3171	7766.6342	0.0000	1.0000
F81G	39	-3881.0805	7840.1611	0.0000	1.0000
F81IG	40	-3880.5276	7841.0551	0.0000	1.0000
F81I	39	-3886.6448	7851.2896	0.0000	1.0000
JC69G	36	-3892.7332	7857.4665	0.0000	1.0000
JC69IG	37	-3892.1397	7858.2795	0.0000	1.0000
JC69I	36	-3898.3996	7868.7991	0.0000	1.0000
F81	38	-3944.9271	7965.8541	0.0000	1.0000
JC69	35	-3955.8910	7981.7820	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-3750.8344	7767.0542	0.9485	0.9485
HKYIG	41	-3750.4397	7772.8993	0.0510	0.9996
HKYI	40	-3758.5341	7782.4535	0.0004	1.0000
K2PG	37	-3772.6405	7790.7624	0.0000	1.0000
GTRG	44	-3752.1415	7796.2069	0.0000	1.0000
K2PIG	38	-3772.2941	7796.7042	0.0000	1.0000
GTRIG	45	-3751.6271	7801.8127	0.0000	1.0000
GTRI	44	-3756.9635	7805.8508	0.0000	1.0000
K2PI	37	-3780.3710	7806.2235	0.0000	1.0000
SYMG	41	-3771.6782	7815.3764	0.0000	1.0000
SYMIG	42	-3771.3532	7821.3610	0.0000	1.0000
SYMI	41	-3779.6313	7831.2826	0.0000	1.0000
HKY	39	-3828.0876	7914.9260	0.0000	1.0000
K2P	36	-3846.6975	7932.2419	0.0000	1.0000
GTR	43	-3824.4203	7934.1299	0.0000	1.0000
SYM	40	-3843.3171	7952.0196	0.0000	1.0000
F81G	39	-3881.0805	8020.9118	0.0000	1.0000
JC69G	36	-3892.7332	8024.3133	0.0000	1.0000
F81IG	40	-3880.5276	8026.4405	0.0000	1.0000
JC69IG	37	-3892.1397	8029.7609	0.0000	1.0000
F81I	39	-3886.6448	8032.0403	0.0000	1.0000
JC69I	36	-3898.3996	8035.6459	0.0000	1.0000
F81	38	-3944.9271	8141.9702	0.0000	1.0000
JC69	35	-3955.8910	8143.9942	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_11854.phylip = 1 - 761;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_11854.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_11854.phylip = 1 - 761;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_11854.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (HKYG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_11854.phylip = 1 - 761;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_11854.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

