Siliceous and Delicious

Had a frustrating morning not knowing quite where to pick up—continue struggling with the morphological/molecular distance plot or move on to the rest of the list. There was some amount of distraction, too, with the lab’s new intern, Sarah, an undergrad from France, arriving and needing some assistance getting settled it. Fortunately, I got an email back from Maude after lunch, introducing me to her labmate Allison, an expert in all things R and phylogeny. She sat down with my problem and was extremely helpful—she very quickly pointed out some pretty major things (retrospectively pretty obvious):

1. I was trying to read an RNA sequence with a function called read.dna(). 
2. DNA and RNA are not the same.
3. The distance measures implemented in the dist.dna() function probably do not apply to RNA.

She pointed me toward the package {seqinr}, which includes a function read.alignment() that reads clustal files like mine at face value, including the U nucleotides. Hallelujah. She also pointed out another fact that was almost immediately obvious retrospectively: pairwise genetic distances computed directly from sequences are not the same as the genetic distances measured along branches of a tree (those are known as patristic distances). It might make more sense to use patristic distances, since they actually reflect the ‘evolutionary’ distances, at least in so far as the phylogeny is correct. This, however, would require obtaining the files with the trees themselves from Sorhannus and Medlin. Perhaps I will send them another email.

Anyway, this was quite an eye-opener of a conversation and a very stark reminder of how much better it is to ask for help rather than to just bash your head against the wall alone. Other people are invariably smarter than oneself. It eventually resulted in this plot, which looks a lot more reasonable:

Now this, I think, is something I can work with. The correlation is still quite weak, but it’s now not ridiculous, and it’s not binned, and there’s still some vague positive association between the two variables, which does make sense.