我尝试使用以下 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")
那么有人知道我做错了什么吗?或者我怎样才能在没有振动的情况下直线下降?我真的不知道我做错了什么