Raman Spectra of Chondrocytes in Cartilage: hyperspec s chondro data set
|
|
- Lucinda Stanley
- 6 years ago
- Views:
Transcription
1 Raman Spectra of Chondrocytes in Cartilage: hyperspec s chondro data set Claudia Beleites <chemometrie@beleites.de> DIA Raman Spectroscopy Group, University of Trieste/Italy (2 28) Spectroscopy Imaging, IPHT, Jena/Germany (28 216) ÖPV, JKI, Berlin/Germany Chemometric Consulting and Chemometrix GmbH, Wölfersheim/Germany October, 217 Reproducing this vignette The data set and source file of this vignette are available at hyperspec s github page at and Renishaw/chondro.txt (ca. 31 MB). The.Rnw looks for rawdata/chondro.txt, so please save the data file in subdirectory rawdata of your working directory. In order to reproduce the examples by typing in the commands, have a look at the definitions of color palettes used in this document are defined in vignettes.defs. Also, the package hyperspec needs to be loaded first via library (hyperspec). Contents 1 Introduction 2 2 The Data Set 2 3 Data Import 3 4 Preprocessing Aligning the Spatial Grid Spectral Smoothing Baseline Correction Normalization Subtracting the Overall Composition Outlier Removal by Principal Component Analysis (PCA) Hierarchical Cluster Analysis (HCA) 7 6 Plotting a False-Colour Map of Certain Spectral Regions 8 1
2 Figure 1 Microphotograph of the cartilage section. The frame indicates the measurement area (3 2 µm). 7 Saving the data set 9 1 Introduction This vignette describes the chondro data set. It shows a complete data analysis work flow on a Raman map demonstrating frequently needed preprocessing methods baseline correction normalization smoothing / interpolating spectra preprocessing the spatial grid and other basic work techniques plotting spectra plotting false color maps cutting the spectral range, selecting (extracting) or deleting spectra, and aggregating spectra (e.g. calculating cluster mean spectra). The chemometric methods used are Principal Component Analysis (PCA) and Hierarchical Cluster Analysis, showing how to use data analysis procedures provided by R and other packages. 2 The Data Set Raman spectra of a cartilage section were measured on each point of a grid, resulting in a so-called Raman map. Figure 1 shows a microscope picture of the measured area and its surroundings. The measurement parameters were: Excitation wavelength: 633 nm Exposure time: s per spectrum Objective:, NA.8 2
3 Measurement grid: 3 2 µm, 1 µm step size Spectrometer: Renishaw InVia 3 Data Import Renishaw provides a converter to export their proprietary data in a so-called long format ASCII file. Raman maps are exported having four columns, y, x, raman shift, and intensity. hyperspec comes with a function to import such files, read.txt.renishaw. The function assumes a map as default, but can also handle single spectra (data = "spc"), time series (data = "ts"), and depth profiles (data = "depth"). In addition, large files may be processed in chunks. In order to speed up the reading read.txt.renishaw does not allow missing values, but it does work with NA. > chondro <- read.txt.renishaw ("rawdata/chondro.txt", data = "xyspc") > chondro hyperspec object 87 spectra 4 data columns 1272 data points / spectrum wavelength: Delta * tilde(nu)/cm^-1 [numeric] data: (87 rows x 4 columns) 1. y: y/(mu * m) [numeric] x: x/(mu * m) [numeric] spc: I / a.u. [matrix1272] NA 4. filename: filename [character] rawdata/chondro.txt rawdata/chondro.txt... rawdata/chondro.txt To get an overview of the spectra (figure 2a): > plot (chondro, "spcprctl") A mean intensity map (figure 2b) is produced by: > plotmap (chondro, func.args = list (na.rm = TRUE), col.regions = seq.palette (2)) plotmap applies a function to squeeze all spectral intensities into a summary characteristic for the whole spectrum. This function defaults to the mean. Further arguments that should be handed to this function can be given in list func.args. As the raw data contains NAs due to deleting cosmic ray spikes, this argument is needed here. I / a.u. y µm ν ~ cm 1 (a) The raw spectra: median, 16 th and 84 th, and th and 9 th percentile spectra. 1 2 x µm (b) The sum intensity of the raw spectra. Figure 2 The raw spectra. 3
4 y y y 4 Preprocessing As usual in Raman spectroscopy of biological tissues, the spectra need some preprocessing. 4.1 Aligning the Spatial Grid The spectra were acquired on a regular square grid of measurement positions over the sample. hyperspec however does not store the data on this grid but rather stores the position where the spectra were taken. This allows more handling of irregular point patterns, and of sparse measurements where not a complete (square) is retained but only few positions on this grid. Occasionally, numeric or rounding errors will lead to seemingly irregular spacing of the data points. This is most obvious in false-colour maps of the sample (fig. fig:raw:b), as plotmap by default uses panel.levelplot.raster which assumes a regular grid is underlying the data. The symptoms are warnings x values are not equispaced; output may be wrong and white stripes in the false colour map (fig 3a) which occur even thought the points are almost at their correct place (fig 3b) > ## disturb a few points > chondro$x [] <- chondro$x [] + rnorm (1, sd =.1) > chondro$y [66] <- chondro$y [66] + rnorm (1, sd =.1) > plotmap (chondro, col.regions = seq.palette (2)) > require ("latticeextra") > plotmap (chondro, col.regions = seq.palette (2), panel = panel.levelplot.points, + col = NA, pch = 22, cex = 1.9) Such slight rounding errors can be corrected by fitraster or makeraster. fitraster needs the step size of the raster and a starting coordinate, whereas makeraster tries to guess these parameters for fitraster. As long as just a few points are affected by rounding errors, makeraster works fine: > chondro$x <- fitraster (chondro$x)$x > chondro$y <- fitraster (chondro$y)$x the result is shown in figure 3c. > plotmap (chondro, col.regions = seq.palette (2)) x (a) Off-grid points cause stripes 1 2 x (b) Slightly wrong coordinages 2 x (c) Corrected grid Figure 3 Rounding erros in the point coordinates of lateral measurement grid. 4.2 Spectral Smoothing As the overview shows that the spectra contain NAs (from cosmic spike removal that was done previously), the first step is to remove these. Also, the wavelength axis of the raw spectra is not 4
5 evenly spaced (the data points are between.8 and 1 cm 1 apart from each other). Furthermore, it would be good to trade some spectral resolution for higher signal to noise ratio. All three of these issues are tackled by interpolating and smoothing of the wavelength axis by spc.loess. The resolution is to be reduced to 8 cm 1, or 4 cm 1 data point spacing. > chondro <- spc.loess (chondro, seq (62, 18, 4)) > chondro hyperspec object 87 spectra 4 data columns 3 data points / spectrum wavelength: Delta * tilde(nu)/cm^-1 [numeric] data: (87 rows x 4 columns) 1. y: y [numeric] x: x [numeric] spc: I / a.u. [matrix3] filename: filename [character] rawdata/chondro.txt rawdata/chondro.txt... rawdata/chondro.txt The spectra are now the same as in the data set chondro. However, the data set also contains the clustering results (see at the very end of this document). They are stored for saving as the distributed demo data: > spectra.to.save <- chondro 4.3 Baseline Correction The next step is a linear baseline correction. spc.fit.poly.below tries to automatically find appropriate support points for polynomial baselines. The default is a linear baseline, which is appropriate in our case: > baselines <- spc.fit.poly.below (chondro) > chondro <- chondro - baselines I / a.u I / a.u ν ~ cm 1 (a) The spectra after smoothing, baseline correction, and normalization ν ~ cm 1 (b) The spectra after subtracting the th percentile spectrum. Figure 4 The preprocessed spectra. 4.4 Normalization As the spectra are quite similar, area normalization should work well:. > chondro <- chondro / rowmeans (chondro) > plot (chondro, "spcprctl")
6 Note that normalization effectively cancels the information of one variate (wavelength), it introduces a collinearity. If needed, this collinearity can be broken by removing one of the variates involved in the normalization. Note that the chondro object shipped with hyperspec set has multiple collinearities as only the first principal components are shipped (see below). For the results of these preprocessing steps, see figure 4a. 4. Subtracting the Overall Composition The spectra are very homogeneous, but I m interested in the differences between the different regions of the sample. Subtracting the minimum spectrum cancels out the matrix composition that is common to all spectra. But the minimum spectrum also picks up a lot of noise. So instead, the th percentile spectrum is subtracted: > chondro <- chondro - quantile (chondro,.) > plot (chondro, "spcprctl") The resulting data set is shown in figure 4b. Some interesting differences start to show up: there are distinct lipid bands in some but not all of the spectra. 4.6 Outlier Removal by Principal Component Analysis (PCA) PCA is a technique that decomposes the data into scores and loadings (virtual spectra). It is known to be quite sensitive to outliers. Thus, I use it for outlier detection. The resulting scores and loadings are put again into hyperspec objects by decomposition: > pca <- prcomp (chondro, center = TRUE) > scores <- decomposition (chondro, pca$x, label.wavelength = "PC", + label.spc = "score / a.u.") > loadings <- decomposition (chondro, t(pca$rotation), scores = FALSE, + label.spc = "loading I / a.u.") Plotting the scores of each PC against all other gives a good idea where to look for outliers. > pairs (scores [[,,1:2]], pch = 19, cex =.) Now the spectra can be found either by plotting two scores against each other (by plot) and identifying with identify, or they can be identified in the score map by map.identify. There is also a function to identify spectra in a spectra plot, spc.identify, which could be used to identify principal components that are heavily influenced e. g. by cosmic ray spikes. > ## omit the first 4 PCs > out <- map.identify (scores [,,]) > out <- c (out, map.identify (scores [,,6])) > out <- c (out, map.identify (scores [,,7])) > out [1] > outcols <- c ("red", "blue", "#88", "orange", "magenta", "brown") > cols <- rep ("black", nrow(chondro)) > cols [out] <- outcols We can check our findings by comparing the spectra to the bulk of spectra (figure a): 6
7 > plot(chondro[1], plot.args = list (ylim = c (1, length (out) +.7)), + lines.args = list( type = "n")) > for (i in seq (along = out)){ + plot(chondro, "spcprctl", yoffset = i, add = TRUE, col = "gray") + plot (chondro [out[i]], yoffset = i, col = outcols[i], add = TRUE, + lines.args = list (lwd = 2)) + text (6, i +.33, out [i]) + } and also by looking where these spectra appear in the scores pairs plot (figure b): > png ("chondro-fig--pca-pairs2.png", width =, height = ) > pch <- rep (46L, nrow (chondro)) > pch [out] <- 19L > pairs (scores [[,,1:7]], pch = pch, cex = 1, col = cols) > dev.off () Finally, the outliers are removed: > chondro <- chondro [- out] 69 I / a.u ν ~ cm 1 (a) The suspected outlier spectra. (b) pairs plot of the first 7 scores. Hierarchical Cluster Analysis (HCA) HCA merges objects according to their (dis)similarity into clusters. The result is a dendrogram, a graph stating at which level two objects are similar and thus grouped together. The first step in HCA is the choice of the distance. The R function dist offers a variety of distance measures to be computed. The so-called Pearson distance D 2 Pear son = 1 COR(X ) 2 is popular in data analysis of vibrational spectra and is provided by hyperspec s pearson.dist function. Also for computing the dendrogram, a number of choices are available. Here we choose Ward s method, and, as it uses Euclidean distance for calculating the dendrogram, Euclidean distance also for the distance matrix : 7
8 > dist <- dist (chondro) > dendrogram <- hclust (dist, method = "ward.d") > plot (dendrogram) In order to get clusters, the dendrogram is cut at a level specified either by height or by the number of clusters. > chondro$clusters <- as.factor (cutree (dendrogram, k = 3)) > cols <- c ("dark blue", "orange", "#C22") The result for k= 3 clusters is plotted as a map (figure b). If the color-coded variate (left hand side of the formula) is a factor, the legend bar does not show intermediate colors, and hyperspec s levelplot method uses the levels of the factor for the legend. Thus meaningful names are assigned > levels (chondro$clusters) <- c ("matrix", "lacuna", "cell") and the cluster membership map is plotted: > print (plotmap (chondro, clusters ~ x * y, col.regions = cols)) The cluster membership can also be marked in the dendrogram: > par (xpd = TRUE) # allow plotting the markers into the margin > plot (dendrogram, labels = FALSE, hang = -1) > mark.dendrogram (dendrogram, chondro$clusters, col = cols) Figure a shows the dendrogram and b the resulting cluster map. The three clusters correspond to the cartilage matrix, the lacuna and the cells. The left cell is destroyed and its contents are leaking into the matrix, while the right cells looks intact. We can plot the cluster mean spectra ± 1 standard deviation using aggregate (see figure 6): > cluster.means <- aggregate (chondro, chondro$clusters, mean_pm_sd) > plot (cluster.means, stacked = ".aggregate", fill = ".aggregate", col = cols) Height Cluster Dendrogram y 1 cell lacuna matrix matrix cell lacuna (a) dist The dendrogram. Figure 2 x (b) The cluster map for k = 3 clusters. Hierarchical cluster analysis. 6 Plotting a False-Colour Map of Certain Spectral Regions hyperspec comes with a sophisticated interface for specifying spectral ranges. Expressing things like cm -1 ± 1 data points is easily possible. Thus, we can have a fast look at the nucleic acid distribution, using the DNA bands at 728, 782, 98, 124, 1482, and 177 cm -1 : 8
9 I / a.u. matrix lacuna cell ν ~ cm 1 Figure 6 The cluster mean ± 1 standard deviation spectra. The blue cluster shows distinct lipid bands, the yellow cluster collagen, and the red cluster proteins and nucleic acids. > plotmap (chondro[,, c( 728, 782, 98, 124, 1482, 177)], + col.regions = colorramppalette (c("white", "gold", "dark green"), space = "Lab") (2)) y.. 2 x Figure 7 False colour map of the DNA band intensities. The result is shown in figure 7. While the nucleus of the right cell shows up nicely, nothing is detected in the remainders of the left cell. 7 Saving the data set Finally, the example data set is put together and saved. In order to keep the package size small, only a PCA-compressed version with PCs is shipped as example data set of the package. > spectra.to.save$clusters <- factor (NA, levels = levels (chondro$clusters)) > spectra.to.save$clusters[- out] <- chondro$clusters > pca <- prcomp (spectra.to.save) >.chondro.scores <- pca$x [, seq_len ()] >.chondro.loadings <- pca$rot [, seq_len ()] >.chondro.center <- pca$center >.chondro.wl <- wl (chondro) >.chondro.labels <- lapply (labels (chondro), as.expression) >.chondro.extra <- spectra.to.save$.. > save (.chondro.scores,.chondro.loadings,.chondro.center, 9
10 +.chondro.wl,.chondro.labels,.chondro.extra, + file = "chondro-internal.rda") This is the file distributed with hyperspec as example data set. Session Info R version ( ) Platform: x86_64-pc-linux-gnu (64-bit) Running under: Ubuntu LTS Matrix products: default BLAS/LAPACK: /usr/lib/openblas/lib/libopenblas_haswellp-r.2.19.so locale: [1] LC_CTYPE=de_DE.UTF-8 LC_NUMERIC=C LC_TIME=de_DE.UTF-8 [4] LC_COLLATE=de_DE.UTF-8 LC_MONETARY=de_DE.UTF-8 LC_MESSAGES=de_DE.UTF-8 [7] LC_PAPER=de_DE.UTF-8 LC_NAME=C LC_ADDRESS=C [] LC_TELEPHONE=C LC_MEASUREMENT=de_DE.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] grid stats graphics grdevices utils datasets methods base other attached packages: [1] latticeextra_.6-28 RColorBrewer_1.1-2 hyperspec_ ggplot2_2.2.1 [] lattice_.2-3 loaded via a namespace (and not attached): [1] Rcpp_ crayon_1.3.4 plyr_1.8.4 R6_2.2.2 gtable_.2. [6] magrittr_1. scales_.. rlang_.1.2 lazyeval_.2. rstudioapi_.7 [11] testthat_1..2 tools_3.4.2 munsell_.4.3 compiler_3.4.2 colorspace_1.3-2 [16] tibble_1.3.4
Raman Spectra of Chondrocytes in Cartilage: hyperspec s chondro data set
Raman Spectra of Chondrocytes in Cartilage: hyperspec s chondro data set Claudia Beleites CENMAT and DI3, University of Trieste Spectroscopy Imaging, IPHT Jena e.v. February 13,
More informationFitting Baselines to Spectra
Fitting Baselines to Spectra Claudia Beleites DIA Raman Spectroscopy Group, University of Trieste/Italy (2005 2008) Spectroscopy Imaging, IPHT, Jena/Germany (2008 2016) ÖPV, JKI,
More informationChemometric Analysis of Bio-Spectroscopic Data in : hyperspec
Chemometric Analysis of Bio-Spectroscopic Data in : hyperspec Claudia Beleites 1,2 (Claudia.Beleites@ipht-jena.de), Christoph Krafft 2, Jürgen Popp 2,3, and Valter Sergo 1 1 CENMAT and Dept. of Industrial
More informationCalibration of Quinine Fluorescence Emission Vignette for the Data Set flu of the R package hyperspec
Calibration of Quinine Fluorescence Emission Vignette for the Data Set flu of the R package hyperspec Claudia Beleites CENMAT and DI3, University of Trieste Spectroscopy Imaging,
More informationhyperspec Introduction
hyperspec Introduction Claudia Beleites DIA Raman Spectroscopy Group, University of Trieste/Italy (2005 2008) Spectroscopy Imaging, IPHT, Jena/Germany (2008 2016) ÖPV, JKI, Berlin/Germany
More informationExamples of implementation of pre-processing method described in paper with R code snippets - Electronic Supplementary Information (ESI)
Electronic Supplementary Material (ESI) for Analyst. This journal is The Royal Society of Chemistry 2015 Examples of implementation of pre-processing method described in paper with R code snippets - Electronic
More informationImage Analysis with beadarray
Mike Smith October 30, 2017 Introduction From version 2.0 beadarray provides more flexibility in the processing of array images and the extraction of bead intensities than its predecessor. In the past
More informationhyperspec Plotting functions
hperspec Plotting functions Claudia Beleites DIA Raman Spectroscop Group, Universit of Trieste/Ital (2 28) Spectroscop Imaging, IPHT, Jena/German (28 216) ÖPV, JKI, Berlin/German
More informationThe analysis of rtpcr data
The analysis of rtpcr data Jitao David Zhang, Markus Ruschhaupt October 30, 2017 With the help of this document, an analysis of rtpcr data can be performed. For this, the user has to specify several parameters
More informationRMassBank for XCMS. Erik Müller. January 4, Introduction 2. 2 Input files LC/MS data Additional Workflow-Methods 2
RMassBank for XCMS Erik Müller January 4, 2019 Contents 1 Introduction 2 2 Input files 2 2.1 LC/MS data........................... 2 3 Additional Workflow-Methods 2 3.1 Options..............................
More informationsegmentseq: methods for detecting methylation loci and differential methylation
segmentseq: methods for detecting methylation loci and differential methylation Thomas J. Hardcastle October 30, 2018 1 Introduction This vignette introduces analysis methods for data from high-throughput
More informationIntroduction to the TSSi package: Identification of Transcription Start Sites
Introduction to the TSSi package: Identification of Transcription Start Sites Julian Gehring, Clemens Kreutz October 30, 2017 Along with the advances in high-throughput sequencing, the detection of transcription
More informationAnalysing Spatial Data in R: Vizualising Spatial Data
Analysing Spatial Data in R: Vizualising Spatial Data Roger Bivand Department of Economics Norwegian School of Economics and Business Administration Bergen, Norway 31 August 2007 Vizualising Spatial Data
More informationThe latest trend of hybrid instrumentation
Multivariate Data Processing of Spectral Images: The Ugly, the Bad, and the True The results of various multivariate data-processing methods of Raman maps recorded with a dispersive Raman microscope are
More informationCQN (Conditional Quantile Normalization)
CQN (Conditional Quantile Normalization) Kasper Daniel Hansen khansen@jhsph.edu Zhijin Wu zhijin_wu@brown.edu Modified: August 8, 2012. Compiled: April 30, 2018 Introduction This package contains the CQN
More informationVisualisation, transformations and arithmetic operations for grouped genomic intervals
## Warning: replacing previous import ggplot2::position by BiocGenerics::Position when loading soggi Visualisation, transformations and arithmetic operations for grouped genomic intervals Thomas Carroll
More informationPackage superheat. February 4, 2017
Type Package Package superheat February 4, 2017 Title A Graphical Tool for Exploring Complex Datasets Using Heatmaps Version 0.1.0 Description A system for generating extendable and customizable heatmaps
More informationsurvsnp: Power and Sample Size Calculations for SNP Association Studies with Censored Time to Event Outcomes
survsnp: Power and Sample Size Calculations for SNP Association Studies with Censored Time to Event Outcomes Kouros Owzar Zhiguo Li Nancy Cox Sin-Ho Jung Chanhee Yi June 29, 2016 1 Introduction This vignette
More informationAdvanced analysis using bayseq; generic distribution definitions
Advanced analysis using bayseq; generic distribution definitions Thomas J Hardcastle October 30, 2017 1 Generic Prior Distributions bayseq now offers complete user-specification of underlying distributions
More informationPreliminary Figures for Renormalizing Illumina SNP Cell Line Data
Preliminary Figures for Renormalizing Illumina SNP Cell Line Data Kevin R. Coombes 17 March 2011 Contents 1 Executive Summary 1 1.1 Introduction......................................... 1 1.1.1 Aims/Objectives..................................
More informationsegmentseq: methods for detecting methylation loci and differential methylation
segmentseq: methods for detecting methylation loci and differential methylation Thomas J. Hardcastle October 13, 2015 1 Introduction This vignette introduces analysis methods for data from high-throughput
More informationTriform: peak finding in ChIP-Seq enrichment profiles for transcription factors
Triform: peak finding in ChIP-Seq enrichment profiles for transcription factors Karl Kornacker * and Tony Håndstad October 30, 2018 A guide for using the Triform algorithm to predict transcription factor
More informationPackage hyperspec. October 6, 2017
Encoding UTF-8 Type Package Package hyperspec October 6, 2017 Title Work with Hyperspectral Data, i.e. Spectra + Meta Information (Spatial, Time, Concentration,...) Version 0.99-20171005 Date 2017-10-05
More informationMultivariate Calibration Quick Guide
Last Updated: 06.06.2007 Table Of Contents 1. HOW TO CREATE CALIBRATION MODELS...1 1.1. Introduction into Multivariate Calibration Modelling... 1 1.1.1. Preparing Data... 1 1.2. Step 1: Calibration Wizard
More informationExperiment 5: Exploring Resolution, Signal, and Noise using an FTIR CH3400: Instrumental Analysis, Plymouth State University, Fall 2013
Experiment 5: Exploring Resolution, Signal, and Noise using an FTIR CH3400: Instrumental Analysis, Plymouth State University, Fall 2013 Adapted from JP Blitz and DG Klarup, "Signal-to-Noise Ratio, Signal
More informationAnalyse RT PCR data with ddct
Analyse RT PCR data with ddct Jitao David Zhang, Rudolf Biczok and Markus Ruschhaupt October 30, 2017 Abstract Quantitative real time PCR (qrt PCR or RT PCR for short) is a laboratory technique based on
More informationAn Overview of the S4Vectors package
Patrick Aboyoun, Michael Lawrence, Hervé Pagès Edited: February 2018; Compiled: June 7, 2018 Contents 1 Introduction.............................. 1 2 Vector-like and list-like objects...................
More informationCardinal design and development
Kylie A. Bemis October 30, 2017 Contents 1 Introduction.............................. 2 2 Design overview........................... 2 3 iset: high-throughput imaging experiments............ 3 3.1 SImageSet:
More informationHow to use CNTools. Overview. Algorithms. Jianhua Zhang. April 14, 2011
How to use CNTools Jianhua Zhang April 14, 2011 Overview Studies have shown that genomic alterations measured as DNA copy number variations invariably occur across chromosomal regions that span over several
More informationIntroduction to Raman spectroscopy measurement data processing using Igor Pro
Introduction to Raman spectroscopy measurement data processing using Igor Pro This introduction is intended to minimally guide beginners to processing Raman spectroscopy measurement data, which includes
More informationPackage arphit. March 28, 2019
Type Package Title RBA-style R Plots Version 0.3.1 Author Angus Moore Package arphit March 28, 2019 Maintainer Angus Moore Easily create RBA-style graphs
More informationA complete analysis of peptide microarray binding data using the pepstat framework
A complete analysis of peptide microarray binding data using the pepstat framework Greg Imholte, Renan Sauteraud, Mike Jiang and Raphael Gottardo April 30, 2018 This document present a full analysis, from
More informationAnalysis of two-way cell-based assays
Analysis of two-way cell-based assays Lígia Brás, Michael Boutros and Wolfgang Huber April 16, 2015 Contents 1 Introduction 1 2 Assembling the data 2 2.1 Reading the raw intensity files..................
More informationMIGSA: Getting pbcmc datasets
MIGSA: Getting pbcmc datasets Juan C Rodriguez Universidad Católica de Córdoba Universidad Nacional de Córdoba Cristóbal Fresno Instituto Nacional de Medicina Genómica Andrea S Llera Fundación Instituto
More informationMMDiff2: Statistical Testing for ChIP-Seq Data Sets
MMDiff2: Statistical Testing for ChIP-Seq Data Sets Gabriele Schwikert 1 and David Kuo 2 [1em] 1 The Informatics Forum, University of Edinburgh and The Wellcome
More informationSCnorm: robust normalization of single-cell RNA-seq data
SCnorm: robust normalization of single-cell RNA-seq data Rhonda Bacher and Christina Kendziorski February 18, 2019 Contents 1 Introduction.............................. 1 2 Run SCnorm.............................
More informationNacoStringQCPro. Dorothee Nickles, Thomas Sandmann, Robert Ziman, Richard Bourgon. Modified: April 10, Compiled: April 24, 2017.
NacoStringQCPro Dorothee Nickles, Thomas Sandmann, Robert Ziman, Richard Bourgon Modified: April 10, 2017. Compiled: April 24, 2017 Contents 1 Scope 1 2 Quick start 1 3 RccSet loaded with NanoString data
More informationriboseqr Introduction Getting Data Workflow Example Thomas J. Hardcastle, Betty Y.W. Chung October 30, 2017
riboseqr Thomas J. Hardcastle, Betty Y.W. Chung October 30, 2017 Introduction Ribosome profiling extracts those parts of a coding sequence currently bound by a ribosome (and thus, are likely to be undergoing
More informationIntroduction to the Codelink package
Introduction to the Codelink package Diego Diez October 30, 2018 1 Introduction This package implements methods to facilitate the preprocessing and analysis of Codelink microarrays. Codelink is a microarray
More informationAnalyzing Genomic Data with NOJAH
Analyzing Genomic Data with NOJAH TAB A) GENOME WIDE ANALYSIS Step 1: Select the example dataset or upload your own. Two example datasets are available. Genome-Wide TCGA-BRCA Expression datasets and CoMMpass
More informationMultipleAlignment Objects
MultipleAlignment Objects Marc Carlson Bioconductor Core Team Fred Hutchinson Cancer Research Center Seattle, WA April 30, 2018 Contents 1 Introduction 1 2 Creation and masking 1 3 Analytic utilities 7
More informationModule 10. Data Visualization. Andrew Jaffe Instructor
Module 10 Data Visualization Andrew Jaffe Instructor Basic Plots We covered some basic plots on Wednesday, but we are going to expand the ability to customize these basic graphics first. 2/37 But first...
More informationRaman Images. Jeremy M. Shaver 1, Eunah Lee 2, Andrew Whitley 2, R. Scott Koch. 1. Eigenvector Research, Inc. 2. HORIBA Jobin Yvon, Inc.
Analyzing and Visualizing Large Raman Images Jeremy M. Shaver 1, Eunah Lee 2, Andrew Whitley 2, R. Scott Koch 1 1. Eigenvector Research, Inc. 2. HORIBA Jobin Yvon, Inc. What is a Large Image? Gone from
More informationHow to use cghmcr. October 30, 2017
How to use cghmcr Jianhua Zhang Bin Feng October 30, 2017 1 Overview Copy number data (arraycgh or SNP) can be used to identify genomic regions (Regions Of Interest or ROI) showing gains or losses that
More informationHierarchical Clustering
Hierarchical Clustering Hierarchical Clustering Produces a set of nested clusters organized as a hierarchical tree Can be visualized as a dendrogram A tree-like diagram that records the sequences of merges
More informationHow to Use pkgdeptools
How to Use pkgdeptools Seth Falcon February 10, 2018 1 Introduction The pkgdeptools package provides tools for computing and analyzing dependency relationships among R packages. With it, you can build
More informationNavigator. The Navigator
TM002-02-A Introduction to WiRE and System start-up Navigator The Navigator Enables control of what data is open and where it is viewed Windows Measurement - Viewers and data housed together (different
More informationPoints Lines Connected points X-Y Scatter. X-Y Matrix Star Plot Histogram Box Plot. Bar Group Bar Stacked H-Bar Grouped H-Bar Stacked
Plotting Menu: QCExpert Plotting Module graphs offers various tools for visualization of uni- and multivariate data. Settings and options in different types of graphs allow for modifications and customizations
More informationHierarchical Clustering
What is clustering Partitioning of a data set into subsets. A cluster is a group of relatively homogeneous cases or observations Hierarchical Clustering Mikhail Dozmorov Fall 2016 2/61 What is clustering
More informationComputing with large data sets
Computing with large data sets Richard Bonneau, spring 2009 Lecture 8(week 5): clustering 1 clustering Clustering: a diverse methods for discovering groupings in unlabeled data Because these methods don
More informationUsing metama for differential gene expression analysis from multiple studies
Using metama for differential gene expression analysis from multiple studies Guillemette Marot and Rémi Bruyère Modified: January 28, 2015. Compiled: January 28, 2015 Abstract This vignette illustrates
More informationGetting started with flowstats
Getting started with flowstats F. Hahne, N. Gopalakrishnan October, 7 Abstract flowstats is a collection of algorithms for the statistical analysis of flow cytometry data. So far, the focus is on automated
More informationReproducible Homerange Analysis
Reproducible Homerange Analysis (Sat Aug 09 15:28:43 2014) based on the rhr package This is an automatically generated file with all parameters and settings, in order to enable later replication of the
More informationCount outlier detection using Cook s distance
Count outlier detection using Cook s distance Michael Love August 9, 2014 1 Run DE analysis with and without outlier removal The following vignette produces the Supplemental Figure of the effect of replacing
More informationIntroduction to R (BaRC Hot Topics)
Introduction to R (BaRC Hot Topics) George Bell September 30, 2011 This document accompanies the slides from BaRC s Introduction to R and shows the use of some simple commands. See the accompanying slides
More informationMethodology for spot quality evaluation
Methodology for spot quality evaluation Semi-automatic pipeline in MAIA The general workflow of the semi-automatic pipeline analysis in MAIA is shown in Figure 1A, Manuscript. In Block 1 raw data, i.e..tif
More informationDupChecker: a bioconductor package for checking high-throughput genomic data redundancy in meta-analysis
DupChecker: a bioconductor package for checking high-throughput genomic data redundancy in meta-analysis Quanhu Sheng, Yu Shyr, Xi Chen Center for Quantitative Sciences, Vanderbilt University, Nashville,
More informationQuick Start Guide: Welcome to OceanView
Quick Start Guide: Welcome to OceanView Contents: Ctrl + Click a Topic Below for More Information Introduction to OceanView... 3 Welcome Screen... 3 Data Window... 3 Schematic View... 3 Persistence...
More informationComputer Graphics Fundamentals. Jon Macey
Computer Graphics Fundamentals Jon Macey jmacey@bournemouth.ac.uk http://nccastaff.bournemouth.ac.uk/jmacey/ 1 1 What is CG Fundamentals Looking at how Images (and Animations) are actually produced in
More informationChapter 6: Cluster Analysis
Chapter 6: Cluster Analysis The major goal of cluster analysis is to separate individual observations, or items, into groups, or clusters, on the basis of the values for the q variables measured on each
More informationUsing the qrqc package to gather information about sequence qualities
Using the qrqc package to gather information about sequence qualities Vince Buffalo Bioinformatics Core UC Davis Genome Center vsbuffalo@ucdavis.edu 2012-02-19 Abstract Many projects in bioinformatics
More informationPROPER: PROspective Power Evaluation for RNAseq
PROPER: PROspective Power Evaluation for RNAseq Hao Wu [1em]Department of Biostatistics and Bioinformatics Emory University Atlanta, GA 303022 [1em] hao.wu@emory.edu October 30, 2017 Contents 1 Introduction..............................
More informationClustering algorithms 6CCS3WSN-7CCSMWAL
Clustering algorithms 6CCS3WSN-7CCSMWAL Contents Introduction: Types of clustering Hierarchical clustering Spatial clustering (k means etc) Community detection (next week) What are we trying to cluster
More informationCluster Analysis for Microarray Data
Cluster Analysis for Microarray Data Seventh International Long Oligonucleotide Microarray Workshop Tucson, Arizona January 7-12, 2007 Dan Nettleton IOWA STATE UNIVERSITY 1 Clustering Group objects that
More informationR-software multims-toolbox. (User Guide)
R-software multims-toolbox (User Guide) Pavel Cejnar 1, Štěpánka Kučková 2 1 Department of Computing and Control Engineering, Institute of Chemical Technology, Technická 3, 166 28 Prague 6, Czech Republic,
More informationTutorial 7: Automated Peak Picking in Skyline
Tutorial 7: Automated Peak Picking in Skyline Skyline now supports the ability to create custom advanced peak picking and scoring models for both selected reaction monitoring (SRM) and data-independent
More informationLD vignette Measures of linkage disequilibrium
LD vignette Measures of linkage disequilibrium David Clayton June 13, 2018 Calculating linkage disequilibrium statistics We shall first load some illustrative data. > data(ld.example) The data are drawn
More informationBar Charts and Frequency Distributions
Bar Charts and Frequency Distributions Use to display the distribution of categorical (nominal or ordinal) variables. For the continuous (numeric) variables, see the page Histograms, Descriptive Stats
More informationPackage ClustGeo. R topics documented: July 14, Type Package
Type Package Package ClustGeo July 14, 2017 Title Hierarchical Clustering with Spatial Constraints Version 2.0 Author Marie Chavent [aut, cre], Vanessa Kuentz [aut], Amaury Labenne [aut], Jerome Saracco
More informationTutorial script for whole-cell MALDI-TOF analysis
Tutorial script for whole-cell MALDI-TOF analysis Julien Textoris June 19, 2013 Contents 1 Required libraries 2 2 Data loading 2 3 Spectrum visualization and pre-processing 4 4 Analysis and comparison
More informationClustering and Dimensionality Reduction
Clustering and Dimensionality Reduction Some material on these is slides borrowed from Andrew Moore's excellent machine learning tutorials located at: Data Mining Automatically extracting meaning from
More informationNature Methods: doi: /nmeth Supplementary Figure 1
Supplementary Figure 1 Schematic representation of the Workflow window in Perseus All data matrices uploaded in the running session of Perseus and all processing steps are displayed in the order of execution.
More informationSimBindProfiles: Similar Binding Profiles, identifies common and unique regions in array genome tiling array data
SimBindProfiles: Similar Binding Profiles, identifies common and unique regions in array genome tiling array data Bettina Fischer October 30, 2018 Contents 1 Introduction 1 2 Reading data and normalisation
More informationHyperspectral Chemical Imaging: principles and Chemometrics.
Hyperspectral Chemical Imaging: principles and Chemometrics aoife.gowen@ucd.ie University College Dublin University College Dublin 1,596 PhD students 6,17 international students 8,54 graduate students
More informationflowchic - Analyze flow cytometric data of complex microbial communities based on histogram images
flowchic - Analyze flow cytometric data of complex microbial communities based on histogram images Schumann, J., Koch, C., Fetzer, I., Müller, S. Helmholtz Centre for Environmental Research - UFZ Department
More informationSpot Analyzer TiVi80 User Manual User Manual 3.2 Version 3.2 October 2013
Spot Analyzer TiVi80 User Manual User Manual 3.2 Version 3.2 October 2013 PIONEERS IN TISSUE VIABILITY IMAGING Dear Valued Customer! Welcome to the WheelsBridge Spot Analyzer TiVi80 system designed for
More information3-D. Here red spheres show the location of gold nanoparticles inside/around a cell nucleus.
3-D The CytoViva 3-D System allows the user can locate objects of interest in a 3-D space. It does this by acquiring multiple Z planes and performing our custom software routines to locate and observe
More informationUNIFIT - Spectrum Processing, Peak Fitting, Analysis and Presentation Software for XPS, AES, XAS and RAMAN Spectroscopy Based on WINDOWS
UNIFIT - Spectrum Processing, Peak Fitting, Analysis and Presentation Software for XPS, AES, XAS and RAMAN Spectroscopy Based on WINDOWS UNIFIT FOR WINDOWS is an universal processing, analysis and presentation
More informationMany Patterns & Many Methods
Many Patterns & Many Methods New methods for visualising & utilising multiple analysis techniques in polymorph and salt screening systems 2009 Gordon Barr, Chris Gilmore & Gordon Cunningham WestCHEM, Chemistry
More informationMultivariate analyses in ecology. Cluster (part 2) Ordination (part 1 & 2)
Multivariate analyses in ecology Cluster (part 2) Ordination (part 1 & 2) 1 Exercise 9B - solut 2 Exercise 9B - solut 3 Exercise 9B - solut 4 Exercise 9B - solut 5 Multivariate analyses in ecology Cluster
More informationMetabolomic Data Analysis with MetaboAnalyst
Metabolomic Data Analysis with MetaboAnalyst User ID: guest6522519400069885256 April 14, 2009 1 Data Processing and Normalization 1.1 Reading and Processing the Raw Data MetaboAnalyst accepts a variety
More informationHESP PIPELINE v. 1.0
HESP PIPELINE v. 1.0 Installation and Usage Arun Surya 20/06/2017 1. INSTALLATION The HESP pipeline installation zip file can be downloaded from, https://www.iiap.res.in/hesp/hesp_pipeline.zip. The zip
More informationPresentation and analysis of multidimensional data sets
Presentation and analysis of multidimensional data sets Overview 1. 3D data visualisation Multidimensional images Data pre-processing Visualisation methods Multidimensional images wavelength time 3D image
More informationData Visualization. Andrew Jaffe Instructor
Module 9 Data Visualization Andrew Jaffe Instructor Basic Plots We covered some basic plots previously, but we are going to expand the ability to customize these basic graphics first. 2/45 Read in Data
More informationClassification of Hyperspectral Breast Images for Cancer Detection. Sander Parawira December 4, 2009
1 Introduction Classification of Hyperspectral Breast Images for Cancer Detection Sander Parawira December 4, 2009 parawira@stanford.edu In 2009 approximately one out of eight women has breast cancer.
More informationMALDIquant: Quantitative Analysis of Mass Spectrometry Data
MALDIquant: Quantitative Analysis of Mass Spectrometry Data Sebastian Gibb November 12, 2017 Abstract MALDIquant provides a complete analysis pipeline for MALDI- TOF and other 2D mass spectrometry data.
More informationRNAseq Differential Expression Analysis Jana Schor and Jörg Hackermüller November, 2017
RNAseq Differential Expression Analysis Jana Schor and Jörg Hackermüller November, 2017 RNA-Seq: Example dataset and download} Dataset taken from Somel et al. 2012: (check publication for details) 3 human
More informationPackage rafalib. R topics documented: August 29, Version 1.0.0
Version 1.0.0 Package rafalib August 29, 2016 Title Convenience Functions for Routine Data Eploration A series of shortcuts for routine tasks originally developed by Rafael A. Irizarry to facilitate data
More informationChemometrics. Description of Pirouette Algorithms. Technical Note. Abstract
19-1214 Chemometrics Technical Note Description of Pirouette Algorithms Abstract This discussion introduces the three analysis realms available in Pirouette and briefly describes each of the algorithms
More informationUser Manual STATISTICS
Version 6 User Manual STATISTICS 2006 BRUKER OPTIK GmbH, Rudolf-Plank-Straße 27, D-76275 Ettlingen, www.brukeroptics.com All rights reserved. No part of this manual may be reproduced or transmitted in
More informationPackage JBTools. R topics documented: June 2, 2015
Package JBTools June 2, 2015 Title Misc Small Tools and Helper Functions for Other Code of J. Buttlar Version 0.7.2.9 Date 2015-05-20 Author Maintainer Collection of several
More informationRaman Spectrometer Installation Manual
RI Raman Spectrometer Instruction Manual Application software (included in standard configuration) Connection RI RAMAN is recognized by Windows operational system as standard HID device. That is why there
More informationUser Manual Release Note. PhaserMatch
User Manual Release Note PhaserMatch Copyright Tektronix, Inc. Printed in the U.S.A. Unpublished rights reserved under the copyright laws of the United States. Contents of this publication may not be reproduced
More informationVery large searches present a number of challenges. These are the topics we will cover during this presentation.
1 Very large searches present a number of challenges. These are the topics we will cover during this presentation. 2 The smartest way to merge files, like fractions from a MudPIT run, is using Mascot Daemon.
More informationSCnorm: robust normalization of single-cell RNA-seq data
SCnorm: robust normalization of single-cell RNA-seq data Rhonda Bacher and Christina Kendziorski September 24, 207 Contents Introduction.............................. 2 Run SCnorm.............................
More informationIntroduction to BatchtoolsParam
Nitesh Turaga 1, Martin Morgan 2 Edited: March 22, 2018; Compiled: January 4, 2019 1 Nitesh.Turaga@ RoswellPark.org 2 Martin.Morgan@ RoswellPark.org Contents 1 Introduction..............................
More informationOperating Procedure for Horiba Raman Microscope
Operating Procedure for Horiba Raman Microscope SAFETY Be aware of Laser radiation at all times! Do not remove the covers of the instrument. Components are supplied with 110V electric source. Do not touch
More informationPackage SC3. September 29, 2018
Type Package Title Single-Cell Consensus Clustering Version 1.8.0 Author Vladimir Kiselev Package SC3 September 29, 2018 Maintainer Vladimir Kiselev A tool for unsupervised
More informationShrinkage of logarithmic fold changes
Shrinkage of logarithmic fold changes Michael Love August 9, 2014 1 Comparing the posterior distribution for two genes First, we run a DE analysis on the Bottomly et al. dataset, once with shrunken LFCs
More information1. ABOUT INSTALLATION COMPATIBILITY SURESIM WORKFLOWS a. Workflow b. Workflow SURESIM TUTORIAL...
SuReSim manual 1. ABOUT... 2 2. INSTALLATION... 2 3. COMPATIBILITY... 2 4. SURESIM WORKFLOWS... 2 a. Workflow 1... 3 b. Workflow 2... 4 5. SURESIM TUTORIAL... 5 a. Import Data... 5 b. Parameter Selection...
More information