NIM : 220605110070

KELAS : A

MATKUL : LINEAR ALGEBRA

DOSEN PENGAMPU : Prof.Dr.Suhartono,M.Kom

JURUSAN : TEKNIK INFORMATIKA

LEMBAGA : UNIVERSITAS ISLAM NEGERI MAULANA MALIK IBRAHIM MALANG

## library need to open a jpg file
library(jpeg)
## Warning: package 'jpeg' was built under R version 4.2.2
## Download from the author’s website
myurl <- "http://polytopes.net/Tora_Sleeping.JPG"
z <- tempfile()
download.file(myurl,z,mode="wb")
Kitty <- readJPEG(z)
## dimension
d <- dim(Kitty)
d
## [1] 2340 4160    3

Kami menggunakan contoh kerja di Bagian 2.2. Seperti pada Bagian 2.2, kami menghasilkan titik data simulasi di bawah distribusi normal bivariat. Ingatlah bahwa a distribusi normal bivariat didefinisikan oleh mean dan standar deviasi, tetapi tidak seperti distribusi normal univariat, rata-rata didefinisikan oleh vektor dua dimensi dan standar deviasi didefinisikan oleh matriks 2 × 2. Sekali lagi, kami menggunakan paket R mvtnorm [18] dan kami menggunakan paket matlib untuk operasi matriks. Pertama, kita unggah semua paket yang dibutuhkan:

library(mvtnorm)
library(ggplot2)
## Warning: package 'ggplot2' was built under R version 4.2.3
library(matlib)
## Warning: package 'matlib' was built under R version 4.2.2
sigma <- matrix(c(4,3,3,4), ncol = 2, nrow = 2)
## Mean
mu <- c(5, 5)
n <- 1000
x <- rmvnorm(n = n, mean = mu, sigma = sigma)
d <- data.frame(x)
p2 <- ggplot(d, aes(x = X1, y = X2)) +
geom_point(alpha = .5) +
geom_density_2d()
p2

di mana µ adalah rata-rata dan σ adalah matriks kovarians yang mendefinisikan deviasi standar untuk distribusi normal bivariat. Perhatikan bahwa matriks kovarians ini matriks kovarians ini adalah simetris. Untuk menentukan standar deviasi untuk distribusi normal multivariat, matriks kovarians harus simetris. Namun, dalam kehidupan nyata, kita tidak mengamati µ dan σ. Kita hanya mengamati set data x yang berisi 1000 titik data. Oleh karena itu kita harus mengestimasi mereka dari titik-titik data. Sekarang kita mengestimasi matriks kovarian. Misalkan x adalah matriks yang baris-barisnya mewakili setiap titik data dalam set data dan yang kolomnya mewakili setiap variabel. Maka kita dapat mengestimasi matriks kovarian σ^ dengan xT - x, di mana xT adalah transpose dari matriks x. ^ σ adalah matriks simetris berukuran k × k simetris dengan k adalah jumlah variabel acak. Dalam R, kita dapat mengestimasi dengan:

y <- t(x) %*% x
eigen(y)
## eigen() decomposition
## $values
## [1] 57892.198   988.564
## 
## $vectors
##           [,1]       [,2]
## [1,] 0.7056294 -0.7085811
## [2,] 0.7085811  0.7056294
prcomp(d, center = FALSE)
## Standard deviations (1, .., p=2):
## [1] 7.6124995 0.9947631
## 
## Rotation (n x k) = (2 x 2):
##           PC1        PC2
## X1 -0.7056294  0.7085811
## X2 -0.7085811 -0.7056294