Latihan 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: a. prediksikan permintaan es teh selama 5 hari ke depan. b. prediksikan permintaan es teh selama 20 hari ke depan.

Jawab :

set.seed(123)
data <- data.frame(
  permintaan = c(50, 60, 70, 80, 90),
  frekuensi = c(10, 15, 30, 25, 10)
)
data
##   permintaan frekuensi
## 1         50        10
## 2         60        15
## 3         70        30
## 4         80        25
## 5         90        10
data$probabilitas <- round(data$frekuensi / sum(data$frekuensi),2)
data
##   permintaan frekuensi probabilitas
## 1         50        10         0.11
## 2         60        15         0.17
## 3         70        30         0.33
## 4         80        25         0.28
## 5         90        10         0.11
data$probcum <- cumsum(data$probabilitas)
data
##   permintaan frekuensi probabilitas probcum
## 1         50        10         0.11    0.11
## 2         60        15         0.17    0.28
## 3         70        30         0.33    0.61
## 4         80        25         0.28    0.89
## 5         90        10         0.11    1.00
expvalue <- round(sum(data$permintaan*data$probabilitas),2)
expvalue
## [1] 71.1

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         70
## 2       2       0.79         80
## 3       3       0.41         70
## 4       4       0.88         80
## 5       5       0.94         90
print(simulasi_a$average)
## [1] 78

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         70
## 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         80
## 14      14       0.57         70
## 15      15       0.10         50
## 16      16       0.90         90
## 17      17       0.25         60
## 18      18       0.04         50
## 19      19       0.33         70
## 20      20       0.95         90
print(simulasi_b$average)
## [1] 72.5

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

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

Latihan 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

Jawab :

set.seed(456)
data <- data.frame(
  permintaan = round(rexp(20, rate = 0.015)),
  frekuensi  = pmax(0, ceiling(rnorm(20, mean = 12, sd = 4))))
print(data)
##    permintaan frekuensi
## 1         167        19
## 2         138        18
## 3          31        13
## 4          15        12
## 5         160         8
## 6          56         7
## 7         110        13
## 8         127        18
## 9          13        16
## 10         20        12
## 11         46        17
## 12         60        14
## 13        122        10
## 14         30        12
## 15         21        14
## 16         27        11
## 17          3        12
## 18         47        13
## 19         10        17
## 20         20         7

Menghitung probabilitas

data$probabilitas <- round(data$frekuensi / sum(data$frekuensi),2)
data
##    permintaan frekuensi probabilitas
## 1         167        19         0.07
## 2         138        18         0.07
## 3          31        13         0.05
## 4          15        12         0.05
## 5         160         8         0.03
## 6          56         7         0.03
## 7         110        13         0.05
## 8         127        18         0.07
## 9          13        16         0.06
## 10         20        12         0.05
## 11         46        17         0.06
## 12         60        14         0.05
## 13        122        10         0.04
## 14         30        12         0.05
## 15         21        14         0.05
## 16         27        11         0.04
## 17          3        12         0.05
## 18         47        13         0.05
## 19         10        17         0.06
## 20         20         7         0.03

Menghitung probabilitas kumulatif

data$probcum <- cumsum(data$probabilitas)
data
##    permintaan frekuensi probabilitas probcum
## 1         167        19         0.07    0.07
## 2         138        18         0.07    0.14
## 3          31        13         0.05    0.19
## 4          15        12         0.05    0.24
## 5         160         8         0.03    0.27
## 6          56         7         0.03    0.30
## 7         110        13         0.05    0.35
## 8         127        18         0.07    0.42
## 9          13        16         0.06    0.48
## 10         20        12         0.05    0.53
## 11         46        17         0.06    0.59
## 12         60        14         0.05    0.64
## 13        122        10         0.04    0.68
## 14         30        12         0.05    0.73
## 15         21        14         0.05    0.78
## 16         27        11         0.04    0.82
## 17          3        12         0.05    0.87
## 18         47        13         0.05    0.92
## 19         10        17         0.06    0.98
## 20         20         7         0.03    1.01

Menghitung nilai ekspektasi

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

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         56
## 2       2       0.79         27
## 3       3       0.41        127
## 4       4       0.88         47
## 5       5       0.94         10
print(simulasi_a$average)
## [1] 53.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         56
## 2        2       0.79         27
## 3        3       0.41        127
## 4        4       0.88         47
## 5        5       0.94         10
## 6        6       0.05        167
## 7        7       0.53         20
## 8        8       0.89         47
## 9        9       0.55         46
## 10      10       0.46         13
## 11      11       0.96         10
## 12      12       0.45         13
## 13      13       0.68        122
## 14      14       0.57         46
## 15      15       0.10        138
## 16      16       0.90         47
## 17      17       0.25        160
## 18      18       0.04        167
## 19      19       0.33        110
## 20      20       0.95         10
print(simulasi_b$average)
## [1] 69.15

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

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