Kod 3.4 - NMSwBE

# Kod 3.4. Generowanie wartości z dwuwymiarowego rozkładu normalnego

library(MASS)

library(ellipse)

Id=matrix(c(1,0,0,1),2,2)

Sigma <- matrix(c(10,8,8,10),2,2)

z=mvrnorm(n=1000, rep(0, 2), Id)

plot(z)

z=z%*%Sigma

plot(z)

plot(z%*%solve(Sigma))