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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	40	-1075.0554	2243.1266	0.6724	0.6724
HKYIG	41	-1075.0188	2245.7587	0.1803	0.8527
HKYG	40	-1076.7247	2246.4654	0.1266	0.9793
GTRI	44	-1073.5066	2250.9809	0.0132	0.9926
GTRG	44	-1074.7373	2253.4424	0.0039	0.9964
GTRIG	45	-1073.4258	2253.6127	0.0036	1.0000
HKY	39	-1091.8765	2274.0849	0.0000	1.0000
F81I	39	-1093.8685	2278.0691	0.0000	1.0000
F81G	39	-1094.7249	2279.7819	0.0000	1.0000
F81IG	40	-1093.8863	2280.7884	0.0000	1.0000
GTR	43	-1090.0428	2281.2824	0.0000	1.0000
SYMI	41	-1098.0996	2291.9204	0.0000	1.0000
SYMG	41	-1098.7592	2293.2395	0.0000	1.0000
K2PI	37	-1104.6833	2294.3941	0.0000	1.0000
SYMIG	42	-1098.7592	2295.9664	0.0000	1.0000
K2PG	37	-1105.5085	2296.0445	0.0000	1.0000
K2PIG	38	-1105.5086	2298.6865	0.0000	1.0000
F81	38	-1109.9558	2307.5808	0.0000	1.0000
JC69I	36	-1120.4446	2323.2954	0.0000	1.0000
JC69G	36	-1121.3504	2325.1071	0.0000	1.0000
JC69IG	37	-1120.8306	2326.6887	0.0000	1.0000
K2P	36	-1124.8834	2332.1730	0.0000	1.0000
SYM	40	-1120.6439	2334.3037	0.0000	1.0000
JC69	35	-1138.5862	2356.9778	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1075.0554	2230.1107	0.5829	0.5829
HKYIG	41	-1075.0188	2232.0376	0.2224	0.8053
HKYG	40	-1076.7247	2233.4495	0.1098	0.9151
GTRI	44	-1073.5066	2235.0132	0.0502	0.9653
GTRIG	45	-1073.4258	2236.8516	0.0200	0.9853
GTRG	44	-1074.7373	2237.4746	0.0147	1.0000
HKY	39	-1091.8765	2261.7529	0.0000	1.0000
F81I	39	-1093.8685	2265.7371	0.0000	1.0000
GTR	43	-1090.0428	2266.0856	0.0000	1.0000
F81G	39	-1094.7249	2267.4498	0.0000	1.0000
F81IG	40	-1093.8863	2267.7725	0.0000	1.0000
SYMI	41	-1098.0996	2278.1992	0.0000	1.0000
SYMG	41	-1098.7592	2279.5184	0.0000	1.0000
SYMIG	42	-1098.7592	2281.5184	0.0000	1.0000
K2PI	37	-1104.6833	2283.3666	0.0000	1.0000
K2PG	37	-1105.5085	2285.0171	0.0000	1.0000
K2PIG	38	-1105.5086	2287.0172	0.0000	1.0000
F81	38	-1109.9558	2295.9115	0.0000	1.0000
JC69I	36	-1120.4446	2312.8892	0.0000	1.0000
JC69G	36	-1121.3504	2314.7009	0.0000	1.0000
JC69IG	37	-1120.8306	2315.6612	0.0000	1.0000
SYM	40	-1120.6439	2321.2879	0.0000	1.0000
K2P	36	-1124.8834	2321.7667	0.0000	1.0000
JC69	35	-1138.5862	2347.1724	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1075.0554	2377.3176	0.8006	0.8006
HKYG	40	-1076.7247	2380.6564	0.1508	0.9514
HKYIG	41	-1075.0188	2382.9247	0.0485	0.9999
GTRI	44	-1073.5066	2396.9408	0.0000	1.0000
GTRG	44	-1074.7373	2399.4022	0.0000	1.0000
GTRIG	45	-1073.4258	2402.4593	0.0000	1.0000
HKY	39	-1091.8765	2405.2797	0.0000	1.0000
F81I	39	-1093.8685	2409.2638	0.0000	1.0000
F81G	39	-1094.7249	2410.9766	0.0000	1.0000
F81IG	40	-1093.8863	2414.9794	0.0000	1.0000
K2PI	37	-1104.6833	2419.5330	0.0000	1.0000
K2PG	37	-1105.5085	2421.1835	0.0000	1.0000
GTR	43	-1090.0428	2424.3330	0.0000	1.0000
K2PIG	38	-1105.5086	2426.8638	0.0000	1.0000
SYMI	41	-1098.0996	2429.0863	0.0000	1.0000
SYMG	41	-1098.7592	2430.4054	0.0000	1.0000
F81	38	-1109.9558	2435.7581	0.0000	1.0000
SYMIG	42	-1098.7592	2436.0856	0.0000	1.0000
JC69I	36	-1120.4446	2445.3754	0.0000	1.0000
JC69G	36	-1121.3504	2447.1871	0.0000	1.0000
JC69IG	37	-1120.8306	2451.8276	0.0000	1.0000
K2P	36	-1124.8834	2454.2530	0.0000	1.0000
SYM	40	-1120.6439	2468.4948	0.0000	1.0000
JC69	35	-1138.5862	2475.9784	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10882.phylip = 1 - 293;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10882.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10882.phylip = 1 - 293;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10882.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10882.phylip = 1 - 293;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10882.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

