
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:30 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/chr5_11293.phylip" (19 taxa, 864 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	-3222.6592	6538.1536	0.7167	0.7167
GTRIG	45	-3222.7350	6540.5311	0.2183	0.9351
GTRI	44	-3225.0605	6542.9561	0.0649	1.0000
HKYG	40	-3246.2647	6576.5148	0.0000	1.0000
HKYIG	41	-3246.2647	6578.7191	0.0000	1.0000
HKYI	40	-3249.8108	6583.6071	0.0000	1.0000
GTR	43	-3269.2167	6629.0480	0.0000	1.0000
SYMG	41	-3274.1402	6634.4702	0.0000	1.0000
SYMIG	42	-3274.0943	6636.5880	0.0000	1.0000
SYMI	41	-3275.3027	6636.7952	0.0000	1.0000
K2PG	37	-3281.1960	6639.7964	0.0000	1.0000
K2PIG	38	-3281.1960	6641.9848	0.0000	1.0000
K2PI	37	-3282.7941	6642.9925	0.0000	1.0000
HKY	39	-3294.2088	6670.2040	0.0000	1.0000
SYM	40	-3323.0619	6730.1093	0.0000	1.0000
K2P	36	-3333.0184	6741.2581	0.0000	1.0000
F81G	39	-3362.7754	6807.3372	0.0000	1.0000
F81IG	40	-3361.7628	6807.5109	0.0000	1.0000
F81I	39	-3363.8624	6809.5112	0.0000	1.0000
JC69G	36	-3389.2124	6853.6461	0.0000	1.0000
JC69IG	37	-3389.2123	6855.8291	0.0000	1.0000
JC69I	36	-3390.3227	6855.8667	0.0000	1.0000
F81	38	-3406.8306	6893.2539	0.0000	1.0000
JC69	35	-3435.4926	6944.0286	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	44	-3222.6592	6533.3185	0.6985	0.6985
GTRIG	45	-3222.7350	6535.4700	0.2382	0.9367
GTRI	44	-3225.0605	6538.1209	0.0633	1.0000
HKYG	40	-3246.2647	6572.5294	0.0000	1.0000
HKYIG	41	-3246.2647	6574.5293	0.0000	1.0000
HKYI	40	-3249.8108	6579.6217	0.0000	1.0000
GTR	43	-3269.2167	6624.4333	0.0000	1.0000
SYMG	41	-3274.1402	6630.2804	0.0000	1.0000
SYMIG	42	-3274.0943	6632.1885	0.0000	1.0000
SYMI	41	-3275.3027	6632.6054	0.0000	1.0000
K2PG	37	-3281.1960	6636.3921	0.0000	1.0000
K2PIG	38	-3281.1960	6638.3921	0.0000	1.0000
K2PI	37	-3282.7941	6639.5881	0.0000	1.0000
HKY	39	-3294.2088	6666.4176	0.0000	1.0000
SYM	40	-3323.0619	6726.1238	0.0000	1.0000
K2P	36	-3333.0184	6738.0368	0.0000	1.0000
F81IG	40	-3361.7628	6803.5255	0.0000	1.0000
F81G	39	-3362.7754	6803.5508	0.0000	1.0000
F81I	39	-3363.8624	6805.7248	0.0000	1.0000
JC69G	36	-3389.2124	6850.4248	0.0000	1.0000
JC69IG	37	-3389.2123	6852.4247	0.0000	1.0000
JC69I	36	-3390.3227	6852.6454	0.0000	1.0000
F81	38	-3406.8306	6889.6612	0.0000	1.0000
JC69	35	-3435.4926	6940.9851	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	44	-3222.6592	6742.8277	0.8911	0.8911
GTRI	44	-3225.0605	6747.6301	0.0807	0.9719
GTRIG	45	-3222.7350	6749.7408	0.0281	1.0000
HKYG	40	-3246.2647	6762.9923	0.0000	1.0000
HKYIG	41	-3246.2647	6769.7538	0.0000	1.0000
HKYI	40	-3249.8108	6770.0846	0.0000	1.0000
K2PG	37	-3281.1960	6812.5703	0.0000	1.0000
K2PI	37	-3282.7941	6815.7663	0.0000	1.0000
K2PIG	38	-3281.1960	6819.3318	0.0000	1.0000
SYMG	41	-3274.1402	6825.5049	0.0000	1.0000
SYMI	41	-3275.3027	6827.8299	0.0000	1.0000
GTR	43	-3269.2167	6829.1809	0.0000	1.0000
SYMIG	42	-3274.0943	6832.1746	0.0000	1.0000
HKY	39	-3294.2088	6852.1189	0.0000	1.0000
K2P	36	-3333.0184	6909.4534	0.0000	1.0000
SYM	40	-3323.0619	6916.5868	0.0000	1.0000
F81G	39	-3362.7754	6989.2521	0.0000	1.0000
F81I	39	-3363.8624	6991.4261	0.0000	1.0000
F81IG	40	-3361.7628	6993.9884	0.0000	1.0000
JC69G	36	-3389.2124	7021.8414	0.0000	1.0000
JC69I	36	-3390.3227	7024.0620	0.0000	1.0000
JC69IG	37	-3389.2123	7028.6029	0.0000	1.0000
F81	38	-3406.8306	7070.6009	0.0000	1.0000
JC69	35	-3435.4926	7107.6402	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/chr5_11293.phylip = 1 - 864;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11293.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/chr5_11293.phylip = 1 - 864;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11293.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/chr5_11293.phylip = 1 - 864;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr5_11293.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

