这可能是直截了当，基本等等，但我无法得到它。我有2个Nx1向量u和w，它们由负值和正值组成。我正在尝试计算w'u u'w，它应该是二次形式。我应该能够像
t(w)%*%u%*%t(u)%*%w
但有时我会得到一个负值，具体取决于两个向量中的值。这是不可能的，因为那是一种二次形式。我试过
crossprod(w, u)%*%crossprod(u, w)
和
crossprod(w, u)*crossprod(u, w)
给出积极和平等的结果。但是，由于我正在处理Nx1向量，我还应该能够将其写为
`sum(w*u)^2`
给出正值但与上述值不同。
所以我想我在某处做错了什么。所以，问题是：我如何表达对于向量和矩阵都有效的w'u u'w？
编辑：这里有一个csv file，其原始向量可以重现完全相同的问题