
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:40 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_22077.phylip" (19 taxa, 970 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	-4069.2582	8222.0470	0.6872	0.6872
HKYIG	41	-4068.9646	8223.6405	0.3098	0.9969
GTRG	44	-4070.6502	8233.5816	0.0021	0.9991
GTRIG	45	-4070.4183	8235.3171	0.0009	1.0000
HKYI	40	-4087.7048	8258.9403	0.0000	1.0000
GTRI	44	-4087.9104	8268.1018	0.0000	1.0000
SYMG	41	-4123.4398	8332.5908	0.0000	1.0000
K2PG	37	-4127.8919	8332.8009	0.0000	1.0000
K2PIG	38	-4127.4867	8334.1570	0.0000	1.0000
SYMIG	42	-4123.1857	8334.2677	0.0000	1.0000
SYMI	41	-4135.4431	8356.5975	0.0000	1.0000
K2PI	37	-4141.7742	8360.5655	0.0000	1.0000
HKY	39	-4179.9145	8441.1837	0.0000	1.0000
GTR	43	-4177.9256	8445.9376	0.0000	1.0000
F81G	39	-4200.1691	8481.6930	0.0000	1.0000
F81IG	40	-4200.1012	8483.7331	0.0000	1.0000
F81I	39	-4211.0442	8503.4432	0.0000	1.0000
K2P	36	-4233.5359	8541.9270	0.0000	1.0000
SYM	40	-4230.0178	8543.5664	0.0000	1.0000
JC69G	36	-4240.3373	8555.5298	0.0000	1.0000
JC69IG	37	-4240.2730	8557.5632	0.0000	1.0000
JC69I	36	-4250.0627	8574.9807	0.0000	1.0000
F81	38	-4296.6705	8672.5247	0.0000	1.0000
JC69	35	-4337.3808	8747.4597	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-4069.2582	8218.5164	0.6666	0.6666
HKYIG	41	-4068.9646	8219.9293	0.3289	0.9956
GTRG	44	-4070.6502	8229.3005	0.0030	0.9986
GTRIG	45	-4070.4183	8230.8366	0.0014	1.0000
HKYI	40	-4087.7048	8255.4096	0.0000	1.0000
GTRI	44	-4087.9104	8263.8207	0.0000	1.0000
SYMG	41	-4123.4398	8328.8796	0.0000	1.0000
K2PG	37	-4127.8919	8329.7838	0.0000	1.0000
SYMIG	42	-4123.1857	8330.3713	0.0000	1.0000
K2PIG	38	-4127.4867	8330.9734	0.0000	1.0000
SYMI	41	-4135.4431	8352.8863	0.0000	1.0000
K2PI	37	-4141.7742	8357.5483	0.0000	1.0000
HKY	39	-4179.9145	8437.8289	0.0000	1.0000
GTR	43	-4177.9256	8441.8512	0.0000	1.0000
F81G	39	-4200.1691	8478.3382	0.0000	1.0000
F81IG	40	-4200.1012	8480.2024	0.0000	1.0000
F81I	39	-4211.0442	8500.0884	0.0000	1.0000
K2P	36	-4233.5359	8539.0717	0.0000	1.0000
SYM	40	-4230.0178	8540.0357	0.0000	1.0000
JC69G	36	-4240.3373	8552.6745	0.0000	1.0000
JC69IG	37	-4240.2730	8554.5460	0.0000	1.0000
JC69I	36	-4250.0627	8572.1254	0.0000	1.0000
F81	38	-4296.6705	8669.3410	0.0000	1.0000
JC69	35	-4337.3808	8744.7616	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-4069.2582	8413.6082	0.9587	0.9587
HKYIG	41	-4068.9646	8419.8984	0.0413	1.0000
GTRG	44	-4070.6502	8443.9015	0.0000	1.0000
GTRIG	45	-4070.4183	8450.3149	0.0000	1.0000
HKYI	40	-4087.7048	8450.5014	0.0000	1.0000
GTRI	44	-4087.9104	8478.4218	0.0000	1.0000
K2PG	37	-4127.8919	8510.2437	0.0000	1.0000
K2PIG	38	-4127.4867	8516.3106	0.0000	1.0000
SYMG	41	-4123.4398	8528.8488	0.0000	1.0000
SYMIG	42	-4123.1857	8535.2177	0.0000	1.0000
K2PI	37	-4141.7742	8538.0083	0.0000	1.0000
SYMI	41	-4135.4431	8552.8554	0.0000	1.0000
HKY	39	-4179.9145	8628.0434	0.0000	1.0000
GTR	43	-4177.9256	8651.5750	0.0000	1.0000
F81G	39	-4200.1691	8668.5527	0.0000	1.0000
F81IG	40	-4200.1012	8675.2943	0.0000	1.0000
F81I	39	-4211.0442	8690.3029	0.0000	1.0000
K2P	36	-4233.5359	8714.6544	0.0000	1.0000
JC69G	36	-4240.3373	8728.2572	0.0000	1.0000
JC69IG	37	-4240.2730	8735.0060	0.0000	1.0000
SYM	40	-4230.0178	8735.1275	0.0000	1.0000
JC69I	36	-4250.0627	8747.7081	0.0000	1.0000
F81	38	-4296.6705	8854.6782	0.0000	1.0000
JC69	35	-4337.3808	8915.4669	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_22077.phylip = 1 - 970;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_22077.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_22077.phylip = 1 - 970;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_22077.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_22077.phylip = 1 - 970;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr2_22077.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

