Comment comparer deux matrices pour voir si elles sont identiques dans R?

Par exemple, j'ai deux matrices et je veux savoir si elles sont identiques dans chaque élément.

mymatrix<-Matrix(rnorm(20),ncol=5)
mysvd<-svd(mymatrix) 
newmatrix<-mysvd$u %*% diag(mysvd$d) %*% t(mysvd$v)

J'ai utilisé les méthodes suivantes pour comparer:

identical(Matrix(newmatrix), mymatrix)
all.equal(Matrix(newmatrix), mymatrix)

Pourquoi le premier ne retourne pas VRAI?
Peu importe, je utiliser la Matrice de la matrice de paquet ou de la matrice à partir de la base de package

OriginalL'auteur Tyler 十三将士归玉门 | 2014-04-12

  1. 8

    Ils ne sont pas exactement égales (par identical) parce que de très petites différences:

    > max(abs(Matrix(newmatrix) - mymatrix))
[1] 1.110223e-15

    mais ces différences sont inférieures à la valeur par défaut tolerance à l'intérieur de all.equal:

    > .Machine$double.eps ^ 0.5
[1] 1.490116e-08

    donc identical sera de retour FALSE et all.equal sera de retour TRUE.

    OriginalL'auteur flodel

  2. 7

    Essayez ce qui suit:

    > dput(Matrix(newmatrix))
new("dgeMatrix"
    , x = c(-0.705193264426589, 0.68023073271425, 0.0726318059033283, -0.111055227906436, 
-0.113940777963113, 0.726463241417717, -0.343435098646076, 0.885225942372688, 
-0.549848405897803, -0.0227469387867766, -0.927524398860002, 
0.58047674424687, 0.521144348439824, 0.279602090928527, -1.31686400403363, 
0.906874499735628, -0.276997805548975, 0.632960950203858, 0.453881309098762, 
-0.00528540521655077)
    , Dim = 4:5
    , Dimnames = list(NULL, NULL)
    , factors = list()
)

> dput(newmatrix)
structure(c(-0.705193264426589, 0.68023073271425, 0.0726318059033283, 
-0.111055227906436, -0.113940777963113, 0.726463241417717, -0.343435098646076, 
0.885225942372688, -0.549848405897803, -0.0227469387867766, -0.927524398860002, 
0.58047674424687, 0.521144348439824, 0.279602090928527, -1.31686400403363, 
0.906874499735628, -0.276997805548975, 0.632960950203858, 0.453881309098762, 
-0.00528540521655077), .Dim = 4:5)

    Clairement ce ne sont pas les mêmes structures de données.

    OriginalL'auteur G. Grothendieck

  3. 3

    Le premier vous de comparer:

    > Matrix(newmatrix)
4 x 5 Matrix of class "dgeMatrix"
           [,1]       [,2]       [,3]        [,4]       [,5]
[1,]  0.5052901 -0.3264201 -0.8576401 -0.62666359  2.1076090
[2,]  0.2356111  0.4911067 -1.2376674  1.11231840  0.8576557
[3,] -0.6244670  1.4423943 -1.2820541 -0.05297437 -2.0458810
[4,] -0.2669079  1.1218459  0.6371571 -0.52168139  0.2163623

    avec:

    > mymatrix
           [,1]       [,2]       [,3]        [,4]       [,5]
[1,]  0.5052901 -0.3264201 -0.8576401 -0.62666359  2.1076090
[2,]  0.2356111  0.4911067 -1.2376674  1.11231840  0.8576557
[3,] -0.6244670  1.4423943 -1.2820541 -0.05297437 -2.0458810
[4,] -0.2669079  1.1218459  0.6371571 -0.52168139  0.2163623

    Et si vous lisez l'aide pour identical vous voir: The safe and reliable way to test two objects for being _exactly_ equal. Vous comparez un Matrix-objet de classe avec une R de la matrice de l'objet. Ils sont très différents. Même deux objets avec des attributs différents ne seront pas identiques:

    > x
a b c d 
1 2 3 4 
> y
[1] 1 2 3 4
> x==y
   a    b    c    d 
TRUE TRUE TRUE TRUE 
> identical(x,y)
[1] FALSE

    OriginalL'auteur Spacedman

  4. -2

    Vous pouvez vérifier si de 2 matrices sont identiques (à l'identique) ou non comme suit.

    Supposons que vous disposez de 2 matrices, newMatrix et oldMatrix, qui pourrait être n'importe quelle dimension.

    which (which (newMatrix == oldMatrix) == FALSE) sera de retour integer (0) si les deux matrices sont identiques.

    OriginalL'auteur Saket Shrivastava