``` r
simpan=read.table("data4.1.csv",header=TRUE, sep=",") #membaca data4.1
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 12 30 5
## 13 13 40 5
## 14 14 40 6
## 15 15 40 6
## 16 16 50 NA
## 17 17 50 NA
## 18 18 50 NA
## 19 NA 50 NA
data_1=simpan$data1 #data_1 menyimpan data1
data_11 = na.omit(data_1) #data_11 menyimpan data1, tanpa NA
data_2=simpan$data2 #data_2 menyimpan data2
data_21=na.omit(data_2) #data_21 menyimpan data2, tanpa NA
data_3=simpan$data3 #data_3 menyimpan data3
data_31=na.omit(data_3) #data_31 menyimpan data3, tanpa NA
data_1 #menampilkan data1 atau data pada Tabel 4.1
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 NA
data_11 #menampilkan data1 atau data pada Tabel 4.1, dengan menghapus NA
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
## attr(,"na.action")
## [1] 19
## attr(,"class")
## [1] "omit"
data_2 #menampilkan data2 atau data pada Tabel 4.3
## [1] 10 10 10 20 20 30 30 30 30 30 30 30 40 40 40 50 50 50 50
data_21 #menampilkan data2 atau data pada Tabel 4.3, dengan menghapus NA
## [1] 10 10 10 20 20 30 30 30 30 30 30 30 40 40 40 50 50 50 50
data_3 #menampilkan data3 atau data pada Tabel 4.9
## [1] 1 2 2 3 3 3 3 3 4 4 4 5 5 6 6 NA NA NA NA
data_31 #menampilkan data3 atau data pada Tabel 4.9, dengan menghapus NA
## [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"
library(psych) #mengaktifkan package "psych" yang di dalamnya menyediakan fungsi describe()
## Warning: package 'psych' was built under R version 4.5.3
describe(data_11) #menyajikan ukuran gejala pusat, letak, pencaran, kemiringan, dan kurtosis untuk data1
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 18 9.5 5.34 9.5 9.5 6.67 1 18 17 0 -1.4 1.26
describe(data_21) #menyajikan ukuran gejala pusat, letak, pencaran, kemiringan, dan kurtosis untuk data2
## 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) #menyajikan ukuran gejala pusat, letak, pencaran, kemiringan, dan kurtosis untuk data3
## 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
# ============================================================
# APLIKASI DALAM R (DATA BERKELOMPOK)
# ============================================================
simpan=read.table("data4.2.csv",header=TRUE, sep=",") #membaca data4.2
simpan
## jurusan.angka jurusan.label jenis.angka jenis.label IQ EQ
## 1 1 matematika 1 laki-laki 101 102
## 2 1 matematika 1 laki-laki 104 104
## 3 1 matematika 1 laki-laki 102 104
## 4 1 matematika 2 perempuan 111 105
## 5 1 matematika 2 perempuan 101 106
## 6 1 matematika 2 perempuan 105 102
## 7 2 statistika 1 laki-laki 104 101
## 8 2 statistika 1 laki-laki 103 112
## 9 2 statistika 1 laki-laki 104 104
## 10 2 statistika 2 perempuan 102 101
## 11 2 statistika 2 perempuan 101 102
## 12 2 statistika 2 perempuan 101 104
library(psych)
describeBy(simpan$IQ, simpan$jurusan.angka)
##
## Descriptive statistics by group
## group: 1
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 6 104 3.79 103 104 2.97 101 111 10 0.86 -0.93 1.55
## ------------------------------------------------------------
## group: 2
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 6 102.5 1.38 102.5 102.5 2.22 101 104 3 0 -2.06 0.56
describeBy(simpan$IQ, simpan$jurusan.label)
##
## Descriptive statistics by group
## group: matematika
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 6 104 3.79 103 104 2.97 101 111 10 0.86 -0.93 1.55
## ------------------------------------------------------------
## group: statistika
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 6 102.5 1.38 102.5 102.5 2.22 101 104 3 0 -2.06 0.56
describeBy(simpan$EQ, simpan$jurusan.angka)
##
## Descriptive statistics by group
## group: 1
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 6 103.83 1.6 104 103.83 2.22 102 106 4 -0.02 -1.82 0.65
## ------------------------------------------------------------
## group: 2
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 6 104 4.15 103 104 2.22 101 112 11 1.05 -0.59 1.69
describeBy(simpan$EQ, simpan$jurusan.label)
##
## Descriptive statistics by group
## group: matematika
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 6 103.83 1.6 104 103.83 2.22 102 106 4 -0.02 -1.82 0.65
## ------------------------------------------------------------
## group: statistika
## vars n mean sd median trimmed mad min max range skew kurtosis se
## X1 1 6 104 4.15 103 104 2.22 101 112 11 1.05 -0.59 1.69
library(doBy)
## Warning: package 'doBy' was built under R version 4.5.3
summaryBy(IQ + EQ ~ jurusan.label + jenis.label, data = simpan,
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 matematika laki-laki 102.3333 1.5275252 307
## 2 matematika perempuan 105.6667 5.0332230 317
## 3 statistika laki-laki 103.6667 0.5773503 311
## 4 statistika perempuan 101.3333 0.5773503 304
## EQ.ratarata EQ.standar_deviasi EQ.jumlah
## 1 103.3333 1.154701 310
## 2 104.3333 2.081666 313
## 3 105.6667 5.686241 317
## 4 102.3333 1.527525 307