R中的极值函数

我正在研究一个极端值问题 根据Pevsner（2003，第103页）给出了一个有趣的极值分布

我试图从标准正态分布生成一个样本（大小为1000），并重复1000次。然后，从这些极大值中减去给定函数，然后除以bn，其中

---

fn只需关闭for循环

}

上的大括号，然后移除

）

。尽管如此，你还是试图在那里记录负数…@timriffe，然后它会弹出很多警告这些警告主要是因为n中有负数，而你正在记录这些，生成

NaN

值。另外，看起来您没有按照该文本中的说明操作：您应该取最大值

rnorm（1000）

，1000次。其余的是给你的you@timriffe哦，我明白了！这就是为什么。谢谢

> my.stat <-NULL 
> for (i in 1:1000) { 
+  n <- rnorm(1000)
+   fn <- exp(-n)*exp(-exp(-n)) 
+   an <- sqrt(2*log(n)) - 0.5*(log(log(n))+log(4*pi))*(2*log(n))^(-1/2)
+   bn <- (2*log(n))^(-1/2)
+ my.stat <- c(my.stat, sum(sum(fn-an)/bn)))
> par(mfrow=c(2,2))
> hist(my.stat,freq=FALSE,main="histogram of 1000 M",xlab="M")

+ my.stat <- c(my.stat, sum(sum(fn-an)/bn)))