## Load Glioblastoma multiforme data
data(GBM)
observations <- SMAPObservations(value=as.numeric(GBM[,2]),
chromosome=as.character(GBM[,3]),
startPosition=as.numeric(GBM[,4]),
endPosition=as.numeric(GBM[,5]),
name="G24460",
reporterId=as.character(GBM[,1]))
## plot observations
plot(observations, ylim=c(0,2))
## plot subset of observations (chromosome 9)
ids <- which(chromosome(observations) == "9")
plot(observations[ids])
Run the code above in your browser using DataLab