Compensation: Basic Principles

Similar documents
1 Some Statistical Basics.

A Fast Algorithm for Cosmic Rays Removal from Single Images

Homework on Properties of Galaxies in the Hubble Deep Field Name: Due: Friday, April 8 30 points Prof. Rieke & TA Melissa Halford

The full, blue supermoon is coming to the night sky near you

Homework #7: Properties of Galaxies in the Hubble Deep Field Name: Due: Friday, October points Profs. Rieke

GS Analysis of Microarray Data

arxiv:astro-ph/ v1 12 Nov 2003

GS Analysis of Microarray Data

Influence of moon light on tip-tilt sensing accuracy using STRAP

Assignment #9 Star Colors & the B-V Index

Measuring Colocalization within Fluorescence Microscopy Images

TEK 8 Test Review. 15. Galaxies are best described as -

Chapter 6 The Standard Deviation as a Ruler and the Normal Model

Earth Science: Second Quarter Grading Rubric Kindergarten

Hertzsprung-Russel Diagrams and Distance to Stars

The full, blue supermoon is coming to the night sky near you

Sun, Moon, and Stars. The Stars

AS 102 Lab The Luminosity of the Sun

Measuring Mitochondrial Membrane Potential with JC-1 Using the Cellometer Vision Image Cytometer

Observing Night 1 DON T FORGET TO KEEP NOTES ON YOUR PROCEDURE AND KEEP AN OBSERVING LOG FOR THE NIGHT!!!!!

Standardized and flexible eight colour flow cytometry panels harmonized between different. laboratories to study human NK cell phenotype and function

An introduction to closure phases

Ground Based Gravitational Microlensing Searches for Extra-Solar Terrestrial Planets Sun Hong Rhie & David Bennett (University of Notre Dame)

Why Don t We See Stars in the Daytime?

Laboratory Exercise. Atomic Spectra A Kirchoff Potpourri

The Kepler Mission: 20% of all Stars in the Milky Way Have Earth like Planets!

WFC3 IR Blobs, IR Sky Flats and the measured IR background levels

Lecture #8. Light-matter interaction. Kirchoff s laws

By Dr. Hany Farid, Dartmouth College

SEQUENCING THE STARS

Observing the Stars. radius: the distance from the center of a sphere to its surface; half its diameter. VY Canis Majoris

PMT Charge Response and Reconstruction of Position and Energy for a Reactor θ 13 Experiment

Parallax: Measuring the distance to Stars

Prelab 7: Sunspots and Solar Rotation

Announcements. Lecture 1 - Data and Data Summaries. Data. Numerical Data. all variables. continuous discrete. Homework 1 - Out 1/15, due 1/22

LAB: Photometry of the Pleiades Cluster

Decomposing white light into its components:

1. Why photons? 2. Photons in a vacuum

2. Discrete means unique, that other states don t overlap it. 3. Electrons in the outer electron shells have greater potential energy.

ACTIVITY 2 Exploring Light Patterns

AstroBITS: Open Cluster Project

Relationships between variables. Visualizing Bivariate Distributions: Scatter Plots

Verification of Continuous Forecasts

MAT Mathematics in Today's World

Stellar Masses: Binary Stars

Physics 6720 Introduction to Statistics April 4, 2017

Photoelectric Photometry of the Pleiades Student Manual

Data Processing in DES

The purpose of this visit is to investigate lunar phases. After this lab, the students will be able to demonstrate and apply these concepts:

Analysis of wavelength shifts reported for X-shooter spectra

SKINAKAS OBSERVATORY Astronomy Projects for University Students COLOUR IN ASTRONOMY

Complex Numbers. April 10, 2015

Apoptosis Detection Using the BD Accuri C6 Flow Cytometer

Spectrometers. Materials: Easy Spectrometer. Old CD Razor Index card Cardboard tube at least 10 inches long

Lecture 8. October 25, 2017 Lab 5

Advanced Spectroscopy Laboratory

Astronomical "color"

Stars and Planets GPS S4E1 A-D: Students will compare and contrast the physical attributes of stars, star patterns, and planets. Ms.

ASTRO 114 Lecture Okay. What we re going to discuss today are what we call radiation laws. We ve

AS 101: Day Lab #2 Summer Spectroscopy

The sun then falls below the horizon, and it is nighttime for our little person:

How hot is the Sun? hydrogen atom energy levels: Answer now (on your own):

Origin, early history, and fate of the Universe Does the Universe have a beginning? An end? What physics processes caused the Universe to be what it

Descriptive Statistics (And a little bit on rounding and significant digits)

X Rays must be viewed from space used for detecting exotic objects such as neutron stars and black holes also observing the Sun.

Astron 104 Laboratory #5 Colors of Stars

Cross-Talk in the ACS WFC Detectors. I: Description of the Effect

Chapter 8. Linear Regression. Copyright 2010 Pearson Education, Inc.

Open Cluster Research Project

VARIABLE STARS. What will you learn in this Lab?

At Home Phases Demo. Astronomy 210. Section 1 MWF Astronomy Building. Geocentric vs. Heliocentric system. The Motion of the Planets

More on Stochastics and the Phenomenon of Line-Edge Roughness

Test ABCDE. 1. What is the oldest era on the geological timescale? A. Precambrian B. Paleozoic C. Mesozoic D. Cenozoic

Using a Quantitative Measure to Achieve a Qualitative Lighting Design Effect! Introduction What is Luminance?

DAY LABORATORY EXERCISE: SPECTROSCOPY

NAME: VARIABLE STARS. What will you learn in this Lab?

Ensemble Data Assimilation and Uncertainty Quantification

High-dimensional regression

Clusters of galaxies

Phases of the Moon. Two perspectives: On Earth, or outside the Moon s orbit. More Phases. What if we zoom out? Phases of the Moon Demo 2/3/17

ACTIVITY 1. Exploring Light from Gases

Lecture 32. Lidar Error and Sensitivity Analysis

Observing Dark Worlds (Final Report)

Rational Numbers and Exponents

5. A particular star has an angle of parallax of 0.2 arcsecond. What is the distance to this star? A) 50 pc B) 2 pc C) 5 pc D) 0.

12.2 Dimensionality Reduction

1 The Preliminary Processing

A blue flame is hotter than a yellow one.

Quest Chapter 28. Remember: We see by color addition. How do we get white with color addition? So, is it a color?

Dark Skies Outreach to Sub-Saharan Africa Program

A Stellar Spectra 3. Stars shine at night (during the day too!). A star is a self-luminous sphere of gas. Stars are held together by gravity.

For instance, for a particular star cluster, these data were derived:

PHYS 352. On Measurement, Systematic and Statistical Errors. Errors in Measurement

Lecture 7: the Local Group and nearby clusters

Correlation. January 11, 2018

Experiment 4 Radiation in the Visible Spectrum

The Hertzsprung-Russell Diagram and Stellar Evolution

ADVANCED ANALYTICAL LAB TECH (Lecture) CHM

appstats27.notebook April 06, 2017

Southern African Large Telescope

Transcription:

Compensation: Basic Principles Bonus topic: Daytime Stargazing Mario Roederer ImmunoTechnology Section Vaccine Research Center, NIAID, NIH

Spectral Spillover Necessitates Compensation FITC PE Argon Laser FL1 FL2 450 500 550 600

Spectral Spillover Necessitates Compensation Total signal detected in FL1 FITC PE Unwanted signal detected in FL2 = roughly 15% Argon Laser FL1 FL2 450 500 550 600 True PE = Total FL2 15% FL1

Compensation Converts From Measurement to Expression Uncompensated Compensated PE - no stain FL2-15%FL1 PE - no stain FITC CD3 Measured Fluorescence FITC CD3 Calculated Expression

Compensation in Linear Domain Uncompensated Lines indicate cells that express equal amounts of an antigen Compensated

Compensation in Log Domain Uncompensated Compensated Lines indicate cells that express equal amounts of an antigen

Compensation is a Transformation Grid Lines Equal Measurement Equal Expression Bivariate Plots Histograms of each measurement

Compensation in 2 colors: Mostly aesthetic? Accurate identification and enumeration of subsets is still possible in two color experiments (but not visual estimation of expression levels ) 10 4 Uncompensated 10 4 Compensated 10 3 10 3 CD4 10 2 10 2 10 1 10 1 10 0 10 0 10 1 10 2 10 3 10 4 10 0 10 0 10 1 10 2 10 3 10 4 CD3 CD3

Compensation with 3 colors Accurate separation of populations can require a curved, two-dimensional surface This gets worse with more colors...

Compensation: Not just aesthetic Accurate discrimination of subsets is possible with uncompensated data However, this is true only when the expression of all antigens is uniform on each subset (e.g., CD45 / CD3 / CD4 / CD8) Separation of populations may require multidimensional surfaces.

Compensation for more colors: It s not just pretty pictures Spillover from unviewed measurement channel can alter event positions without obvious visual evidence (no diagnostic diagonals!) Thus, gate positions may depend on unviewed measurement channels and be different for various tubes in a panel Gates could also be highly dependant on which subset is being analyzed!

Compensation in the Polychromatic World 10 4 Uncompensated Compensated? APC (Unstained) e 10 3 10 2 10 1 10 3 10 2 10 1 10 0 10 0 10 1 10 2 10 3 10 4 Cy5.5APC CD57 10 0 10 0 10 1 10 2 10 3

Compensation in the Polychromatic World But is it correctly compensated? Compensated? 10 3 10 2 10 1 10 0 10 0 10 1 10 2 10 3

Compensation in the Polychromatic World But is it correctly compensated? Compensated? No! 10 3 10 2 10 1 10 0 10 0 10 1 10 2 10 3

Compensation in the Polychromatic World But is it correctly compensated? Compensated? No! 10 3 10 2 10 1 10 0 10 0 10 1 10 2 10 3

Compensation in the Polychromatic World But is it correctly compensated? Compensated? No! Yes! 10 3 10 2 10 1 10 0 10 0 10 1 10 2 10 3

Compensation in the Polychromatic World But is it correctly compensated? Compensated? No! Yes! 10 3 10 2 10 1 10 0 10 0 10 1 10 2 10 3

Compensation in the Polychromatic World But is it correctly compensated? Compensated? No! Yes! We can t tell just by looking! 10 3 10 2 10 1 10 0 10 0 10 1 10 2 10 3

What s Missing From This Picture?

What s Missing From This Picture? Where are the stars?

Why can t we see stars during the day? On cloudless nights, we see lots of stars. On cloudless days, we rarely see stars. We call this the Blue Sky Effect. But.. Why? Nowhere on the web is this commonly-asked question answered correctly.

The Blue Sky Effect The stars are easily visible at night (well, in most places). When the sun rises the sky becomes blue and the stars are no longer detectable.

The Blue Sky Effect The stars are easily visible at night (well, in most places). When the sun rises the sky becomes blue and the stars are no longer detectable.

Why Doesn t it Look Like This? But light is additive. The white light from the stars should add to the blue scattered light from the sun why don t we see white pinpricks against the blue background?

Signal Quantification Intensity 100 0

How Bright Are the Stars? Intensity 100 0 A star s intensity is the area under the histogram (i.e., counting the number of photons coming from the star).

How Bright Is The Sky? 1500 Intensity Y 1000 500 0 Intensity 100 0

Light Signals Are Additive 1500 Intensity 1000 500 0

Star Light Brightness Is Unchanged Intensity 1700 1650 1600 1550 1500 200 150 100 50 0 Background-subtracted Intensity We determine the intensity of a star by subtracting the background. Here, the sky (background) is 1500 units, so subtract 1500 units from the measurements. A star s intensity is now the area under the histogram (# of photons only from the star)

But Measurements Are Not Perfect Every measurement we make has an error associated with it. The variance in a measurement is the distribution of values that you would get if you made the same measurement on exactly the same object time and again (and the object did not change over time). Variances can arise from a variety of sources: e.g., inherent variability, measurement error. The final variance is the sum of all of these. Measurement error can arise from many sources. When counting photons, for example, the minimum error is governed by Poisson statistics.

Measurement Errors Light intensity measurements are, essentially, photon counting. The minimum error will be ± n, where n = # of photons counted. Number of counts Error Precision 100 ±10 10% 10,000 ±100 1% 1,000,000 ±1,000 0.1% As we increase the signal level, while the precision (CV) is getting better (smaller), the absolute error is increasing.

Measuring Sky Intensity 1500 Intensity 1000 500 0

Star Light is Still Additive 1500 Intensity 1000 500 0

Variation in Intensity Intensity 1600 1500 1400

Variation in Intensity Intensity 1600 1500 1400 We estimate the intensity of a star by subtracting the background. Here, the sky (background) is, on average, 1500 units, so subtract 1500 units from the measurements.

Estimating Star Light Intensity 100 Intensity 0-100

Estimating Star Light Intensity 100 Intensity 0-100 A star s intensity is, as before, the area under the histogram (# of photons only from the star)

Background Variance Ideal Measured Intensity 100 0-100 Intensity 100 0-100

Background Variance Ideal Measured Intensity 100 0-100 Intensity 100 0-100

Background Variance Is this negative light intensity? No. It is a negative estimated light intensity! And arbitrarily setting it to zero would introduce a systematic bias. Intensity 100 0-100

Why Can t We See Stars During the Day? Because the intensity of the star light is less than the variance of the intensity of the day light. Stars that are brighter than this variance can still be seen. The error on estimating the intensity of the star light will be at least as great as the variance of the day light... and can sometimes be <0.

Is This Related to Compensation? Completely. 0 100 Time since illumination Position in sky 29

Is This Related to Compensation? Completely. Compensation is the process by which we estimate a probe s contribution to a signal by subtracting estimated contributions from other sources. All the same principles apply: errors propagate; you cannot distinguish dim objects against high variance; you can get negative estimated signal intensities. 29

Imperfect Measurement Leads to Apparent Spread in Compensation 10 4 Uncompensated Compensated Spillover Fluorescence 10 3 10 2 1100 700 10 3 10 2 200 10 1 10 1 10 0 10 0 10 1 10 2 10 3 10 4 Primary Fluorescence 10 0 10 0 10 1 10 2 10 3 (-200)

Imperfect Measurement Leads to Apparent Spread in Compensation 10 4 Uncompensated Compensated 10 3 10 3 Spillover Fluorescence 10 2 10 2 10 1 10 1 10 0 10 0 10 1 10 2 10 3 10 4 Primary Fluorescence 10 0 10 0 10 1 10 2 10 3 Spreading manifests with a 1:2 Log:Log slope: square root relationship! Undercompensation manifests with a 1:1 Log:Log slope: linear relationship

The Actual Spread Background from PE-CD8 10 5 Lymphocytes <Cy7PE-A>: CD20 10 4 10 3 10 2 10 1 Variance contributed from background 10 1 10 2 10 3 10 4 10 5 <PE-A>: CD8

Compensation Does NOT Introduce nor Increase Error: Compensation Only Reveals It! The measurement error is already present. Compensation does not increase this error, it does not change it, it does not introduce any more error. Compensation simply makes the error more apparent by shifting it to the low end of the logscale. Staining controls are necessary to define gate placement!

Staining Controls Staining controls are necessary to identify cells which do or do not express a given antigen. The threshold for positivity may depend on the amount of fluorescence in other channels!

FMO Controls FMO controls are a much better way to identify positive vs. negative cells FMO controls can also help identify problems in compensation that are not immediately visible FMO controls should be used whenever accurate discrimination is essential or when antigen expression is relatively low

Compensation Controls The purpose of a compensation control is to allow software to estimate how much spillover is happening. This needs to be as precise as possible! Precise and accurate compensation requires determining the slope as precisely as possible. 10 5 10 5 B560-A PE 10 4 10 3 10 2 0 Cy7PE B710-A 10 4 10 3 10 2 0 0 10 2 10 3 10 4 10 5 B530-A 0 10 2 10 3 10 4 10 5 B530-A 0 10 2 10 3 10 4 10 5 B530-A FITC FITC FITC Compensation = 7.4% = 0.56%

Compensation Controls The purpose of a compensation control is to allow software to estimate how much spillover is happening. This needs to be as precise as possible! Precise and accurate compensation requires determining the slope as precisely as possible. Small absolute errors can be large relative errors in compensation and dramatically impact data 10 5 10 5 PE 10 4 10 3 10 2 0 Cy7PE 10 4 10 3 10 2 0 0 10 2 10 3 10 4 10 5 0 10 2 10 3 10 4 10 5 0 102 10 3 10 4 10 5 FITC FITC FITC Compensation = 0.56% Compensation = 0.86%

Compensation Controls There are only three rules to follow to ensure accurate compensation. All errors in compensation can be traced to a problem with one of these. 1. The experimental samples and control samples must be matched. Identical fluorochrome Identical instrument sensitivity (not settings) 2. The background (autofluorescence) for any given control s positive and negative events must be identical. 3. Each control s positive must be at least as bright as any experimental sample. To ensure precise compensation, collect enough events to precisely estimate median fluorescences in all channels: Beads: 5-10K Positive cells: 20-50K