我想使用Ryacas
进行符号矩阵运算,使用将基本R
矩阵转换为Ryacas
格式的函数。函数的结果似乎与Ryacas
格式匹配。但是当我尝试将矩阵相乘时,错误
# Error in aa %*% aa : requires numeric/complex matrix/vector arguments
抛出。下面的代码是显示这种情况的最小示例。
有什么建议吗?
library(Ryacas)
conv.mat <- function(x) {
conv <- lapply(1:nrow(x), function(i) paste0(x[i, ], collapse = ", "))
conv <- paste0("List(", paste0("List(", unlist(conv), ")", collapse = ", "), ")")
noquote(conv)
}
# Writing a matrix manually for Ryacas format
a <- List(List(1, 2), List(3, 7))
a * a
# expression(list(list(7, 16), list(24, 55)))
# Writing a matrix in R and convert it to Ryacas format by the function conv.mat
aa <- matrix(c(1, 2, 3, 7), 2, byrow = TRUE)
aa <- conv.mat(aa)
# [1] List(List(1, 2), List(3, 7))
aa * aa
# Error in aa * aa : non-numeric argument to binary operator
答案 0 :(得分:0)
首先,将Ryacas
而不是aa * aa
乘以aa %*% aa
矩阵。但这本身对您的情况无济于事,因为conv.mat
并不能完全满足我们的需要(一个表达式)。
我们可能会使用,例如
conv.mat <- function(x)
do.call(List, lapply(1:nrow(x), function(r) do.call(List, as.list(x[r, ]))))
然后
M <- matrix(c(1, 2, 3, 7), 2, byrow = TRUE)
M %*% M
# [,1] [,2]
# [1,] 7 16
# [2,] 24 55
M <- conv.mat(M)
M * M
# expression(list(list(7, 16), list(24, 55)))