simpan=read.table("data4.2.csv.", header=TRUE, sep=",")
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)
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