求不沉！

setwd("C:/Users/ll123/Desktop/R")
jsyh <- read.csv("建设银行.csv",head=T,stringsAsFactors = F)
jsyh$时间 <- as.Date(jsyh$时间)
n <- nrow(jsyh)
re <- log(jsyh$收盘价[2:n]/jsyh$收盘价[1:(n-1)])
par (mfrow = c(1,2))
plot(jsyh$时间,jsyh$收盘价,type='l',col='blue',xlab='日期',ylab='收盘价')
plot(jsyh$时间[2:n],re,type='l',col="red",xlab='日期',ylab='收益率')
redata <- exp(re)
lgwbull <- function (params)  {
  lambda <- params [1]
  k <- params [2]
  n <- length(redata)
  lgw <- 0
  for(i in 1:n) {
    lgw <- lgw + log((k/lambda) * ((redata [i] / lambda)^(k-1))*exp(-(redata [i] / lambda)^k))
    }
   return(lgw)  
}
xmle <- optim(c(1,1),lgwbull,method="BFGS",control=list(fnscale=-1))
lambda <- xmle$par[1]
k <- xmle$par[2]
alpha <- 0.05
dVaR <- log(qweibull(1-alpha,k,lambda))


这个接下来要怎么做呀

是全球某工商的同学吗

你好，请问你的问题解决了吗

哈哈哈哈  是的