Как я могу наложить произвольное параметрическое распределение на гистограмму, используя ggplot?

Этот вопрос был перенесен из Cross Validated, поскольку на него можно ответить в Stack Overflow. Перенесено7 лет назад.

Я предпринял попытку на основе примера Quick -R , но я не понимаю, откуда берется коэффициент масштабирования. Является ли этот метод разумным? Как я могу изменить его для использования ggplot?

Ниже приведен пример наложения нормального и логнормального распределения с использованием этого метода :

## Get a log-normalish data set: the number of characters per word in "Alice in Wonderland"
alice.raw <- readLines(con = "http://www.gutenberg.org/cache/epub/11/pg11.txt", 
                       n = -1L, ok = TRUE, warn = TRUE,
                       encoding = "UTF-8")

alice.long <- paste(alice.raw, collapse=" ")
alice.long.noboilerplate <- strsplit(alice.long, split="\\*\\*\\*")[[1]][3]
alice.words <- strsplit(alice.long.noboilerplate, "[[:space:]]+")[[1]]
alice.nchar <- nchar(alice.words)
alice.nchar <- alice.nchar[alice.nchar > 0]

# Now we want to plot both the histogram and then log-normal probability dist
require(MASS)
h <- hist(alice.nchar, breaks=1:50, xlab="Characters in word", main="Count")
xfit <- seq(1, 50, 0.1)

# Plot a normal curve
yfit<-dnorm(xfit,mean=mean(alice.nchar),sd=sd(alice.nchar))
yfit <- yfit * diff(h$mids[1:2]) * length(alice.nchar) 
lines(xfit, yfit, col="blue", lwd=2) 

# Now plot a log-normal curve
params <- fitdistr(alice.nchar, densfun="lognormal")
yfit <- dlnorm(xfit, meanlog=params$estimate[1], sdlog=params$estimate[1])
yfit <- yfit * diff(h$mids[1:2]) * length(alice.nchar)
lines(xfit, yfit, col="red", lwd=2)

. Это дает следующий график: Plot produced by the code above, showing a histogram of word length superimposed with a normal distribution curve and a log-normal distribution curve

Чтобы уточнить, я хотел бы иметь подсчет по оси y -, а не оценку плотности.

задан fmark 2 July 2012 в 00:12

0 ответов

Другие вопросы по тегам:

data-visualization ggplot2 histogram r

Как я могу наложить произвольное параметрическое распределение на гистограмму, используя ggplot?

0 ответов

Похожие вопросы: