R: Наложение распределения Пуассона по гистограмме данных - PullRequest
0 голосов
/ 09 января 2019

У меня есть некоторые дискретные данные, которые я нанес на гистограмму. Я хотел бы наложить распределение Пуассона, чтобы показать, что данные примерно распределены по Пуассону. Представьте, что два графика из приведенного ниже кода объединяются в один график, и это то, чего я хотел бы достичь.

# Read data
data <- read.csv("data.csv")

# Plot data
hist(data, prob=TRUE)

# Plot Poisson
c <- c(0:7)
plot(c, dpois(c, mean(data)), type="l")

Я попробовал функцию кривой:

curve(c, dpois(x=c, lambda=mean(data)), add=T)

Но все, что я получаю, это:

Кривая Пуассона, похоже, внезапно останавливается, но я ожидаю, что она будет следовать форме гистограммы.

Attempt at overlaying Poisson distribution

Мне бы хотелось, чтобы это выглядело так (необязательно с цветами или несколькими наборами данных): Correctly overlayed Poisson distribution

1 Ответ

0 голосов
/ 09 января 2019

Код ниже делает то, что вы хотите.

set.seed(12111978)
vec <- rpois(50, 3)
hist(vec, prob=TRUE, ylim = c(0, .25)) # may need to tweak the y axis.
lines(0:max(vec), dpois(0:max(vec), mean(vec)), col = 'red')

Poisson distribution represented with lines

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