#KELOMPOK 2#
#1. Maulida Syifa (2502056010), 2. Wildia Nuratna (2502056045), 3. Destya Natalyne Obehetan (2502056022), 4. Kanaya Shafa Az Zahra (2502056023)
#Program Studi : Ilmu Komunikasi A (Angkatan 2025)
#Mata Kuliah : Statistika

#Ukuran Gejala Pusat, Letak,Pencaran, Kemiringan dan Keruncingan#

#Aplikasi dalam R#
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
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"
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
##    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) #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
library(pastecs)
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)
##                       x
## nbr.val      15.0000000
## nbr.null      0.0000000
## nbr.na        0.0000000
## min           1.0000000
## max           6.0000000
## range         5.0000000
## sum          54.0000000
## median        3.0000000
## mean          3.6000000
## SE.mean       0.3754363
## CI.mean.0.95  0.8052307
## var           2.1142857
## std.dev       1.4540584
## coef.var      0.4039051
#Aplikasi dalam R (Data Berkelompok)#
simpan = read.table ("data4.2.csv", header = TRUE, sep = ",") 
simpan
##    jurusan.angka jurusan.label jenis.angka  IQ  EQ
## 1              1             1   laki-laki 101 102
## 2              1             1   laki-laki 104 104
## 3              1             1   laki-laki 102 104
## 4              1             2   perempuan 111 105
## 5              1             2   perempuan 101 106
## 6              1             2   perempuan 105 102
## 7              2             1   laki-laki 104 101
## 8              2             1   laki-laki 103 112
## 9              2             1   laki-laki 104 104
## 10             2             2   perempuan 102 101
## 11             2             2   perempuan 101 102
## 12             2             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: 1
##    vars n mean   sd median trimmed  mad min max range  skew kurtosis   se
## X1    1 6  103 1.26  103.5     103 0.74 101 104     3 -0.49     -1.7 0.52
## ------------------------------------------------------------ 
## group: 2
##    vars n  mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 6 103.5 3.99  101.5   103.5 0.74 101 111    10 0.99    -0.83 1.63
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: 1
##    vars n  mean   sd median trimmed  mad min max range skew kurtosis   se
## X1    1 6 104.5 3.89    104   104.5 1.48 101 112    11 1.03    -0.55 1.59
## ------------------------------------------------------------ 
## group: 2
##    vars n   mean   sd median trimmed  mad min max range skew kurtosis  se
## X1    1 6 103.33 1.97    103  103.33 2.22 101 106     5 0.14    -1.95 0.8
library(doBy)
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 IQ.ratarata IQ.standar_deviasi IQ.jumlah EQ.ratarata
## 1             1       103.0           1.264911       618    104.5000
## 2             2       103.5           3.987480       621    103.3333
##   EQ.standar_deviasi EQ.jumlah
## 1           3.885872       627
## 2           1.966384       620