Забегая наперед хочу сказать что данное задание учебное, поэтому не обессудьте. Стоит задача построить график нормального распределения. На этом графике должно быть около 100000 попыток замеров и сам график. Попытки должны отображаться в виде гистограммы. Соответственно будет очень много столбиков, поверх которых должна идти линия распределения. Вот нашел график без гистограммы:
С виду то что нужно, но не хватает еще гистограммы. Далее нашел формулы плотности и функции распределения:
Вот мой код для отрисовки:
xseq<-seq(-4,4,.01)
y<-2*xseq + rnorm(length(xseq),0,5.5)
hist(y, prob=TRUE, ylim=c(0,.06), breaks=20)
curve(dnorm(x, mean(y), sd(y)), add=TRUE, col="darkblue", lwd=2)
par(new=TRUE)
plot(xseq, cumulative , col="red",type="l",lwd=2, cex=2, cex.axis=.8)
И вот тот ужас который получился:
И соответственно мои проблемы:
- Как добавить больше столбиков, потому что сейчас их больше 20 где-то, а нужно 100000.
- Как починить маркировку осей графика, а то там что-то не очень понятное происходит.
И главный вопрос- я вообще правильно выполняю поставленную задачу, или есть способ получше?
break=20
, судя по всему. Но 100000 столбцов будет полная ерунда. Нужно разделять на интервалы и считать количество исходов в каждом из них.