并排打印矩阵和向量

Question

出于教学目的,我希望能够并排打印或显示矩阵和向量,通常用于说明矩阵方程的结果,如$ A x = b $.

我可以使用SAS/IML执行此操作,其中print语句采用(空格分隔)表达式的任意集合,计算它们并打印结果,例如,

print A ' * ' x '=' (A * x) '=' b;

                    A                     X   #TEM1001       B
                    1     1    -4  *  0.733 =        2 =     2
                    1    -2     1     -0.33          1       1
                    1     1     1      -0.4          0       0

请注意,引用的字符串按原样打印.

我已经搜索过了,但在R中却找不到这样的东西.我想这样的事情可以通过一个函数来完成,showObj(object, ...)它取一个参数列表,将每一个格式化为一个字符块,并将它们并排加入.

这种方法的另一个用途是将3D阵列显示为其切片的并排集合的紧凑方式.

这是响铃还是有人有建议开始？

Answer 1

我创建了一个非常简单的函数，可以打印中间带有任意字符串（通常是运算符）的矩阵和向量。它允许具有不同行数的矩阵，并将向量视为列矩阵。它不是很详细，所以我担心有很多失败的例子。但对于像你问题中的例子一样简单的例子，应该足够了。

format()用于将数字转换为字符。这样做的优点是矩阵的所有行都具有相同的宽度，因此在打印时可以很好地对齐。如果需要，您可以添加一些参数作为format()参数mat_op_print()以进行配置。作为示例，我添加了width可用于控制列的最小宽度的参数。

如果矩阵和向量是函数调用中的名称，则这些名称将作为标题打印在第一行中。否则，仅打印数字。

所以，这是这个函数：

mat_op_print <- function(..., width = 0) {

  # get arguments
  args <- list(...)
  chars <- sapply(args, is.character)

  # auxilliary function to create character of n spaces
  spaces <- function(n) paste(rep(" ", n), collapse = "")

  # convert vectors to row matrix
  vecs <- sapply(args, is.vector)
  args[vecs & !chars] <- lapply(args[vecs & !chars], function(v) matrix(v, ncol = 1))

  # convert all non-characters to character with format
  args[!chars] <- lapply(args[!chars], format, width = width)

  # print names as the first line, if present
  arg_names <- names(args)
  if (!is.null(arg_names)) {
    get_title <- function(x, name) {
      if (is.matrix(x)) {
        paste0(name, spaces(sum(nchar(x[1, ])) + ncol(x) - 1 - nchar(name)))
      } else {
        spaces(nchar(x))
      }
    }
  cat(mapply(get_title, args, arg_names), "\n")
  }

  # auxiliary function to create the lines
  get_line <- function(x, n) {
    if (is.matrix(x)) {
      if (nrow(x) < n) {
       spaces(sum(nchar(x[1, ])) + ncol(x) - 1)
      } else {
        paste(x[n, ], collapse = " ")
      }
    } else if (n == 1) {
      x
    } else {
      spaces(nchar(x))
    }
  }

  # print as many lines as needed for the matrix with most rows
  N <- max(sapply(args[!chars], nrow))
  for (n in 1:N) {
    cat(sapply(args, get_line, n), "\n")
  }
}

这是它如何工作的一个例子：

A = matrix(c(0.5, 1, 3, 0.75, 2.8, 4), nrow = 2)
x = c(0.5, 3.7, 2.3)
y = c(0.7, -1.2)
b = A %*% x - y
mat_op_print(A = A, " * ", x = x, " - ", y = y, " = ", b = b, width = 6)
## A                        x          y          b      
##   0.50   3.00   2.80  *     0.5  -     0.7  =  17.090 
##   1.00   0.75   4.00        3.7       -1.2     13.675 
##                             2.3    

还可以并排打印 3 维数组的切片：

A <- array(1:12, dim = c(2, 2, 3))
mat_op_print(A1 = A[, , 1], " | ", A2 = A[, , 2], " | ", A3 = A[, , 3])
## A1      A2      A3    
## 1 3  |  5 7  |   9 11 
## 2 4     6 8     10 12