# Reproducing VaR from Table 2.1 in page 47 of
# McNeal A., Frey R. and Embrechts P (2005).
alpha <- c(.90, .95, .975, .99, .995)
VaR(variance=(10000*0.2/sqrt(250))^2, alpha=alpha, model='both', df=4)
# only normal VaR results
VaR(variance=(10000*0.2/sqrt(250))^2, alpha=alpha)
Run the code above in your browser using DataLab