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

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


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	-528.5481	1162.5225	0.4821	0.4821
HKYI	40	-529.1835	1163.7934	0.2554	0.7375
HKY	39	-531.4559	1164.9118	0.1460	0.8835
HKYIG	41	-528.5479	1166.0021	0.0846	0.9682
K2PG	37	-537.3813	1170.0657	0.0111	0.9793
K2PI	37	-538.0139	1171.3309	0.0059	0.9852
K2P	36	-539.8823	1171.7946	0.0047	0.9898
F81G	39	-535.2073	1172.4147	0.0034	0.9933
K2PIG	38	-537.3809	1173.3877	0.0021	0.9954
F81I	39	-536.2316	1174.4632	0.0012	0.9966
F81	38	-538.2159	1175.0578	0.0009	0.9975
GTRG	44	-527.9137	1175.5073	0.0007	0.9983
F81IG	40	-535.2071	1175.8405	0.0006	0.9989
GTRI	44	-528.8256	1177.3312	0.0003	0.9992
JC69G	36	-542.8368	1177.7038	0.0002	0.9994
GTR	43	-530.9022	1177.8362	0.0002	0.9996
GTRIG	45	-527.9131	1179.2134	0.0001	0.9997
JC69I	36	-543.8641	1179.7583	0.0001	0.9998
JC69	35	-545.7232	1180.2524	0.0001	0.9999
JC69IG	37	-542.8365	1180.9760	0.0000	1.0000
SYMG	41	-536.5750	1182.0562	0.0000	1.0000
SYM	40	-539.4221	1184.2707	0.0000	1.0000
SYMI	41	-537.7871	1184.4804	0.0000	1.0000
SYMIG	42	-536.7787	1185.9983	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	40	-528.5481	1137.0961	0.4706	0.4706
HKYI	40	-529.1835	1138.3670	0.2493	0.7199
HKYIG	41	-528.5479	1139.0958	0.1732	0.8931
HKY	39	-531.4559	1140.9118	0.0698	0.9629
GTRG	44	-527.9137	1143.8273	0.0163	0.9792
GTRI	44	-528.8256	1145.6512	0.0065	0.9857
GTRIG	45	-527.9131	1145.8263	0.0060	0.9917
GTR	43	-530.9022	1147.8045	0.0022	0.9939
F81G	39	-535.2073	1148.4147	0.0016	0.9955
K2PG	37	-537.3813	1148.7626	0.0014	0.9969
K2PI	37	-538.0139	1150.0279	0.0007	0.9976
F81IG	40	-535.2071	1150.4142	0.0006	0.9982
F81I	39	-536.2316	1150.4632	0.0006	0.9988
K2PIG	38	-537.3809	1150.7617	0.0005	0.9993
K2P	36	-539.8823	1151.7646	0.0003	0.9997
F81	38	-538.2159	1152.4318	0.0002	0.9999
SYMG	41	-536.5750	1155.1500	0.0001	0.9999
SYMIG	42	-536.7787	1157.5574	0.0000	0.9999
SYMI	41	-537.7871	1157.5741	0.0000	1.0000
JC69G	36	-542.8368	1157.6737	0.0000	1.0000
SYM	40	-539.4221	1158.8443	0.0000	1.0000
JC69IG	37	-542.8365	1159.6730	0.0000	1.0000
JC69I	36	-543.8641	1159.7282	0.0000	1.0000
JC69	35	-545.7232	1161.4464	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	40	-528.5481	1262.5281	0.3070	0.3070
HKY	39	-531.4559	1263.2079	0.2185	0.5255
HKYI	40	-529.1835	1263.7990	0.1626	0.6882
K2P	36	-539.8823	1264.6533	0.1061	0.7942
K2PG	37	-537.3813	1264.7872	0.0992	0.8935
K2PI	37	-538.0139	1266.0524	0.0527	0.9462
HKYIG	41	-528.5479	1267.6635	0.0235	0.9697
K2PIG	38	-537.3809	1269.9220	0.0076	0.9773
JC69G	36	-542.8368	1270.5624	0.0055	0.9828
F81G	39	-535.2073	1270.7108	0.0051	0.9880
JC69	35	-545.7232	1271.1993	0.0040	0.9920
F81	38	-538.2159	1271.5921	0.0033	0.9953
JC69I	36	-543.8641	1272.6169	0.0020	0.9973
F81I	39	-536.2316	1272.7593	0.0018	0.9991
JC69IG	37	-542.8365	1275.6975	0.0004	0.9995
F81IG	40	-535.2071	1275.8461	0.0004	0.9999
GTRG	44	-527.9137	1281.8025	0.0000	1.0000
GTR	43	-530.9022	1282.6438	0.0000	1.0000
GTRI	44	-528.8256	1283.6263	0.0000	1.0000
SYMG	41	-536.5750	1283.7177	0.0000	1.0000
SYM	40	-539.4221	1284.2762	0.0000	1.0000
SYMI	41	-537.7871	1286.1419	0.0000	1.0000
GTRIG	45	-527.9131	1286.9372	0.0000	1.0000
SYMIG	42	-536.7787	1289.2609	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/chr3_18256.phylip = 1 - 170;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_18256.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/chr3_18256.phylip = 1 - 170;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr3_18256.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 (HKYG)]

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


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

