
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:20 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/chr4_16026.phylip" (19 taxa, 336 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: K2PG


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	-1149.1987	2389.5160	0.4941	0.4941
HKYI	40	-1149.9028	2390.9243	0.2443	0.7384
HKYIG	41	-1149.2008	2392.1158	0.1347	0.8731
K2PG	37	-1155.3881	2394.2124	0.0472	0.9203
GTRG	44	-1146.8991	2395.4065	0.0260	0.9463
K2PI	37	-1156.1815	2395.7993	0.0214	0.9676
K2PIG	38	-1155.3881	2396.7559	0.0132	0.9809
GTRI	44	-1147.7939	2397.1961	0.0106	0.9915
GTRIG	45	-1146.8994	2398.0747	0.0068	0.9983
SYMG	41	-1154.2323	2402.1790	0.0009	0.9992
SYMI	41	-1155.1743	2404.0630	0.0003	0.9995
SYMIG	42	-1154.2353	2404.7982	0.0002	0.9998
F81G	39	-1159.3321	2407.2047	0.0001	0.9999
HKY	39	-1159.4480	2407.4365	0.0001	0.9999
F81I	39	-1160.0838	2408.7081	0.0000	0.9999
F81IG	40	-1159.3318	2409.7823	0.0000	1.0000
JC69G	36	-1164.8196	2410.5489	0.0000	1.0000
K2P	36	-1165.4108	2411.7312	0.0000	1.0000
JC69I	36	-1165.6379	2412.1855	0.0000	1.0000
JC69IG	37	-1164.8193	2413.0748	0.0000	1.0000
GTR	43	-1157.8208	2414.6004	0.0000	1.0000
SYM	40	-1164.2780	2419.6747	0.0000	1.0000
F81	38	-1169.1782	2424.3362	0.0000	1.0000
JC69	35	-1174.4760	2427.3521	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-1149.1987	2378.3974	0.4414	0.4414
HKYI	40	-1149.9028	2379.8057	0.2183	0.6597
HKYIG	41	-1149.2008	2380.4016	0.1620	0.8217
GTRG	44	-1146.8991	2381.7983	0.0806	0.9023
GTRI	44	-1147.7939	2383.5879	0.0329	0.9353
GTRIG	45	-1146.8994	2383.7988	0.0296	0.9649
K2PG	37	-1155.3881	2384.7762	0.0182	0.9831
K2PI	37	-1156.1815	2386.3631	0.0082	0.9913
K2PIG	38	-1155.3881	2386.7761	0.0067	0.9980
SYMG	41	-1154.2323	2390.4647	0.0011	0.9991
SYMI	41	-1155.1743	2392.3487	0.0004	0.9995
SYMIG	42	-1154.2353	2392.4706	0.0004	0.9999
F81G	39	-1159.3321	2396.6642	0.0000	0.9999
HKY	39	-1159.4480	2396.8959	0.0000	0.9999
F81I	39	-1160.0838	2398.1675	0.0000	1.0000
F81IG	40	-1159.3318	2398.6636	0.0000	1.0000
JC69G	36	-1164.8196	2401.6392	0.0000	1.0000
GTR	43	-1157.8208	2401.6415	0.0000	1.0000
K2P	36	-1165.4108	2402.8215	0.0000	1.0000
JC69I	36	-1165.6379	2403.2758	0.0000	1.0000
JC69IG	37	-1164.8193	2403.6386	0.0000	1.0000
SYM	40	-1164.2780	2408.5561	0.0000	1.0000
F81	38	-1169.1782	2414.3564	0.0000	1.0000
JC69	35	-1174.4760	2418.9521	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1155.3881	2526.0093	0.6125	0.6125
K2PI	37	-1156.1815	2527.5962	0.2770	0.8895
HKYG	40	-1149.1987	2531.0818	0.0485	0.9380
K2PIG	38	-1155.3881	2531.8264	0.0334	0.9714
HKYI	40	-1149.9028	2532.4901	0.0240	0.9954
HKYIG	41	-1149.2008	2536.9031	0.0026	0.9980
JC69G	36	-1164.8196	2539.0552	0.0009	0.9989
K2P	36	-1165.4108	2540.2375	0.0005	0.9994
JC69I	36	-1165.6379	2540.6918	0.0004	0.9998
JC69IG	37	-1164.8193	2544.8717	0.0000	0.9999
F81G	39	-1159.3321	2545.5315	0.0000	0.9999
HKY	39	-1159.4480	2545.7633	0.0000	0.9999
SYMG	41	-1154.2323	2546.9663	0.0000	1.0000
F81I	39	-1160.0838	2547.0349	0.0000	1.0000
SYMI	41	-1155.1743	2548.8502	0.0000	1.0000
GTRG	44	-1146.8991	2549.7512	0.0000	1.0000
F81IG	40	-1159.3318	2551.3481	0.0000	1.0000
GTRI	44	-1147.7939	2551.5408	0.0000	1.0000
JC69	35	-1174.4760	2552.5510	0.0000	1.0000
SYMIG	42	-1154.2353	2552.7892	0.0000	1.0000
GTRIG	45	-1146.8994	2555.5688	0.0000	1.0000
F81	38	-1169.1782	2559.4066	0.0000	1.0000
SYM	40	-1164.2780	2561.2405	0.0000	1.0000
GTR	43	-1157.8208	2565.7773	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/chr4_16026.phylip = 1 - 336;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr4_16026.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/chr4_16026.phylip = 1 - 336;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr4_16026.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 (K2PG)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr4_16026.phylip = 1 - 336;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr4_16026.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

