

library(PerformanceAnalytics)
n_shares<-100
confident<-0.99
x<-read.csv("http://canisius.edu/~yany/data/ibmDaily.csv")
n<-nrow(x)
p<-x$Adj.Close
ret<-p[2:n]/p[1:(n-1)]-1
n<-nrow(x)
S<-skewness(ret) 
S
K<-kurtosis(ret) 
K

