==============================================================================

1. ANALISIS UNIVARIAT

==============================================================================

simpan = read.table("4.1komstat.csv", header = TRUE, sep=";")
simpan
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 11 12 13 14 15 16 17 18
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 11 12 13 14 15 16 17 18
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"
suppressWarnings(library(psych)) #mengaktifkan package "psych" yang di dalamnya menyediakan fungsi describe()
describe(data_11) #menyajikan ukuran gejala pusat, letak, pencaran, kemiringan, dan kurtosis untuk data1
describe(data_21) #menyajikan ukuran gejala pusat, letak, pencaran, kemiringan, dan kurtosis untuk data2
describe(data_31) #menyajikan ukuran gejala pusat, letak, pencaran, kemiringan, dan kurtosis untuk data3
library(pastecs)
## Warning: package 'pastecs' was built under R version 4.5.3
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
stat.desc(data_21)
##      nbr.val     nbr.null       nbr.na          min          max        range 
##   19.0000000    0.0000000    0.0000000   10.0000000   50.0000000   40.0000000 
##          sum       median         mean      SE.mean CI.mean.0.95          var 
##  600.0000000   30.0000000   31.5789474    3.0839291    6.4790946  180.7017544 
##      std.dev     coef.var 
##   13.4425353    0.4256803
stat.desc(data_31)

==============================================================================

2. ANALISIS BIVARIAT & MULTIVARIAT

==============================================================================

simpan=read.table("4.2komstat.csv",header=TRUE, sep = ";") #membaca data4.1
simpan
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)   ) } )