
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:46 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_10298.phylip" (19 taxa, 438 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	-1112.6502	2313.5623	0.4701	0.4701
HKYG	40	-1112.9436	2314.1492	0.3506	0.8207
HKYIG	41	-1112.9855	2316.6679	0.0995	0.9202
GTRI	44	-1110.3068	2318.6900	0.0362	0.9564
GTRG	44	-1110.4662	2319.0088	0.0309	0.9873
GTRIG	45	-1110.4876	2321.5363	0.0087	0.9960
HKY	39	-1118.7081	2323.2554	0.0037	0.9997
GTR	43	-1116.3259	2328.2558	0.0003	1.0000
F81I	39	-1132.3387	2350.5166	0.0000	1.0000
F81G	39	-1132.6751	2351.1895	0.0000	1.0000
F81IG	40	-1132.7118	2353.6855	0.0000	1.0000
F81	38	-1137.9421	2359.3128	0.0000	1.0000
K2PI	37	-1146.4644	2373.9589	0.0000	1.0000
K2PG	37	-1146.8120	2374.6539	0.0000	1.0000
K2PIG	38	-1146.8703	2377.1691	0.0000	1.0000
SYMI	41	-1145.2132	2381.1233	0.0000	1.0000
SYMG	41	-1145.5557	2381.8083	0.0000	1.0000
K2P	36	-1152.4429	2383.5292	0.0000	1.0000
SYMIG	42	-1145.6047	2384.3537	0.0000	1.0000
SYM	40	-1151.0474	2390.3568	0.0000	1.0000
JC69I	36	-1162.0577	2402.7588	0.0000	1.0000
JC69G	36	-1162.4029	2403.4491	0.0000	1.0000
JC69IG	37	-1162.4485	2405.9270	0.0000	1.0000
JC69	35	-1167.7629	2411.7944	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-1112.6502	2305.3004	0.4118	0.4118
HKYG	40	-1112.9436	2305.8873	0.3071	0.7189
HKYIG	41	-1112.9855	2307.9709	0.1083	0.8272
GTRI	44	-1110.3068	2308.6136	0.0786	0.9058
GTRG	44	-1110.4662	2308.9324	0.0670	0.9727
GTRIG	45	-1110.4876	2310.9751	0.0241	0.9969
HKY	39	-1118.7081	2315.4162	0.0026	0.9995
GTR	43	-1116.3259	2318.6518	0.0005	1.0000
F81I	39	-1132.3387	2342.6774	0.0000	1.0000
F81G	39	-1132.6751	2343.3503	0.0000	1.0000
F81IG	40	-1132.7118	2345.4235	0.0000	1.0000
F81	38	-1137.9421	2351.8842	0.0000	1.0000
K2PI	37	-1146.4644	2366.9289	0.0000	1.0000
K2PG	37	-1146.8120	2367.6239	0.0000	1.0000
K2PIG	38	-1146.8703	2369.7405	0.0000	1.0000
SYMI	41	-1145.2132	2372.4264	0.0000	1.0000
SYMG	41	-1145.5557	2373.1113	0.0000	1.0000
SYMIG	42	-1145.6047	2375.2094	0.0000	1.0000
K2P	36	-1152.4429	2376.8858	0.0000	1.0000
SYM	40	-1151.0474	2382.0948	0.0000	1.0000
JC69I	36	-1162.0577	2396.1154	0.0000	1.0000
JC69G	36	-1162.4029	2396.8057	0.0000	1.0000
JC69IG	37	-1162.4485	2398.8970	0.0000	1.0000
JC69	35	-1167.7629	2405.5257	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	40	-1112.6502	2468.5891	0.5468	0.5468
HKYG	40	-1112.9436	2469.1760	0.4077	0.9545
HKY	39	-1118.7081	2474.6227	0.0268	0.9813
HKYIG	41	-1112.9855	2475.3419	0.0187	0.9999
GTRI	44	-1110.3068	2488.2313	0.0000	1.0000
GTRG	44	-1110.4662	2488.5501	0.0000	1.0000
GTR	43	-1116.3259	2494.1872	0.0000	1.0000
GTRIG	45	-1110.4876	2494.6750	0.0000	1.0000
F81I	39	-1132.3387	2501.8840	0.0000	1.0000
F81G	39	-1132.6751	2502.5568	0.0000	1.0000
F81	38	-1137.9421	2507.0085	0.0000	1.0000
F81IG	40	-1132.7118	2508.7123	0.0000	1.0000
K2PI	37	-1146.4644	2517.9710	0.0000	1.0000
K2PG	37	-1146.8120	2518.6660	0.0000	1.0000
K2P	36	-1152.4429	2523.8456	0.0000	1.0000
K2PIG	38	-1146.8703	2524.8649	0.0000	1.0000
SYMI	41	-1145.2132	2539.7973	0.0000	1.0000
SYMG	41	-1145.5557	2540.4823	0.0000	1.0000
JC69I	36	-1162.0577	2543.0753	0.0000	1.0000
JC69G	36	-1162.4029	2543.7656	0.0000	1.0000
SYM	40	-1151.0474	2545.3836	0.0000	1.0000
SYMIG	42	-1145.6047	2546.6626	0.0000	1.0000
JC69	35	-1167.7629	2548.4034	0.0000	1.0000
JC69IG	37	-1162.4485	2549.9391	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_10298.phylip = 1 - 438;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10298.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_10298.phylip = 1 - 438;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10298.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_10298.phylip = 1 - 438;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr7_10298.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

