This post is one of a series (see dental project here or in the sidebar).
After getting a set of sequences and removing chimeras, the next step is almost anticlimactic. We just copy a modified version of the shell script (from here, without the
cdcalls) or paste in the commands working from the
dentaldirectory (either individually, or all at once):
It's all over in a few seconds.
The heatmap QIIME produced is at the top of the post. It is truly a remarkable html page, with a graphic where you can reorder the columns or rows by drag and drop, and redo the map at different threshholds for the OTUs, etc. I've never seen anything quite like it. But (and this is just me), it's not pretty enough.
So what I'd like to do from here is to show you how I currently make heatmaps with matplotlib, and we'll get into that next time.
First, I have to extract the data from QIIME. The script is complicated a bit by an additional job: I'm going to organize the rows and columns. (QIIME can do this too---see the tutorial).
The columns will be in the order they appear in
sample_names.txtand the rows as they appear in
genera_and_colors.txt. These files are in the same directory. The second one starts like this:
I just do this:
Here's the first part of
The leading comma on line 1 is so the column headers line up properly in a spreadsheet. Speaking of spreadsheets, here is a screenshot after dropping the
data.csvfile onto Numbers (you could use Excel, of course):
That was painless! Zipped project files in Dropbox (here).