时间戳

首页 > 解决方案 > 在R中乘法矩阵时结果不正确？

r - 在R中乘法矩阵时结果不正确？

问题描述

在 R 中将这两个矩阵相乘时，我得到了一些奇怪的结果：

> matrix3

    [1,]  3.19747172 -2.806e-05 -0.00579284 -0.00948720 -0.01054026  0.17575719
    [2,] -0.00002806  2.000e-08  0.00000057  0.00000006 -0.00000009 -0.00000358
    [3,] -0.00579284  5.700e-07  0.00054269  0.00001793 -0.00002686 -0.00310465
    [4,] -0.00948720  6.000e-08  0.00001793  0.00003089  0.00002527 -0.00066290
    [5,] -0.01054026 -9.000e-08 -0.00002686  0.00002527  0.00023776 -0.00100898
    [6,]  0.17575719 -3.580e-06 -0.00310465 -0.00066290 -0.00100898  0.03725362

> matrix4
              [,1]
    x0    2428.711
    x1 1115178.561
    x2   74411.013
    x3  925700.445
    x4   74727.396
    x5   13342.182

> matrix3%*%matrix4
                  [,1]
    [1,] 78.4244581753
    [2,] -0.0023802299
    [3,]  0.1164568885
    [4,] -0.0018504732
    [5,] -0.0006493249
    [6,] -0.1497822396

问题是，如果您尝试在 excel 中将这两个矩阵相乘，您会得到：

>78.4824494081686
>-0.0000419022486847151
>0.112430295996347
>-0.000379343461780479
>0.000340414687578061
>-0.14454024116344

并且使用在线矩阵我也得到了 excel 的结果。希望您能帮助理解如何在 R 中获得相同的结果。

标签： rmatrix

解决方案

由于使用了库（matlib）中的函数 inv() 导致出现问题。matrix3 是使用 inv() 函数求逆的结果。不知道为什么当我使用solve() 求逆然后正常继续时，我得到了正确的矩阵。也许在 inv() 函数中有某种舍入。

推荐阅读