data(chromosome17)
## pSimCCA model on chromosome 17p
models17ppSimCCA <- screen.cgh.mrna(geneExp, geneCopyNum, 10, 17, 'p')
plot(models17ppSimCCA,
hilightGenes=c("ENSG00000108342", "ENSG00000108298"), showDensity = TRUE)
## Dependency model around 50th gene
model <- models17ppSimCCA[[50]]
## example annnotation types
ann.types <- factor(c(rep("Samples 1 - 10", 10), rep("Samples 11 - 51", 41)))
plot(model, geneExp, geneCopyNum, ann.types, legend.x = 40, legend.y = -4,
order = TRUE)
Run the code above in your browser using DataLab