Я только что понял, что в этом посте, который я собрал , гистограмма имеет только значения индекса по оси x вместо предполагаемого диапазона от 0 до 5:
Here is the code for convenience:
x1 = 5 # Maximum value
x0 = 0.1 # It can't be zero; otherwise X^0^(neg) is 1/0.
alpha = -2.5 # It has to be negative.
y = runif(1e5) # Number of samples
x = ((x1^(alpha+1) - x0^(alpha+1))*y + x0^(alpha+1))^(1/(alpha+1))
h = hist(x, prob=T, breaks=40, plot=F)
plot(h$count, log="xy", type='l', lwd=1, lend=2,
xlab="", ylab="", main="Density in logarithmic scale")
По логике я нашел и нашел несколько подходящих сообщений, таких как this и this , и пробовал такие вещи, как
x1 = 5 # Maximum value
x0 = 0.1 # It can't be zero; otherwise X^0^(neg) is 1/0.
alpha = -2.5 # It has to be negative.
y = runif(1e5) # Number of samples
x = ((x1^(alpha+1) - x0^(alpha+1))*y + x0^(alpha+1))^(1/(alpha+1))
h = hist(x, prob=T, breaks=40, plot=F)
plot(h$count, log="xy", type='l', lwd=1, lend=2,
xlab="", ylab="", xaxt = 'n', main="Density in logarithmic scale")
axis(side=1, at=seq(0, 5, .2), labels=seq(0, 5, .2))
получение бессмысленных результатов: