# \dontshow{
data(Hedenfalk)
X <- Hedenfalk
k <- dim(X)[1]
s <- apply(X[1:100, 1:7], 1, density)
# }
# \donttest{
data(Hedenfalk)
X <- Hedenfalk
k <- dim(X)[1]
### Estimated densities of logged gene expression levels for patients with BRCA1 mutations.
s <- apply(X[, 1:7], 1, density)
### Plot of estimated densities of 6 randomly selected genes.
set.seed(375)
rs <- sample(1:k, 6)
plot(s[[rs[1]]], main = "Kernel estimates for 6 randomly selected genes",
xlab = "x", ylab = "density", xlim = c(-1, 4.3), ylim = c(0, 2))
for (i in 2:6){
lines(s[[rs[i]]], col = i)
}
# }
Run the code above in your browser using DataLab