
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:40 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/chr6_9631.phylip" (19 taxa, 558 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	-1856.1535	3810.3930	0.8009	0.8009
GTRG	44	-1858.7412	3813.2017	0.1966	0.9975
GTRI	44	-1863.3605	3822.4403	0.0019	0.9994
HKYIG	41	-1868.2490	3825.1724	0.0005	0.9999
HKYG	40	-1871.4393	3829.2229	0.0001	1.0000
HKYI	40	-1875.2755	3836.8953	0.0000	1.0000
K2PIG	38	-1878.5357	3838.7825	0.0000	1.0000
SYMIG	42	-1873.9587	3838.9310	0.0000	1.0000
SYMG	41	-1876.7311	3842.1367	0.0000	1.0000
K2PG	37	-1881.7360	3842.8797	0.0000	1.0000
SYMI	41	-1880.2967	3849.2678	0.0000	1.0000
K2PI	37	-1885.3375	3850.0828	0.0000	1.0000
F81IG	40	-1907.0270	3900.3982	0.0000	1.0000
F81G	39	-1909.6380	3903.2991	0.0000	1.0000
JC69IG	37	-1914.6600	3908.7278	0.0000	1.0000
F81I	39	-1913.1296	3910.2823	0.0000	1.0000
JC69G	36	-1917.3672	3911.8476	0.0000	1.0000
JC69I	36	-1920.8483	3918.8099	0.0000	1.0000
GTR	43	-1921.4094	3936.1806	0.0000	1.0000
HKY	39	-1932.5432	3949.1095	0.0000	1.0000
SYM	40	-1936.5299	3959.4041	0.0000	1.0000
K2P	36	-1942.5909	3962.2950	0.0000	1.0000
F81	38	-1968.3079	4018.3267	0.0000	1.0000
JC69	35	-1976.2357	4027.2990	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-1856.1535	3802.3070	0.8287	0.8287
GTRG	44	-1858.7412	3805.4824	0.1694	0.9980
GTRI	44	-1863.3605	3814.7210	0.0017	0.9997
HKYIG	41	-1868.2490	3818.4980	0.0003	1.0000
HKYG	40	-1871.4393	3822.8786	0.0000	1.0000
HKYI	40	-1875.2755	3830.5510	0.0000	1.0000
SYMIG	42	-1873.9587	3831.9174	0.0000	1.0000
K2PIG	38	-1878.5357	3833.0715	0.0000	1.0000
SYMG	41	-1876.7311	3835.4622	0.0000	1.0000
K2PG	37	-1881.7360	3837.4720	0.0000	1.0000
SYMI	41	-1880.2967	3842.5933	0.0000	1.0000
K2PI	37	-1885.3375	3844.6751	0.0000	1.0000
F81IG	40	-1907.0270	3894.0539	0.0000	1.0000
F81G	39	-1909.6380	3897.2759	0.0000	1.0000
JC69IG	37	-1914.6600	3903.3201	0.0000	1.0000
F81I	39	-1913.1296	3904.2591	0.0000	1.0000
JC69G	36	-1917.3672	3906.7344	0.0000	1.0000
JC69I	36	-1920.8483	3913.6967	0.0000	1.0000
GTR	43	-1921.4094	3928.8188	0.0000	1.0000
HKY	39	-1932.5432	3943.0864	0.0000	1.0000
SYM	40	-1936.5299	3953.0598	0.0000	1.0000
K2P	36	-1942.5909	3957.1818	0.0000	1.0000
F81	38	-1968.3079	4012.6157	0.0000	1.0000
JC69	35	-1976.2357	4022.4714	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-1858.7412	3995.7542	0.2273	0.2273
HKYIG	41	-1868.2490	3995.7967	0.2225	0.4498
HKYG	40	-1871.4393	3995.8530	0.2164	0.6662
GTRIG	45	-1856.1535	3996.9032	0.1280	0.7942
K2PIG	38	-1878.5357	3997.3971	0.1000	0.8942
K2PG	37	-1881.7360	3997.4733	0.0962	0.9904
HKYI	40	-1875.2755	4003.5254	0.0047	0.9951
K2PI	37	-1885.3375	4004.6764	0.0026	0.9977
GTRI	44	-1863.3605	4004.9928	0.0022	0.9999
SYMG	41	-1876.7311	4012.7610	0.0000	1.0000
SYMIG	42	-1873.9587	4013.5405	0.0000	1.0000
SYMI	41	-1880.2967	4019.8921	0.0000	1.0000
JC69G	36	-1917.3672	4062.4113	0.0000	1.0000
JC69IG	37	-1914.6600	4063.3214	0.0000	1.0000
F81G	39	-1909.6380	4065.9259	0.0000	1.0000
F81IG	40	-1907.0270	4067.0283	0.0000	1.0000
JC69I	36	-1920.8483	4069.3736	0.0000	1.0000
F81I	39	-1913.1296	4072.9091	0.0000	1.0000
HKY	39	-1932.5432	4111.7364	0.0000	1.0000
K2P	36	-1942.5909	4112.8587	0.0000	1.0000
GTR	43	-1921.4094	4114.7662	0.0000	1.0000
SYM	40	-1936.5299	4126.0342	0.0000	1.0000
JC69	35	-1976.2357	4173.8240	0.0000	1.0000
F81	38	-1968.3079	4176.9414	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/chr6_9631.phylip = 1 - 558;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_9631.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/chr6_9631.phylip = 1 - 558;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_9631.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/chr6_9631.phylip = 1 - 558;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_9631.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

