1

I am trying to append the "matrix" class and in turn overwrite the default behaviour of "[". Code examples below:

annMatrix <- function(mat=NULL, rowAnn=NULL, colAnn=NULL) {
  if(is.null(mat)) mat <- matrix(nrow=0, ncol=0)
  mat <- as.matrix(mat)
  if(is.null(rowAnn)) rowAnn <- data.frame(row.names=seq_len(nrow(mat)))
  if(is.null(colAnn)) colAnn <- data.frame(row.names=seq_len(ncol(mat)))
  rowAnn <- data.frame(rowAnn, stringsAsFactors=FALSE)
  colAnn <- data.frame(colAnn, stringsAsFactors=FALSE)
  stopifnot(nrow(mat)==nrow(rowAnn) & ncol(mat)==nrow(colAnn))
  attr(mat, "colAnn") <- colAnn
  attr(mat, "rowAnn") <- rowAnn
  class(mat) <- append(class(mat), "annMatrix")
  mat
}

`[.annMatrix` <- function(annMat, rowExpr=NULL, colExpr=NULL) {
  stopifnot(is.valid.annMatrix(annMat))
  rowExpr <- eval(substitute(list(rowExpr)), attr(annMat, "rowAnn"), parent.frame())
  colExpr <- eval(substitute(list(colExpr)), attr(annMat, "colAnn"), parent.frame())
  indsR   <- unlist(rowExpr)
  indsC   <- unlist(colExpr)
  if(is.null(indsR)) indsR <- seq_len(nrow(annMat))
  if(is.null(indsC)) indsC <- seq_len(ncol(annMat))
  attr(annMat, "rowAnn") <- attr(annMat, "rowAnn")[indsR,,drop=FALSE]
  attr(annMat, "colAnn") <- attr(annMat, "colAnn")[indsC,,drop=FALSE]
  annMat <- unclass(annMat)
  annMat <- annMat[indsR,indsC,drop=FALSE]
  class(annMat) <- append(class(annMat), "annMatrix")
  annMat
}

The basic idea is to make matrix preserve it's specific attributes after subsetting.

However I am running into a problem:

How to write "[" function in such a way that it behaves differently when called with and without a comma:

annMat[i]
annMat[i,]

as the default "[" for matrices seems to do.

I was thinking to set second argument to some value by default, but the value will not change because of an added comma.

Karolis Koncevičius
  • 9,417
  • 9
  • 56
  • 89

0 Answers0