Learn R Programming

charm (version 2.18.0)

qcReport: Microarray quality report

Description

Calculate microarray quality scores and produce an optional pdf report

Usage

qcReport(dat, file = NULL, utRange = c(30, 100), enRange = c(8, 12), numProbes = 5e+05, blockSize)

Arguments

dat
a TilingFeatureSet
file
name of output pdf file
utRange
color-scale range for the untreated channel plots
enRange
color-scale range for the methyl-depleted channel plots
numProbes
maximum number of probes to use for plots. If smaller than the number of probes on the array numProbes are chosen at random, speeding up calculations for high-density arrays with several million probes.
blockSize
The array is divided into a series of blockSize x blockSize rectangular blocks and the average signal level calculated for each. If blockSize is unspecified a size is chosen that gives about 1250 probes per block.

Value

a matrix with a row for each sample. The 3 columns contain array signal strength score, untreated channel standard deviation and methyl-depleted channel standard deviation.

Details

This function calculates microarray quality scores and produces an optional pdf report. Three quality metrics are calculated for each array:
Average signal strength.
The average percentile rank of untreated channel signal probes among the background (anti-genomic) probes. Since the untreated channel contains total DNA a successful hybridization would have strong signal for all untreated channel genomic probes.

Untreated channel signal standard deviation.
The array is divided into a series of rectangular blocks and the average signal level calculated for each. Since probes are arranged randomly on the array there should be no large differences between blocks. Arrays with spatial artifacts have a larger standard deviation between blocks.

Methyl-depleted channel signal standard deviation

Examples

Run this code
	if (require(charmData)) {
		phenodataDir <- system.file("extdata", package="charmData")
		pd <- read.delim(file.path(phenodataDir, "phenodata.txt"))
		dataDir <- system.file("data", package="charmData")
		setwd(dataDir)
		rawData <- readCharm(files=pd$filename, sampleKey=pd)
                ## Not run:
		#qcReport(rawData, file="qcReport.pdf")
	}

Run the code above in your browser using DataLab