r - 使用 R 进行重尾参数估计

Question

我需要估计类帕累托分布尾部的 alpha 参数。我想用 R 来做。经过一番谷歌搜索后，我找到了这个参考。我尝试直接安装带有函数的 asw 库install.packages('aws')。R回答了我

In getDependencies(pkgs, dependencies, available, lib) :
package ‘aws’ is not available

所以我尝试了sos库

findFn('awstindex')
found 0 matches
x has zero rows;  nothing to display.
Warning message:
In findFn("awstindex") : HIT not found in HTML;  processing one page only.

为什么 R 找不到这个库？是否有不同的函数可以用 R 测量尾部参数？

Answer 1

您对其他实现持开放态度吗？这是 fitdistrplus 包文档中提到的函数的使用，该函数在估计 Pareto 参数中具有特殊用途。

require(fitdistrplus)
require(actuar)  # to get a random Pareto function
rpar1000 <- rpareto(1000, 2.5, 3)
mgedist(rpar1000,"pareto",gof="ADR", start=list(shape=1, scale=1))
#------ the ADR measure weights the right tail observations ----------
$estimate
   shape    scale 
2.356235 2.871350 

$convergence
[1] 0

$value
[1] 0.2150991

$hessian
          shape     scale
shape  90.15040 -55.60725
scale -55.60725  35.04152

$gof
[1] "ADR"

$optim.function
[1] "optim"

$loglik
[1] -1612.233