<p><div><div><div><div data-width="100" data-show-faces="false" data-send="true" data-layout="button_count" data-height="21" data-href="http://www.r-bloggers.com/bioc-2016-conference-overview-and-few-ways-of-downloading-tcga-data/"></div></div><div><a data-url="http://www.r-bloggers.com/bioc-2016-conference-overview-and-few-ways-of-downloading-tcga-data/" data-counturl="http://www.r-bloggers.com/bioc-2016-conference-overview-and-few-ways-of-downloading-tcga-data/" data-count="horizontal" data-text="BioC 2016 Conference Overview and Few Ways of Downloading TCGA Data" class="twitter-share-button" data-via="rbloggers" href="https://twitter.com/share"></a></div><div></div></div><div></div></div>
<div>
(This article was first published on <strong><a href="http://r-addict.com/2016/07/22/BioC2016-RTCGA.html"> http://r-addict.com</a></strong>, and kindly contributed to <a rel="nofollow" href="http://www.r-bloggers.com/">R-bloggers)</a> </div>
<p><img src="http://i2.wp.com/r-addict.com/images/fulls/intro_rtcga_presentation.PNG?w=456" data-recalc-dims="1"> Few weeks ago I have a great pleasure of attending <a target="_blank" href="http://bioconductor.org/help/course-materials/2016/BioC2016/">BioC 2016: Where Software and Biology Connect</a> Conference at Stanford, where I have learned a lot! It wouldn’t be possible without the scholarship that I received from Bioconductor (organizers), which I deeply appreciate. It was an excellent place for software developers, statisticians and biologists to exchange their experiences and to better explain their work, as the understanding between collaborators in interdisciplinary teams is essential. In this post I present my thoughts and feelings about the event and I share the knowledge that I have learned during the event, i.e. about many ways of downloading <a target="_blank" href="http://cancergenome.nih.gov/">The Cancer Genome Atlas</a> data.</p>
<ul>
<li><a target="_blank" href="http://r-addict.com/2016/07/22/BioC2016-RTCGA.html#conference-overview">Conference overview</a>
<ul>
<li><a target="_blank" href="http://r-addict.com/2016/07/22/BioC2016-RTCGA.html#developer-day">Developer Day</a></li>
<li><a target="_blank" href="http://r-addict.com/2016/07/22/BioC2016-RTCGA.html#conference-talks">Conference Talks</a></li>
<li><a target="_blank" href="http://r-addict.com/2016/07/22/BioC2016-RTCGA.html#big-names">Big names</a></li>
<li><a target="_blank" href="http://r-addict.com/2016/07/22/BioC2016-RTCGA.html#workshops">Workshops</a></li>
</ul>
</li>
<li><a target="_blank" href="http://r-addict.com/2016/07/22/BioC2016-RTCGA.html#get-tcga-data">Get TCGA Data!</a>
<ul>
<li><a target="_blank" href="http://r-addict.com/2016/07/22/BioC2016-RTCGA.html#rtcga-factory-of-r-packages">RTCGA factory of R packages</a></li>
<li><a target="_blank" href="http://r-addict.com/2016/07/22/BioC2016-RTCGA.html#rtcgatoolbox-package">RTCGAToolbox package</a></li>
<li><a target="_blank" href="http://r-addict.com/2016/07/22/BioC2016-RTCGA.html#bigrquery---r-interface-to-google-bigquery">bigrquery – R interface to Google BigQuery</a></li>
</ul>
</li>
</ul>
<h1>Conference overview</h1>
<h3>Developer Day</h3>
<p>The day before the conference a Developer Day was organized. It was a good opportunity to find out what are currents projects about in Bioconductor and what are future plans for Bioconductor project. I didn’t know that most of Bioconductor developers are so keen on object oriented programming (I normally create tools for <a target="_blank" href="http://r-addict.com/2016/05/23/Informative-Survival-Plots.html">visualizations</a> and <a target="_blank" href="http://r-addict.com/2016/06/13/RHero-Saves-Backup-City-With-archivist-github.html">reproducible research</a>) and most of them were explaining improvements in Bioconductor S4 classes. Developer Day was held by <a target="_blank" href="https://www.roswellpark.edu/martin-morgan">Martin Morgan </a> who let everyone to pronounce themselves which was a spark reviving the networking! The funniest idea during this day was splitting attenders to 2 separate rooms (even though we could fit in 1 room), so most of the talks (and jokes!) have to be said twice, once in each room. This was cool idea besides that I, as a volunteer for lightning talk about RTCGA project (presentation – <a target="_blank" href="http://rtcga.github.io/RTCGA/BioC2016/">rtcga.github.io/RTCGA/BioC2016/</a>), couldn’t hear the presentations that were said in the other room as I spoke in mine room.</p>
<h3>Conference Talks</h3>
<p>The conference was divided into 2 days. Each day started with talks form invited and community speakers – you can find all presentations at the <a target="_blank" href="http://bioconductor.org/help/course-materials/2016/BioC2016/">conference’s website</a>. I really enjoyed Sandrine Dudoit speech about various approaches to normalization and clustering and metrics to select the best one (new package <a target="_blank" href="https://bioconductor.org/packages/devel/bioc/html/clusterExperiment.html">clusterExperiment</a>). Also interesting talk had Susan Holmes about <strong>Multicomponent data integration for the Human Microbiome</strong> but my favourite talk gave Jenny Bryan about <strong>Spreadsheets: the Data Format we Love to Hate</strong>. I’ve been overwhelmed by statistics she presented about Excel popularity again python and R. It seems that we’ll need to learn how to live with people using Excel because there are to many of them an R users are in minority!</p>
<h3>Big names</h3>
<p>Moreover, organizers have invited 2 big names: Rob Tibshirani and Dr. Robert Gentleman (Developer Day Keynote Address). Rob Tibshirani tried to convince everyone that most of the inference done from GLMnet feature selection was not right and proposed a new method, that is covered in his new package: <a target="_blank" href="https://cran.r-project.org/web/packages/selectiveInference/index.html">selectiveInference</a>. To my big surprise we also had a chance to listen to Ramnath Vaidya from Alteryx about simple and quick interface for creating a <a target="_blank" href="http://www.htmlwidgets.org/">htmlwidget</a> – interactive visualizations that does not have to be powered by <a target="_blank" href="http://shiny.rstudio.com/">shiny</a> applications. Ramnath created a new widget during his 15 min talk!</p>
<p>I would add Jim Hester to big names section, as his presentation sparked heated discussion about the future of installation packages from Bioconductor. New function <a target="_blank" href="https://github.com/hadley/devtools/blob/master/R/install-bioc.r"><code>devtools::install_bioc()</code></a> might be a competition for <code>BiocInstaller::biocLite()</code>!</p>
<h3>Workshops</h3>
<p>After lunch attenders had a chance to join workshops. I would like to recommend the <strong>Introduction to ImmuneSpaceR</strong> workshop by Renan Sauteraud from Fred Hutchinson Cancer Research Center for intermediate R users. I’ve learned a bit about R5 classes that Renan uses in his package and discovered many database from which I can drain as a data vampire with the help of his <a target="_blank" href="http://bioconductor.org/packages/release/bioc/html/ImmuneSpaceR.html">ImmuneSpaceR</a> package which is already in Bioconductor’s release branch. <a target="_blank" href="http://bioconductor.org/help/course-materials/2016/BioC2016/ConcurrentWorkshops1/Sauteraud.html">slides</a> <a target="_blank" href="https://github.com/SRenan/bioc2016">github</a></p>
<h1>Get TCGA Data!</h1>
<p>The Cancer Genome Atlas project is an effort to provide publicly available, high quality clinical and genomic information collected from tumour (and healthy) samples gathered from patients suffering from over 30 types of cancers. If you haven’t hear about this project you can check slides from <a target="_blank" href="https://rtcga.github.io/RTCGA/BioC2016/">my BioC2016 lightning talk</a> where you will find more helpful definitions and links.</p>
<h2>RTCGA factory of R packages</h2>
<p>I have started <a target="_blank" href="http://rtcga.github.io/RTCGA/">RTCGA</a> project almost year and a half ago. It provides software package called <a target="_blank" href="https://bioconductor.org/packages/release/bioc/html/RTCGA.html">RTCGA, which is available to download from Biocondcutor</a>. This package enables you to download any dataset from TCGA for any cohort type and for every release date (TCGA releases it’s dataset over time) and also provides datasets in the whole <a target="_blank" href="http://r-addict.com/2016/05/04/RTCGA-Quick-Guide.html">family of R data packages</a>.</p>
<p>To get information about available datasets in RTCGA see <code>?datasetsTCGA</code>, or check <code>?checkTCGA</code> to find out possible parameters’ values for <code>?downloadTCGA</code> to download any dataset if those already prepared for R do not satisfy your needs.</p>
<p>Today I’ve got a message from Dario Strbenac (author of <a target="_blank" href="http://bioconductor.org/packages/ClassifyR/">ClassifyR</a> package) who wrote that</p>
<blockquote>
<p>TCGA Data Portal has closed down and the data moved to Genomic Data Commons</p>
</blockquote>
<p>so you might encounter some issues with <code>downloadTCGA</code> for few days, but this is good news as I can close RTCGA project with putting final TCGA data snapshot into <code>ExperimentHub</code>. <a target="_blank" href="https://github.com/RTCGA/RTCGA/issues/85">This process</a> is already moderated by <a target="_blank" href="https://github.com/vobencha">Valerie Obenchain</a> on project’s <a target="_blank" href="https://github.com/RTCGA/">GitHub repository</a>.</p>
<h2>RTCGAToolbox package</h2>
<p>Almost at the same time as RTCGA reached Bioconductor, a similar package called <a target="_blank" href="https://www.bioconductor.org/packages/3.3/bioc/html/RTCGAToolbox.html">RTCGAToolbox</a> was published with suchlike functionality. It is as popular as RTCGA, if you look at the <a target="_blank" href="http://bioconductor.org/packages/stats/">packages’ download stats</a></p>
<p>To download data from <a target="_blank" href="http://gdac.broadinstitute.org/">Firehose</a> you could use <code>getFirehoseData</code> function but it allows only to downloaded 15 datatypes and the highest possible snapshot date is <code>2015-11-01</code> (<code>RTCGAToolbox::getFirehoseRunningDates(1)</code>) where <code>RTCGA::downloadTCGA</code> could download datasets even from <code>2016-01-28</code> (<code>tail(RTCGA::checkTCGA('Dates'),1)</code>).</p>
<h2>bigrquery – R interface to Google BigQuery</h2>
<p>The approach using <a target="_blank" href="https://cran.r-project.org/web/packages/bigrquery/index.html">bigrquery</a> was new for me. We have used this package during <strong>Facilitated Discussion: Approaches to Data Modeling</strong> at the Developer Day.</p>
<blockquote>
<p>Developer Day precedes the main conference on June 24, providing developers and would-be developers an opportunity to gain insights into project direction and software development best practice</p>
</blockquote>
<p>You can have a look at basic way of extracting clinical information with bigrquery</p>
<figure>
<pre><code data-lang="r"><span>library</span><span>(</span><span>devtools</span><span>)</span><span>
</span><span>devtools</span><span>::</span><span>install_github</span><span>(</span><span>&quot;rstats-db/bigrquery&quot;</span><span>,</span><span>
</span><span>upgrade_dependencies</span><span>=</span><span>FALSE</span><span>)</span><span>
</span><span>library</span><span>(</span><span>bigrquery</span><span>)</span><span>
</span><span>bq</span><span> </span><span>=</span><span> </span><span>src_bigquery</span><span>(</span><span>&quot;isb-cgc&quot;</span><span>,</span><span>
</span><span>&quot;tcga_201510_alpha&quot;</span><span>,</span><span>
</span><span>&quot;cgc-05-0011&quot;</span><span>)</span><span>
</span><span>bq</span><span>
</span><span>library</span><span>(</span><span>dplyr</span><span>)</span><span>
</span><span>bq</span><span> </span><span>%>%</span><span> </span><span>tbl</span><span>(</span><span>&quot;Clinical_data&quot;</span><span>)</span><span> </span><span>%>%</span><span>
</span><span>glimpse</span><span>()</span></code></pre>
</figure>
<p>This package looks to work faster than regular <code>downloadTCGA</code> but is not as fast as</p>
<figure>
<pre><code data-lang="r"><span>library</span><span>(</span><span>RTCGA.clicnial</span><span>)</span><span>
</span><span>dim</span><span>(</span><span>BRCA.clinical</span><span>)</span><span>
</span><span>?</span><span>clinical</span></code></pre>
</figure>
<p>for clinical observations.</p>
<p>To compare, bigrquery returns only 65 variables for clinical information where RTCGA returns all of them and bigrquery allows to extract only specific datatypes (approximately 10) where RTCGA allows downloading any datatype (for BRCA there are 43 possible datasets – <code>nrow(checkTCGA(what = 'DataSets', 'BRCA'))</code>).</p>
<p>If you are aware of any other method of downloading TCGA into R please let me know. You can write messages in the Disqus panel that is below. I would be happy to be able to compare them to RTCGA package.</p>
<p><div><div><div><div data-width="55" data-show-faces="false" data-send="true" data-layout="box_count" data-height="62" data-href="http://www.r-bloggers.com/bioc-2016-conference-overview-and-few-ways-of-downloading-tcga-data/"></div></div><div><a data-url="http://www.r-bloggers.com/bioc-2016-conference-overview-and-few-ways-of-downloading-tcga-data/" data-counturl="http://www.r-bloggers.com/bioc-2016-conference-overview-and-few-ways-of-downloading-tcga-data/" data-count="vertical" data-text="BioC 2016 Conference Overview and Few Ways of Downloading TCGA Data" class="twitter-share-button" data-via="rbloggers" href="https://twitter.com/share"></a></div><div></div></div><div></div></div>
<div>
<div>To <strong>leave a comment</strong> for the author, please follow the link and comment on their blog: <strong><a href="http://r-addict.com/2016/07/22/BioC2016-RTCGA.html"> http://r-addict.com</a></strong>.</div>
<hr>
<a rel="nofollow" href="http://www.r-bloggers.com/">R-bloggers.com</a> offers <strong><a rel="nofollow" href="http://feedburner.google.com/fb/a/mailverify?uri=RBloggers">daily e-mail updates</a></strong> about <a title="The R Project for Statistical Computing" rel="nofollow" href="http://www.r-project.org/">R</a> news and <a title="R tutorials" rel="nofollow" href="http://www.r-bloggers.com/?s=tutorial">tutorials</a> on topics such as: <a title="Data science" rel="nofollow" href="http://www.r-bloggers.com/?s=data%20science">Data science</a>, <a title="Big Data" rel="nofollow" href="http://www.r-bloggers.com/?s=Big%20Data">Big Data, <a title="R jobs" rel="nofollow" href="http://www.r-users.com/">R jobs</a>, visualization (<a title="ggplot and ggplot2 tutorials" rel="nofollow" href="http://www.r-bloggers.com/?s=ggplot2">ggplot2</a>, <a title="Boxplots using lattice and ggplot2 tutorials" rel="nofollow" href="http://www.r-bloggers.com/?s=boxplot">Boxplots</a>, <a title="Maps and gis" rel="nofollow" href="http://www.r-bloggers.com/?s=map">maps</a>, <a title="Animation in R" rel="nofollow" href="http://www.r-bloggers.com/?s=animation">animation</a>), programming (<a title="RStudio IDE for R" rel="nofollow" href="http://www.r-bloggers.com/?s=RStudio">RStudio</a>, <a title="Sweave and literate programming" rel="nofollow" href="http://www.r-bloggers.com/?s=sweave">Sweave</a>, <a title="LaTeX in R" rel="nofollow" href="http://www.r-bloggers.com/?s=LaTeX">LaTeX</a>, <a title="SQL and databases" rel="nofollow" href="http://www.r-bloggers.com/?s=SQL">SQL</a>, <a title="Eclipse IDE for R" rel="nofollow" href="http://www.r-bloggers.com/?s=eclipse">Eclipse</a>, <a title="git and github, Version Control System" rel="nofollow" href="http://www.r-bloggers.com/?s=git">git</a>, <a title="Large data in R using Hadoop" rel="nofollow" href="http://www.r-bloggers.com/?s=hadoop">hadoop</a>, <a title="Web Scraping of google, facebook, yahoo, twitter and more using R" rel="nofollow" href="http://www.r-bloggers.com/?s=Web+Scraping">Web Scraping</a>) statistics (<a title="Regressions and ANOVA analysis tutorials" rel="nofollow" href="http://www.r-bloggers.com/?s=regression">regression</a>, <a title="principal component analysis tutorial" rel="nofollow" href="http://www.r-bloggers.com/?s=PCA">PCA</a>, <a title="Time series" rel="nofollow" href="http://www.r-bloggers.com/?s=time+series">time series</a>, <a title="finance trading" rel="nofollow" href="http://www.r-bloggers.com/?s=trading">trading</a>) and more...
<img height="1" alt="" width="1" src="http://feeds.feedburner.com/~r/RBloggers/~4/-bZ_uW8hxHg"></a></div>