This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Forum: Reproducible Science Short Example

I want to give a 20-30 minute talk on the subject of Reproducible Science in biology. What I think would be thrilling for the audience is that I do a self-contained small experiment live while coding, documenting, and commenting it.

I would like to start from an existing template. I have seen nice examples posted on blogs in the past that featured both code and results. Is there such an example that you would suggest for my talk?

Here are some ideas I have for the talk:

  1. Do analyses that are interesting for biologists
  2. Produce cool figures
  3. Keep it simple for biologists (e.g., avoid complicated analyses or packages)
  4. Start from an empty directory
  5. Get the data with wget
  6. Use R (most familiar language for biologists where I work)
  7. Use RStudio to make it more palatable than my green on black terminal :)
  8. Use a bash script (or a VERY simple make script) to run everything
  9. Create all the output dynamically (figures, tables...)
  10. Use GitHub to version control everything, including the data (smallish text file)
  11. Test the whole pipeline on another data set to show how much time this approach can save

Do you have any suggestions for such a talk?

Thanks in advance for your input

reproducible-science computing open-science
• 3,974 views • link
updated by Ram Ram • written by Eric Normandeau Eric Normandeau

3 answers

Here is a simple and short illustration using a Makefile: it downloads some DNA, extracts the protein accessions, downloads the protein, transforms it to SVG, convert to PNG and merge everything in a latex document:

• 0 views • link
updated by Ram Ram • written by Pierre Lindenbaum Pierre Lindenbaum

Thank you Pierre. It's definitely a cool example, but from the point of view of a biologist with basic notions of R, I don't think we can pretend this is even remotely simple. I am trying to lure them in, not scare them off ;)

• 0 views • link
written by Eric Normandeau Eric Normandeau

I would definitely suggest using RStudio with RMarkdown (and knitr) for your reproducible analysis example.

Do data download (within R), compute, explanation, figures and tables all within the same document. Some sources you could borrow from:

• 1 views • link
updated by Ram Ram • written by David Fredman David Fredman

I like the idea! Do you have a link to a blog post that would have done something similar to what I am planning?

• 1 views • link
written by Eric Normandeau Eric Normandeau

How about running through a differential expression analysis based on the DESeq2 vignette, with a read count matrix as input?

To showcase reuseability, re-run the analysis a second time using a different count matrix and design matrix but identical analysis/plotting code.

Read count matrices for various datasets are either built into R packages (e.g. Pasilla), or you could download them from e.g. the recount project

Edit: Actually, there's ongoing work in the Software Carpentry project to create a "capstone" showcase just like this to round off the novice lessons. See here for the Rmd and data (and, when more finished, the Software Carpentry main bc repo under novice/r/capstones).

• 0 views • link
updated by Ram Ram • written by David Fredman David Fredman

One example might be analysis to find the coverage over a gene or a transcript with RNA-seq data.

The analysis and the question are relatively simple yet slightly differing methods may easily lead to drastically different results.

• 0 views • link
written by donfreed donfreed

Log in to answer this question.

More posts like this