Loading phylogenetic data into R

Martin R. Smith

2019-06-03

Loading phylogenetic data into R

It can be a bit fiddly to get a phylogenetic dataset into R, particularly if you are not used to working with files in the NEXUS format.

Filesystem navigation

To open a file, you’re going to have to tell R where it is. You can do this by providing the full path to the file on your system. Be careful to use forward slashes (/, not \, which you’ll get if you copy file paths in Windows).

You can save typing by giving R a working directory. You can think of R as having a file explorer window open invisibly in the background. You can see the folder that’s open at the moment by typing getwd() at the console. setwd() tells R to open a different folder instead. `setwd(‘../’) tells R to go up to a parent directory. (You can do this using the Graphical User Interface in RStudio).

By setting the directory that your files are in as the working directory, you only need to specify the filename:

Getting raw data:

From an Excel spreadsheet

If your data is in an Excel spreadsheet, one way to load it into R is using the xlsx package. First you’ll have to install it:

Then you should prepare your Excel spreadsheet such that each row corresponds to a taxon, and each column to a character.

Then you can read the data from the Excel file by telling R which sheet, rows and columns contain your data:

From a Nexus file

TreeSearch contains an inbuilt Nexus parser:

This will extract character names and codings from a dataset. It’s been written to work with datasets downloaded from MorphoBank, but my aim is for this function to handle most valid (and many invalid) NEXUS files. If you find a file that this function can’t handle, please let me know and I’ll try to fix it.

In the meantime, alternative Nexus parsers are available: try

Non-standard elements of a Nexus file might be beyond the capabilities of ape’s parser. In particular, you will need to replace spaces in taxon names with an underscore, and to arrange all data into a single block starting BEGIN DATA. You’ll need to strip out comments, character definitions and separate taxon blocks.

The function readNexus in package phylobase uses the NCL library and promises to be more powerful, but I’ve not been able to get it to work.

From a TNT file

A TNT format dataset downloaded from MorphoBank can be parsed with ReadTntCharacters, which might also handle other TNT-compatible files. If there’s a file that’s not being read correctly, please let me know and I’ll try to fix it.

Processing raw data

The next stage is to get the raw data into a format that TreeSearch can understand. If you’ve used the ReadAsPhyDat or ReadTntAsPhyDat functions, then you can skip this step – you’re already there.

Otherwise, you can try

or if that doesn’t work,

These functions haven’t been exhaustively tested – if they don’t work on your dataset, please let me know.

Failing that, you can enlist the help of the phangorn package, which was installed when you installed TreeSearch:

type='USER' tells the parser to expect morphological data.

The levels parameter simply lists all the states that any character might take. 0:9 includes all the integer digits from 0 to 9. If you have inapplicable data in your matrix, you should list - as a separate level as it represents an additional state (as handled by the Morphy implementation of (Brazeau, Guillerme, & Smith, 2019)).

If you have more complicated ambiguities, you may need to use a contrast matrix to decode your matrix.

A contrast matrix translates the tokens used in your dataset to the character states to which they correspond: for example decoding ‘A’ to {01}. For more details, see the ‘phangorn-specials’ vignette in the phangorn package, accessible by typing ‘?phangorn’ in the R prompt and navigating to index > package vignettes.

##   0 1 -
## 0 1 0 0
## 1 0 1 0
## - 0 0 1
## A 1 1 0
## + 1 1 0
## ? 1 1 1

If you need to use a contrast matrix, convert the data using

What next?

You might want to:

References

Brazeau, M. D., Guillerme, T., & Smith, M. R. (2019). An algorithm for morphological phylogenetic analysis with inapplicable data. Systematic Biology, in press. doi:10.1093/sysbio/syy083