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

Minimum AIC  model: JC69I
Minimum AICc model: JC69I
Minimum BIC  model: JC69I


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
JC69I	36	-463.4828	1018.6990	0.4027	0.4027
JC69G	36	-463.6772	1019.0877	0.3316	0.7343
K2PI	37	-463.3748	1021.7347	0.0883	0.8225
K2PG	37	-463.5553	1022.0958	0.0737	0.8962
JC69IG	37	-463.6769	1022.3388	0.0653	0.9615
K2PIG	38	-463.5550	1025.3957	0.0142	0.9756
F81I	39	-462.7375	1027.1114	0.0060	0.9816
F81G	39	-462.9435	1027.5233	0.0049	0.9865
SYMI	41	-459.6408	1027.7740	0.0043	0.9908
SYMG	41	-459.8875	1028.2674	0.0034	0.9942
JC69	35	-470.7081	1029.9456	0.0015	0.9956
HKYI	40	-462.6323	1030.3028	0.0012	0.9969
HKYG	40	-462.8240	1030.6863	0.0010	0.9979
F81IG	40	-462.9431	1030.9244	0.0009	0.9987
SYMIG	42	-459.8872	1031.7744	0.0006	0.9993
K2P	36	-470.5651	1032.8635	0.0003	0.9997
HKYIG	41	-462.8236	1034.1396	0.0002	0.9998
GTRI	44	-458.6599	1036.5009	0.0001	0.9999
GTRG	44	-458.9332	1037.0475	0.0000	0.9999
F81	38	-469.8875	1038.0607	0.0000	1.0000
SYM	40	-466.8219	1038.6820	0.0000	1.0000
GTRIG	45	-458.9329	1040.7229	0.0000	1.0000
HKY	39	-469.7541	1041.1445	0.0000	1.0000
GTR	43	-465.5272	1046.6170	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
JC69I	36	-463.4828	998.9657	0.2465	0.2465
JC69G	36	-463.6772	999.3544	0.2030	0.4495
K2PI	37	-463.3748	1000.7497	0.1010	0.5506
K2PG	37	-463.5553	1001.1107	0.0844	0.6349
SYMI	41	-459.6408	1001.2817	0.0774	0.7124
JC69IG	37	-463.6769	1001.3538	0.0747	0.7871
SYMG	41	-459.8875	1001.7751	0.0605	0.8476
K2PIG	38	-463.5550	1003.1100	0.0310	0.8786
F81I	39	-462.7375	1003.4751	0.0259	0.9045
SYMIG	42	-459.8872	1003.7744	0.0223	0.9267
F81G	39	-462.9435	1003.8869	0.0210	0.9478
HKYI	40	-462.6323	1005.2647	0.0106	0.9584
GTRI	44	-458.6599	1005.3198	0.0103	0.9687
HKYG	40	-462.8240	1005.6481	0.0087	0.9774
GTRG	44	-458.9332	1005.8664	0.0078	0.9852
F81IG	40	-462.9431	1005.8862	0.0077	0.9929
HKYIG	41	-462.8236	1007.6473	0.0032	0.9962
GTRIG	45	-458.9329	1007.8657	0.0029	0.9990
JC69	35	-470.7081	1011.4161	0.0005	0.9995
K2P	36	-470.5651	1013.1302	0.0002	0.9997
SYM	40	-466.8219	1013.6438	0.0002	0.9999
F81	38	-469.8875	1015.7750	0.0001	0.9999
GTR	43	-465.5272	1017.0545	0.0000	1.0000
HKY	39	-469.7541	1017.5081	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
JC69I	36	-463.4828	1112.2755	0.4854	0.4854
JC69G	36	-463.6772	1112.6642	0.3997	0.8851
K2PI	37	-463.3748	1117.2070	0.0412	0.9264
K2PG	37	-463.5553	1117.5680	0.0344	0.9608
JC69IG	37	-463.6769	1117.8111	0.0305	0.9913
JC69	35	-470.7081	1121.5784	0.0046	0.9959
K2PIG	38	-463.5550	1122.7148	0.0026	0.9986
F81I	39	-462.7375	1126.2273	0.0005	0.9990
K2P	36	-470.5651	1126.4400	0.0004	0.9994
F81G	39	-462.9435	1126.6392	0.0004	0.9998
SYMI	41	-459.6408	1130.3289	0.0001	0.9998
SYMG	41	-459.8875	1130.8224	0.0000	0.9999
HKYI	40	-462.6323	1131.1645	0.0000	0.9999
HKYG	40	-462.8240	1131.5479	0.0000	1.0000
F81IG	40	-462.9431	1131.7860	0.0000	1.0000
F81	38	-469.8875	1135.3798	0.0000	1.0000
SYMIG	42	-459.8872	1135.9692	0.0000	1.0000
HKYIG	41	-462.8236	1136.6946	0.0000	1.0000
SYM	40	-466.8219	1139.5436	0.0000	1.0000
HKY	39	-469.7541	1140.2604	0.0000	1.0000
GTRI	44	-458.6599	1143.8096	0.0000	1.0000
GTRG	44	-458.9332	1144.3562	0.0000	1.0000
GTRIG	45	-458.9329	1149.5030	0.0000	1.0000
GTR	43	-465.5272	1152.3967	0.0000	1.0000



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

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

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


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

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


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

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


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

