R/bicluster.r

In biclust: BiCluster Algorithms

### Function to extract bicluster

bicluster <- function(x, BicRes, number=1:BicRes@Number)
{
    res <- list()
    for(i in 1:length(number))
    {
        res[[i]] <- x[BicRes@RowxNumber[,number[i]],BicRes@NumberxCol[number[i],]]
    }
    names(res) <- paste("Bicluster",number,sep="")
    return(res)
}

biclusternumber <- function(BicRes, number=1:BicRes@Number)
{
    res <- list()
    for(i in 1:length(number))
    {
        res[[i]] <- list()
        res[[i]][[1]] <- which(BicRes@RowxNumber[,number[i]])
        res[[i]][[2]] <- which(BicRes@NumberxCol[number[i],])
        names(res[[i]])<-c("Rows","Cols")
    }
    names(res) <- paste("Bicluster",number,sep="")
    return(res)
}

biclusterES <- function(bicres, number, plot=FALSE, ...)
{
    res <- bicres@info$Rowvalues[,number] %*% t(bicres@info$Colvalues[number,])
    if(plot==1)
    {
        #biclust:::heatmapBC(res, bicres, number=number, ...)
		heatmapBC(res, bicres, number=number, ...)
		
    }
    else
    {
        if(plot==2)
        {
			#biclust:::drawHeatmap(res, bicres, number=number, ...)
			drawHeatmap(res, bicres, number=number, ...)
        }
    }


    return(res)
}