У меня есть программа, которая имитирует честные кости 100 раз.
При использовании этой программы мне нужно иметь кривую колокольчика над гистограммой, показывающую нормальное / гауссовское распределение.
В моем коде переменная, представляющая результат эксперимента - это z.
Определение кривой колокола содержит два параметра: ave определяет, где на оси z кривая должна иметь максимум, а sig определяет ширину пика.
set.seed(123)
x <- sample(1:6, size=100, replace = TRUE)
hist(x,
main="100 fair rolls",
xlab = "Dice Result",
ylab = "Probability",
xlim=c(0.5,6.5),
breaks=-1:100+.5,
prob=TRUE )
Уравнение распределения Гаусса у меня пока:
bell<-function(z,ave,sig){
exp(-(z-ave)^2/(2*sig))/sqrt(2*pi*sig)
}
points<-0:600
lines(points,bell(points,350,292))