
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:29 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/chrZ_8159.phylip" (19 taxa, 436 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	-1318.2422	2724.7882	0.4467	0.4467
HKYIG	41	-1317.2578	2725.2568	0.3534	0.8001
HKYI	40	-1319.3499	2727.0037	0.1476	0.9477
GTRG	44	-1315.9383	2730.0045	0.0329	0.9806
GTRIG	45	-1315.8908	2732.3969	0.0100	0.9906
GTRI	44	-1317.1905	2732.5089	0.0094	1.0000
HKY	39	-1329.1861	2744.2511	0.0000	1.0000
GTR	43	-1328.4341	2752.5213	0.0000	1.0000
K2PG	37	-1337.7751	2756.6155	0.0000	1.0000
K2PI	37	-1338.6836	2758.4326	0.0000	1.0000
K2PIG	38	-1337.7394	2758.9448	0.0000	1.0000
SYMG	41	-1336.8209	2764.3829	0.0000	1.0000
SYMI	41	-1337.7279	2766.1969	0.0000	1.0000
SYMIG	42	-1336.7805	2766.7519	0.0000	1.0000
K2P	36	-1349.0434	2776.7634	0.0000	1.0000
SYM	40	-1348.2741	2784.8519	0.0000	1.0000
F81G	39	-1360.2294	2806.3376	0.0000	1.0000
F81IG	40	-1360.2033	2808.7103	0.0000	1.0000
F81I	39	-1363.1689	2812.2165	0.0000	1.0000
F81	38	-1372.4519	2828.3699	0.0000	1.0000
JC69G	36	-1378.4408	2835.5583	0.0000	1.0000
JC69I	36	-1378.9354	2836.5475	0.0000	1.0000
JC69IG	37	-1378.4361	2837.9376	0.0000	1.0000
JC69	35	-1388.6123	2853.5245	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1318.2422	2716.4844	0.3815	0.3815
HKYIG	41	-1317.2578	2716.5157	0.3755	0.7570
HKYI	40	-1319.3499	2718.6999	0.1260	0.8830
GTRG	44	-1315.9383	2719.8766	0.0700	0.9530
GTRIG	45	-1315.8908	2721.7815	0.0270	0.9800
GTRI	44	-1317.1905	2722.3811	0.0200	1.0000
HKY	39	-1329.1861	2736.3723	0.0000	1.0000
GTR	43	-1328.4341	2742.8682	0.0000	1.0000
K2PG	37	-1337.7751	2749.5501	0.0000	1.0000
K2PI	37	-1338.6836	2751.3672	0.0000	1.0000
K2PIG	38	-1337.7394	2751.4788	0.0000	1.0000
SYMG	41	-1336.8209	2755.6417	0.0000	1.0000
SYMI	41	-1337.7279	2757.4558	0.0000	1.0000
SYMIG	42	-1336.7805	2757.5611	0.0000	1.0000
K2P	36	-1349.0434	2770.0867	0.0000	1.0000
SYM	40	-1348.2741	2776.5481	0.0000	1.0000
F81G	39	-1360.2294	2798.4588	0.0000	1.0000
F81IG	40	-1360.2033	2800.4065	0.0000	1.0000
F81I	39	-1363.1689	2804.3377	0.0000	1.0000
F81	38	-1372.4519	2820.9039	0.0000	1.0000
JC69G	36	-1378.4408	2828.8816	0.0000	1.0000
JC69I	36	-1378.9354	2829.8708	0.0000	1.0000
JC69IG	37	-1378.4361	2830.8722	0.0000	1.0000
JC69	35	-1388.6123	2847.2245	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-1318.2422	2879.5901	0.6854	0.6854
HKYI	40	-1319.3499	2881.8055	0.2264	0.9118
HKYIG	41	-1317.2578	2883.6990	0.0878	0.9997
HKY	39	-1329.1861	2895.4003	0.0003	0.9999
GTRG	44	-1315.9383	2899.2929	0.0000	1.0000
K2PG	37	-1337.7751	2900.4229	0.0000	1.0000
GTRI	44	-1317.1905	2901.7973	0.0000	1.0000
K2PI	37	-1338.6836	2902.2400	0.0000	1.0000
GTRIG	45	-1315.8908	2905.2754	0.0000	1.0000
K2PIG	38	-1337.7394	2906.4292	0.0000	1.0000
K2P	36	-1349.0434	2916.8819	0.0000	1.0000
GTR	43	-1328.4341	2918.2068	0.0000	1.0000
SYMG	41	-1336.8209	2922.8251	0.0000	1.0000
SYMI	41	-1337.7279	2924.6392	0.0000	1.0000
SYMIG	42	-1336.7805	2928.8221	0.0000	1.0000
SYM	40	-1348.2741	2939.6538	0.0000	1.0000
F81G	39	-1360.2294	2957.4869	0.0000	1.0000
F81I	39	-1363.1689	2963.3658	0.0000	1.0000
F81IG	40	-1360.2033	2963.5122	0.0000	1.0000
JC69G	36	-1378.4408	2975.6767	0.0000	1.0000
F81	38	-1372.4519	2975.8543	0.0000	1.0000
JC69I	36	-1378.9354	2976.6659	0.0000	1.0000
JC69IG	37	-1378.4361	2981.7450	0.0000	1.0000
JC69	35	-1388.6123	2989.9420	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/chrZ_8159.phylip = 1 - 436;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_8159.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/chrZ_8159.phylip = 1 - 436;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_8159.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/chrZ_8159.phylip = 1 - 436;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_8159.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

