
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13: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/chr7_5167.phylip" (19 taxa, 664 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	-1836.2408	3766.8791	0.4910	0.4910
GTRI	44	-1836.5698	3767.5369	0.3534	0.8444
GTRIG	45	-1836.2407	3769.1805	0.1554	0.9998
SYMG	41	-1848.1593	3783.8555	0.0001	0.9999
SYMI	41	-1848.3803	3784.2975	0.0001	1.0000
SYMIG	42	-1848.1617	3786.1398	0.0000	1.0000
HKYG	40	-1853.3511	3791.9671	0.0000	1.0000
HKYI	40	-1853.8228	3792.9105	0.0000	1.0000
HKYIG	41	-1853.2317	3794.0004	0.0000	1.0000
K2PG	37	-1861.9669	3802.4257	0.0000	1.0000
K2PI	37	-1862.3844	3803.2608	0.0000	1.0000
K2PIG	38	-1861.8687	3804.4797	0.0000	1.0000
GTR	43	-1866.1173	3824.3378	0.0000	1.0000
SYM	40	-1877.8272	3840.9193	0.0000	1.0000
HKY	39	-1884.9195	3852.8390	0.0000	1.0000
K2P	36	-1893.7988	3863.8464	0.0000	1.0000
F81G	39	-1892.6074	3868.2149	0.0000	1.0000
F81I	39	-1892.8670	3868.7339	0.0000	1.0000
F81IG	40	-1892.6016	3870.4680	0.0000	1.0000
JC69G	36	-1899.8287	3875.9062	0.0000	1.0000
JC69I	36	-1900.0864	3876.4217	0.0000	1.0000
JC69IG	37	-1899.8313	3878.1547	0.0000	1.0000
F81	38	-1922.3883	3925.5190	0.0000	1.0000
JC69	35	-1930.0816	3934.1759	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1836.2408	3760.4817	0.4789	0.4789
GTRI	44	-1836.5698	3761.1395	0.3447	0.8236
GTRIG	45	-1836.2407	3762.4814	0.1762	0.9999
SYMG	41	-1848.1593	3778.3186	0.0001	0.9999
SYMI	41	-1848.3803	3778.7605	0.0001	1.0000
SYMIG	42	-1848.1617	3780.3233	0.0000	1.0000
HKYG	40	-1853.3511	3786.7023	0.0000	1.0000
HKYI	40	-1853.8228	3787.6456	0.0000	1.0000
HKYIG	41	-1853.2317	3788.4634	0.0000	1.0000
K2PG	37	-1861.9669	3797.9337	0.0000	1.0000
K2PI	37	-1862.3844	3798.7688	0.0000	1.0000
K2PIG	38	-1861.8687	3799.7373	0.0000	1.0000
GTR	43	-1866.1173	3818.2346	0.0000	1.0000
SYM	40	-1877.8272	3835.6544	0.0000	1.0000
HKY	39	-1884.9195	3847.8390	0.0000	1.0000
K2P	36	-1893.7988	3859.5976	0.0000	1.0000
F81G	39	-1892.6074	3863.2149	0.0000	1.0000
F81I	39	-1892.8670	3863.7339	0.0000	1.0000
F81IG	40	-1892.6016	3865.2032	0.0000	1.0000
JC69G	36	-1899.8287	3871.6574	0.0000	1.0000
JC69I	36	-1900.0864	3872.1729	0.0000	1.0000
JC69IG	37	-1899.8313	3873.6627	0.0000	1.0000
F81	38	-1922.3883	3920.7766	0.0000	1.0000
JC69	35	-1930.0816	3930.1631	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-1836.2408	3958.4061	0.4803	0.4803
GTRI	44	-1836.5698	3959.0639	0.3457	0.8260
SYMG	41	-1848.1593	3962.7481	0.0548	0.8808
SYMI	41	-1848.3803	3963.1901	0.0439	0.9247
K2PG	37	-1861.9669	3964.3701	0.0243	0.9491
GTRIG	45	-1836.2407	3964.9041	0.0186	0.9677
K2PI	37	-1862.3844	3965.2052	0.0160	0.9837
HKYG	40	-1853.3511	3966.6336	0.0079	0.9916
HKYI	40	-1853.8228	3967.5769	0.0049	0.9965
SYMIG	42	-1848.1617	3969.2512	0.0021	0.9986
K2PIG	38	-1861.8687	3970.6721	0.0010	0.9997
HKYIG	41	-1853.2317	3972.8929	0.0003	1.0000
GTR	43	-1866.1173	4011.6607	0.0000	1.0000
SYM	40	-1877.8272	4015.5857	0.0000	1.0000
K2P	36	-1893.7988	4021.5357	0.0000	1.0000
HKY	39	-1884.9195	4023.2720	0.0000	1.0000
JC69G	36	-1899.8287	4033.5956	0.0000	1.0000
JC69I	36	-1900.0864	4034.1110	0.0000	1.0000
F81G	39	-1892.6074	4038.6479	0.0000	1.0000
F81I	39	-1892.8670	4039.1669	0.0000	1.0000
JC69IG	37	-1899.8313	4040.0991	0.0000	1.0000
F81IG	40	-1892.6016	4045.1345	0.0000	1.0000
JC69	35	-1930.0816	4087.6030	0.0000	1.0000
F81	38	-1922.3883	4091.7113	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_5167.phylip = 1 - 664;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_5167.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_5167.phylip = 1 - 664;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_5167.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_5167.phylip = 1 - 664;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_5167.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

