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

Minimum AIC  model: K2PIG
Minimum AICc model: K2PG
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
K2PG	37	-1074.1862	2230.0977	0.3054	0.3054
K2PIG	38	-1073.1572	2230.4797	0.2523	0.5577
K2PI	37	-1074.3893	2230.5039	0.2493	0.8069
HKYG	40	-1072.0512	2233.1883	0.0651	0.8721
HKYIG	41	-1071.0135	2233.5938	0.0532	0.9252
HKYI	40	-1072.2788	2233.6435	0.0519	0.9771
SYMG	41	-1072.9524	2237.4715	0.0077	0.9848
SYMI	41	-1073.1777	2237.9220	0.0061	0.9909
SYMIG	42	-1072.0542	2238.1698	0.0054	0.9963
GTRG	44	-1070.8192	2240.7308	0.0015	0.9978
GTRI	44	-1071.0696	2241.2317	0.0012	0.9989
GTRIG	45	-1069.8848	2241.3987	0.0011	1.0000
JC69G	36	-1089.5559	2258.4105	0.0000	1.0000
JC69I	36	-1089.7923	2258.8833	0.0000	1.0000
JC69IG	37	-1088.5835	2258.8923	0.0000	1.0000
F81G	39	-1087.5047	2261.6283	0.0000	1.0000
F81IG	40	-1086.5070	2262.0999	0.0000	1.0000
F81I	39	-1087.7570	2262.1327	0.0000	1.0000
K2P	36	-1097.2623	2273.8232	0.0000	1.0000
HKY	39	-1094.9188	2276.4564	0.0000	1.0000
SYM	40	-1095.8174	2280.7206	0.0000	1.0000
GTR	43	-1093.5221	2283.6141	0.0000	1.0000
JC69	35	-1112.0940	2301.0733	0.0000	1.0000
F81	38	-1109.9501	2304.0655	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PIG	38	-1073.1572	2222.3144	0.2405	0.2405
K2PG	37	-1074.1862	2222.3725	0.2336	0.4741
K2PI	37	-1074.3893	2222.7787	0.1907	0.6647
HKYIG	41	-1071.0135	2224.0271	0.1021	0.7668
HKYG	40	-1072.0512	2224.1024	0.0984	0.8652
HKYI	40	-1072.2788	2224.5576	0.0783	0.9435
SYMG	41	-1072.9524	2227.9048	0.0147	0.9582
SYMIG	42	-1072.0542	2228.1085	0.0133	0.9715
SYMI	41	-1073.1777	2228.3553	0.0117	0.9832
GTRG	44	-1070.8192	2229.6383	0.0062	0.9894
GTRIG	45	-1069.8848	2229.7695	0.0058	0.9952
GTRI	44	-1071.0696	2230.1392	0.0048	1.0000
JC69G	36	-1089.5559	2251.1118	0.0000	1.0000
JC69IG	37	-1088.5835	2251.1670	0.0000	1.0000
JC69I	36	-1089.7923	2251.5846	0.0000	1.0000
F81G	39	-1087.5047	2253.0095	0.0000	1.0000
F81IG	40	-1086.5070	2253.0140	0.0000	1.0000
F81I	39	-1087.7570	2253.5139	0.0000	1.0000
K2P	36	-1097.2623	2266.5245	0.0000	1.0000
HKY	39	-1094.9188	2267.8376	0.0000	1.0000
SYM	40	-1095.8174	2271.6347	0.0000	1.0000
GTR	43	-1093.5221	2273.0443	0.0000	1.0000
JC69	35	-1112.0940	2294.1881	0.0000	1.0000
F81	38	-1109.9501	2295.9002	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1074.1862	2370.2412	0.5108	0.5108
K2PI	37	-1074.3893	2370.6474	0.4169	0.9276
K2PIG	38	-1073.1572	2374.1796	0.0713	0.9989
HKYG	40	-1072.0512	2383.9605	0.0005	0.9995
HKYI	40	-1072.2788	2384.4157	0.0004	0.9999
HKYIG	41	-1071.0135	2387.8816	0.0001	1.0000
SYMG	41	-1072.9524	2391.7593	0.0000	1.0000
SYMI	41	-1073.1777	2392.2099	0.0000	1.0000
JC69G	36	-1089.5559	2394.9841	0.0000	1.0000
JC69I	36	-1089.7923	2395.4569	0.0000	1.0000
SYMIG	42	-1072.0542	2395.9595	0.0000	1.0000
JC69IG	37	-1088.5835	2399.0358	0.0000	1.0000
GTRG	44	-1070.8192	2405.4822	0.0000	1.0000
GTRI	44	-1071.0696	2405.9831	0.0000	1.0000
F81G	39	-1087.5047	2408.8711	0.0000	1.0000
F81I	39	-1087.7570	2409.3756	0.0000	1.0000
GTRIG	45	-1069.8848	2409.6099	0.0000	1.0000
K2P	36	-1097.2623	2410.3968	0.0000	1.0000
F81IG	40	-1086.5070	2412.8721	0.0000	1.0000
HKY	39	-1094.9188	2423.6992	0.0000	1.0000
SYM	40	-1095.8174	2431.4928	0.0000	1.0000
JC69	35	-1112.0940	2434.0639	0.0000	1.0000
GTR	43	-1093.5221	2444.8917	0.0000	1.0000
F81	38	-1109.9501	2447.7654	0.0000	1.0000



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

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

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


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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_4772.phylip = 1 - 402;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_4772.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;


[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/chrZ_4772.phylip = 1 - 402;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chrZ_4772.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

