1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
|
\name{normalize.quantiles.in.blocks}
\alias{normalize.quantiles.in.blocks}
\title{Quantile Normalization carried out separately within blocks of rows}
\description{
Using a normalization based upon quantiles this function
normalizes the columns of a matrix such that different subsets of rows
get normalized together.
}
\usage{
normalize.quantiles.in.blocks(x,blocks,copy=TRUE)
}
\arguments{
\item{x}{A matrix of intensities where each column corresponds to a
chip and each row is a probe.}
\item{copy}{Make a copy of matrix before normalizing. Usually safer to
work with a copy}
\item{blocks}{A vector giving block membership for each each row}
}
\details{This method is based upon the concept of a quantile-quantile
plot extended to n dimensions. No special allowances are made for
outliers. If you make use of quantile normalization either through
\code{\link[affy]{rma}} or \code{\link[affy]{expresso}}
please cite Bolstad et al, Bioinformatics (2003).
}
\value{
From \code{normalize.quantiles.use.target} a normalized \code{matrix}.
}
\references{
Bolstad, B (2001) \emph{Probe Level Quantile Normalization of High Density
Oligonucleotide Array Data}. Unpublished manuscript
\url{http://bmbolstad.com/stuff/qnorm.pdf}
Bolstad, B. M., Irizarry R. A., Astrand, M, and Speed, T. P. (2003)
\emph{A Comparison of Normalization Methods for High Density
Oligonucleotide Array Data Based on Bias and Variance.}
Bioinformatics 19(2) ,pp 185-193. \url{http://bmbolstad.com/misc/normalize/normalize.html}
}
\examples{
### setup the data
blocks <- c(rep(1,5),rep(2,5),rep(3,5))
par(mfrow=c(3,2))
x <- matrix(c(rexp(5,0.05),rnorm(5),rnorm(5,10)))
boxplot(x ~ blocks)
y <- matrix(c(-rexp(5,0.05),rnorm(5,10),rnorm(5)))
boxplot(y ~ blocks)
pre.norm <- cbind(x,y)
### the in.blocks version
post.norm <- normalize.quantiles.in.blocks(pre.norm,blocks)
boxplot(post.norm[,1] ~ blocks)
boxplot(post.norm[,2] ~ blocks)
### the usual version
post.norm <- normalize.quantiles(pre.norm)
boxplot(post.norm[,1] ~ blocks)
boxplot(post.norm[,2] ~ blocks)
}
\author{Ben Bolstad, \email{bmb@bmbolstad.com}}
\seealso{\code{\link{normalize.quantiles}}}
\keyword{manip}
|
