n<-3e2 # number of realizations
p<-.5*n # number of assets
gamma<-1
x <- matrix(data = rnorm(n*p), nrow = p, ncol = n)
# Simple MV portfolio
cov_mtrx <- Sigma_sample_estimator(x)
means <- rowMeans(x)
cust_port_simp <- new_MeanVar_portfolio(mean_vec=means,
cov_mtrx=cov_mtrx, gamma=2)
str(cust_port_simp)
# Portfolio with Bayes-Stein shrunk means
# and a Ledoit and Wolf estimator for covariance matrix
TM <- matrix(0, p, p)
diag(TM) <- 1
cov_mtrx <- CovarEstim(x, type="LW20", TM=TM)
means <- mean_bs(x)
cust_port_BS_LW <- new_MeanVar_portfolio(mean_vec=means$means,
cov_mtrx=cov_mtrx, gamma=2)
str(cust_port_BS_LW)
Run the code above in your browser using DataLab