急求！！！关于R语言算极值问题。

a<-rnorm(200)
gpd.fitrange(a, 0, 1, nint = 100)

是不是你上下限设置有问题

你看看你的data数据结构，summary(data)

我看你的三分位数才到0.01，那么改为0.002到0.01

数据量估计不够

function (data, umin, umax, nint = 10, show = FALSE, ...)
{
m <- s <- up <- ul <- matrix(0, nrow = nint, ncol = 2)
u <- seq(umin, umax, length = nint)
for (i in 1:nint) {
z <- gpd.fit(data, u[i], show = show, ...)
m[i, ] <- z$mle
m[i, 1] <- m[i, 1] - m[i, 2] * u[i]
d <- matrix(c(1, -u[i]), ncol = 1)
v <- t(d) %*% z$cov %*% d
s[i, ] <- z$se
s[i, 1] <- sqrt(v)
up[i, ] <- m[i, ] + 1.96 * s[i, ]
ul[i, ] <- m[i, ] - 1.96 * s[i, ]
}
names <- c("Modified Scale", "Shape")
oldpar <- par(mfrow = c(2, 1))
for (i in 1:2) {
um <- max(up[, i])
ud <- min(ul[, i])
plot(u, m[, i], ylim = c(ud, um), xlab = "Threshold",
ylab = names[i], type = "b")
for (j in 1:nint) lines(c(u[j], u[j]), c(ul[j, i], up[j,
i]))
}
par(oldpar)
invisible(list(thresholds = u, mle = m, se = s, ci.low = ul,
ci.up = up))
}
<environment: namespace:ismev>

@米粒cheer
找到了min(y)的位置了
在gpd.fit函数里

u就是threshold，GPD的分布参数有两个，ζ和δ，分别对应下图的第二个图和第一个图，横坐标就是u的不同取值，函数gpd.fitrange就是寻找使gpd分布最小的那对参数ζ和δ（不知理解对不对）：

@米粒cheer
终于知道你的问题在哪里了
你自己看看length（data）是多少，把你的data.frame转换成向量吧