NO 1

Sebuah warung es teh di pinggir jalan melakukan pengamatan dan mencatat frekuensi harian permintaan es teh selama 100 hari terakhir. Data pengamatan adalah sebagai berikut: Permintaan (gelas) (50,60,70,80,90) Frekuensi Observasi (10,20,40,20,10) a. prediksikan permintaan es teh selama 5 hari ke depan. b. prediksikan permintaan es teh selama 20 hari ke depan.

# Data permintaan es teh dan frekuensi observasi
data_es_teh <- data.frame(
  permintaan = c(50, 60, 70, 80, 90),
  frekuensi = c(10, 20, 40, 20, 10)
)
data_es_teh
##   permintaan frekuensi
## 1         50        10
## 2         60        20
## 3         70        40
## 4         80        20
## 5         90        10

Menghitung Probabilitas

# Menghitung probabilitas berdasarkan frekuensi
data_es_teh$probabilitas <- round(data_es_teh$frekuensi / sum(data_es_teh$frekuensi), 2)
data_es_teh
##   permintaan frekuensi probabilitas
## 1         50        10          0.1
## 2         60        20          0.2
## 3         70        40          0.4
## 4         80        20          0.2
## 5         90        10          0.1

Mengitung Probabilitas

# Menghitung probabilitas berdasarkan frekuensi
data_es_teh$probabilitas <- round(data_es_teh$frekuensi / sum(data_es_teh$frekuensi), 2)
data_es_teh
##   permintaan frekuensi probabilitas
## 1         50        10          0.1
## 2         60        20          0.2
## 3         70        40          0.4
## 4         80        20          0.2
## 5         90        10          0.1

Menghitung Probabilitas Kumulatif

# Menghitung probabilitas kumulatif
data_es_teh$probcum <- cumsum(data_es_teh$probabilitas)
data_es_teh
##   permintaan frekuensi probabilitas probcum
## 1         50        10          0.1     0.1
## 2         60        20          0.2     0.3
## 3         70        40          0.4     0.7
## 4         80        20          0.2     0.9
## 5         90        10          0.1     1.0

Menghitung nilai ekspektasi

# Menghitung nilai ekspektasi
expvalue_es_teh <- round(sum(data_es_teh$permintaan * data_es_teh$probabilitas), 2)
expvalue_es_teh
## [1] 70

Membuat fungsi simulasi

simulasi_es_teh <- function(a, permintaan, kumulatif){
  # Fungsi untuk mensimulasikan permintaan es teh selama 'a' hari
  simulasi <- data.frame(
    hari_ke = 1:a,
    rand_value = round(runif(a), 2)  # Menghasilkan angka acak antara 0 dan 1
  )
  simulasi$permintaan <- 0  # Menambahkan kolom permintaan yang diinisialisasi dengan 0
  
  for (i in 1:nrow(simulasi)){
    for (j in 1:length(kumulatif)){
      if (simulasi$rand_value[i] <= kumulatif[j]){  # Membandingkan angka acak dengan nilai kumulatif
        simulasi$permintaan[i] = permintaan[j]
        break
      }
    }
  }
  average <- mean(simulasi$permintaan)  # Menghitung rata-rata permintaan
  return(list(data = simulasi, average = average))  # Mengembalikan hasil simulasi dan rata-rata
}

Prediksi 5 hari

set.seed(123)
simulasi_5_hari <- simulasi_es_teh(5, data_es_teh$permintaan, data_es_teh$probcum)

print(simulasi_5_hari$data)  # Menampilkan hasil simulasi untuk 5 hari
##   hari_ke rand_value permintaan
## 1       1       0.29         60
## 2       2       0.79         80
## 3       3       0.41         70
## 4       4       0.88         80
## 5       5       0.94         90
print(simulasi_5_hari$average)  # Menampilkan rata-rata permintaan untuk 5 hari
## [1] 76

Prediksi 20hari kedepan

set.seed(123)
simulasi_20_hari <- simulasi_es_teh(20, data_es_teh$permintaan, data_es_teh$probcum)

print(simulasi_20_hari$data)  # Menampilkan hasil simulasi untuk 20 hari
##    hari_ke rand_value permintaan
## 1        1       0.29         60
## 2        2       0.79         80
## 3        3       0.41         70
## 4        4       0.88         80
## 5        5       0.94         90
## 6        6       0.05         50
## 7        7       0.53         70
## 8        8       0.89         80
## 9        9       0.55         70
## 10      10       0.46         70
## 11      11       0.96         90
## 12      12       0.45         70
## 13      13       0.68         70
## 14      14       0.57         70
## 15      15       0.10         50
## 16      16       0.90         80
## 17      17       0.25         60
## 18      18       0.04         50
## 19      19       0.33         70
## 20      20       0.95         90
print(simulasi_20_hari$average)  # Menampilkan rata-rata permintaan untuk 20 hari
## [1] 71

NO 2

Bangkitkanlah data dengan distribusi eksponensial untuk variabel permintaan, dan data dengan distribusi normal untuk variabel frekuensi a. prediksikan permintaan selama 5 hari ke depan b. prediksikan permintaan selama 20 hari ke depan c. prediksikan permintaan selama 100 hari ke depan d. prediksikan permintaan selama 1000 hari ke depan ## Membangkitkan sata vaiabel permintaan dan freakuensi

set.seed(123)
data <- data.frame( 
  permintaan = round(rexp(10,0.01234)),
  frekuensi = ceiling(rnorm(10,10,5))
  )
data
##    permintaan frekuensi
## 1          68         8
## 2          47        17
## 3         108        12
## 4           3        13
## 5           5        11
## 6          26         8
## 7          25        19
## 8          12        13
## 9         221         1
## 10          2        14

Menghitung Probabilitas

data$probabilitas <- round(data$frekuensi / sum(data$frekuensi),2)
data
##    permintaan frekuensi probabilitas
## 1          68         8         0.07
## 2          47        17         0.15
## 3         108        12         0.10
## 4           3        13         0.11
## 5           5        11         0.09
## 6          26         8         0.07
## 7          25        19         0.16
## 8          12        13         0.11
## 9         221         1         0.01
## 10          2        14         0.12

Menghitung Probabilitas Kumulatif

data$probcum <- cumsum(data$probabilitas)
data
##    permintaan frekuensi probabilitas probcum
## 1          68         8         0.07    0.07
## 2          47        17         0.15    0.22
## 3         108        12         0.10    0.32
## 4           3        13         0.11    0.43
## 5           5        11         0.09    0.52
## 6          26         8         0.07    0.59
## 7          25        19         0.16    0.75
## 8          12        13         0.11    0.86
## 9         221         1         0.01    0.87
## 10          2        14         0.12    0.99

Menghitung nilai ekspektasi

expvalue <- round(sum(data$permintaan*data$probabilitas),2)
expvalue
## [1] 32.98

Membuat fungsi simulasi

simulasi <- function(a, permintaan, kumulatif){
  simulasi <- data.frame(
  hari_ke = 1:a,
  rand_value = round(runif(a),2)
  )
  simulasi$permintaan = 0
  
  for (i in 1:nrow(simulasi)){
    for (j in 1:length(kumulatif)){
      if (simulasi$rand_value[i] <= kumulatif[j]){
        simulasi$permintaan[i] = permintaan[j]
        break
      }
    }
  }
  average <- mean(simulasi$permintaan)
  return(list(data = simulasi, average = average))

}

Prediksi Permintaan 5 hari kedepan

set.seed(123)
simulasi_a <- simulasi(5, data$permintaan, data$probcum)

print(simulasi_a$data)
##   hari_ke rand_value permintaan
## 1       1       0.29        108
## 2       2       0.79         12
## 3       3       0.41          3
## 4       4       0.88          2
## 5       5       0.94          2
print(simulasi_a$average)
## [1] 25.4

Prediksi permintaan 20 hari kedepan

set.seed(123)
simulasi_b <- simulasi(20, data$permintaan, data$probcum)

print(simulasi_b$data)
##    hari_ke rand_value permintaan
## 1        1       0.29        108
## 2        2       0.79         12
## 3        3       0.41          3
## 4        4       0.88          2
## 5        5       0.94          2
## 6        6       0.05         68
## 7        7       0.53         26
## 8        8       0.89          2
## 9        9       0.55         26
## 10      10       0.46          5
## 11      11       0.96          2
## 12      12       0.45          5
## 13      13       0.68         25
## 14      14       0.57         26
## 15      15       0.10         47
## 16      16       0.90          2
## 17      17       0.25        108
## 18      18       0.04         68
## 19      19       0.33          3
## 20      20       0.95          2
print(simulasi_b$average)
## [1] 27.1

Prediksi permintaan 100 hari kedepan

set.seed(123)
simulasi_c <- simulasi(100, data$permintaan, data$probcum)

print(simulasi_c$data)
##     hari_ke rand_value permintaan
## 1         1       0.29        108
## 2         2       0.79         12
## 3         3       0.41          3
## 4         4       0.88          2
## 5         5       0.94          2
## 6         6       0.05         68
## 7         7       0.53         26
## 8         8       0.89          2
## 9         9       0.55         26
## 10       10       0.46          5
## 11       11       0.96          2
## 12       12       0.45          5
## 13       13       0.68         25
## 14       14       0.57         26
## 15       15       0.10         47
## 16       16       0.90          2
## 17       17       0.25        108
## 18       18       0.04         68
## 19       19       0.33          3
## 20       20       0.95          2
## 21       21       0.89          2
## 22       22       0.69         25
## 23       23       0.64         25
## 24       24       0.99          2
## 25       25       0.66         25
## 26       26       0.71         25
## 27       27       0.54         26
## 28       28       0.59         26
## 29       29       0.29        108
## 30       30       0.15         47
## 31       31       0.96          2
## 32       32       0.90          2
## 33       33       0.69         25
## 34       34       0.80         12
## 35       35       0.02         68
## 36       36       0.48          5
## 37       37       0.76         12
## 38       38       0.22         47
## 39       39       0.32        108
## 40       40       0.23        108
## 41       41       0.14         47
## 42       42       0.41          3
## 43       43       0.41          3
## 44       44       0.37          3
## 45       45       0.15         47
## 46       46       0.14         47
## 47       47       0.23        108
## 48       48       0.47          5
## 49       49       0.27        108
## 50       50       0.86         12
## 51       51       0.05         68
## 52       52       0.44          5
## 53       53       0.80         12
## 54       54       0.12         47
## 55       55       0.56         26
## 56       56       0.21         47
## 57       57       0.13         47
## 58       58       0.75         25
## 59       59       0.90          2
## 60       60       0.37          3
## 61       61       0.67         25
## 62       62       0.09         47
## 63       63       0.38          3
## 64       64       0.27        108
## 65       65       0.81         12
## 66       66       0.45          5
## 67       67       0.81         12
## 68       68       0.81         12
## 69       69       0.79         12
## 70       70       0.44          5
## 71       71       0.75         25
## 72       72       0.63         25
## 73       73       0.71         25
## 74       74       0.00         68
## 75       75       0.48          5
## 76       76       0.22         47
## 77       77       0.38          3
## 78       78       0.61         25
## 79       79       0.35          3
## 80       80       0.11         47
## 81       81       0.24        108
## 82       82       0.67         25
## 83       83       0.42          3
## 84       84       0.79         12
## 85       85       0.10         47
## 86       86       0.43          3
## 87       87       0.98          2
## 88       88       0.89          2
## 89       89       0.89          2
## 90       90       0.18         47
## 91       91       0.13         47
## 92       92       0.65         25
## 93       93       0.34          3
## 94       94       0.66         25
## 95       95       0.32        108
## 96       96       0.19         47
## 97       97       0.78         12
## 98       98       0.09         47
## 99       99       0.47          5
## 100     100       0.51          5
print(simulasi_c$average)
## [1] 29.96

Prediksi permintaan 1000 hari kedepan

set.seed(123)
simulasi_d <- simulasi(1000, data$permintaan, data$probcum)

print(simulasi_d$data)
##      hari_ke rand_value permintaan
## 1          1       0.29        108
## 2          2       0.79         12
## 3          3       0.41          3
## 4          4       0.88          2
## 5          5       0.94          2
## 6          6       0.05         68
## 7          7       0.53         26
## 8          8       0.89          2
## 9          9       0.55         26
## 10        10       0.46          5
## 11        11       0.96          2
## 12        12       0.45          5
## 13        13       0.68         25
## 14        14       0.57         26
## 15        15       0.10         47
## 16        16       0.90          2
## 17        17       0.25        108
## 18        18       0.04         68
## 19        19       0.33          3
## 20        20       0.95          2
## 21        21       0.89          2
## 22        22       0.69         25
## 23        23       0.64         25
## 24        24       0.99          2
## 25        25       0.66         25
## 26        26       0.71         25
## 27        27       0.54         26
## 28        28       0.59         26
## 29        29       0.29        108
## 30        30       0.15         47
## 31        31       0.96          2
## 32        32       0.90          2
## 33        33       0.69         25
## 34        34       0.80         12
## 35        35       0.02         68
## 36        36       0.48          5
## 37        37       0.76         12
## 38        38       0.22         47
## 39        39       0.32        108
## 40        40       0.23        108
## 41        41       0.14         47
## 42        42       0.41          3
## 43        43       0.41          3
## 44        44       0.37          3
## 45        45       0.15         47
## 46        46       0.14         47
## 47        47       0.23        108
## 48        48       0.47          5
## 49        49       0.27        108
## 50        50       0.86         12
## 51        51       0.05         68
## 52        52       0.44          5
## 53        53       0.80         12
## 54        54       0.12         47
## 55        55       0.56         26
## 56        56       0.21         47
## 57        57       0.13         47
## 58        58       0.75         25
## 59        59       0.90          2
## 60        60       0.37          3
## 61        61       0.67         25
## 62        62       0.09         47
## 63        63       0.38          3
## 64        64       0.27        108
## 65        65       0.81         12
## 66        66       0.45          5
## 67        67       0.81         12
## 68        68       0.81         12
## 69        69       0.79         12
## 70        70       0.44          5
## 71        71       0.75         25
## 72        72       0.63         25
## 73        73       0.71         25
## 74        74       0.00         68
## 75        75       0.48          5
## 76        76       0.22         47
## 77        77       0.38          3
## 78        78       0.61         25
## 79        79       0.35          3
## 80        80       0.11         47
## 81        81       0.24        108
## 82        82       0.67         25
## 83        83       0.42          3
## 84        84       0.79         12
## 85        85       0.10         47
## 86        86       0.43          3
## 87        87       0.98          2
## 88        88       0.89          2
## 89        89       0.89          2
## 90        90       0.18         47
## 91        91       0.13         47
## 92        92       0.65         25
## 93        93       0.34          3
## 94        94       0.66         25
## 95        95       0.32        108
## 96        96       0.19         47
## 97        97       0.78         12
## 98        98       0.09         47
## 99        99       0.47          5
## 100      100       0.51          5
## 101      101       0.60         25
## 102      102       0.33          3
## 103      103       0.49          5
## 104      104       0.95          2
## 105      105       0.48          5
## 106      106       0.89          2
## 107      107       0.91          2
## 108      108       0.61         25
## 109      109       0.41          3
## 110      110       0.15         47
## 111      111       0.94          2
## 112      112       0.30        108
## 113      113       0.06         68
## 114      114       0.95          2
## 115      115       0.72         25
## 116      116       0.14         47
## 117      117       0.55         26
## 118      118       0.95          2
## 119      119       0.59         26
## 120      120       0.40          3
## 121      121       0.65         25
## 122      122       0.32        108
## 123      123       0.31        108
## 124      124       0.22         47
## 125      125       0.37          3
## 126      126       0.98          2
## 127      127       0.15         47
## 128      128       0.09         47
## 129      129       0.14         47
## 130      130       0.69         25
## 131      131       0.62         25
## 132      132       0.89          2
## 133      133       0.67         25
## 134      134       0.74         25
## 135      135       0.52          5
## 136      136       0.66         25
## 137      137       0.82         12
## 138      138       0.79         12
## 139      139       0.98          2
## 140      140       0.44          5
## 141      141       0.31        108
## 142      142       0.41          3
## 143      143       0.01         68
## 144      144       0.18         47
## 145      145       0.84         12
## 146      146       0.23        108
## 147      147       0.24        108
## 148      148       0.08         47
## 149      149       0.25        108
## 150      150       0.73         25
## 151      151       0.85         12
## 152      152       0.50          5
## 153      153       0.39          3
## 154      154       0.25        108
## 155      155       0.11         47
## 156      156       0.39          3
## 157      157       0.57         26
## 158      158       0.22         47
## 159      159       0.44          5
## 160      160       0.22         47
## 161      161       0.50          5
## 162      162       0.35          3
## 163      163       0.65         25
## 164      164       0.37          3
## 165      165       0.36          3
## 166      166       0.53         26
## 167      167       0.74         25
## 168      168       0.22         47
## 169      169       0.41          3
## 170      170       0.27        108
## 171      171       0.63         25
## 172      172       0.18         47
## 173      173       0.86         12
## 174      174       0.75         25
## 175      175       0.67         25
## 176      176       0.62         25
## 177      177       0.37          3
## 178      178       0.53         26
## 179      179       0.87        221
## 180      180       0.58         26
## 181      181       0.84         12
## 182      182       0.31        108
## 183      183       0.71         25
## 184      184       0.27        108
## 185      185       0.59         26
## 186      186       0.48          5
## 187      187       0.27        108
## 188      188       0.56         26
## 189      189       0.91          2
## 190      190       0.90          2
## 191      191       0.27        108
## 192      192       0.32        108
## 193      193       0.99          2
## 194      194       0.62         25
## 195      195       0.94          2
## 196      196       0.47          5
## 197      197       0.41          3
## 198      198       0.66         25
## 199      199       0.15         47
## 200      200       0.57         26
## 201      201       0.24        108
## 202      202       0.96          2
## 203      203       0.60         25
## 204      204       0.52          5
## 205      205       0.40          3
## 206      206       0.88          2
## 207      207       0.36          3
## 208      208       0.29        108
## 209      209       0.17         47
## 210      210       0.17         47
## 211      211       0.48          5
## 212      212       0.25        108
## 213      213       0.22         47
## 214      214       0.67         25
## 215      215       0.05         68
## 216      216       0.70         25
## 217      217       0.35          3
## 218      218       0.41          3
## 219      219       0.82         12
## 220      220       0.92          2
## 221      221       0.28        108
## 222      222       0.96          2
## 223      223       0.73         25
## 224      224       0.69         25
## 225      225       0.05         68
## 226      226       0.40          3
## 227      227       0.48          5
## 228      228       0.56         26
## 229      229       0.70         25
## 230      230       0.92          2
## 231      231       0.62         25
## 232      232       0.43          3
## 233      233       0.54         26
## 234      234       0.06         68
## 235      235       0.26        108
## 236      236       0.40          3
## 237      237       0.20         47
## 238      238       0.83         12
## 239      239       0.15         47
## 240      240       0.80         12
## 241      241       0.55         26
## 242      242       0.66         25
## 243      243       0.17         47
## 244      244       0.63         25
## 245      245       0.31        108
## 246      246       0.72         25
## 247      247       0.40          3
## 248      248       0.97          2
## 249      249       0.97          2
## 250      250       0.73         25
## 251      251       0.26        108
## 252      252       0.22         47
## 253      253       0.59         26
## 254      254       0.27        108
## 255      255       0.53         26
## 256      256       0.79         12
## 257      257       0.17         47
## 258      258       0.40          3
## 259      259       0.47          5
## 260      260       0.87        221
## 261      261       0.93          2
## 262      262       0.88          2
## 263      263       0.67         25
## 264      264       0.95          2
## 265      265       0.52          5
## 266      266       0.58         26
## 267      267       0.34          3
## 268      268       0.35          3
## 269      269       0.02         68
## 270      270       0.50          5
## 271      271       0.87        221
## 272      272       0.01         68
## 273      273       0.07         68
## 274      274       0.16         47
## 275      275       0.77         12
## 276      276       0.74         25
## 277      277       0.97          2
## 278      278       0.47          5
## 279      279       0.07         68
## 280      280       0.65         25
## 281      281       0.76         12
## 282      282       0.14         47
## 283      283       0.40          3
## 284      284       0.22         47
## 285      285       0.06         68
## 286      286       0.40          3
## 287      287       0.06         68
## 288      288       0.23        108
## 289      289       0.05         68
## 290      290       0.67         25
## 291      291       0.30        108
## 292      292       0.10         47
## 293      293       0.07         68
## 294      294       0.88          2
## 295      295       0.75         25
## 296      296       0.82         12
## 297      297       0.98          2
## 298      298       0.10         47
## 299      299       0.10         47
## 300      300       0.80         12
## 301      301       0.78         12
## 302      302       0.01         68
## 303      303       0.78         12
## 304      304       0.73         25
## 305      305       0.63         25
## 306      306       0.48          5
## 307      307       0.16         47
## 308      308       0.01         68
## 309      309       0.45          5
## 310      310       0.49          5
## 311      311       0.39          3
## 312      312       0.46          5
## 313      313       0.71         25
## 314      314       0.06         68
## 315      315       0.35          3
## 316      316       0.80         12
## 317      317       0.84         12
## 318      318       0.24        108
## 319      319       0.35          3
## 320      320       0.86         12
## 321      321       0.85         12
## 322      322       0.30        108
## 323      323       0.15         47
## 324      324       0.70         25
## 325      325       0.10         47
## 326      326       0.03         68
## 327      327       1.00          0
## 328      328       0.03         68
## 329      329       0.34          3
## 330      330       0.92          2
## 331      331       0.62         25
## 332      332       0.29        108
## 333      333       0.74         25
## 334      334       0.83         12
## 335      335       0.31        108
## 336      336       0.49          5
## 337      337       0.70         25
## 338      338       0.64         25
## 339      339       0.64         25
## 340      340       0.98          2
## 341      341       0.41          3
## 342      342       0.12         47
## 343      343       0.53         26
## 344      344       0.23        108
## 345      345       0.49          5
## 346      346       0.37          3
## 347      347       0.98          2
## 348      348       0.39          3
## 349      349       0.23        108
## 350      350       0.62         25
## 351      351       0.14         47
## 352      352       0.97          2
## 353      353       0.52          5
## 354      354       0.16         47
## 355      355       0.62         25
## 356      356       0.99          2
## 357      357       0.67         25
## 358      358       0.42          3
## 359      359       0.32        108
## 360      360       0.84         12
## 361      361       0.14         47
## 362      362       0.19         47
## 363      363       0.90          2
## 364      364       0.31        108
## 365      365       0.36          3
## 366      366       0.78         12
## 367      367       0.19         47
## 368      368       0.02         68
## 369      369       0.41          3
## 370      370       0.48          5
## 371      371       0.42          3
## 372      372       0.34          3
## 373      373       0.87        221
## 374      374       0.46          5
## 375      375       0.53         26
## 376      376       0.96          2
## 377      377       0.77         12
## 378      378       0.21         47
## 379      379       0.31        108
## 380      380       0.97          2
## 381      381       0.58         26
## 382      382       0.76         12
## 383      383       0.37          3
## 384      384       0.77         12
## 385      385       0.54         26
## 386      386       0.91          2
## 387      387       0.19         47
## 388      388       0.28        108
## 389      389       0.09         47
## 390      390       0.21         47
## 391      391       0.98          2
## 392      392       0.30        108
## 393      393       0.73         25
## 394      394       0.79         12
## 395      395       0.11         47
## 396      396       0.24        108
## 397      397       0.27        108
## 398      398       0.10         47
## 399      399       0.12         47
## 400      400       0.99          2
## 401      401       0.99          2
## 402      402       0.14         47
## 403      403       0.91          2
## 404      404       0.58         26
## 405      405       0.40          3
## 406      406       0.45          5
## 407      407       0.71         25
## 408      408       0.08         47
## 409      409       0.34          3
## 410      410       0.68         25
## 411      411       0.32        108
## 412      412       0.83         12
## 413      413       0.22         47
## 414      414       0.50          5
## 415      415       0.28        108
## 416      416       0.19         47
## 417      417       0.95          2
## 418      418       0.32        108
## 419      419       0.48          5
## 420      420       0.03         68
## 421      421       0.55         26
## 422      422       0.64         25
## 423      423       0.60         25
## 424      424       0.32        108
## 425      425       0.89          2
## 426      426       0.63         25
## 427      427       0.30        108
## 428      428       0.39          3
## 429      429       0.16         47
## 430      430       0.86         12
## 431      431       0.95          2
## 432      432       0.56         26
## 433      433       0.33          3
## 434      434       1.00          0
## 435      435       0.23        108
## 436      436       0.61         25
## 437      437       0.11         47
## 438      438       0.49          5
## 439      439       0.10         47
## 440      440       0.16         47
## 441      441       0.28        108
## 442      442       0.58         26
## 443      443       0.73         25
## 444      444       0.17         47
## 445      445       0.87        221
## 446      446       0.71         25
## 447      447       0.76         12
## 448      448       0.15         47
## 449      449       0.36          3
## 450      450       0.67         25
## 451      451       0.52          5
## 452      452       0.35          3
## 453      453       0.24        108
## 454      454       0.06         68
## 455      455       0.24        108
## 456      456       0.89          2
## 457      457       0.81         12
## 458      458       0.75         25
## 459      459       0.15         47
## 460      460       0.12         47
## 461      461       0.97          2
## 462      462       0.44          5
## 463      463       0.46          5
## 464      464       0.17         47
## 465      465       0.58         26
## 466      466       0.27        108
## 467      467       0.23        108
## 468      468       0.69         25
## 469      469       0.28        108
## 470      470       0.81         12
## 471      471       0.09         47
## 472      472       0.82         12
## 473      473       0.43          3
## 474      474       0.76         12
## 475      475       0.66         25
## 476      476       0.44          5
## 477      477       0.63         25
## 478      478       0.00         68
## 479      479       0.22         47
## 480      480       0.70         25
## 481      481       0.22         47
## 482      482       0.81         12
## 483      483       0.31        108
## 484      484       0.69         25
## 485      485       0.93          2
## 486      486       0.12         47
## 487      487       0.13         47
## 488      488       0.68         25
## 489      489       0.43          3
## 490      490       0.83         12
## 491      491       0.97          2
## 492      492       0.07         68
## 493      493       0.46          5
## 494      494       0.70         25
## 495      495       0.09         47
## 496      496       0.99          2
## 497      497       0.25        108
## 498      498       0.05         68
## 499      499       0.69         25
## 500      500       0.79         12
## 501      501       0.35          3
## 502      502       0.37          3
## 503      503       0.29        108
## 504      504       0.08         47
## 505      505       0.37          3
## 506      506       0.18         47
## 507      507       0.54         26
## 508      508       0.50          5
## 509      509       0.95          2
## 510      510       0.34          3
## 511      511       0.46          5
## 512      512       0.08         47
## 513      513       0.86         12
## 514      514       0.40          3
## 515      515       0.74         25
## 516      516       0.17         47
## 517      517       0.45          5
## 518      518       0.77         12
## 519      519       0.06         68
## 520      520       0.82         12
## 521      521       0.30        108
## 522      522       0.36          3
## 523      523       0.31        108
## 524      524       0.04         68
## 525      525       0.52          5
## 526      526       0.68         25
## 527      527       0.90          2
## 528      528       0.03         68
## 529      529       0.99          2
## 530      530       0.30        108
## 531      531       0.94          2
## 532      532       0.69         25
## 533      533       0.45          5
## 534      534       0.82         12
## 535      535       0.04         68
## 536      536       0.74         25
## 537      537       0.35          3
## 538      538       0.83         12
## 539      539       0.54         26
## 540      540       0.27        108
## 541      541       0.80         12
## 542      542       0.09         47
## 543      543       0.83         12
## 544      544       0.28        108
## 545      545       0.75         25
## 546      546       0.96          2
## 547      547       0.08         47
## 548      548       0.85         12
## 549      549       0.80         12
## 550      550       0.39          3
## 551      551       0.33          3
## 552      552       0.20         47
## 553      553       0.57         26
## 554      554       0.89          2
## 555      555       0.53         26
## 556      556       0.59         26
## 557      557       0.67         25
## 558      558       0.53         26
## 559      559       0.51          5
## 560      560       0.02         68
## 561      561       0.05         68
## 562      562       0.93          2
## 563      563       0.77         12
## 564      564       0.20         47
## 565      565       0.65         25
## 566      566       0.65         25
## 567      567       0.40          3
## 568      568       0.81         12
## 569      569       0.55         26
## 570      570       0.89          2
## 571      571       0.55         26
## 572      572       0.91          2
## 573      573       0.59         26
## 574      574       0.42          3
## 575      575       0.95          2
## 576      576       0.71         25
## 577      577       0.41          3
## 578      578       0.02         68
## 579      579       0.57         26
## 580      580       0.49          5
## 581      581       0.88          2
## 582      582       0.81         12
## 583      583       0.85         12
## 584      584       0.37          3
## 585      585       0.87        221
## 586      586       0.15         47
## 587      587       0.28        108
## 588      588       0.67         25
## 589      589       0.98          2
## 590      590       0.58         26
## 591      591       0.53         26
## 592      592       0.06         68
## 593      593       0.97          2
## 594      594       0.12         47
## 595      595       0.09         47
## 596      596       0.88          2
## 597      597       0.51          5
## 598      598       0.34          3
## 599      599       0.89          2
## 600      600       0.03         68
## 601      601       0.24        108
## 602      602       0.69         25
## 603      603       0.23        108
## 604      604       0.32        108
## 605      605       0.17         47
## 606      606       0.80         12
## 607      607       0.15         47
## 608      608       0.82         12
## 609      609       0.33          3
## 610      610       0.37          3
## 611      611       0.63         25
## 612      612       0.10         47
## 613      613       0.02         68
## 614      614       0.99          2
## 615      615       0.58         26
## 616      616       0.78         12
## 617      617       0.89          2
## 618      618       0.75         25
## 619      619       0.98          2
## 620      620       0.04         68
## 621      621       0.90          2
## 622      622       0.87        221
## 623      623       0.78         12
## 624      624       0.38          3
## 625      625       0.04         68
## 626      626       0.36          3
## 627      627       0.27        108
## 628      628       0.85         12
## 629      629       0.36          3
## 630      630       0.30        108
## 631      631       0.76         12
## 632      632       0.84         12
## 633      633       0.46          5
## 634      634       0.73         25
## 635      635       0.10         47
## 636      636       0.22         47
## 637      637       0.95          2
## 638      638       0.75         25
## 639      639       0.82         12
## 640      640       0.42          3
## 641      641       0.59         26
## 642      642       0.80         12
## 643      643       0.89          2
## 644      644       0.39          3
## 645      645       0.09         47
## 646      646       0.63         25
## 647      647       0.75         25
## 648      648       0.09         47
## 649      649       0.30        108
## 650      650       0.61         25
## 651      651       0.75         25
## 652      652       0.92          2
## 653      653       0.48          5
## 654      654       0.57         26
## 655      655       0.74         25
## 656      656       0.86         12
## 657      657       0.91          2
## 658      658       0.06         68
## 659      659       0.50          5
## 660      660       0.35          3
## 661      661       0.85         12
## 662      662       0.81         12
## 663      663       0.12         47
## 664      664       0.71         25
## 665      665       0.24        108
## 666      666       0.07         68
## 667      667       0.94          2
## 668      668       0.16         47
## 669      669       0.65         25
## 670      670       0.17         47
## 671      671       0.02         68
## 672      672       0.52          5
## 673      673       0.09         47
## 674      674       0.28        108
## 675      675       0.42          3
## 676      676       0.59         26
## 677      677       0.81         12
## 678      678       0.20         47
## 679      679       0.46          5
## 680      680       0.45          5
## 681      681       0.73         25
## 682      682       0.71         25
## 683      683       0.83         12
## 684      684       0.89          2
## 685      685       0.95          2
## 686      686       0.55         26
## 687      687       0.52          5
## 688      688       0.17         47
## 689      689       0.48          5
## 690      690       0.25        108
## 691      691       0.04         68
## 692      692       0.63         25
## 693      693       0.54         26
## 694      694       0.14         47
## 695      695       0.28        108
## 696      696       0.58         26
## 697      697       0.17         47
## 698      698       0.10         47
## 699      699       0.43          3
## 700      700       0.36          3
## 701      701       0.84         12
## 702      702       0.26        108
## 703      703       0.02         68
## 704      704       0.86         12
## 705      705       0.33          3
## 706      706       0.63         25
## 707      707       0.55         26
## 708      708       0.38          3
## 709      709       0.19         47
## 710      710       0.43          3
## 711      711       0.63         25
## 712      712       0.52          5
## 713      713       0.66         25
## 714      714       0.73         25
## 715      715       0.49          5
## 716      716       0.38          3
## 717      717       0.01         68
## 718      718       0.00         68
## 719      719       0.99          2
## 720      720       0.11         47
## 721      721       0.42          3
## 722      722       0.72         25
## 723      723       0.74         25
## 724      724       0.87        221
## 725      725       0.61         25
## 726      726       0.76         12
## 727      727       0.85         12
## 728      728       0.61         25
## 729      729       0.79         12
## 730      730       0.02         68
## 731      731       0.42          3
## 732      732       0.87        221
## 733      733       0.65         25
## 734      734       0.92          2
## 735      735       0.17         47
## 736      736       0.32        108
## 737      737       0.80         12
## 738      738       0.99          2
## 739      739       0.32        108
## 740      740       0.00         68
## 741      741       0.99          2
## 742      742       0.15         47
## 743      743       0.05         68
## 744      744       0.60         25
## 745      745       0.32        108
## 746      746       0.53         26
## 747      747       0.80         12
## 748      748       0.07         68
## 749      749       0.70         25
## 750      750       0.95          2
## 751      751       0.28        108
## 752      752       0.25        108
## 753      753       0.16         47
## 754      754       0.07         68
## 755      755       0.56         26
## 756      756       0.49          5
## 757      757       0.49          5
## 758      758       0.89          2
## 759      759       0.04         68
## 760      760       0.20         47
## 761      761       0.51          5
## 762      762       0.24        108
## 763      763       0.58         26
## 764      764       0.48          5
## 765      765       0.57         26
## 766      766       0.14         47
## 767      767       0.15         47
## 768      768       0.41          3
## 769      769       0.68         25
## 770      770       0.65         25
## 771      771       0.92          2
## 772      772       0.81         12
## 773      773       0.68         25
## 774      774       0.81         12
## 775      775       0.13         47
## 776      776       0.25        108
## 777      777       0.33          3
## 778      778       0.41          3
## 779      779       0.64         25
## 780      780       0.81         12
## 781      781       0.26        108
## 782      782       0.82         12
## 783      783       0.02         68
## 784      784       0.65         25
## 785      785       0.81         12
## 786      786       0.46          5
## 787      787       0.20         47
## 788      788       0.02         68
## 789      789       0.28        108
## 790      790       0.92          2
## 791      791       0.93          2
## 792      792       0.54         26
## 793      793       0.22         47
## 794      794       0.49          5
## 795      795       0.80         12
## 796      796       0.41          3
## 797      797       0.10         47
## 798      798       0.28        108
## 799      799       0.36          3
## 800      800       0.26        108
## 801      801       0.47          5
## 802      802       0.37          3
## 803      803       0.12         47
## 804      804       0.05         68
## 805      805       0.26        108
## 806      806       0.97          2
## 807      807       0.49          5
## 808      808       0.48          5
## 809      809       0.75         25
## 810      810       0.67         25
## 811      811       0.05         68
## 812      812       0.70         25
## 813      813       0.36          3
## 814      814       0.88          2
## 815      815       0.78         12
## 816      816       0.14         47
## 817      817       0.30        108
## 818      818       0.13         47
## 819      819       0.59         26
## 820      820       0.56         26
## 821      821       0.69         25
## 822      822       0.31        108
## 823      823       0.61         25
## 824      824       0.99          2
## 825      825       0.74         25
## 826      826       0.08         47
## 827      827       0.45          5
## 828      828       0.05         68
## 829      829       0.34          3
## 830      830       0.73         25
## 831      831       0.00         68
## 832      832       0.77         12
## 833      833       0.46          5
## 834      834       0.72         25
## 835      835       0.67         25
## 836      836       0.57         26
## 837      837       0.70         25
## 838      838       0.66         25
## 839      839       0.29        108
## 840      840       0.10         47
## 841      841       0.96          2
## 842      842       0.74         25
## 843      843       0.61         25
## 844      844       0.12         47
## 845      845       0.55         26
## 846      846       0.26        108
## 847      847       0.90          2
## 848      848       0.01         68
## 849      849       0.24        108
## 850      850       0.13         47
## 851      851       0.33          3
## 852      852       0.73         25
## 853      853       0.99          2
## 854      854       0.72         25
## 855      855       0.50          5
## 856      856       0.44          5
## 857      857       0.95          2
## 858      858       0.12         47
## 859      859       0.08         47
## 860      860       0.89          2
## 861      861       0.42          3
## 862      862       0.04         68
## 863      863       0.65         25
## 864      864       0.47          5
## 865      865       0.62         25
## 866      866       0.27        108
## 867      867       0.16         47
## 868      868       0.11         47
## 869      869       0.51          5
## 870      870       0.55         26
## 871      871       0.14         47
## 872      872       0.17         47
## 873      873       0.76         12
## 874      874       0.53         26
## 875      875       0.86         12
## 876      876       0.67         25
## 877      877       0.01         68
## 878      878       0.69         25
## 879      879       0.89          2
## 880      880       0.63         25
## 881      881       0.11         47
## 882      882       0.92          2
## 883      883       0.68         25
## 884      884       0.15         47
## 885      885       0.75         25
## 886      886       0.94          2
## 887      887       0.42          3
## 888      888       0.30        108
## 889      889       0.26        108
## 890      890       0.22         47
## 891      891       0.57         26
## 892      892       0.76         12
## 893      893       0.67         25
## 894      894       0.55         26
## 895      895       0.81         12
## 896      896       0.76         12
## 897      897       0.02         68
## 898      898       0.38          3
## 899      899       0.05         68
## 900      900       0.80         12
## 901      901       0.92          2
## 902      902       0.54         26
## 903      903       0.85         12
## 904      904       0.58         26
## 905      905       0.67         25
## 906      906       0.51          5
## 907      907       0.76         12
## 908      908       0.90          2
## 909      909       0.82         12
## 910      910       0.07         68
## 911      911       0.00         68
## 912      912       0.05         68
## 913      913       0.87        221
## 914      914       0.58         26
## 915      915       0.31        108
## 916      916       0.96          2
## 917      917       0.59         26
## 918      918       0.53         26
## 919      919       0.38          3
## 920      920       0.32        108
## 921      921       0.81         12
## 922      922       0.04         68
## 923      923       0.36          3
## 924      924       0.86         12
## 925      925       0.70         25
## 926      926       0.68         25
## 927      927       0.35          3
## 928      928       0.55         26
## 929      929       0.14         47
## 930      930       0.78         12
## 931      931       0.89          2
## 932      932       0.20         47
## 933      933       0.77         12
## 934      934       0.60         25
## 935      935       0.96          2
## 936      936       0.16         47
## 937      937       0.53         26
## 938      938       0.87        221
## 939      939       0.87        221
## 940      940       0.02         68
## 941      941       0.98          2
## 942      942       0.49          5
## 943      943       0.39          3
## 944      944       0.42          3
## 945      945       0.09         47
## 946      946       0.16         47
## 947      947       0.41          3
## 948      948       0.34          3
## 949      949       0.42          3
## 950      950       0.30        108
## 951      951       0.56         26
## 952      952       0.16         47
## 953      953       0.96          2
## 954      954       0.04         68
## 955      955       0.37          3
## 956      956       0.96          2
## 957      957       0.65         25
## 958      958       0.06         68
## 959      959       0.41          3
## 960      960       0.43          3
## 961      961       0.51          5
## 962      962       0.45          5
## 963      963       0.62         25
## 964      964       0.14         47
## 965      965       0.91          2
## 966      966       0.57         26
## 967      967       0.55         26
## 968      968       0.12         47
## 969      969       0.76         12
## 970      970       0.48          5
## 971      971       0.78         12
## 972      972       0.05         68
## 973      973       0.82         12
## 974      974       0.27        108
## 975      975       0.28        108
## 976      976       0.64         25
## 977      977       0.95          2
## 978      978       0.01         68
## 979      979       0.35          3
## 980      980       0.42          3
## 981      981       0.46          5
## 982      982       0.71         25
## 983      983       0.92          2
## 984      984       0.63         25
## 985      985       0.90          2
## 986      986       0.76         12
## 987      987       0.14         47
## 988      988       0.15         47
## 989      989       0.19         47
## 990      990       0.43          3
## 991      991       0.09         47
## 992      992       0.22         47
## 993      993       0.57         26
## 994      994       0.40          3
## 995      995       0.57         26
## 996      996       0.83         12
## 997      997       0.64         25
## 998      998       0.39          3
## 999      999       0.71         25
## 1000    1000       0.11         47
print(simulasi_d$average)
## [1] 34.072