如何使用 R SPREDA 包提取 Fréchet 拟合模型的参数 xi 值？

如何解决如何使用 R SPREDA 包提取 Fréchet 拟合模型的参数 xi 值？

我正在尝试遵循 this post，它使 Frechet 分布适合某些风测量，如下所示：

require('SPREDA')
require(e1071)
require(extRemes)
require(survival)

#The R functions for performing the extreme value analysis (for maxima)
#can be downloaded from: 
source("http://www.datall-analyse.nl/R/eva_max.R")

#data from Castillo et al.,Table 1.1,p. 9-10
#yearly maximum wind speed (in miles/hour)
evobs <- scan("http://www.datall-analyse.nl/blog_data/extremes_Table1-1.txt")

#explore data visually
options(repr.plot.width=7,repr.plot.height=7)
hist(evobs)

#Fréchet
frechetmod <- Lifedata.MLE(Surv(evobs) ~ 1,dist="frechet")
frechetmod

在某些时候，作者使用输出或函数 Lifedata.MLE() 来计算以下返回级别：

#note：对于wind数据，返回周期为20表示每 #20 年风速预计（平均）大于 #muG+qlev(1-1/20)*sigmaG=49.4 mph（在 Gumbel 分布的情况下）， #或 exp(muF+qlev(1-1/20)*sigmaF)=51.4 mph（在 Fréchet 分布的情况下） #（这个预期的风速也称为回归水平）

我想这需要 xi 形状参数来计算 51.4 mph。但是，此参数未包含在 frechetmod 的输出中：

Call:
model.frame(formula = Surv(evobs) ~ 1)

Coefficients:
(Intercept)    logsigma 
  3.3579367   0.1957559 

Loglikelihod:
-173.4035 (df=2)

我怎样才能得到这个形状参数估计？