Book Image

R Bioinformatics Cookbook - Second Edition

By : Dan MacLean
Book Image

R Bioinformatics Cookbook - Second Edition

By: Dan MacLean

Overview of this book

The updated second edition of R Bioinformatics Cookbook takes a recipe-based approach to show you how to conduct practical research and analysis in computational biology with R. You’ll learn how to create a useful and modular R working environment, along with loading, cleaning, and analyzing data using the most up-to-date Bioconductor, ggplot2, and tidyverse tools. This book will walk you through the Bioconductor tools necessary for you to understand and carry out protocols in RNA-seq and ChIP-seq, phylogenetics, genomics, gene search, gene annotation, statistical analysis, and sequence analysis. As you advance, you'll find out how to use Quarto to create data-rich reports, presentations, and websites, as well as get a clear understanding of how machine learning techniques can be applied in the bioinformatics domain. The concluding chapters will help you develop proficiency in key skills, such as gene annotation analysis and functional programming in purrr and base R. Finally, you'll discover how to use the latest AI tools, including ChatGPT, to generate, edit, and understand R code and draft workflows for complex analyses. By the end of this book, you'll have gained a solid understanding of the skills and techniques needed to become a bioinformatics specialist and efficiently work with large and complex bioinformatics datasets.
Table of Contents (16 chapters)

Visualizing trees of many genes quickly with ggtree

Once you have computed a tree, the first thing you will want to do with it is take a look. That’s possible in many programs, but R has an extremely powerful, flexible, and fast system in the form of the ggtree package. In this recipe, we’ll learn how to get into ggtree and re-layout, highlight, and annotate tree images in just a few commands.

Getting ready

You’ll need the ggplot2, ggtree, and ape packages. You’ll also require the itol.nwk file from the rbioinfcookbook package. The file is a Newick tree of 191 species from the Interactive Tree of Life online tool’s public dataset. At the time of writing, there is an issue with an upstream dependency that causes this code to fail, though it is correct. We hope this will have gone away by the time you read this. If it hasn’t, a workaround is to install the source version of ggtree from Biocmanager, like this:

BiocManager::install(&quot...