Cara melapisi kurva normal pada histogram di r (2 contoh)

Seringkali Anda mungkin ingin melapisi kurva normal pada histogram di R.

Contoh berikut menunjukkan cara melakukan ini di basis R dan ggplot2 .

Contoh 1: menempatkan kurva normal pada histogram dasar R

Kita dapat menggunakan kode berikut untuk membuat histogram di basis R dan menempatkan kurva normal pada histogram:

 #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 )

Kurva hitam pada plot mewakili kurva normal.

Jangan ragu untuk menggunakan argumen col , lwd dan lty untuk mengubah warna, lebar, dan jenis garis:

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

Contoh 2: Letakkan kurva normal pada histogram di ggplot2

Kita dapat menggunakan kode berikut untuk membuat histogram di ggplot2 dan melapisi kurva normal pada histogram:

 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)))

Kurva hitam pada plot mewakili kurva normal.

Jangan ragu untuk menggunakan argumen col , lwd dan lty untuk mengubah warna, lebar, dan jenis garis:

 #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 '))

Catatan : Anda dapat menemukan dokumentasi lengkap untuk stat_functiondi sini .

Sumber daya tambahan

Tutorial berikut menjelaskan cara melakukan operasi umum lainnya di R:

Cara membuat histogram frekuensi relatif di R
Cara menentukan jeda histogram di R