График распределения значений вероятности - PullRequest
0 голосов
/ 29 сентября 2011

у меня

значения вероятности: 0.06,0.06,0.1,0.08,0.12,0.16,0.14,0.14,0.08,0.02,0.04, суммирование до 1

соответствующие интервалы, в которых стохастическая переменная может принимать свое значение с соответствующей вероятностью из приведенного выше списка:

126,162,233,304,375,446,517,588,659,730,801,839

Как я могу построить распределение вероятностей?

На оси x значения интервалов, между гистограммами интервалов со значением вероятности?

Спасибо.

1 Ответ

2 голосов
/ 29 сентября 2011

Как насчет

x <- c(126,162,233,304,375,446,517,588,659,730,801,839)
p <- c(0.06,0.06,0.1,0.08,0.12,0.16,0.14,0.14,0.08,0.02,0.04)
plot(x,c(p,0),type="s")
lines(x,c(0,p),type="S")
rect(x[-1],0,x[-length(x)],p,col="lightblue")

для быстрого ответа?(С включенным rect вам может не понадобиться вызов lines, и вы могли бы изменить его на plot(x,p,type="n"). Как обычно, я бы порекомендовал par(bty="l",lty=1) для моих предпочтительных графических значений по умолчанию ...)

(Объяснение: "s" и "S" - это два разных типа ступеней (см. Подробности в ?plot): я использовал их оба, чтобы получить как левую, так и правую границы распределения.)

edit : В ваших комментариях вы говорите "(это) не похоже на гистограмму".Не совсем понятно, чего вы хотите.Я добавил прямоугольники в приведенном выше примере - может быть, это так?Или вы могли бы сделать

b <- barplot(p,width=diff(x),space=0)

, но получить правильные метки оси X - это боль.

...