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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	45	-2501.1362	5098.4607	0.6177	0.6177
GTRG	44	-2502.7732	5099.4568	0.3754	0.9931
SYMIG	42	-2509.6891	5108.7532	0.0036	0.9967
SYMG	41	-2511.2241	5109.5656	0.0024	0.9991
HKYG	40	-2513.4692	5111.8048	0.0008	0.9998
HKYIG	41	-2514.0272	5115.1717	0.0001	1.0000
GTRI	44	-2513.5291	5120.9686	0.0000	1.0000
SYMI	41	-2521.6318	5130.3810	0.0000	1.0000
K2PIG	38	-2525.6754	5131.7355	0.0000	1.0000
HKYI	40	-2523.6411	5132.1487	0.0000	1.0000
K2PG	37	-2527.2295	5132.6126	0.0000	1.0000
K2PI	37	-2535.2578	5148.6692	0.0000	1.0000
F81G	39	-2551.2715	5185.1651	0.0000	1.0000
F81IG	40	-2555.8674	5196.6012	0.0000	1.0000
F81I	39	-2557.7548	5198.1318	0.0000	1.0000
JC69IG	37	-2573.5459	5225.2454	0.0000	1.0000
JC69G	36	-2574.8209	5225.5711	0.0000	1.0000
JC69I	36	-2581.9595	5239.8482	0.0000	1.0000
GTR	43	-2583.4878	5258.6149	0.0000	1.0000
SYM	40	-2591.1529	5267.1722	0.0000	1.0000
HKY	39	-2596.2223	5275.0668	0.0000	1.0000
K2P	36	-2607.7656	5291.4603	0.0000	1.0000
F81	38	-2622.4122	5325.2090	0.0000	1.0000
JC69	35	-2649.6929	5373.0971	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	45	-2501.1362	5092.2724	0.6511	0.6511
GTRG	44	-2502.7732	5093.5464	0.3444	0.9955
SYMIG	42	-2509.6891	5103.3782	0.0025	0.9980
SYMG	41	-2511.2241	5104.4482	0.0015	0.9995
HKYG	40	-2513.4692	5106.9384	0.0004	0.9999
HKYIG	41	-2514.0272	5110.0543	0.0001	1.0000
GTRI	44	-2513.5291	5115.0582	0.0000	1.0000
SYMI	41	-2521.6318	5125.2636	0.0000	1.0000
HKYI	40	-2523.6411	5127.2822	0.0000	1.0000
K2PIG	38	-2525.6754	5127.3509	0.0000	1.0000
K2PG	37	-2527.2295	5128.4590	0.0000	1.0000
K2PI	37	-2535.2578	5144.5156	0.0000	1.0000
F81G	39	-2551.2715	5180.5429	0.0000	1.0000
F81IG	40	-2555.8674	5191.7348	0.0000	1.0000
F81I	39	-2557.7548	5193.5096	0.0000	1.0000
JC69IG	37	-2573.5459	5221.0918	0.0000	1.0000
JC69G	36	-2574.8209	5221.6419	0.0000	1.0000
JC69I	36	-2581.9595	5235.9190	0.0000	1.0000
GTR	43	-2583.4878	5252.9755	0.0000	1.0000
SYM	40	-2591.1529	5262.3058	0.0000	1.0000
HKY	39	-2596.2223	5270.4446	0.0000	1.0000
K2P	36	-2607.7656	5287.5311	0.0000	1.0000
F81	38	-2622.4122	5320.8244	0.0000	1.0000
JC69	35	-2649.6929	5369.3857	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-2513.4692	5289.8297	0.6399	0.6399
SYMG	41	-2511.2241	5291.9118	0.2260	0.8659
GTRG	44	-2502.7732	5294.7268	0.0553	0.9212
SYMIG	42	-2509.6891	5295.4141	0.0392	0.9604
HKYIG	41	-2514.0272	5297.5179	0.0137	0.9741
K2PG	37	-2527.2295	5297.6334	0.0129	0.9870
GTRIG	45	-2501.1362	5298.0251	0.0106	0.9977
K2PIG	38	-2525.6754	5301.0976	0.0023	1.0000
HKYI	40	-2523.6411	5310.1735	0.0000	1.0000
SYMI	41	-2521.6318	5312.7272	0.0000	1.0000
K2PI	37	-2535.2578	5313.6900	0.0000	1.0000
GTRI	44	-2513.5291	5316.2386	0.0000	1.0000
F81G	39	-2551.2715	5358.8619	0.0000	1.0000
F81I	39	-2557.7548	5371.8286	0.0000	1.0000
F81IG	40	-2555.8674	5374.6261	0.0000	1.0000
JC69G	36	-2574.8209	5386.2440	0.0000	1.0000
JC69IG	37	-2573.5459	5390.2663	0.0000	1.0000
JC69I	36	-2581.9595	5400.5212	0.0000	1.0000
SYM	40	-2591.1529	5445.1971	0.0000	1.0000
HKY	39	-2596.2223	5448.7636	0.0000	1.0000
GTR	43	-2583.4878	5449.5837	0.0000	1.0000
K2P	36	-2607.7656	5452.1333	0.0000	1.0000
F81	38	-2622.4122	5494.5711	0.0000	1.0000
JC69	35	-2649.6929	5529.4156	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6735.phylip = 1 - 715;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6735.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6735.phylip = 1 - 715;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6735.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.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_6735.phylip = 1 - 715;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_6735.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

