search

hackseq / 2017_project_5

Developing advanced R tutorials for genomic data analysis
https://hackseq.github.io/2017_project_5/
MIT License
1 stars 2 forks source link

Showing how to use R to make a fully reproducible pipeline #5

Open privefl opened 6 years ago

privefl commented 6 years ago

A biological analysis is sometimes more appropriately called a pipeline. This is because it generally consists of many steps, using many different software and data formats. Yet, these analysis pipelines are becoming very complex and usually makes use of many bash/perl scripts. For people like me who don't really know that much bash or perl, it can be really hard to understand those scripts.

What is important in these pipelines? To list what comes to my mind:

I think we can do each of these operations in R. And I think we should.

The main reason would be to put all your analysis in a single notebook where you have all your code, results and possibly some writing. Using notebooks is good practice and makes it possible to have a fully reproducible analysis, which will a standard in years to come. Another reason is simply that it's easier!

In this tutorial, I'll show an example of a moderately complex analysis of the 1000 Genomes data, all in R.

You can find the first version of the tuto there.

zhenyisong commented 6 years ago

I cann't agree u more. I am now trying to use this style to maintain my code, which will be sent to my collaborators.