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

Minimum AIC  model: K2PG
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	-1036.8875	2156.2195	0.6808	0.6808
K2PI	37	-1038.7260	2159.8965	0.1083	0.7891
K2PIG	38	-1037.6150	2160.1577	0.0950	0.8841
HKYG	40	-1036.0628	2162.0649	0.0366	0.9207
K2P	36	-1041.4103	2162.7966	0.0254	0.9461
HKYI	40	-1037.1415	2164.2224	0.0125	0.9586
HKYIG	41	-1035.9008	2164.2696	0.0122	0.9707
SYM	40	-1037.5314	2165.0022	0.0084	0.9792
HKY	39	-1038.9420	2165.3101	0.0072	0.9864
GTR	43	-1034.2617	2166.0953	0.0049	0.9913
GTRG	44	-1033.4700	2167.0873	0.0030	0.9943
SYMG	41	-1037.5691	2167.6063	0.0023	0.9965
GTRI	44	-1034.5421	2169.2314	0.0010	0.9976
GTRIG	45	-1033.2822	2169.3029	0.0010	0.9985
SYMIG	42	-1037.4449	2169.9020	0.0007	0.9993
SYMI	41	-1038.7170	2169.9021	0.0007	1.0000
JC69G	36	-1055.9825	2191.9411	0.0000	1.0000
F81G	39	-1053.2903	2194.0066	0.0000	1.0000
JC69I	36	-1057.0621	2194.1003	0.0000	1.0000
JC69IG	37	-1055.8713	2194.1871	0.0000	1.0000
F81I	39	-1054.3424	2196.1108	0.0000	1.0000
F81IG	40	-1053.1669	2196.2731	0.0000	1.0000
JC69	35	-1062.8500	2203.2224	0.0000	1.0000
F81	38	-1060.3694	2205.6665	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	37	-1036.8875	2147.7751	0.5844	0.5844
K2PIG	38	-1037.6150	2151.2300	0.1039	0.6883
K2PI	37	-1038.7260	2151.4520	0.0930	0.7812
HKYG	40	-1036.0628	2152.1255	0.0664	0.8476
HKYIG	41	-1035.9008	2153.8015	0.0287	0.8763
HKYI	40	-1037.1415	2154.2830	0.0226	0.8989
GTR	43	-1034.2617	2154.5234	0.0200	0.9189
K2P	36	-1041.4103	2154.8205	0.0173	0.9362
GTRG	44	-1033.4700	2154.9400	0.0163	0.9524
SYM	40	-1037.5314	2155.0628	0.0153	0.9677
HKY	39	-1038.9420	2155.8841	0.0101	0.9778
GTRIG	45	-1033.2822	2156.5644	0.0072	0.9850
GTRI	44	-1034.5421	2157.0842	0.0056	0.9906
SYMG	41	-1037.5691	2157.1382	0.0054	0.9960
SYMIG	42	-1037.4449	2158.8898	0.0023	0.9983
SYMI	41	-1038.7170	2159.4340	0.0017	1.0000
JC69G	36	-1055.9825	2183.9650	0.0000	1.0000
F81G	39	-1053.2903	2184.5806	0.0000	1.0000
JC69IG	37	-1055.8713	2185.7427	0.0000	1.0000
JC69I	36	-1057.0621	2186.1243	0.0000	1.0000
F81IG	40	-1053.1669	2186.3337	0.0000	1.0000
F81I	39	-1054.3424	2186.6849	0.0000	1.0000
JC69	35	-1062.8500	2195.7000	0.0000	1.0000
F81	38	-1060.3694	2196.7387	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	37	-1036.8875	2292.6746	0.7173	0.7173
K2P	36	-1041.4103	2295.8038	0.1500	0.8673
K2PI	37	-1038.7260	2296.3515	0.1141	0.9814
K2PIG	38	-1037.6150	2300.0456	0.0180	0.9994
HKY	39	-1038.9420	2308.6160	0.0002	0.9996
HKYG	40	-1036.0628	2308.7736	0.0002	0.9999
HKYI	40	-1037.1415	2310.9311	0.0001	0.9999
SYM	40	-1037.5314	2311.7109	0.0001	1.0000
HKYIG	41	-1035.9008	2314.3658	0.0000	1.0000
SYMG	41	-1037.5691	2317.7025	0.0000	1.0000
SYMI	41	-1038.7170	2319.9983	0.0000	1.0000
GTR	43	-1034.2617	2322.9201	0.0000	1.0000
SYMIG	42	-1037.4449	2323.3703	0.0000	1.0000
JC69G	36	-1055.9825	2324.9483	0.0000	1.0000
JC69I	36	-1057.0621	2327.1076	0.0000	1.0000
GTRG	44	-1033.4700	2327.2529	0.0000	1.0000
GTRI	44	-1034.5421	2329.3971	0.0000	1.0000
JC69IG	37	-1055.8713	2330.6422	0.0000	1.0000
JC69	35	-1062.8500	2332.7671	0.0000	1.0000
GTRIG	45	-1033.2822	2332.7935	0.0000	1.0000
F81G	39	-1053.2903	2337.3125	0.0000	1.0000
F81I	39	-1054.3424	2339.4167	0.0000	1.0000
F81IG	40	-1053.1669	2342.9818	0.0000	1.0000
F81	38	-1060.3694	2345.5544	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset /Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_6119.phylip = 1 - 371;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_6119.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 AICc model (K2PG)]

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

