
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:45 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/chr22_189.phylip" (19 taxa, 503 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	-1667.5426	3431.7314	0.6161	0.6161
GTRIG	45	-1667.0971	3433.2533	0.2879	0.9039
GTRI	44	-1669.4062	3435.4588	0.0956	0.9995
SYMG	41	-1678.7443	3446.9593	0.0003	0.9998
SYMIG	42	-1678.3709	3448.5939	0.0001	0.9999
SYMI	41	-1680.2423	3449.9554	0.0001	1.0000
HKYG	40	-1683.3283	3453.7562	0.0000	1.0000
HKYIG	41	-1682.6456	3454.7620	0.0000	1.0000
HKYI	40	-1685.7050	3458.5095	0.0000	1.0000
K2PG	37	-1689.6061	3459.2596	0.0000	1.0000
K2PIG	38	-1688.9832	3460.3543	0.0000	1.0000
K2PI	37	-1691.8781	3463.8034	0.0000	1.0000
GTR	43	-1703.7121	3501.6681	0.0000	1.0000
SYM	40	-1712.4401	3511.9797	0.0000	1.0000
F81G	39	-1716.6280	3517.9946	0.0000	1.0000
F81IG	40	-1716.0394	3519.1783	0.0000	1.0000
JC69G	36	-1721.1432	3520.0031	0.0000	1.0000
HKY	39	-1718.0902	3520.9190	0.0000	1.0000
JC69IG	37	-1720.5410	3521.1294	0.0000	1.0000
F81I	39	-1718.8640	3522.4667	0.0000	1.0000
JC69I	36	-1723.4485	3524.6138	0.0000	1.0000
K2P	36	-1723.5610	3524.8387	0.0000	1.0000
F81	38	-1749.5666	3581.5211	0.0000	1.0000
JC69	35	-1753.9693	3583.3347	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-1667.5426	3423.0852	0.5781	0.5781
GTRIG	45	-1667.0971	3424.1942	0.3320	0.9101
GTRI	44	-1669.4062	3426.8125	0.0897	0.9997
SYMG	41	-1678.7443	3439.4886	0.0002	0.9999
SYMIG	42	-1678.3709	3440.7417	0.0001	1.0000
SYMI	41	-1680.2423	3442.4847	0.0000	1.0000
HKYG	40	-1683.3283	3446.6566	0.0000	1.0000
HKYIG	41	-1682.6456	3447.2913	0.0000	1.0000
HKYI	40	-1685.7050	3451.4100	0.0000	1.0000
K2PG	37	-1689.6061	3453.2123	0.0000	1.0000
K2PIG	38	-1688.9832	3453.9664	0.0000	1.0000
K2PI	37	-1691.8781	3457.7561	0.0000	1.0000
GTR	43	-1703.7121	3493.4241	0.0000	1.0000
SYM	40	-1712.4401	3504.8802	0.0000	1.0000
F81G	39	-1716.6280	3511.2560	0.0000	1.0000
F81IG	40	-1716.0394	3512.0787	0.0000	1.0000
HKY	39	-1718.0902	3514.1803	0.0000	1.0000
JC69G	36	-1721.1432	3514.2863	0.0000	1.0000
JC69IG	37	-1720.5410	3515.0821	0.0000	1.0000
F81I	39	-1718.8640	3515.7281	0.0000	1.0000
JC69I	36	-1723.4485	3518.8970	0.0000	1.0000
K2P	36	-1723.5610	3519.1220	0.0000	1.0000
F81	38	-1749.5666	3575.1332	0.0000	1.0000
JC69	35	-1753.9693	3577.9386	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-1667.5426	3608.7911	0.4253	0.4253
K2PG	37	-1689.6061	3609.3741	0.3178	0.7431
GTRI	44	-1669.4062	3612.5184	0.0660	0.8091
SYMG	41	-1678.7443	3612.5328	0.0655	0.8746
K2PI	37	-1691.8781	3613.9179	0.0328	0.9074
GTRIG	45	-1667.0971	3614.1207	0.0296	0.9370
K2PIG	38	-1688.9832	3614.3488	0.0264	0.9634
HKYG	40	-1683.3283	3615.4802	0.0150	0.9784
SYMI	41	-1680.2423	3615.5289	0.0146	0.9930
SYMIG	42	-1678.3709	3618.0065	0.0042	0.9973
HKYI	40	-1685.7050	3620.2336	0.0014	0.9987
HKYIG	41	-1682.6456	3620.3355	0.0013	1.0000
JC69G	36	-1721.1432	3666.2276	0.0000	1.0000
JC69I	36	-1723.4485	3670.8383	0.0000	1.0000
K2P	36	-1723.5610	3671.0632	0.0000	1.0000
JC69IG	37	-1720.5410	3671.2439	0.0000	1.0000
SYM	40	-1712.4401	3673.7038	0.0000	1.0000
GTR	43	-1703.7121	3674.9095	0.0000	1.0000
F81G	39	-1716.6280	3675.8590	0.0000	1.0000
HKY	39	-1718.0902	3678.7834	0.0000	1.0000
F81I	39	-1718.8640	3680.3311	0.0000	1.0000
F81IG	40	-1716.0394	3680.9023	0.0000	1.0000
JC69	35	-1753.9693	3725.6593	0.0000	1.0000
F81	38	-1749.5666	3735.5156	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/chr22_189.phylip = 1 - 503;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr22_189.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/chr22_189.phylip = 1 - 503;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr22_189.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/chr22_189.phylip = 1 - 503;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr22_189.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

