Membaca dan Menyiapkan Data
# Membaca data4.1
simpan = read.table("data4.1.csv", header=TRUE, sep=",")
simpan
## data1 data2 data3
## 1 1 10 1
## 2 2 10 2
## 3 3 10 2
## 4 4 20 3
## 5 5 20 3
## 6 6 30 3
## 7 7 30 3
## 8 8 30 3
## 9 9 30 4
## 10 10 30 4
## 11 11 30 4
## 12 11 30 5
## 13 12 40 5
## 14 13 40 6
## 15 14 40 6
## 16 15 50 NA
## 17 16 50 NA
## 18 17 50 NA
## 19 18 50 NA
# Memisahkan data dan menghilangkan nilai NA (Missing Values)
data_1 = simpan$data1
data_11 = na.omit(data_1)
data_2 = simpan$data2
data_21 = na.omit(data_2)
data_3 = simpan$data3
data_31 = na.omit(data_3)
# Menampilkan data asli dan data tanpa NA
data_1
## [1] 1 2 3 4 5 6 7 8 9 10 11 11 12 13 14 15 16 17 18
data_11
## [1] 1 2 3 4 5 6 7 8 9 10 11 11 12 13 14 15 16 17 18
data_2
## [1] 10 10 10 20 20 30 30 30 30 30 30 30 40 40 40 50 50 50 50
data_21
## [1] 10 10 10 20 20 30 30 30 30 30 30 30 40 40 40 50 50 50 50
data_3
## [1] 1 2 2 3 3 3 3 3 4 4 4 5 5 6 6 NA NA NA NA
data_31
## [1] 1 2 2 3 3 3 3 3 4 4 4 5 5 6 6
## attr(,"na.action")
## [1] 16 17 18 19
## attr(,"class")
## [1] "omit"
Menggunakan Package ‘psych’ untuk Ukuran Gejala Pusat dan
Pencaran
library(psych)
## Warning: package 'psych' was built under R version 4.6.1
# Menyajikan ukuran gejala pusat, letak, pencaran, kemiringan, dan kurtosis
describe(data_11)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 19 9.58 5.2 10 9.59 5.93 1 18 17 -0.04 -1.32 1.19
describe(data_21)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 19 31.58 13.44 30 31.76 14.83 10 50 40 -0.14 -1.13 3.08
describe(data_31)
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 15 3.6 1.45 3 3.62 1.48 1 6 5 0.14 -1.01 0.38
Menggunakan Package ‘pastecs’ untuk Statistik Tambahan
library(pastecs)
## Warning: package 'pastecs' was built under R version 4.6.1
# Menyajikan nilai statistik tambahan seperti variance, sum, SE mean, dan koefisien variasi
stat.desc(data_11)
## nbr.val nbr.null nbr.na min max range
## 19.0000000 0.0000000 0.0000000 1.0000000 18.0000000 17.0000000
## sum median mean SE.mean CI.mean.0.95 var
## 182.0000000 10.0000000 9.5789474 1.1928535 2.5060921 27.0350877
## std.dev coef.var
## 5.1995276 0.5428078
Analisis Data Berkelompok (Grouped Data)
# Membaca data4.2
simpan_kelompok = read.table("data4.2.csv", header=TRUE, sep=",")
simpan_kelompok
## jurusan.angka jurusan.label jenis.angka jenis.label IQ EQ
## 1 1 Matematika 1 Laki-Laki 110 105
## 2 1 Matematika 2 Perempuan 115 110
## 3 1 Matematika 1 Laki-Laki 112 108
## 4 1 Matematika 2 Perempuan 118 112
## 5 2 Biologi 1 Laki-Laki 105 100
## 6 2 Biologi 2 Perempuan 108 102
## 7 2 Biologi 1 Laki-Laki 106 101
## 8 2 Biologi 2 Perempuan 110 105
library(psych)
# Menyajikan nilai statistik berdasarkan kelompok variabel tertentu
describeBy(simpan_kelompok$IQ, simpan_kelompok$jurusan.angka)
##
## Descriptive statistics by group
## group: 1
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 4 113.75 3.5 113.5 113.75 3.71 110 118 8 0.12 -2.11 1.75
## ------------------------------------------------------------
## group: 2
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 4 107.25 2.22 107 107.25 2.22 105 110 5 0.18 -2.11 1.11
describeBy(simpan_kelompok$IQ, simpan_kelompok$jurusan.label)
##
## Descriptive statistics by group
## group: Biologi
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 4 107.25 2.22 107 107.25 2.22 105 110 5 0.18 -2.11 1.11
## ------------------------------------------------------------
## group: Matematika
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 4 113.75 3.5 113.5 113.75 3.71 110 118 8 0.12 -2.11 1.75
describeBy(simpan_kelompok$EQ, simpan_kelompok$jurusan.angka)
##
## Descriptive statistics by group
## group: 1
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 4 108.75 2.99 109 108.75 2.97 105 112 7 -0.16 -2.02 1.49
## ------------------------------------------------------------
## group: 2
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 4 102 2.16 101.5 102 1.48 100 105 5 0.45 -1.88 1.08
describeBy(simpan_kelompok$EQ, simpan_kelompok$jurusan.label)
##
## Descriptive statistics by group
## group: Biologi
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 4 102 2.16 101.5 102 1.48 100 105 5 0.45 -1.88 1.08
## ------------------------------------------------------------
## group: Matematika
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 4 108.75 2.99 109 108.75 2.97 105 112 7 -0.16 -2.02 1.49
Menggunakan Package ‘doBy’ untuk Ringkasan Multi-Variabel
library(doBy)
## Warning: package 'doBy' was built under R version 4.6.1
# Menyajikan ringkasan statistik (rata-rata, standar deviasi, jumlah) untuk kombinasi kategori
summaryBy(IQ + EQ ~ jurusan.label + jenis.label, data = simpan_kelompok,
FUN = function(x) c(ratarata = mean(x), standar_deviasi = sd(x), jumlah = sum(x)))
## jurusan.label jenis.label IQ.ratarata IQ.standar_deviasi IQ.jumlah
## 1 Biologi Laki-Laki 105.5 0.7071068 211
## 2 Biologi Perempuan 109.0 1.4142136 218
## 3 Matematika Laki-Laki 111.0 1.4142136 222
## 4 Matematika Perempuan 116.5 2.1213203 233
## EQ.ratarata EQ.standar_deviasi EQ.jumlah
## 1 100.5 0.7071068 201
## 2 103.5 2.1213203 207
## 3 106.5 2.1213203 213
## 4 111.0 1.4142136 222