Question

我需要在特定矩阵中获取每对列的内积，我通过计算

来实现

t(M) %*% M

然而，这自然产生了对称的结果，只做了两倍以上的必要工作（我也不需要对角线）。显然，我可以将乘法分解为单个内部产品操作，但有没有更好的方法来计算产品的上三角部分？

Answer 1

来自help("crossprod")中的说明：

给定矩阵x和y作为参数，返回矩阵交叉积。这正式等同于（但通常略快）致电t(x) %*% y（crossprod）或x %*% t(y)（tcrossprod）。

因此，请使用crossprod(M)。