Jak nałożyć krzywą normalną na histogram w r (2 przykłady)

Często możesz chcieć nałożyć krzywą normalną na histogram w R.

Poniższe przykłady pokazują, jak to zrobić w podstawowym R i ggplot2 .

Przykład 1: nałóż krzywą normalną na bazowy histogram R

Możemy użyć poniższego kodu, aby utworzyć histogram o podstawie R i nałożyć na histogram krzywą normalną:

 #make this example reproducible
set. seed ( 0 )

#define data
data <-rnorm( 1000 )

#create histogram
hist_data <- hist(data)

#define x and y values to use for normal curve
x_values <- seq(min(data), max(data), length = 100 )
y_values <- dnorm(x_values, mean = mean(data), sd = sd(data)) 
y_values <- y_values * diff(hist_data$mids[1:2]) * length(data) 

#overlay normal curve on histogram
lines(x_values, y_values, lwd = 2 )

Czarna krzywa na wykresie przedstawia krzywą normalną.

Możesz użyć argumentów col , lwd i lty , aby zmodyfikować odpowiednio kolor, szerokość i rodzaj linii:

 #overlay normal curve with custom aesthetics
lines(x_values, y_values, col=' red ', lwd= 5 , lty=' dashed ')

Przykład 2: Nałóż krzywą normalną na histogram w ggplot2

Możemy użyć następującego kodu, aby utworzyć histogram w ggplot2 i nałożyć na histogram krzywą normalną:

 library (ggplot2)

#make this example reproducible
set. seed ( 0 )

#define data
data <- data. frame (x=rnorm( 1000 ))

#create histogram and overlay normal curve
ggplot(data, aes(x)) +
  geom_histogram(aes(y = ..density..), fill=' lightgray ', col=' black ') +
  stat_function(fun = dnorm, args = list(mean=mean(data$x), sd=sd(data$x)))

Czarna krzywa na wykresie przedstawia krzywą normalną.

Możesz użyć argumentów col , lwd i lty , aby zmodyfikować odpowiednio kolor, szerokość i rodzaj linii:

 #overlay normal curve with custom aesthetics
ggplot(data, aes(x)) +
  geom_histogram(aes(y = ..density..), fill=' lightgray ', col=' black ') +
  stat_function(fun = dnorm, args = list(mean=mean(data$x), sd=sd(data$x)),
                col=' red ', lwd= 2 , lty=' dashed '))

Uwaga : Pełną dokumentację funkcji stat_function znajdziesztutaj .

Dodatkowe zasoby

Poniższe samouczki wyjaśniają, jak wykonywać inne typowe operacje w języku R:

Jak utworzyć histogram częstotliwości względnej w R
Jak określić przerwy w histogramie w R