La matrice de la Multiplication dans r

Je suis en train de multiplier chaque ligne d'une matrice par la colonne d'une autre matrice. Par exemple:

mat1 <- matrix(rnorm(10), nrow=5, ncol=2)
mat2 <- matrix(rnorm(5), nrow=5)

Je veux multiplier chaque ligne de mat1 par mat2. La forme souhaitée de la sortie est de 5*2.

Si vous voulez multiplication de matrice, comme le suggère le titre, c'est à dire, vous voulez multiplier les lignes de la première matrice par des colonnes de la deuxième, c'est mat1 %*% mat2. Mais cela exige de la première matrice d'avoir autant de colonnes que le deuxième a des lignes -- qui n'est pas le cas dans ton exemple.

OriginalL'auteur mike | 2012-03-15

  1. 5

    Vous pouvez simplement utiliser apply() à multiplier chaque colonne de mat1 par mat2. (Le "*" effectuera R est généralement vectorisé élément-sage multiplication des deux égale à la longueur des vecteurs). 

    apply(mat1, 2, "*", mat2)
           [,1]       [,2]
[1,]  0.1785476  0.4175557
[2,]  0.2644247 -0.3745997
[3,] -0.5328542  0.8945527
[4,] -2.7351502 -0.7715341
[5,] -0.9719129 -0.1346929

    Ou, mieux encore, de les convertir mat1 à un vecteur pour prendre avantage de la R de règles en matière de recyclage:

    mat2 <- matrix(1:10, ncol=2)
mat1 <- matrix(1:5, ncol=1)

as.vector(mat1)*mat2
     [,1] [,2]
[1,]    1    6
[2,]    4   14
[3,]    9   24
[4,]   16   36
[5,]   25   50

    OriginalL'auteur Josh O'Brien

  2. 1

    Votre première matrice a cinq lignes et deux colonnes; votre deuxième matrice a cinq lignes et une colonne. Si ils ont le même nombre de lignes et la seconde a toujours une colonne vous pouvez faire 

    mat1 * rep(mat2,ncol(mat1))
           [,1]        [,2]
[1,] -0.2327958  0.76093047
[2,] -0.3636661 -0.18991299
[3,] -0.8729468  0.58214118
[4,]  0.8017349 -0.59781909
[5,] -0.2230380 -0.08296606

    Si mat1 effectivement eu autant d'éléments dans ses lignes que mat2avait dans sa colonne unique (tel que vos mots suggèrent) vous permettra d'adapter un peu

    mat1 <- matrix(rnorm(10), nrow=2, ncol=5)
mat2 <- matrix(rnorm(5), nrow=5, ncol=1)
mat1 * rep(mat2,nrow(mat1))
            [,1]        [,2]       [,3]        [,4]       [,5]
[1,] -0.19818805 -0.05938007 -1.7792597  0.06937307 -0.7193403
[2,] -0.05087793  0.10781853  0.2243285 -0.11416273  2.4063926

    ou dans de sarah version

    mat1 <- matrix(rnorm(10), nrow=5, ncol=2)
mat2 <- matrix(rnorm(2), nrow=2, ncol=1)
mat1 * rep(mat2,nrow(mat1))
           [,1]        [,2]
[1,]  0.1528393  0.68646359
[2,]  0.2420454  0.22987250
[3,] -0.2592124 -0.07626098
[4,]  0.4431273  0.27320838
[5,] -0.1698307  0.47578667

    OriginalL'auteur Henry