Learn R Programming

hyperSpec (version 0.98-20140523)

chondro: Raman spectra of 2 Chondrocytes in Cartilage A Raman-map (laterally resolved Raman spectra) of chondrocytes in cartilage.

Description

See the vignette.

Arguments

format

The data set has 796 Raman spectra measured on a 23 $\times$ 35 grid with 1 micron step size. Spatial information is in chondro$x and chondro$y. Each spectrum has 300 data points in the range of ca. 600 - 1800 cm$^{-1}$.

References

The raw data is available at http://hyperspec.r-forge.r-project.org/blob/chondro.zip

Examples

Run this code
chondro

## do baseline correction
baselines <- spc.fit.poly.below (chondro)
chondro <- chondro - baselines

## area normalization
chondro <- chondro / colMeans (chondro)

## substact common composition
chondro <- chondro - quantile (chondro, 0.05)

## PCA
pca <- prcomp (~ spc, data = chondro, center = TRUE)
scores <- decomposition (chondro, pca$x, label.wavelength = "PC", label.spc = "score / a.u.")
loadings <- decomposition (chondro, pca$rotation, scores = FALSE, label.spc = "loading I / a.u.")
center <- decomposition (chondro, pca$center, scores = FALSE)

# remove outliers
out <- c(105, 140, 216, 289, 75, 69)
chondro <- chondro [- out]

# Hierarchical cluster analysis
# calculation omitted for speed - see vignette for details

cols <- c ("dark blue", "orange", "#C02020")
plotmap (chondro, clusters ~ x * y, col.regions = cols)

cluster.means <- aggregate (chondro, chondro$clusters, mean_pm_sd)
plot (cluster.means, stacked = ".aggregate", fill = ".aggregate", col = cols)

## plot nucleic acids
plotmap (chondro[, , c( 728, 782, 1098, 1240, 1482, 1577)],
       col.regions = colorRampPalette (c("white", "gold", "dark green"), space = "Lab") (20))

vignette ("chondro", package = "hyperSpec")

Run the code above in your browser using DataLab