Statistika Deskriptif dengan R

#install.packages("dplyr")
library(dplyr)

## Warning: package 'dplyr' was built under R version 4.3.3

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

Persiapan Data

#Nama Mahasiswa
nama <- c("Chevroline", "Lutfiatin", "Dean", "Nathalia", "Futhy", "Najwa", "Lia", "Oca", "Adyun", "Iyyah")

#Berat Badan Mahasiswa
bb <- c(40, 56, 64, 49, 47, 62, 69, 57, 50, 53)

#Tinggi Badan Mahasiswa
tb <- c(157, 174, 179, 165, 172, 167, 158, 153, 178, 177)

kelas <- data.frame(nama, bb, tb)

Ukuran Pemusatan Data

mean_bb <- mean(bb)
cat("Rata-rata berat badan mahasiswa:", mean_bb, "kg")

## Rata-rata berat badan mahasiswa: 54.7 kg

median_bb <- median(bb)
cat("Median berat badan mahasiswa:", median_bb, "kg")

## Median berat badan mahasiswa: 54.5 kg

modus <- function(x) {
  ux <- unique(x)
  ux[which.max(tabulate(match(x, ux)))]
}

modus_bb <- modus(bb)
cat("Modus berat badan mahasiswa:", modus_bb, "kg")

## Modus berat badan mahasiswa: 40 kg

mean_tb <- mean(tb)
cat("Rata-rata tinggi badan mahasiswa:", mean_tb, "cm")

## Rata-rata tinggi badan mahasiswa: 168 cm

median_tb <- median(tb)
cat("Median tinggi badan mahasiswa:", median_tb, "cm")

## Median tinggi badan mahasiswa: 169.5 cm

modus <- function(x) {
  ux <- unique(x)
  ux[which.max(tabulate(match(x, ux)))]
}

modus_tb <- modus(tb)
cat("Modus tinggi badan mahasiswa:", modus_tb, "cm")

## Modus tinggi badan mahasiswa: 157 cm

Ukuran Penyebaran Data

range_bb <- range(bb)
cat("Range berat badan mahasiswa:", range_bb, "\n")

## Range berat badan mahasiswa: 40 69

variance_bb <- var(bb)
cat("Variansi berat badan mahasiswa:", variance_bb, "\n")

## Variansi berat badan mahasiswa: 76.01111

sd_bb <- sd(bb)
cat("Simpangan baku berat badan mahasiswa:", sd_bb, "\n")

## Simpangan baku berat badan mahasiswa: 8.718435

Visualisasi Data

hist(kelas$bb, col = "pink",
     main = "Histogram Berat Badan Mahasiswa",
     xlab = "Berat badan (kg)")

hist(kelas$tb, col = "skyblue",
     main = "Histogram Tinggi Badan Mahasiswa",
     xlab = "Tinggi badan (cm)")

boxplot(kelas$bb, col = "red",
        main = "Boxplot Berat Badan Mahasiswa",
        horizontal = TRUE)

boxplot(kelas$tb, col = "purple",
        main = "Boxplot Tinggi Badan Mahasiswa",
        horizontal = TRUE)