我最近想出了一个问题，cross prod 似乎表现得很奇怪。
> is.symmetric.matrix(crossprod(M1))
[1] TRUE
> is.symmetric.matrix(crossprod(M1,M1))
[1] FALSE
> dput(M1)
structure(c(0.662345679012346, -0.0537037037037037, -0.0537037037037037, -0.0537037037037037, -0.0537037037037037, -0.0537037037037037, -0.0537037037037037, -0.0537037037037037, -0.0555533333333333, -0.0555533333333333, -0.0555533333333333, -0.0555533333333333, -0.0555533333333333, -0.0555533333333333, -0.0555533333333333, -0.0555533333333333 , -0.082849523255722, 0.209560558823297, 0.794380722981334, 0.794380722981334, -0.082849523255722, 0.209560558823297, -0.082849523255722, -0.082849523255722, -0.0333333333333333, 0.05, -0.0333333333333333, 0.133333333333333, -0.0333333333333333, -0.0333333333333333, -0.0333333333333333, 0.05, 0.00722926923433098, 1.91575634709772, 0.00722926923433098, -0.513278115637501 ，-0.686780577261445，-0.339775654013557，-0.686780577261445，-0.686780577261445，0.138472834067548，-0。0585636096649313, 0.0727940194900547, -0.124242424242424, 0.0727940194900547, 0.00711520491256175, -0.124242424242424, -0.124242424242424, -0.665219916932936, -1.33980912847056, -0.271709543535988, 0.0655850622328247, -2.29547717814886, -2.29547717814886, -1.73331950186751, 1.58341078819248, -0.0829545454545454, -0.0829545454545454, -0.0829545454545454, -0.0829545454545454 , 0.0725852272727273, -0.0829545454545454, -0.0829545454545454, -0.0829545454545454, -0.865647461169836, -1.73543677622087, 4.41521623735429, -1.48692554334914, -0.555008420080178, -4.46906033780982, -3.72352663919465, 2.92414884012395, -0.00937499999999999, -0.00937499999999999, -0.00937499999999999, 0.69375, -0.15, -0.00937499999999999 , -0.15, -0.15, -8.09053494329047, -6.04111906819099, 0.852370693507286, 0.666060159407333, -3.24646105669169,-2.31490838619192, -10.13995081839, -3.24646105669169), .Dim = c(8L, 11L), .Dimnames = list(NULL, NULL))
以前任何人都遇到过这样的问题。我已输入输出以重现此问题。它仅对少数类型的矩阵失败。虽然有很多选择，但我想知道在执行 crossprod(X,X) 而不是 crossprod(X) 时我是否遗漏了什么
$version.string [1] "R 版本 3.4.0 (2017-04-21)"