# c<1
n <- 3e2 # number of realizations
p <- .5*n # number of assets
b <- rep(1/p,p)
# Assets with a diagonal covariance matrix
x <- matrix(data = rnorm(n*p), nrow = p, ncol = n)
test <- new_GMV_portfolio_weights_BDPS19(x=x, b=b, beta=0.05)
str(test)
# Assets with a non-diagonal covariance matrix
Mtrx <- RandCovMtrx(p=p)
x <- t(MASS::mvrnorm(n=n , mu=rep(0,p), Sigma=Mtrx))
test <- new_GMV_portfolio_weights_BDPS19(x=x, b=b, beta=0.05)
summary(test)
# c>1
p <- 1.3*n # number of assets
b <- rep(1/p,p)
# Assets with a diagonal covariance matrix
x <- matrix(data = rnorm(n*p), nrow = p, ncol = n)
test <- new_GMV_portfolio_weights_BDPS19_pgn(x=x, b=b, beta=0.05)
str(test)
Run the code above in your browser using DataLab