
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:42 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/chr2_18608.phylip" (19 taxa, 768 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	-2322.7562	4730.0240	0.4488	0.4488
HKYI	40	-2323.2445	4731.0007	0.2754	0.7241
HKYIG	41	-2322.7457	4732.2352	0.1485	0.8727
GTRG	44	-2320.1540	4733.7851	0.0684	0.9411
GTRI	44	-2320.7908	4735.0587	0.0362	0.9773
GTRIG	45	-2320.1305	4735.9950	0.0227	1.0000
K2PG	37	-2345.6202	4769.0925	0.0000	1.0000
K2PI	37	-2346.0196	4769.8913	0.0000	1.0000
K2PIG	38	-2345.6206	4771.3070	0.0000	1.0000
SYMG	41	-2344.6136	4775.9710	0.0000	1.0000
SYMI	41	-2344.9281	4776.5999	0.0000	1.0000
SYMIG	42	-2344.6199	4778.2219	0.0000	1.0000
HKY	39	-2349.1509	4780.5874	0.0000	1.0000
GTR	43	-2347.1342	4785.4949	0.0000	1.0000
F81G	39	-2365.4513	4813.1883	0.0000	1.0000
F81I	39	-2365.8014	4813.8885	0.0000	1.0000
F81IG	40	-2365.4457	4815.4031	0.0000	1.0000
K2P	36	-2371.7515	4819.1473	0.0000	1.0000
SYM	40	-2370.6709	4825.8535	0.0000	1.0000
JC69G	36	-2384.0360	4843.7162	0.0000	1.0000
JC69I	36	-2384.4183	4844.4808	0.0000	1.0000
JC69IG	37	-2384.0352	4845.9224	0.0000	1.0000
F81	38	-2390.3351	4860.7360	0.0000	1.0000
JC69	35	-2409.0149	4891.4724	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-2322.7562	4725.5123	0.4071	0.4071
HKYI	40	-2323.2445	4726.4890	0.2498	0.6569
HKYIG	41	-2322.7457	4727.4914	0.1513	0.8083
GTRG	44	-2320.1540	4728.3080	0.1006	0.9089
GTRI	44	-2320.7908	4729.5815	0.0532	0.9621
GTRIG	45	-2320.1305	4730.2609	0.0379	1.0000
K2PG	37	-2345.6202	4765.2404	0.0000	1.0000
K2PI	37	-2346.0196	4766.0392	0.0000	1.0000
K2PIG	38	-2345.6206	4767.2412	0.0000	1.0000
SYMG	41	-2344.6136	4771.2272	0.0000	1.0000
SYMI	41	-2344.9281	4771.8561	0.0000	1.0000
SYMIG	42	-2344.6199	4773.2398	0.0000	1.0000
HKY	39	-2349.1509	4776.3017	0.0000	1.0000
GTR	43	-2347.1342	4780.2684	0.0000	1.0000
F81G	39	-2365.4513	4808.9025	0.0000	1.0000
F81I	39	-2365.8014	4809.6027	0.0000	1.0000
F81IG	40	-2365.4457	4810.8914	0.0000	1.0000
K2P	36	-2371.7515	4815.5030	0.0000	1.0000
SYM	40	-2370.6709	4821.3418	0.0000	1.0000
JC69G	36	-2384.0360	4840.0719	0.0000	1.0000
JC69I	36	-2384.4183	4840.8365	0.0000	1.0000
JC69IG	37	-2384.0352	4842.0704	0.0000	1.0000
F81	38	-2390.3351	4856.6701	0.0000	1.0000
JC69	35	-2409.0149	4888.0298	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2322.7562	4911.2639	0.6060	0.6060
HKYI	40	-2323.2445	4912.2405	0.3719	0.9779
HKYIG	41	-2322.7457	4917.8867	0.0221	1.0000
GTRG	44	-2320.1540	4932.6347	0.0000	1.0000
GTRI	44	-2320.7908	4933.9082	0.0000	1.0000
K2PG	37	-2345.6202	4937.0607	0.0000	1.0000
K2PI	37	-2346.0196	4937.8595	0.0000	1.0000
GTRIG	45	-2320.1305	4939.2314	0.0000	1.0000
K2PIG	38	-2345.6206	4943.7052	0.0000	1.0000
HKY	39	-2349.1509	4957.4095	0.0000	1.0000
SYMG	41	-2344.6136	4961.6226	0.0000	1.0000
SYMI	41	-2344.9281	4962.2515	0.0000	1.0000
SYMIG	42	-2344.6199	4968.2790	0.0000	1.0000
GTR	43	-2347.1342	4979.9513	0.0000	1.0000
K2P	36	-2371.7515	4982.6794	0.0000	1.0000
F81G	39	-2365.4513	4990.0103	0.0000	1.0000
F81I	39	-2365.8014	4990.7105	0.0000	1.0000
F81IG	40	-2365.4457	4996.6430	0.0000	1.0000
SYM	40	-2370.6709	5007.0934	0.0000	1.0000
JC69G	36	-2384.0360	5007.2484	0.0000	1.0000
JC69I	36	-2384.4183	5008.0129	0.0000	1.0000
JC69IG	37	-2384.0352	5013.8906	0.0000	1.0000
F81	38	-2390.3351	5033.1341	0.0000	1.0000
JC69	35	-2409.0149	5050.5624	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/chr2_18608.phylip = 1 - 768;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_18608.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/chr2_18608.phylip = 1 - 768;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_18608.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/chr2_18608.phylip = 1 - 768;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_18608.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

