Day 3
10
Lecture 5 - Data integration and batch effect correction
Single-cell RNAseq analysis with R/Bioconductor
Welcome
Program
RStudio
Prerequisites
Day 1
1
Lecture 1 - Introduction to scRNAseq analysis
2
Lecture 2 - From sequencing reads to expression matrices
3
Lab 1: Familiarizing yourself with the course AWS instance
4
Lab 2: From .bcl to count matrix
Day 2
5
Lecture 3 - Quality control for scRNA-Seq data
6
Lab 3: Introduction to R/Bioconductor
7
Lab 4 - Single-cell RNA-seq data wrangling
Day 3
8
Lecture 4 - Identifying cell populations
9
Lab 5: Dimension reduction, clustering and annotation
10
Lecture 5 - Data integration and batch effect correction
11
Lab 6: Batch correction
Day 4
12
Lecture 6 - Advances in single-cell genomics: the epigenome
13
Lab 7: Single-cell ATAC-seq analysis workflow
14
Lecture 7 - Trajectories and pseudotimes
15
Lab 8: Pseudotime analyses
Day 5
16
Lecture 8 - Advances in single-cell genomics: spatial transcriptomics
Extra resources
10
Lecture 5 - Data integration and batch effect correction
Slides here
<div id="quarto-navigation-envelope" class="hidden"> <p><span class="hidden" data-render-id="quarto-int-sidebar-title">Single-cell RNAseq analysis with R/Bioconductor</span> <span class="hidden" data-render-id="quarto-int-navbar-title">Single-cell RNAseq analysis with R/Bioconductor</span> <span class="hidden" data-render-id="quarto-int-next"><span class="chapter-number">11</span> <span class="chapter-title">Lab 6: Batch correction</span></span> <span class="hidden" data-render-id="quarto-int-prev"><span class="chapter-number">9</span> <span class="chapter-title">Lab 5: Dimension reduction, clustering and annotation</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:/index.html">Welcome</span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/program.html">Program</span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/Rstudio.html">RStudio</span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/prerequisites.html">Prerequisites</span> <span class="hidden" data-render-id="quarto-int-sidebar:quarto-sidebar-section-1">Day 1</span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day1/Lecture1_introduction.html"><span class="chapter-number">1</span> <span class="chapter-title">Lecture 1 - Introduction to scRNAseq analysis</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day1/Lecture2_Bcl2matrix.html"><span class="chapter-number">2</span> <span class="chapter-title">Lecture 2 - From sequencing reads to expression matrices</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day1/Lab1_Intro.html"><span class="chapter-number">3</span> <span class="chapter-title">Lab 1: Familiarizing yourself with the course AWS instance</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day1/Lab2_processingreads.html"><span class="chapter-number">4</span> <span class="chapter-title">Lab 2: From .bcl to count matrix</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:quarto-sidebar-section-2">Day 2</span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day2/Lecture3_qualitycontrol.html"><span class="chapter-number">5</span> <span class="chapter-title">Lecture 3 - Quality control for scRNA-Seq data</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day2/Lab3_Rbioc.html"><span class="chapter-number">6</span> <span class="chapter-title">Lab 3: Introduction to R/Bioconductor</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day2/Lab4_data_wrangling_scRNAseq.html"><span class="chapter-number">7</span> <span class="chapter-title">Lab 4 - Single-cell RNA-seq data wrangling</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:quarto-sidebar-section-3">Day 3</span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day3/Lecture4_clustering.html"><span class="chapter-number">8</span> <span class="chapter-title">Lecture 4 - Identifying cell populations</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day3/Lab5_clustering.html"><span class="chapter-number">9</span> <span class="chapter-title">Lab 5: Dimension reduction, clustering and annotation</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day3/Lecture5_batchcorrection.html"><span class="chapter-number">10</span> <span class="chapter-title">Lecture 5 - Data integration and batch effect correction</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day3/Lab6_batch_correction.html"><span class="chapter-number">11</span> <span class="chapter-title">Lab 6: Batch correction</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:quarto-sidebar-section-4">Day 4</span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day4/Lecture6_ATAC.html"><span class="chapter-number">12</span> <span class="chapter-title">Lecture 6 - Advances in single-cell genomics: the epigenome</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day4/Lab7_atac-seq.html"><span class="chapter-number">13</span> <span class="chapter-title">Lab 7: Single-cell ATAC-seq analysis workflow</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day4/Lecture7_pseudotime.html"><span class="chapter-number">14</span> <span class="chapter-title">Lecture 7 - Trajectories and pseudotimes</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day4/Lab8_pseudotime.html"><span class="chapter-number">15</span> <span class="chapter-title">Lab 8: Pseudotime analyses</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:quarto-sidebar-section-5">Day 5</span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/day5/Lecture8_spatial-transcriptomics.html"><span class="chapter-number">16</span> <span class="chapter-title">Lecture 8 - Advances in single-cell genomics: spatial transcriptomics</span></span> <span class="hidden" data-render-id="quarto-int-sidebar:/content/extra.html">Extra resources</span> <span class="hidden" data-render-id="footer-left">Single-cell RNAseq analysis with R/Bioconductor |<br /> J. Serizay, O. Ashenberg, F. Almeida-Silva</span> <span class="hidden" data-render-id="footer-right">This book was built with <a href="https://quarto.org/">Quarto</a>.</span> <span class="hidden" data-render-id="quarto-breadcrumbs-223e2965264ed5a7335008666fb8bf98">Day 3</span> <span class="hidden" data-render-id="quarto-breadcrumbs-4e5f6a1e7edd6ec4395d31fd7681eae5"><span class="chapter-number">10</span> <span class="chapter-title">Lecture 5 - Data integration and batch effect correction</span></span></p> </div> <div id="quarto-meta-markdown" class="hidden"> <p><span class="hidden" data-render-id="quarto-metatitle">Single-cell RNAseq analysis with R/Bioconductor - <span class="chapter-number">10</span> <span class="chapter-title">Lecture 5 - Data integration and batch effect correction</span></span> <span class="hidden" data-render-id="quarto-twittercardtitle">Single-cell RNAseq analysis with R/Bioconductor - <span class="chapter-number">10</span> <span class="chapter-title">Lecture 5 - Data integration and batch effect correction</span></span> <span class="hidden" data-render-id="quarto-ogcardtitle">Single-cell RNAseq analysis with R/Bioconductor - <span class="chapter-number">10</span> <span class="chapter-title">Lecture 5 - Data integration and batch effect correction</span></span> <span class="hidden" data-render-id="quarto-metasitename">Single-cell RNAseq analysis with R/Bioconductor</span> <span class="hidden" data-render-id="quarto-twittercarddesc"></span> <span class="hidden" data-render-id="quarto-ogcardddesc"></span></p> </div> </main> <!-- /main --> <script id = "quarto-html-after-body" type="application/javascript"> window.document.addEventListener("DOMContentLoaded", function (event) { const toggleBodyColorMode = (bsSheetEl) => { const mode = bsSheetEl.getAttribute("data-mode"); const bodyEl = window.document.querySelector("body"); if (mode === "dark") { bodyEl.classList.add("quarto-dark"); bodyEl.classList.remove("quarto-light"); } else { bodyEl.classList.add("quarto-light"); bodyEl.classList.remove("quarto-dark"); } } const toggleBodyColorPrimary = () => { const bsSheetEl = window.document.querySelector("link#quarto-bootstrap"); if (bsSheetEl) { toggleBodyColorMode(bsSheetEl); } } toggleBodyColorPrimary(); const icon = ""; const anchorJS = new window.AnchorJS(); anchorJS.options = { placement: 'right', icon: icon }; anchorJS.add('.anchored'); const isCodeAnnotation = (el) => { for (const clz of el.classList) { if (clz.startsWith('code-annotation-')) { return true; } } return false; } const clipboard = new window.ClipboardJS('.code-copy-button', { text: function(trigger) { const codeEl = trigger.previousElementSibling.cloneNode(true); for (const childEl of codeEl.children) { if (isCodeAnnotation(childEl)) { childEl.remove(); } } return codeEl.innerText; } }); clipboard.on('success', function(e) { // button target const button = e.trigger; // don't keep focus button.blur(); // flash "checked" button.classList.add('code-copy-button-checked'); var currentTitle = button.getAttribute("title"); button.setAttribute("title", "Copied!"); let tooltip; if (window.bootstrap) { button.setAttribute("data-bs-toggle", "tooltip"); button.setAttribute("data-bs-placement", "left"); button.setAttribute("data-bs-title", "Copied!"); tooltip = new bootstrap.Tooltip(button, { trigger: "manual", customClass: "code-copy-button-tooltip", offset: [0, -8]}); tooltip.show(); } setTimeout(function() { if (tooltip) { tooltip.hide(); button.removeAttribute("data-bs-title"); button.removeAttribute("data-bs-toggle"); button.removeAttribute("data-bs-placement"); } button.setAttribute("title", currentTitle); button.classList.remove('code-copy-button-checked'); }, 1000); // clear code selection e.clearSelection(); }); function tippyHover(el, contentFn) { const config = { allowHTML: true, content: contentFn, maxWidth: 500, delay: 100, arrow: false, appendTo: function(el) { return el.parentElement; }, interactive: true, interactiveBorder: 10, theme: 'quarto', placement: 'bottom-start' }; window.tippy(el, config); } const noterefs = window.document.querySelectorAll('a[role="doc-noteref"]'); for (var i=0; i<noterefs.length; i++) { const ref = noterefs[i]; tippyHover(ref, function() { // use id or data attribute instead here let href = ref.getAttribute('data-footnote-href') || ref.getAttribute('href'); try { href = new URL(href).hash; } catch {} const id = href.replace(/^#\/?/, ""); const note = window.document.getElementById(id); return note.innerHTML; }); } let selectedAnnoteEl; const selectorForAnnotation = ( cell, annotation) => { let cellAttr = 'data-code-cell="' + cell + '"'; let lineAttr = 'data-code-annotation="' + annotation + '"'; const selector = 'span[' + cellAttr + '][' + lineAttr + ']'; return selector; } const selectCodeLines = (annoteEl) => { const doc = window.document; const targetCell = annoteEl.getAttribute("data-target-cell"); const targetAnnotation = annoteEl.getAttribute("data-target-annotation"); const annoteSpan = window.document.querySelector(selectorForAnnotation(targetCell, targetAnnotation)); const lines = annoteSpan.getAttribute("data-code-lines").split(","); const lineIds = lines.map((line) => { return targetCell + "-" + line; }) let top = null; let height = null; let parent = null; if (lineIds.length > 0) { //compute the position of the single el (top and bottom and make a div) const el = window.document.getElementById(lineIds[0]); top = el.offsetTop; height = el.offsetHeight; parent = el.parentElement.parentElement; if (lineIds.length > 1) { const lastEl = window.document.getElementById(lineIds[lineIds.length - 1]); const bottom = lastEl.offsetTop + lastEl.offsetHeight; height = bottom - top; } if (top !== null && height !== null && parent !== null) { // cook up a div (if necessary) and position it let div = window.document.getElementById("code-annotation-line-highlight"); if (div === null) { div = window.document.createElement("div"); div.setAttribute("id", "code-annotation-line-highlight"); div.style.position = 'absolute'; parent.appendChild(div); } div.style.top = top - 2 + "px"; div.style.height = height + 4 + "px"; let gutterDiv = window.document.getElementById("code-annotation-line-highlight-gutter"); if (gutterDiv === null) { gutterDiv = window.document.createElement("div"); gutterDiv.setAttribute("id", "code-annotation-line-highlight-gutter"); gutterDiv.style.position = 'absolute'; const codeCell = window.document.getElementById(targetCell); const gutter = codeCell.querySelector('.code-annotation-gutter'); gutter.appendChild(gutterDiv); } gutterDiv.style.top = top - 2 + "px"; gutterDiv.style.height = height + 4 + "px"; } selectedAnnoteEl = annoteEl; } }; const unselectCodeLines = () => { const elementsIds = ["code-annotation-line-highlight", "code-annotation-line-highlight-gutter"]; elementsIds.forEach((elId) => { const div = window.document.getElementById(elId); if (div) { div.remove(); } }); selectedAnnoteEl = undefined; }; // Attach click handler to the DT const annoteDls = window.document.querySelectorAll('dt[data-target-cell]'); for (const annoteDlNode of annoteDls) { annoteDlNode.addEventListener('click', (event) => { const clickedEl = event.target; if (clickedEl !== selectedAnnoteEl) { unselectCodeLines(); const activeEl = window.document.querySelector('dt[data-target-cell].code-annotation-active'); if (activeEl) { activeEl.classList.remove('code-annotation-active'); } selectCodeLines(clickedEl); clickedEl.classList.add('code-annotation-active'); } else { // Unselect the line unselectCodeLines(); clickedEl.classList.remove('code-annotation-active'); } }); } const findCites = (el) => { const parentEl = el.parentElement; if (parentEl) { const cites = parentEl.dataset.cites; if (cites) { return { el, cites: cites.split(' ') }; } else { return findCites(el.parentElement) } } else { return undefined; } }; var bibliorefs = window.document.querySelectorAll('a[role="doc-biblioref"]'); for (var i=0; i<bibliorefs.length; i++) { const ref = bibliorefs[i]; const citeInfo = findCites(ref); if (citeInfo) { tippyHover(citeInfo.el, function() { var popup = window.document.createElement('div'); citeInfo.cites.forEach(function(cite) { var citeDiv = window.document.createElement('div'); citeDiv.classList.add('hanging-indent'); citeDiv.classList.add('csl-entry'); var biblioDiv = window.document.getElementById('ref-' + cite); if (biblioDiv) { citeDiv.innerHTML = biblioDiv.innerHTML; } popup.appendChild(citeDiv); }); return popup.innerHTML; }); } } }); </script> <nav class="page-navigation"> <div class="nav-page nav-page-previous"> <a href="/content/day3/Lab5_clustering.html" class="pagination-link"> <i class="bi bi-arrow-left-short"></i> <span class="nav-page-text"><span class='chapter-number'>9</span> <span class='chapter-title'>Lab 5: Dimension reduction, clustering and annotation</span></span> </a> </div> <div class="nav-page nav-page-next"> <a href="/content/day3/Lab6_batch_correction.html" class="pagination-link"> <span class="nav-page-text"><span class='chapter-number'>11</span> <span class='chapter-title'>Lab 6: Batch correction</span></span> <i class="bi bi-arrow-right-short"></i> </a> </div> </nav> </div> <!-- /content --> <footer class="footer"> <div class="nav-footer"> <div class="nav-footer-left"> <div class='footer-contents'>Single-cell RNAseq analysis with R/Bioconductor | J. Serizay, O. Ashenberg, F. Almeida-Silva </div> </div> <div class="nav-footer-center"> </div> <div class="nav-footer-right"> <div class='footer-contents'>This book was built with <a href="https://quarto.org/">Quarto</a>. </div> </div> </div> </footer> </body> </html>