如何用R中bootstrap_p获取Hill图得到的自定义xmin对应的p值？

要获取Hill图得到的xmin=0.44和alpha=2.12对应的Bootstrap p值，只需手动更新幂律模型对象的参数后再调用bootstrap_p函数即可，具体步骤如下：

1. 初始化模型并设置自定义参数
   基于你的数据创建幂律模型对象后，直接用setXmin()和setPars()方法替换为Hill图得到的参数值：

   # 初始化幂律模型
   pl_m <- conpl$new(data)
   # 设置Hill图得到的xmin
   pl_m$setXmin(0.44)
   # 设置Hill图得到的alpha
   pl_m$setPars(2.12)
   

2. 调用bootstrap_p计算p值
   此时调用bootstrap_p函数，会基于你手动指定的参数进行幂律假设检验：

   # 计算指定参数下的Bootstrap p值
   boot_result <- bootstrap_p(pl_m)
   # 输出p值
   print(boot_result$p)
   

额外说明 #

• 若需要更精确的结果，可通过no_of_sims参数调整Bootstrap迭代次数，比如bootstrap_p(pl_m, no_of_sims = 5000)
• 确保你的powerlaw包为最新版本，避免参数兼容性问题

内容的提问来源于stack exchange，提问作者Neda Fathi