## R-code for assignment 4

## 4.6
Mu <- c(1,-1,2)
Sigma <- matrix(c(4,0,-1,0,5,0,-1,0,2), 3,3, byrow=T)
Mu
Sigma

library(mvtnorm)
pairs(rmvnorm(100, Mu, Sigma)) 


## 4.7 
par(mfrow=c(2,1))
plot(density(rnorm(100, 1-0.5*(2-2), sqrt(3.5))), main="Density of x1|x3")
plot(density(rnorm(100, 1-0.5*(-2-2), sqrt(3.5))), main="Density of x1|x2,x3")