
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
13:16 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_11272.phylip" (19 taxa, 754 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	40	-2298.4218	4681.4439	0.5464	0.5464
HKYG	40	-2299.0302	4682.6608	0.2974	0.8438
HKYIG	41	-2299.0303	4684.8978	0.0972	0.9409
GTRI	44	-2296.6976	4686.9806	0.0343	0.9752
GTRG	44	-2297.3028	4688.1910	0.0187	0.9939
GTRIG	45	-2297.3028	4690.4531	0.0060	1.0000
K2PG	37	-2313.0545	4704.0363	0.0000	1.0000
K2PI	37	-2313.9068	4705.7410	0.0000	1.0000
K2PIG	38	-2314.4149	4708.9753	0.0000	1.0000
SYMG	41	-2312.7495	4712.3360	0.0000	1.0000
SYMI	41	-2313.4707	4713.7785	0.0000	1.0000
SYMIG	42	-2314.0592	4717.1985	0.0000	1.0000
HKY	39	-2348.6311	4779.6319	0.0000	1.0000
GTR	43	-2347.7441	4786.8177	0.0000	1.0000
K2P	36	-2362.3779	4800.4712	0.0000	1.0000
F81G	39	-2362.2830	4806.9357	0.0000	1.0000
SYM	40	-2361.6027	4807.8057	0.0000	1.0000
F81I	39	-2362.8109	4807.9916	0.0000	1.0000
F81IG	40	-2362.2831	4809.1665	0.0000	1.0000
JC69G	36	-2373.3118	4822.3391	0.0000	1.0000
JC69I	36	-2373.4158	4822.5470	0.0000	1.0000
JC69IG	37	-2374.2418	4826.4109	0.0000	1.0000
F81	38	-2409.7767	4899.6988	0.0000	1.0000
JC69	35	-2420.8500	4915.2098	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-2298.4218	4676.8436	0.5198	0.5198
HKYG	40	-2299.0302	4678.0605	0.2829	0.8027
HKYIG	41	-2299.0303	4680.0607	0.1041	0.9067
GTRI	44	-2296.6976	4681.3953	0.0534	0.9601
GTRG	44	-2297.3028	4682.6056	0.0291	0.9893
GTRIG	45	-2297.3028	4684.6057	0.0107	1.0000
K2PG	37	-2313.0545	4700.1089	0.0000	1.0000
K2PI	37	-2313.9068	4701.8137	0.0000	1.0000
K2PIG	38	-2314.4149	4704.8298	0.0000	1.0000
SYMG	41	-2312.7495	4707.4989	0.0000	1.0000
SYMI	41	-2313.4707	4708.9414	0.0000	1.0000
SYMIG	42	-2314.0592	4712.1183	0.0000	1.0000
HKY	39	-2348.6311	4775.2622	0.0000	1.0000
GTR	43	-2347.7441	4781.4881	0.0000	1.0000
K2P	36	-2362.3779	4796.7557	0.0000	1.0000
F81G	39	-2362.2830	4802.5660	0.0000	1.0000
SYM	40	-2361.6027	4803.2054	0.0000	1.0000
F81I	39	-2362.8109	4803.6219	0.0000	1.0000
F81IG	40	-2362.2831	4804.5663	0.0000	1.0000
JC69G	36	-2373.3118	4818.6236	0.0000	1.0000
JC69I	36	-2373.4158	4818.8315	0.0000	1.0000
JC69IG	37	-2374.2418	4822.4835	0.0000	1.0000
F81	38	-2409.7767	4895.5533	0.0000	1.0000
JC69	35	-2420.8500	4911.7000	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-2298.4218	4861.8593	0.6341	0.6341
HKYG	40	-2299.0302	4863.0762	0.3450	0.9791
HKYIG	41	-2299.0303	4869.7018	0.0126	0.9917
K2PG	37	-2313.0545	4871.2485	0.0058	0.9975
K2PI	37	-2313.9068	4872.9532	0.0025	0.9999
K2PIG	38	-2314.4149	4880.5947	0.0001	1.0000
GTRI	44	-2296.6976	4884.9125	0.0000	1.0000
GTRG	44	-2297.3028	4886.1229	0.0000	1.0000
GTRIG	45	-2297.3028	4892.7483	0.0000	1.0000
SYMG	41	-2312.7495	4897.1400	0.0000	1.0000
SYMI	41	-2313.4707	4898.5825	0.0000	1.0000
SYMIG	42	-2314.0592	4906.3848	0.0000	1.0000
HKY	39	-2348.6311	4955.6525	0.0000	1.0000
K2P	36	-2362.3779	4963.2698	0.0000	1.0000
GTR	43	-2347.7441	4980.3800	0.0000	1.0000
F81G	39	-2362.2830	4982.9563	0.0000	1.0000
F81I	39	-2362.8109	4984.0122	0.0000	1.0000
JC69G	36	-2373.3118	4985.1377	0.0000	1.0000
JC69I	36	-2373.4158	4985.3456	0.0000	1.0000
SYM	40	-2361.6027	4988.2211	0.0000	1.0000
F81IG	40	-2362.2831	4989.5820	0.0000	1.0000
JC69IG	37	-2374.2418	4993.6230	0.0000	1.0000
F81	38	-2409.7767	5071.3182	0.0000	1.0000
JC69	35	-2420.8500	5073.5888	0.0000	1.0000



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

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

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


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

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


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

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


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

