###################################### #Script TP2 parte 2 Ejercicio 7 ###################################### library(MASS) #c) vect.T <- function(k,p) { media <- rep(0,p) sigma <- diag(1,p) x <- (sqrt(k/rchisq(1,k)))*mvrnorm(1,media,sigma) return(x) } vect.T(7,4) vect.T(10,7) ###################################### #e) vectMT <- function(mu,sigma,k,p) { C <- t(chol(sigma)) z <- vect.T(k,p) x <- mu+C%*%z return(x) } ###################################### #Prueba k <- 10 p<- 5 sigma <- diag(5,p) mu <- c(1,2,1,3,4) m <- 20 w <- matrix(0,p,m) for(i in 1:m) { set.seed(100+i) w[,i] <- vectMT(mu,sigma,k,p) }