Sekenario

Y: Status kelulusan tepat waktu (0: Tidak tepat waktu, 1: Tepat waktu)

X1: IPK mahasiswa (0–4) X2: Jumlah SKS yang diambil per semester X3: Status kerja mahasiswa (0: Tidak bekerja, 1: Bekerja) X4: Jam belajar per minggu

library(dplyr)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
n <- 100

x1 <- round(rnorm(n,3,0.4),2)

x1
##   [1] 3.19 3.05 2.67 2.73 3.46 2.93 3.28 3.53 2.69 2.29 3.41 3.32 2.09 3.12 3.70
##  [16] 3.18 2.74 2.84 2.20 2.95 2.82 2.75 2.63 3.05 2.70 3.25 2.77 3.11 3.16 2.57
##  [31] 3.14 3.10 3.25 3.08 3.34 3.11 2.36 2.69 3.75 2.44 3.08 3.01 2.52 3.10 3.23
##  [46] 3.50 3.13 3.07 2.78 2.95 3.34 3.37 2.57 3.30 3.07 2.75 3.57 3.28 3.62 3.25
##  [61] 3.17 3.05 2.89 2.82 2.38 2.25 2.51 2.91 3.19 2.58 3.04 3.31 3.24 2.72 2.85
##  [76] 2.78 2.97 2.90 2.89 3.26 2.48 2.93 2.93 2.88 3.65 3.07 3.53 3.32 2.53 2.77
##  [91] 2.56 3.33 2.53 2.14 2.87 3.40 3.64 2.91 3.23 3.02
x2 <- floor(rnorm(n,20,2))

x2
##   [1] 22 19 19 20 20 17 15 17 19 21 21 18 22 20 20 23 21 16 21 20 22 19 19 19 22
##  [26] 17 19 17 19 20 24 21 20 18 17 18 21 20 19 17 16 20 22 20 19 18 20 20 19 20
##  [51] 17 21 19 18 22 19 19 19 21 22 19 19 20 21 20 22 19 19 18 20 17 19 21 19 20
##  [76] 19 24 16 21 21 16 22 21 21 18 19 23 21 17 18 22 19 22 17 21 18 16 18 17 19
x3 <- round(runif(n))

x3
##   [1] 1 0 0 1 0 0 1 1 1 1 0 0 0 0 1 0 0 0 1 1 1 1 0 0 1 1 0 1 1 0 1 0 0 1 1 0 1
##  [38] 0 1 1 0 1 0 0 1 0 1 0 0 0 1 0 0 1 1 0 1 0 0 0 1 0 1 0 1 0 1 0 1 1 1 1 1 0
##  [75] 1 0 1 0 0 0 0 0 1 1 1 0 1 0 0 0 0 0 0 1 1 1 1 1 1 1
x4 <- round(rnorm(n,18,5))

x4
##   [1] 20 19 10 20 20 27 21 14 16 22 21 13 17 18 14 28 21 19 18 24 16 20 24 15 18
##  [26] 22 24 21 15 15 17 10 16 17 27 12 18 18 13 18 17  8 17 12 26 19 19 14 13 16
##  [51] 20 16 22 22 21 10 17 19 19 13  8 20 13 15 21 15 10 20 22 21 19 14 10 18 29
##  [76] 22 20 18 19 15 15 18 14 10 17 14 12 24 13 18 23 22 18 25 13  5 13 20 14 20
y <- rbinom(n, 1, 0.6)

y
##   [1] 0 1 0 0 1 0 0 1 0 1 0 1 1 0 0 0 1 0 1 1 0 1 1 1 1 1 1 1 0 1 1 1 1 0 0 0 1
##  [38] 1 1 0 1 0 0 0 0 0 1 1 0 1 0 1 1 0 1 1 0 1 1 1 1 0 1 1 0 1 1 0 1 0 1 0 0 1
##  [75] 1 1 1 1 1 0 1 1 1 0 0 0 0 1 0 1 0 0 0 1 1 1 1 1 1 0