Search icon CANCEL
Subscription
0
Cart icon
Your Cart (0 item)
Close icon
You have no products in your basket yet
Arrow left icon
Explore Products
Best Sellers
New Releases
Books
Videos
Audiobooks
Learning Hub
Free Learning
Arrow right icon
Arrow up icon
GO TO TOP
Bioinformatics with Python Cookbook

You're reading from   Bioinformatics with Python Cookbook Learn how to use modern Python bioinformatics libraries and applications to do cutting-edge research in computational biology

Arrow left icon
Product type Paperback
Published in Nov 2018
Publisher Packt
ISBN-13 9781789344691
Length 360 pages
Edition 2nd Edition
Languages
Arrow right icon
Author (1):
Arrow left icon
Tiago Antao Tiago Antao
Author Profile Icon Tiago Antao
Tiago Antao
Arrow right icon
View More author details
Toc

Table of Contents (12) Chapters Close

Preface 1. Python and the Surrounding Software Ecology FREE CHAPTER 2. Next-Generation Sequencing 3. Working with Genomes 4. Population Genetics 5. Population Genetics Simulation 6. Phylogenetics 7. Using the Protein Data Bank 8. Bioinformatics Pipelines 9. Python for Big Genomics Datasets 10. Other Topics in Bioinformatics 11. Advanced NGS Processing

Using decision trees to explore the data


We are now ready to start exploring the data with the objective of finding some rules on how to filter it. Because we have a lot of annotations to explore (in our case, we reduced them, but generally, that would be the case), we need to find a place to start. It can be daunting to go out on a blind fishing expedition. My personal preference for a first approach is using a machine learning technique called decision trees. Decision trees will suggest what the fundamental annotations segregating the data in correct and error calls are. Another advantage of decision trees is that they barely need any data preparation, as opposed to many other machine learning techniques.

How to do it…

  1. We start with a few imports, most notably of scikit-learn:
import gzip
import pickle

import numpy as np
import graphviz
from sklearn import tree
  1. Let's load the data and split it into inputs and outputs:
balanced_fit = np.load(gzip.open('balanced_fit.npy.gz', 'rb'))
ordered_features...
lock icon The rest of the chapter is locked
Register for a free Packt account to unlock a world of extra content!
A free Packt account unlocks extra newsletters, articles, discounted offers, and much more. Start advancing your knowledge today.
Unlock this book and the full library FREE for 7 days
Get unlimited access to 7000+ expert-authored eBooks and videos courses covering every tech area you can think of
Renews at $19.99/month. Cancel anytime
Banner background image