Computational biology and genomics
I'm going to include my blog as a form of community engagement in my fellowship application. The only objective way I can measure its impact is by the number of views. Luckily the WordPress API makes it easy to download all my web traffic and I have made this data available online. I wrote this...
As promised two years ago, here's a short blog post on Snakemake. I have been using Bpipe to manage my workflows/pipelines but Snakemake has been mentioned to me on more than one occasion; in particular: @davetang31 seems useful if testing many diff types of pipelines. snakemake seems a bit more practical if analyzing many samples...
What's the difference between deleteriousness and pathogenicity? Box 1 from the article "Guidelines for investigating causality of sequence variants in human disease" has the following definitions:
The American College of Medical Genetics and Genomics (ACMG) have recommended that genetic variants that may be pathogenic or likely pathogenic in certain genes should be reported back to the patient. The latest list of genes can be found here. How do I assess whether a variant is pathogenic or likely pathogenic? Use this tool,...
I'm a fan of GEMINI and have been using it for for a year and a half for various exome projects. I have written two scripts that can generate variant summaries from a GEMINI database. I prefer bar plots over the pie charts created by VEP. A summary pie chart created by VEP.
The diagram package makes it easy to create flowcharts in R. In this post I'll show an example of creating a simple flowchart. The most important part is to understand how the coordinate systems works; once you understand that, it's just a matter of placing your arrows and boxes accordingly to create your flowchart. To...
A couple of weeks ago, I wrote a post on identifying OMIM phenotypes that are associated with a gene of interest. I thought I solved the problem by using one of my favourite R packages (biomaRt) but alas. For example, I could not find any OMIM IDs associated with the TTN gene using biomaRt. In...
An adjacency list is simply an unordered list that describes connections between vertices. It's a commonly used input format for graphs. In this post, I use the melt() function from the reshape2 package to create an adjacency list from a correlation matrix. I use the geneData dataset, which consists of real but anonymised microarray expression...