R/pplot.R

In lgatto/MLInterfaces: Uniform interfaces to R machine learning procedures for data in Bioconductor containers

#
# code to provide planar displays of classifier predictions
#
# currently assumes use of two raw expression vectors to define
# the plane
#

setGeneric("getGrid", function(x) standardGeneric("getGrid"))

setMethod("getGrid", "ExpressionSet", function(x) {
    lims <- apply(exprs(x), 1, range)
    if (ncol(lims) > 2) 
        warning("grid based only on first 2 genes")
    nn <- colnames(lims)
    v1 <- seq(lims[1, 1], lims[2, 1], len = 200)
    v2 <- seq(lims[1, 2], lims[2, 2], len = 200)
    og <- expand.grid(v1, v2)
    names(og) <- make.names(nn)
    og
})

setMethod("getGrid", "data.frame", function(x) {
    lims <- apply(x, 2, range)
    if (ncol(lims) > 2) 
        warning("grid based only on first 2 genes")
    nn <- colnames(lims)
    v1 <- seq(lims[1, 1], lims[2, 1], len = 200)
    v2 <- seq(lims[1, 2], lims[2, 2], len = 200)
    og <- expand.grid(v1, v2)
    names(og) <- make.names(nn)
    og
})



setGeneric("planarPlot", function(clo, eset, classifLab ) standardGeneric("planarPlot"))

setMethod("planarPlot", c("classifierOutput", "ExpressionSet", "character"), 
	function(clo, eset, classifLab) {
 requireNamespace("RColorBrewer")
 pal <- brewer.pal("Set2", n=8)
 ff <- getGrid(eset)
 if (clo@learnerSchema@mlFunName %in% c("nnet", "rpart"))
    ps <- predict( RObject(clo), newdata=ff, type="class")
 else if (clo@learnerSchema@mlFunName == "randomForest")
    {
    names(ff) <- rownames(exprs(eset))
    ps <- predict( RObject(clo), newdata=ff )
    }
 else
    ps <- predict( RObject(clo), newdata=ff )
 if (clo@learnerSchema@mlFunName %in% c("lda", "qda")) ps <- ps[[1]]
 plot( ff[,1], ff[,2], col=pal[as.numeric(factor(ps))], pch=19,
 xlab = names(ff)[1], ylab=names(ff)[2] )
 #text( exprs(eset)[1,], exprs(eset)[2,], lab=pData(eset)[[classifLab]])
})

#setMethod("planarPlot", c("classifOutput", "data.frame", "character"), 
#	function(clo, eset, classifLab) {
# require(RColorBrewer)
# pal <- brewer.pal("Set2", n=8)
# ff <- getGrid(eset)
# if (clo@method %in% c("nnet", "rpart"))
#    ps <- predict( clo@RObject, newdata=ff, type="class")
# else if (clo@method == "randomForest")
#    {
#    names(ff) <- rownames(exprs(eset))
#    ps <- predict( clo@RObject, newdata=ff )
#    }
# else
#    ps <- predict( clo@RObject, newdata=ff )
# if (clo@method %in% c("lda", "qda")) ps <- ps[[1]]
# plot( ff[,1], ff[,2], col=pal[as.numeric(factor(ps))], pch=19,
# xlab = names(ff)[1], ylab=names(ff)[2] )
# #text( exprs(eset)[1,], exprs(eset)[2,], lab=pData(eset)[[classifLab]])
#})