#Ejercicio6

#Es necesario instalar el paquete "mvtnorm", el cual contiene la función rmvnorm, con la cual se puede simular lo pedido en el ejercicio.
library("mvtnorm")
m<-c(3,4)
r<-c(10,0.9,0.9,4)
dim(r)<-c(2,2)
Q<-mvtnorm::rmvnorm(1000,m,r)
head(Q)

##           [,1]     [,2]
## [1,] 0.2740138 3.502702
## [2,] 0.6686197 3.946295
## [3,] 5.4256669 4.130718
## [4,] 2.1114204 6.632934
## [5,] 9.3246362 2.028236
## [6,] 2.6081605 5.513332

###Ahora procederemos a graficar la dispersión de la muestra

plot(Q, pch=21, xlim=c(-20, 20), ylim=c(-20,20), xlab="Marginal No. 1", ylab="Marginal No. 2", bg = "lightblue")

###Ahora se obtendrá la media muestral y se comparará con la media poblacional

colMeans(Q)

## [1] 3.177082 4.104100

m

## [1] 3 4

m-colMeans(Q)

## [1] -0.177082 -0.104100