# Kod 3.16 Generowanie wartości losowych z dwuwymiarowego rozkładu normalnego o rozkładach brzegowych=====
require(mvtnorm)
S <- matrix(c(1,.5,.5,1),2,2)
A <- rmvnorm(mean=c(0,0),sig=S,n=1000)
U <- pnorm(A)
x <- qgamma(U[,1],2) #x is gamma distributed
y <- qbeta(U[,2],1,2) #y is beta distributed
plot(x,y) #They correlate!
