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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: K2PI


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	-568.2387	1234.9044	0.4308	0.4308
HKYG	40	-568.3262	1235.0794	0.3947	0.8256
HKYIG	41	-568.3259	1238.1094	0.0868	0.9123
HKY	39	-572.1250	1239.6801	0.0396	0.9519
K2PI	37	-575.9174	1241.3708	0.0170	0.9689
K2PG	37	-575.9903	1241.5165	0.0158	0.9847
K2PIG	38	-575.9900	1244.4467	0.0036	0.9883
GTRI	44	-567.2739	1245.3063	0.0024	0.9907
F81I	39	-575.0133	1245.4568	0.0022	0.9929
GTRG	44	-567.3613	1245.4812	0.0022	0.9951
F81G	39	-575.1006	1245.6313	0.0020	0.9971
K2P	36	-580.0133	1246.6640	0.0012	0.9983
F81IG	40	-575.1003	1248.6275	0.0005	0.9987
GTRIG	45	-567.3610	1248.6527	0.0004	0.9992
GTR	43	-571.0543	1249.7315	0.0003	0.9994
F81	38	-579.0432	1250.5531	0.0002	0.9996
JC69I	36	-582.5494	1251.7362	0.0001	0.9997
JC69G	36	-582.6351	1251.9075	0.0001	0.9998
SYMI	41	-575.3248	1252.1072	0.0001	0.9999
SYMG	41	-575.3867	1252.2309	0.0001	1.0000
JC69IG	37	-582.6348	1254.8056	0.0000	1.0000
SYMIG	42	-575.3978	1255.3184	0.0000	1.0000
SYM	40	-579.2597	1256.9464	0.0000	1.0000
JC69	35	-586.6729	1257.1163	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	40	-568.2387	1216.4775	0.4082	0.4082
HKYG	40	-568.3262	1216.6524	0.3740	0.7821
HKYIG	41	-568.3259	1218.6518	0.1376	0.9198
HKY	39	-572.1250	1222.2499	0.0228	0.9425
GTRI	44	-567.2739	1222.5477	0.0196	0.9622
GTRG	44	-567.3613	1222.7225	0.0180	0.9801
GTRIG	45	-567.3610	1224.7220	0.0066	0.9868
K2PI	37	-575.9174	1225.8349	0.0038	0.9905
K2PG	37	-575.9903	1225.9805	0.0035	0.9941
K2PIG	38	-575.9900	1227.9800	0.0013	0.9954
F81I	39	-575.0133	1228.0267	0.0013	0.9966
GTR	43	-571.0543	1228.1086	0.0012	0.9979
F81G	39	-575.1006	1228.2011	0.0012	0.9990
F81IG	40	-575.1003	1230.2006	0.0004	0.9994
K2P	36	-580.0133	1232.0266	0.0002	0.9996
SYMI	41	-575.3248	1232.6496	0.0001	0.9997
SYMG	41	-575.3867	1232.7733	0.0001	0.9999
F81	38	-579.0432	1234.0865	0.0001	0.9999
SYMIG	42	-575.3978	1234.7956	0.0000	1.0000
JC69I	36	-582.5494	1237.0988	0.0000	1.0000
JC69G	36	-582.6351	1237.2702	0.0000	1.0000
SYM	40	-579.2597	1238.5194	0.0000	1.0000
JC69IG	37	-582.6348	1239.2696	0.0000	1.0000
JC69	35	-586.6729	1243.3458	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	37	-575.9174	1351.2305	0.2611	0.2611
K2PG	37	-575.9903	1351.3762	0.2428	0.5039
HKYI	40	-568.2387	1352.0403	0.1742	0.6781
HKYG	40	-568.3262	1352.2153	0.1596	0.8377
K2P	36	-580.0133	1354.0332	0.0643	0.9020
HKY	39	-572.1250	1354.4237	0.0529	0.9549
K2PIG	38	-575.9900	1356.7647	0.0164	0.9713
HKYIG	41	-568.3259	1357.6038	0.0108	0.9821
JC69I	36	-582.5494	1359.1054	0.0051	0.9871
JC69G	36	-582.6351	1359.2767	0.0047	0.9918
F81I	39	-575.0133	1360.2005	0.0029	0.9948
F81G	39	-575.1006	1360.3749	0.0027	0.9975
JC69	35	-586.6729	1361.9633	0.0012	0.9987
F81	38	-579.0432	1362.8712	0.0008	0.9995
JC69IG	37	-582.6348	1364.6653	0.0003	0.9998
F81IG	40	-575.1003	1365.7634	0.0002	1.0000
SYMI	41	-575.3248	1371.6015	0.0000	1.0000
GTRI	44	-567.2739	1371.6669	0.0000	1.0000
SYMG	41	-575.3867	1371.7252	0.0000	1.0000
GTRG	44	-567.3613	1371.8417	0.0000	1.0000
GTR	43	-571.0543	1373.8387	0.0000	1.0000
SYM	40	-579.2597	1374.0823	0.0000	1.0000
SYMIG	42	-575.3978	1377.1367	0.0000	1.0000
GTRIG	45	-567.3610	1377.2302	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/chr6_3116.phylip = 1 - 219;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3116.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/chr6_3116.phylip = 1 - 219;
 Partition Dummy = 1:/Users/bcf/Git/brant/seqcap/Manuscripts/Tetrapods/Nature/Data/917Loci_19Species/alignments_from_loci/917Loci_19Species_alignments_phylip/chr6_3116.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 (K2PI)]

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


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

