Трудности с установкой оси X гистограммы в R для представления истинной области случайной величины - PullRequest
0 голосов
/ 25 октября 2019

Я только что понял, что в этом посте, который я собрал , гистограмма имеет только значения индекса по оси x вместо предполагаемого диапазона от 0 до 5:

enter image description here

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))

получение бессмысленных результатов:

enter image description here

1 Ответ

1 голос
/ 25 октября 2019

Возможно, вы ищете density().

plot(density(x), log="y", ylab="log density x", col=2)
legend("topright", "x", lty=1, col=2)

enter image description here

Примечание: set.seed(42).

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...