我做错了什么?
> crossprod(1:3,4:6)
[,1]
[1,] 32
根据这个网站:http://onlinemschool.com/math/assistance/vector/multiply1/
它应该给:
{-3; 6; -3}
答案 0 :(得分:6)
这是一个广义的交叉产品:
xprod <- function(...) {
args <- list(...)
# Check for valid arguments
if (length(args) == 0) {
stop("No data supplied")
}
len <- unique(sapply(args, FUN=length))
if (length(len) > 1) {
stop("All vectors must be the same length")
}
if (len != length(args) + 1) {
stop("Must supply N-1 vectors of length N")
}
# Compute generalized cross product by taking the determinant of sub-matricies
m <- do.call(rbind, args)
sapply(seq(len),
FUN=function(i) {
det(m[,-i,drop=FALSE]) * (-1)^(i+1)
})
}
对于你的例子:
> xprod(1:3, 4:6)
[1] -3 6 -3
这适用于任何维度:
> xprod(c(0,1)) # 2d
[1] 1 0
> xprod(c(1,0,0), c(0,1,0)) # 3d
[1] 0 0 1
> xprod(c(1,0,0,0), c(0,1,0,0), c(0,0,1,0)) # 4d
[1] 0 0 0 -1
答案 1 :(得分:2)
crossprod执行以下操作:t(1:3) %*% 4:6
因此,它是1x3向量乘以3x1向量 - &gt;标量
答案 2 :(得分:2)
crossprod计算Matrix Product。要执行Cross Product,请编写您的函数,或者:
> install.packages("pracma")
> require("pracma")
> cross(v1,v2)
如果上面的第一行不起作用,请尝试:
> install.packages("pracma", repos="https://cran.r-project.org/web/packages/pracma/index.html”)
答案 3 :(得分:0)
crossProduct <- function(ab,ac){
abci = ab[2] * ac[3] - ac[2] * ab[3];
abcj = ac[1] * ab[3] - ab[1] * ac[3];
abck = ab[1] * ac[2] - ac[1] * ab[2];
return (c(abci, abcj, abck))
}
答案 4 :(得分:0)
您可以尝试expand.grid
expand.grid(LETTERS[1:3],letters[1:3])
输出:
Var1 Var2
1 A a
2 B a
3 C a
4 A b
5 B b
6 C b
7 A c
8 B c
9 C c