Maticové operácie

A <- matrix(c(1,2,3,4), nrow = 2)
B <- matrix(c(5,6,7,8), nrow = 2)

A * B
     [,1] [,2]
[1,]    5   21
[2,]   12   32
(2*A) %*% (2*B)
     [,1] [,2]
[1,]   92  124
[2,]  136  184
det(A)
[1] -2
det(B)
[1] -2
t(A) + B
     [,1] [,2]
[1,]    6    9
[2,]    9   12

Zlučovanie vektorov do matíc

v1 <- c(1, 2, 3)
v2 <- c(4, 5, 6)
C <- cbind(v1, v2)
D <- rbind(v1, v2)

Vypočítanie zvolenej štatistiky po riadkoch (stĺpcoch) matice

M <- matrix(1:9, nrow = 3)
M
     [,1] [,2] [,3]
[1,]    1    4    7
[2,]    2    5    8
[3,]    3    6    9
apply(M, 1, sum)
[1] 12 15 18
apply(M, 1, mean)
[1] 4 5 6
apply(M, 2, sum)
[1]  6 15 24
apply(M, 1, mean)
[1] 4 5 6

Malé cvičenie

Vytvorte maticu 5x5 s hodnotami po riadkoch 1..25, vypočítajte stĺpcové sumy a súčin matíc \(M^t M\).

M5 <- matrix(1:25, nrow = 5, byrow = TRUE)
colSums(M5)
[1] 55 60 65 70 75
t(M5) %*% M5
     [,1] [,2] [,3] [,4] [,5]
[1,]  855  910  965 1020 1075
[2,]  910  970 1030 1090 1150
[3,]  965 1030 1095 1160 1225
[4,] 1020 1090 1160 1230 1300
[5,] 1075 1150 1225 1300 1375

LS0tCnRpdGxlOiAiRG9tw6FjYSDDumxvaGEgMiAobWF0aWNlKSIKYXV0aG9yOiAiSWdvciBaZWxlbmF5ICA8YnI+CihzIHZ5dcW+aXTDrW0gQ2hhdEdQVCkiCmRhdGU6ICJTZXB0ZW1iZXIgMjAyNSIKb3V0cHV0OiAKICBodG1sX25vdGVib29rOgogICAgdG9jOiB0cnVlCiAgICB0b2NfZmxvYXQ6IHRydWUKICAgIHRoZW1lOiBjb3NtbwogICAgaGlnaGxpZ2h0OiBicmVlemVkYXJrCmVkaXRvcl9vcHRpb25zOiAKICBtYXJrZG93bjogCiAgICB3cmFwOiA3MgotLS0KCiMjIE1hdGljb3bDqSBvcGVyw6FjaWUKYGBge3J9CkEgPC0gbWF0cml4KGMoMSwyLDMsNCksIG5yb3cgPSAyKQpCIDwtIG1hdHJpeChjKDUsNiw3LDgpLCBucm93ID0gMikKCkEgKiBCCigyKkEpICUqJSAoMipCKQpkZXQoQSkKZGV0KEIpCnQoQSkgKyBCCgpgYGAKIyMgWmx1xI1vdmFuaWUgdmVrdG9yb3YgZG8gbWF0w61jIAoKYGBge3J9CnYxIDwtIGMoMSwgMiwgMykKdjIgPC0gYyg0LCA1LCA2KQpDIDwtIGNiaW5kKHYxLCB2MikKRCA8LSByYmluZCh2MSwgdjIpCmBgYAoKIyMgVnlwb8SNw610YW5pZSB6dm9sZW5laiDFoXRhdGlzdGlreSBwbyByaWFka29jaCAoc3TEunBjb2NoKSBtYXRpY2UKCmBgYHtyfQpNIDwtIG1hdHJpeCgxOjksIG5yb3cgPSAzKQpNCmFwcGx5KE0sIDEsIHN1bSkKYXBwbHkoTSwgMSwgbWVhbikKYXBwbHkoTSwgMiwgc3VtKQphcHBseShNLCAxLCBtZWFuKQpgYGAKCiMjIE1hbMOpIGN2acSNZW5pZQoKPiBWeXR2b3J0ZSBtYXRpY3UgNXg1IHMgaG9kbm90YW1pIHBvIHJpYWRrb2NoIDEuLjI1LCB2eXBvxI3DrXRhanRlIHN0xLpwY292w6kgc3VteSBhIHPDusSNaW4gbWF0w61jIFwoTV50IE1cKS4KCmBgYHtyfQpNNSA8LSBtYXRyaXgoMToyNSwgbnJvdyA9IDUsIGJ5cm93ID0gVFJVRSkKY29sU3VtcyhNNSkKdChNNSkgJSolIE01CmBgYAotLS0KCgo=