我正在使用poweRlaw包来拟合一些数据并计算p值.我注意到,在某些情况下,数据没有太多的变化,再加上相对较少的观测,一旦模拟次数超过某个阈值,我就无法计算p值。
考虑以下数据:
library(poweRlaw)
vec<-c(1,1,1,1,1,2,1,1,4,1,1,1,1,1,1,1,1,2,1,1,1,2,2,1,2,1,
1,2,1,1,2,1,1,1,1,1,1,2,2,2,2,1,1,1,1,2)
可以拟合幂律并计算p值:
pl=displ$new(vec)
est=estimate_xmin(pl)
bts=bootstrap_p(pl,no_of_si