具有R的幂律分布

具有R的幂律分布,r,power-law,R,Power Law,我试图用下面的R代码来形象化幂律p(x)=x^(-2.5)。当你在最后使用对数刻度时,你会得到很多振动,这是可以看到的 但是我知道,这是我的问题,,我读了一篇文章,作者说我必须使用累积分布函数来消除最后的振动。但对我来说,这不起作用,从这里可以看出 library(ggplot2) chol_r <- read.table("C:\\Users\\me\\Desktop\\1M_just_random_py.txt", header = FALSE

我试图用下面的R代码来形象化幂律p(x)=x^(-2.5)。当你在最后使用对数刻度时,你会得到很多振动,这是可以看到的

但是我知道,这是我的问题,,我读了一篇文章,作者说我必须使用累积分布函数来消除最后的振动。但对我来说,这不起作用,从这里可以看出

library(ggplot2)

chol_r <- read.table("C:\\Users\\me\\Desktop\\1M_just_random_py.txt", 
                   header = FALSE)


chol <- (chol_r)**(-2.5) #this p(x)
chol2 = (1/1.5)*chol_r**(-1.5) # the cumulative  distribution function

qplot(chol2,
      geom="histogram",
      binwidth = 0.001, #0.001 oder 0.38
      main = "Histogram",
      xlab = "Numbers",
      fill=I("blue"),
      col=I("red"),
      log="xy")
库(ggplot2)

chol_r文本文件中有什么,您能生成代码来创建这个1M_just_random_py文件吗?我在for循环中使用Python命令random()来创建数字。这篇文章可以在第3页和第4页的:[link]()中找到。描述了该过程