1 Bài tập về nhà tuần 1

library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.0     ✔ readr     2.1.4
## ✔ forcats   1.0.0     ✔ stringr   1.5.0
## ✔ ggplot2   3.4.1     ✔ tibble    3.1.8
## ✔ lubridate 1.9.2     ✔ tidyr     1.3.0
## ✔ purrr     1.0.1     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the ]8;;http://conflicted.r-lib.org/conflicted package]8;; to force all conflicts to become errors
library(ggplot2)

1.1 Mô phỏng ngẫu nhiên các biến

1.1.1 Mô phỏng ngẫu nhiên phân phối đều

Số phần bánh căn bán được mỗi sáng trong vòng 1000 ngày

bc <- runif(n=1000,min=45,max=50)
bc
##    [1] 48.17172 46.47293 46.91364 45.96604 47.15818 46.73214 49.82236 45.61467
##    [9] 49.70294 47.74680 49.84103 49.31941 49.68380 48.34924 45.09450 46.16904
##   [17] 45.25343 49.23505 45.68107 45.75778 47.12507 48.33839 48.00711 47.73238
##   [25] 47.88661 46.81893 48.11934 49.38977 49.97373 49.78542 47.35312 45.03052
##   [33] 47.60056 46.26701 45.41744 49.66814 49.64730 46.89350 48.60348 45.58803
##   [41] 45.31645 48.90212 47.17307 47.87675 46.67355 47.77493 46.27518 46.54407
##   [49] 45.50909 46.86035 49.80177 49.03402 47.01553 48.39046 47.40571 45.19106
##   [57] 48.02647 49.09744 48.24869 47.58794 45.86063 46.29926 47.70724 46.34076
##   [65] 45.54983 47.92973 45.08416 49.61199 48.35057 46.14975 45.82769 49.72618
##   [73] 46.44111 46.24579 45.46338 49.87211 49.19816 45.82040 46.67331 46.17758
##   [81] 48.87964 47.65258 47.73922 45.38377 46.37295 48.20347 48.21489 46.94603
##   [89] 45.04724 49.80502 45.07566 49.14473 48.05662 48.26662 49.72329 48.42697
##   [97] 49.61222 48.18529 45.95306 45.28809 48.09786 49.48052 47.38048 46.31435
##  [105] 49.80922 46.41718 49.29673 47.84563 47.48177 46.56419 49.36585 46.53262
##  [113] 48.47406 45.02026 47.64308 49.98507 47.46566 47.90106 45.16946 49.07228
##  [121] 49.33111 45.32904 47.46124 49.77562 47.21551 49.26163 47.30487 48.57247
##  [129] 48.40243 45.40309 48.94963 49.02634 48.63046 45.32984 47.83744 49.28518
##  [137] 47.96411 49.66250 49.61478 45.40087 48.19853 46.80716 47.03463 46.93921
##  [145] 46.31178 48.30937 46.78917 49.76869 49.11812 47.24412 45.41382 49.29656
##  [153] 45.66539 48.19293 49.71538 48.25005 48.67823 49.55110 47.49469 48.25975
##  [161] 48.86862 47.76867 47.93598 47.65069 45.50641 48.68490 45.23007 47.64245
##  [169] 45.94585 45.17748 46.38291 45.68434 49.71689 47.82969 46.51023 47.50102
##  [177] 45.53220 47.84713 46.85936 47.95358 47.54561 48.46347 45.08930 48.86052
##  [185] 45.50122 46.70629 48.59250 49.66539 49.01250 46.51189 48.48178 46.10353
##  [193] 49.80602 45.93927 45.06830 46.78180 45.66443 46.14977 48.02254 47.16571
##  [201] 47.33826 47.29140 49.00702 48.82637 49.15287 48.23330 49.99447 46.16232
##  [209] 49.08835 47.13538 49.13112 48.39841 48.32564 46.96916 48.67797 49.72504
##  [217] 45.80476 48.34385 48.15872 46.56716 49.90354 48.54307 48.18028 49.30314
##  [225] 49.93280 48.16857 46.30328 46.47701 45.72570 49.64580 45.12606 48.48370
##  [233] 45.06025 45.33693 45.11782 49.69504 49.05957 47.69774 46.20491 48.95076
##  [241] 46.76586 46.55711 46.65929 47.92278 46.25024 49.02858 45.41802 48.57012
##  [249] 48.28157 48.91178 47.95757 45.26507 45.49841 45.62847 48.69869 47.70223
##  [257] 45.57349 48.16387 48.83574 45.20186 45.86676 46.23641 45.23954 46.85093
##  [265] 47.25704 48.16151 49.31488 45.82817 49.27874 46.48541 46.43506 45.02145
##  [273] 46.58263 46.65378 48.83061 48.36388 49.53816 49.85775 49.62665 47.84452
##  [281] 45.48711 49.65714 46.86416 45.77179 45.34572 48.04091 48.81036 46.29368
##  [289] 48.63984 48.63583 45.64945 49.19609 48.85641 49.48888 46.00257 48.55675
##  [297] 47.16300 49.79015 48.13159 48.31218 46.66545 46.96919 46.70957 46.05251
##  [305] 48.76189 45.39419 45.81883 48.78549 46.61429 45.13062 47.60104 46.28610
##  [313] 45.48998 47.86473 47.60762 46.56983 47.00572 49.75558 45.65879 49.68164
##  [321] 47.63166 48.45087 49.45225 45.53083 46.24147 45.81163 47.98594 48.87128
##  [329] 45.17756 45.73935 49.70952 48.91457 45.81630 47.29898 48.85465 45.80204
##  [337] 46.55456 46.31083 47.17007 49.85333 49.55822 46.61171 49.45414 45.74820
##  [345] 47.64739 47.14157 49.91663 45.16215 47.34822 49.09389 47.26108 49.22066
##  [353] 45.03473 49.34060 46.91029 46.64127 49.90567 47.45478 45.14138 48.21157
##  [361] 47.32741 47.59519 45.21660 48.20284 48.76827 47.21878 48.18906 46.71139
##  [369] 48.47361 46.72909 49.40793 49.27930 48.60141 48.95549 46.34736 48.23482
##  [377] 48.99981 45.01414 45.36620 45.76675 45.93999 49.28964 49.57143 47.42062
##  [385] 48.86333 46.47595 48.46941 49.95614 49.41120 45.00242 47.64089 45.53692
##  [393] 47.55362 47.77865 49.07959 46.03878 48.31775 45.47507 45.73392 47.18679
##  [401] 49.79980 45.71846 47.21824 45.28002 48.34813 49.65083 47.82554 49.81695
##  [409] 48.25728 47.61898 49.79798 49.96941 49.96846 46.38270 48.44638 48.83837
##  [417] 46.73937 46.80452 46.24069 45.23816 48.10304 48.56880 45.38806 49.69123
##  [425] 47.77679 46.74429 47.49236 45.34701 45.77711 49.53534 46.82884 48.47844
##  [433] 46.53249 46.92141 47.31393 47.02248 45.75426 49.58450 48.50730 49.30222
##  [441] 48.80556 49.23608 49.89267 49.92392 47.11389 46.23309 48.56711 47.71815
##  [449] 46.85568 48.89235 49.27604 49.97635 45.96180 48.66546 45.83893 47.54675
##  [457] 47.47671 48.06201 48.91858 48.04577 45.06111 45.69743 47.68163 46.86916
##  [465] 46.65577 47.12956 48.80717 46.83641 45.29006 47.37255 49.46329 46.71638
##  [473] 46.19898 47.31267 46.88386 45.38112 49.31493 47.42512 46.64790 45.31115
##  [481] 45.52332 47.77819 49.70998 49.03707 47.17543 49.89720 49.39051 48.05174
##  [489] 47.37661 48.24392 49.45934 47.41728 49.80970 47.32578 49.72958 45.84316
##  [497] 46.47645 45.72654 46.40591 47.03888 46.44559 48.84184 49.23883 45.24646
##  [505] 46.48209 48.16749 48.24714 49.17239 47.98585 47.77438 47.58363 49.65095
##  [513] 46.38771 45.00345 47.84202 49.21292 46.56135 46.33184 49.80326 49.93605
##  [521] 49.92381 46.84958 45.11146 49.92055 48.58177 47.50886 48.90509 46.08976
##  [529] 47.88143 46.87516 46.09968 45.50625 48.03748 45.21255 48.77973 47.89591
##  [537] 45.59958 46.20924 49.34475 49.08416 49.97209 47.70742 46.87624 47.68547
##  [545] 47.28681 48.61493 49.48490 49.08776 49.45875 47.98489 48.20264 48.06767
##  [553] 49.82783 48.58430 48.57566 46.92420 47.54939 49.98066 47.55414 48.78005
##  [561] 49.16920 46.23195 49.04452 47.41980 46.21087 45.94453 48.92802 47.85371
##  [569] 45.58234 45.95909 48.41486 45.41238 47.20181 45.35198 48.91646 49.94955
##  [577] 45.08296 47.00536 47.61945 45.88545 49.60149 49.16511 48.53555 47.18636
##  [585] 49.53609 47.83577 46.38737 49.51051 47.98856 46.73886 49.15697 46.53481
##  [593] 49.45521 46.97602 45.28062 47.40467 47.58556 49.88238 45.04804 49.52795
##  [601] 49.56406 47.68981 47.34574 49.95049 47.89939 49.94999 45.90175 47.57059
##  [609] 48.41553 46.46747 46.73022 46.67357 46.66956 46.58247 47.46103 45.76207
##  [617] 45.39000 45.63486 49.05475 46.79490 45.94779 49.26340 46.89789 47.08071
##  [625] 49.77692 45.96332 47.36389 48.86774 46.54661 45.60376 45.94611 49.20996
##  [633] 49.46779 46.34568 45.80230 48.86468 49.51976 46.22825 46.24992 49.64885
##  [641] 45.09650 49.04906 47.28651 48.38841 45.37699 47.29265 46.10174 45.40220
##  [649] 49.85313 47.90292 47.92501 49.89930 45.78570 45.89343 46.24920 49.81394
##  [657] 49.13162 45.55510 49.21792 45.21235 48.11724 45.47058 48.60834 45.37986
##  [665] 46.81824 45.08133 48.07973 46.26349 48.32395 47.39256 45.63904 48.53564
##  [673] 46.81048 49.71839 49.39802 45.32581 48.91207 48.10013 46.85114 45.08741
##  [681] 48.86804 45.94476 49.73787 48.99764 49.33456 45.10232 48.73824 45.87295
##  [689] 45.15217 47.87243 46.50265 45.20005 46.24166 49.80796 48.01785 49.12932
##  [697] 45.87808 49.63871 47.62282 47.65259 45.63294 45.73028 49.53116 47.33232
##  [705] 47.27047 49.70720 45.16349 49.98426 47.91088 48.13134 45.08591 46.76797
##  [713] 49.93536 48.07551 48.22232 45.14063 45.20652 45.82100 48.68797 48.03300
##  [721] 49.10402 47.56069 47.59441 46.30632 47.97312 48.69039 46.18785 48.27858
##  [729] 48.59561 45.83035 49.18838 46.89511 45.38627 46.15258 46.02370 46.67066
##  [737] 47.41292 46.51972 48.13157 45.21884 46.40864 47.46899 45.71062 47.37233
##  [745] 46.82349 46.89006 49.56685 47.65770 45.94542 47.92313 49.06989 45.95954
##  [753] 49.50256 49.16689 49.74605 45.55534 49.23438 46.61615 45.96622 48.24052
##  [761] 45.05308 45.92807 49.08339 47.16105 49.09085 46.18523 49.96831 49.96168
##  [769] 48.24742 49.53101 47.39361 48.11075 47.88961 47.78010 48.86453 49.43933
##  [777] 49.17103 49.62477 47.34435 48.45091 45.06394 46.87366 45.35089 48.52013
##  [785] 47.24576 47.90932 49.63712 48.84580 45.97246 48.74499 49.63803 46.79449
##  [793] 45.02466 45.23924 47.54108 48.10953 45.70431 49.06936 47.29700 45.70247
##  [801] 49.13204 48.22661 47.38300 46.95579 48.22736 45.88525 45.55550 45.64898
##  [809] 47.26330 47.68784 49.56542 48.45470 47.59255 48.58980 47.21432 48.26592
##  [817] 48.67091 46.80459 48.34653 47.17270 46.81337 46.86691 49.57075 48.47087
##  [825] 46.39020 49.10982 45.19402 45.67516 46.26907 47.06629 45.36320 46.19916
##  [833] 49.78098 48.47329 49.37798 45.57636 45.26513 49.70447 48.97577 46.97025
##  [841] 46.83172 45.98193 45.27234 48.59292 46.12675 46.66669 48.50812 47.11464
##  [849] 47.63267 48.81749 49.62880 47.31932 48.56571 46.35623 49.94901 46.30739
##  [857] 47.20741 45.45621 47.79835 49.68915 45.37161 48.33425 48.29023 49.65032
##  [865] 47.44737 49.09696 48.93926 45.88252 47.51894 46.02090 48.28046 48.55829
##  [873] 46.72503 45.34251 49.97511 47.78823 47.18434 49.17093 47.05150 49.69843
##  [881] 46.48831 49.65403 45.96542 48.47152 45.63753 49.41686 46.08906 46.67584
##  [889] 49.98170 45.01192 45.76142 47.45049 45.06551 45.86893 46.80338 46.05554
##  [897] 46.23653 46.35752 45.01209 49.12548 46.58892 47.90554 47.14325 47.21393
##  [905] 45.84321 45.88637 47.90425 47.32940 49.58945 46.54233 49.42519 47.62977
##  [913] 49.17571 47.10594 47.54156 47.34109 45.17522 49.03014 47.47423 48.05533
##  [921] 49.46329 49.55963 49.64440 46.27507 49.82784 47.46247 47.38211 45.31173
##  [929] 46.24783 47.50365 49.05965 46.47551 47.11414 48.58595 48.39450 47.09053
##  [937] 45.97857 48.38392 48.03714 49.19276 49.30940 48.84343 47.21084 48.22086
##  [945] 45.09703 46.06009 49.78361 45.34326 47.18948 49.89352 48.27337 48.81225
##  [953] 49.73670 45.75053 49.84344 45.39174 49.44201 49.75305 48.65104 46.06500
##  [961] 45.48771 47.81595 47.55151 49.76694 46.85545 47.57962 46.07326 46.00463
##  [969] 47.52432 48.18510 48.07315 48.14390 47.72819 46.02485 47.24895 47.03769
##  [977] 48.52834 47.75633 45.12610 48.78656 49.08765 49.70992 49.91333 48.22816
##  [985] 48.34729 49.67361 49.18544 45.12446 49.06338 49.01307 45.34542 48.29249
##  [993] 47.27163 46.59008 48.96167 49.55313 49.66394 45.93210 49.51430 46.40945
round(runif(n=1000,min=45,max=50),0)
##    [1] 45 50 46 47 49 49 49 48 48 46 47 49 46 48 47 47 50 48 49 47 45 48 47 47
##   [25] 49 50 46 45 46 49 47 49 45 48 47 46 49 49 46 50 46 48 48 47 48 47 50 50
##   [49] 48 48 45 48 45 47 49 49 50 46 50 46 49 49 47 48 45 48 47 47 47 49 48 48
##   [73] 49 50 45 49 49 48 46 45 48 47 49 47 46 49 49 47 47 49 48 47 48 47 50 48
##   [97] 46 46 50 48 46 49 50 49 45 47 46 47 48 49 45 48 49 46 46 46 49 49 46 46
##  [121] 48 49 49 47 46 48 47 50 48 49 46 50 47 45 48 48 49 47 45 49 50 49 47 46
##  [145] 45 48 49 50 49 48 48 50 46 46 45 49 47 50 45 48 45 49 47 49 47 49 48 45
##  [169] 45 49 49 49 46 48 49 48 46 48 47 48 46 48 46 45 46 49 50 47 49 50 45 47
##  [193] 46 50 46 49 47 50 49 48 45 46 45 48 48 47 49 48 48 46 49 48 49 45 48 49
##  [217] 50 49 49 49 47 47 48 47 49 49 47 48 49 47 48 49 48 50 46 45 49 48 47 50
##  [241] 50 46 46 49 46 45 50 45 47 46 50 46 46 49 49 47 47 48 49 47 49 48 46 47
##  [265] 48 47 47 48 50 48 49 46 49 49 49 48 49 46 50 46 50 47 47 48 46 49 47 49
##  [289] 47 46 50 48 48 46 46 48 48 48 47 50 47 49 49 49 48 50 48 50 48 45 46 50
##  [313] 49 49 47 50 48 50 47 49 48 47 45 49 49 48 47 47 48 46 45 49 47 50 50 49
##  [337] 48 46 46 46 49 47 49 45 48 47 49 46 49 48 46 48 46 48 49 50 50 50 47 46
##  [361] 49 49 45 48 50 49 50 50 48 48 48 49 46 49 47 48 47 49 49 48 45 49 48 49
##  [385] 48 47 50 47 47 50 47 45 46 45 46 50 50 50 48 47 49 46 46 46 49 47 50 49
##  [409] 49 46 47 47 50 45 47 49 46 46 47 47 47 47 46 45 48 45 45 45 48 48 48 50
##  [433] 46 47 48 48 48 48 48 47 48 49 47 49 49 46 49 48 45 46 50 46 50 49 47 48
##  [457] 49 46 47 50 45 47 49 48 48 45 47 46 47 50 48 46 45 45 49 45 48 47 46 46
##  [481] 46 48 49 46 47 50 50 46 46 47 46 46 46 46 48 48 46 50 46 46 47 48 47 49
##  [505] 50 45 45 48 50 50 48 46 49 48 45 49 47 49 45 45 47 47 46 47 49 50 46 50
##  [529] 47 45 47 47 46 48 47 47 47 48 47 46 48 47 47 45 49 46 47 48 49 45 46 48
##  [553] 47 48 48 46 47 46 46 45 49 46 46 45 48 46 47 49 49 50 48 45 46 45 46 50
##  [577] 45 47 48 47 46 48 47 48 46 47 49 45 48 50 48 48 48 48 46 48 48 49 47 46
##  [601] 48 49 50 50 46 48 45 49 45 46 47 46 48 48 45 46 50 48 47 46 47 48 47 48
##  [625] 49 49 49 49 47 46 49 48 48 48 49 48 45 46 50 50 48 46 50 47 47 49 50 47
##  [649] 47 49 46 50 49 49 46 47 47 49 46 49 50 48 46 46 45 45 49 45 47 50 48 46
##  [673] 46 47 50 49 45 50 46 48 47 46 46 50 46 49 47 47 47 46 46 50 48 47 45 47
##  [697] 49 45 48 49 48 48 49 47 49 50 46 45 46 50 46 47 48 50 48 50 49 48 50 46
##  [721] 50 45 49 46 48 49 50 47 48 48 47 48 50 48 50 49 46 47 45 45 47 48 47 49
##  [745] 47 50 49 49 45 48 46 47 45 50 47 45 50 50 46 46 46 50 47 47 47 46 46 47
##  [769] 45 49 47 45 50 50 49 47 49 50 46 47 48 48 46 49 47 49 47 49 46 47 46 49
##  [793] 46 49 46 47 49 47 45 48 47 48 48 50 48 46 46 49 47 46 46 45 47 47 46 46
##  [817] 47 48 48 45 45 45 48 45 46 47 48 46 47 49 47 50 47 49 47 46 47 48 50 47
##  [841] 45 46 50 46 47 49 48 49 46 49 45 48 49 46 46 46 46 46 48 48 48 49 46 45
##  [865] 48 48 47 46 48 46 49 47 46 49 48 47 46 46 49 47 45 49 45 47 46 49 50 49
##  [889] 47 45 46 46 49 50 47 49 45 47 49 47 48 47 49 48 46 50 49 49 50 48 46 50
##  [913] 46 47 45 47 46 49 47 46 49 49 49 47 49 47 50 48 47 47 46 49 47 48 48 48
##  [937] 49 48 49 47 49 49 49 46 48 49 45 48 49 48 47 45 45 49 48 49 48 45 46 46
##  [961] 47 47 49 48 48 49 46 49 48 50 46 48 45 46 49 47 45 50 48 46 46 46 48 50
##  [985] 47 47 47 46 48 50 46 46 46 45 48 49 48 48 48 50
hist(bc)

Qua đồ thị ta thấy bánh căn quán bán được mỗi sáng trong vòng 1000 ngày có phân phối đều

1.1.2 Mô phỏng ngẫu nhiên phân phối Poisson

Số hủ sữa chua phô mai bán được mỗi sáng trong vòng 1000 ngày, có tham số trung bình lambda là 10 (mỗi ngày bán được 10 hũ)

scpm <- rpois (n=1000,lambda=10)
scpm
##    [1] 12  9 15 10 11  8 12 18 12 12  9  7 10 11  9 13 10  9 12  8 11  9 12 16
##   [25] 10  7 14  8 15 11 11  9 12 15 15  9 10 12  8 10 13  6  6 10 11  5 12  7
##   [49] 11 12 13 13  7 10 10 13  6  7 10 14 13 10 15 11 14  8 11 12  9  8 11 11
##   [73]  9 11 12 10 12 11  9 18 12  6  6 13 15  9 11  4  8  7 18 11 16 12 11 11
##   [97]  5 15 11  6 14 11 14 12  8 10 11 13  8 10  9  5  9 13 12 13  9  4 10 12
##  [121] 10  9 14  8  5 11  7 12 12  7 13  8 11  8  7  9 17  8 11  6 11 13 12  6
##  [145] 13 12 10 14 10  6 10 12 16  8  9  8 13 15  8 13 15 10 12 10  9 12 11 13
##  [169] 11  9  7  7 14  5 11  6 11 11 10  9  8 14 18 10  6 14  9 13 20 13 15 11
##  [193]  9  5  9 14 11  9  8  7  8 10  9 10  8  7  8  8 10 10 15 11  9 12 11  8
##  [217]  9  7 13 13  9 12  8 10  7 10  7 16  7  3 16 10  7  7 11 12 10  8 15  6
##  [241]  4 15 11 13  9 11 11 11 11  8 16 11  9  7  7 18 11  9 13 10  5 10  9 10
##  [265] 17  8 14  6 10  8  7  9 13 10  7 11  7 16  9 13 18  4 10  9  9 10  7 11
##  [289]  9 16  8  7  8 11 13 12  6  6 11 15  4 11 11 10 11  4  3 11  3  8 16  8
##  [313]  8  3 12 15  8  8  9 12  6 10 15 13  4  6 10 12 10 11  5 15  9  7 11 13
##  [337]  8  9 13  9 10  8 14 13  8  6 11  8 12 13  9 15 14 15 12 12  6  9  8  9
##  [361]  8 10 11  4 12  7 14  7  7 14 11 13 11 15 12  5  8  8 14  9  8  7  8 12
##  [385] 10 11 20  8  9 12 12 12  4 12  9 11 14 15  3 11  8 14  6 15 13  7 16 12
##  [409] 10 13 20 12  9 14 10 11  9 10 14  7  6  7 13 10 14  8 10  9  7  5 13  8
##  [433]  6 10  6 11 13  8 11 12  8 15 12  7 11 11  7 11  8 12 17 11 12 14  8 10
##  [457]  9 12 17 12 11  3  5  8  6 10 15  9 13 12  4 10 12  5 10  7 11  9 11 14
##  [481]  6  9 12 14 13  8 14 10 12  5 10 10 10  8 16  5 13  5 14  7 12 12 12  7
##  [505]  9 13  6 12  8 10  8 11 11  8 16 14 13 12 12 12 15 13  6 11 13 11 10  9
##  [529] 10 10 12 10 11 12  5 10 12  9  8  9  8 11 14 16  6 13  8  9  6  6 14 13
##  [553] 12  9 10 13  8 14 11 15  4  6 10 10  8  6 12 10 14  6  9  6  7 12  8  5
##  [577] 11  8 11  8 11 13 10 12 13  3 13  8  9 13 12 13  9 16  7 16  7  9 14  7
##  [601] 11 12 13  5  5  4 15 12 14  7 12  5  6 13 11  6  9 10  6 10 12 10  9 13
##  [625] 13 13  8 12  8 10 13 12 13 11  9 12  7 13 15 12 10  5 13  8  5 11  8 12
##  [649] 12 10  8 19  9  8 15  8 12  5 12 12 13  5  9 10 11 14 13  6 11  7 10  6
##  [673]  7  3  8  9 11  7  9  6  7 14 14 10 15  7 12 11  7  5 11  8 16 12 10  9
##  [697]  6 12  7 11  7 10 10  6 11 15  5  9 19  8 11 12 10  7  4  8  9 17  9  7
##  [721] 12 11 12 14 12 13 11 11  7  9  8 13  8 11 11  7  3  9 14 12  8  7  8 12
##  [745]  8  7 10 10 10  8 10 14 13  8  9  8 10  9 12  7 10 10 10  7  7 16 11 17
##  [769]  9 12 11 12  7  6 11  7  8 11 11 11 14  9  8  9  9 11  7  6  8 13 11 11
##  [793] 13 11 10 12 15  9 11  9 10  8  7 13  9  3 13 11  5  5 11  6  6 11  7  6
##  [817] 11 14 11  9 10 10  8 12 16  6  6  8 11 19  9  9  7 10 11 18 10 12  6 12
##  [841] 12 11  6  9  3  7  5 14 10  9  5  7 12  6  6 17  9  8  8  5  7  7 12 15
##  [865] 13 12  7 12 11 12 11  9 14  3 12 12  7 19  9 13 10  6 12  7  5 14 16  7
##  [889] 11 17 10 15 11 10  7  9  4 12 11 10 14  7 15 10  7 12 16  5  9 13  9  8
##  [913]  7  9 11 10 10  6 13  8  7 12  6 10 11 11  7 10  8  9 16  5 21  9 12  8
##  [937] 10  8 13  8 11 18  9 17 13 11 12 13  9 13 14 13  8 11 12 14  6  9 11  4
##  [961]  8 13  9 16  7 18  8 12  9 10  9 14  8  9  9 11 12  1  3  7  7  9  5 15
##  [985] 10  5  5 10 12  7 13 13 16 10 19  8 11 10  9 10
scpm <- table(scpm)
scpm <- as.data.frame(scpm)
ggplot(scpm,aes(scpm,Freq)) + geom_col(fill='red') + geom_point()+ geom_line(aes(as.integer(scpm),Freq), color='green')

Qua đồ thị ta thấy số hủ sữa chua phô mai mà quán bán được mỗi sáng trong vòng 1000 ngày có phân phối đều với lambda là 10

1.1.3 Mô phỏng ngẫu nhiên phân phối chuẩn

Số ly sữa đậu nành bán được mỗi sáng trong vòng 1000 ngày, có trung bình =10 và phương sai = 5

sdn <- rnorm (n=1000, mean=10,sd=5)
sdn
##    [1] 11.60984915 12.38039371 15.58601460  8.47991719 14.44480130  8.40241351
##    [7]  1.37815521 20.32817608  4.72729858 15.95838422  3.82115992  5.31541709
##   [13] 18.58358892  8.41980546 13.26273309  7.35502680 -3.60710398 12.15742374
##   [19]  8.74049950  9.25364859 11.42899729  9.26330213 16.74170916 10.16653434
##   [25] 14.28564456  5.91982020 15.28775025  8.78605728  7.76174855  4.23705015
##   [31]  3.24321011 16.03617973  6.25337461 11.03696675  2.66824329 18.09335916
##   [37] 11.22871576 17.70652119 18.04806053  8.27426732  7.67246530 19.16562336
##   [43]  3.53419601  8.31063345  5.44034342 10.98523485 11.40268471  2.59252453
##   [49] 14.90940732 15.86762136  6.66042588  6.68958043  9.31770072  1.84049629
##   [55] 10.51683568  4.57325689 12.42523668  9.49100815 11.25176643  3.03173001
##   [61]  1.48547626 10.14425994 18.77449121  8.25521400 10.81530412 14.69278715
##   [67]  3.74354095  4.92334566 15.52925325  9.17499170  4.62802015  7.63962707
##   [73]  6.09990439 11.41061371  6.78559330 11.16248942 21.45558153  5.59378696
##   [79] 11.66758284 10.62977846  9.55847469  4.81422495 11.92057314  7.55077437
##   [85] 12.79512993 20.61813058 15.49703322 14.11545715  2.02563763 10.49875812
##   [91]  5.53762207  8.37659235 10.63095888 13.51647319 11.74193718 12.74619015
##   [97] 17.97959331  5.85441347  9.53268134 12.84699430 20.38239967  5.12780526
##  [103] 13.70175746 12.91293702 17.04692037  7.11598267  5.63317502  2.43234883
##  [109] 16.46362738  3.42300340  1.44884956  5.12361533 13.98573533 12.03365068
##  [115] 10.16055668  8.57384461 14.62081071  2.16754722  9.26505929  5.94114413
##  [121]  6.34526298 13.66784125 16.59080987 10.19177955  7.01638924 15.96357181
##  [127] 10.87768101  8.52865217  4.74874449  8.98022956 10.33220777 10.24863628
##  [133] 16.99081162 14.12551867 12.40021468  3.19839569 -3.24048778  7.00562827
##  [139]  8.23960378  5.29815726 12.75841757  8.36837637  5.32632633 11.76696163
##  [145]  5.26676011  3.06398362 10.18559224  2.02815735 10.66223840 14.34135783
##  [151] 19.18676944  1.32802218  8.60188024 21.83928116  6.26041358  8.41380094
##  [157]  8.48339167  8.42373880  2.06744998  3.42673291 15.85772226 15.21870131
##  [163] 10.02819161 14.54332512 13.64309202  7.40779620  8.09786126  7.57775526
##  [169] 12.49207007 -0.67014779  7.38838637  9.88060521  8.17742198  2.84294346
##  [175] 10.60347966  6.60720937  8.19400342 15.72546688  9.10212469  5.17971079
##  [181]  3.74330081  7.81945025 -1.78161972 11.41974381 13.17300154 16.43094284
##  [187] 10.90424280  2.47165123 18.46457064 13.58388704  2.53416281 20.28038045
##  [193] 10.54401491 21.10296428 12.20676371  6.78661911 11.43798269 16.48553490
##  [199] 15.47039946  8.61283372  0.02135244  8.52039907 -1.15723469 16.80666551
##  [205] 13.78433918 13.48266991  8.35859931 13.37621337  8.29421767 17.48299030
##  [211] 15.82044994 15.75054727  8.17183338  6.21918938 -0.74691570  0.67031502
##  [217]  6.88411622 12.78053437 15.30884202 10.94416354  5.55807589 14.67629748
##  [223] 14.86974353 12.64501429  8.69703684 10.35662986 12.47078412 12.23586086
##  [229]  7.51267904 15.32722100  0.73065835  5.01516291 13.55439000  7.20512019
##  [235]  7.16720661 15.18528735 11.43539890  2.48795688 12.81358766 12.76846682
##  [241] 11.42392725 17.91682212  8.44901679 13.43072337 10.91480121 11.39303994
##  [247] 20.56882028  8.87966233  1.81373886 14.40243383 13.62210029 10.09541501
##  [253] 17.16284415  3.75612121  0.56701347 11.33723542  8.13516413 12.20638564
##  [259]  3.47647785  0.99235647 17.53800172  7.65357199 10.32669338 -1.10282972
##  [265]  4.02277374 12.51249991  5.80183708 11.22280542 15.41461134 15.16204889
##  [271] 17.02269646 12.66243098  5.91645305 12.34619858  5.74546071 10.41728328
##  [277]  7.28641634 12.46024570 15.23398424  3.65971643 10.23446166 10.74324076
##  [283]  5.01853060  7.63630388  8.95445361  7.21648397 10.16621649 10.87551694
##  [289]  8.69679332  6.52498573  2.52421915 17.87960249  7.12725526 11.83996710
##  [295] 15.47463766  7.81290020 19.91725621 12.37386459 17.21900533 12.38189601
##  [301] 12.60447985  8.38708724 10.92042699  9.45830248  4.68426337  4.67855993
##  [307]  4.57760389  9.31710381  3.61771002 15.31139466  9.50850613 11.53413812
##  [313] 11.31508880 11.75466171 11.09242167 17.97399116 10.21244314 20.03180917
##  [319] 11.70741245  9.12755130 12.77926013  1.28268045  2.84211240  7.04436292
##  [325] 21.11020223 10.81926104  8.33641472 14.51299769  8.15005861 12.19673287
##  [331] 11.70793101 15.51372984  6.09500724 12.14170170  3.84049496 14.97483802
##  [337] 10.57245518 23.67106115 11.63439808  9.07002050  3.53096179 11.56295533
##  [343] 14.68656316  7.59648911 14.26893917  4.90210560 11.29013875 15.35436767
##  [349]  5.62596061  9.74427034 12.35261673  4.87385696  6.66388050 15.49963164
##  [355]  9.99122432  3.67509748 11.00982392 11.99368946 11.14995507 13.97013220
##  [361] 16.85287930  8.27152998  5.09984323 10.38944961  7.96534050  6.74360530
##  [367]  8.05259477  4.23578818 15.76812944 10.79939069  8.57710600  8.79857246
##  [373]  1.59761974  6.59234794  9.43522724 20.13207079  1.98606010  6.97634734
##  [379] 17.91563735 13.76061642 10.53628886 10.62633813 16.52125528 12.51275930
##  [385] 12.40796125  6.63913084 17.92855712 10.23922508 11.21720532  7.84216017
##  [391] 10.06013605 11.40540768 13.67580348 10.93488318  1.17058829 14.11187489
##  [397] 11.61294963  5.21838603  8.81019902 18.66369482 11.56146533  6.48296151
##  [403]  8.75421671  6.22804928 10.07746489 18.00703388 12.85154981 14.57722984
##  [409]  9.64637174  9.29449474 13.13572607  7.93453028 15.34115835 10.58438505
##  [415] 12.39879391 20.75676773  7.87137505 14.05862016  7.31082785  0.93707857
##  [421] 13.96366456  8.81141939  9.26564770  9.16764777 15.90278089 13.72903676
##  [427]  5.63019686  3.27226523 10.69315989  5.44982673 11.64275469  7.31473261
##  [433]  7.99619934  7.79320772  4.47285444  8.59210173  1.95842710 11.13664894
##  [439]  9.96928878 12.95146673 11.62333819  0.87369741  7.92238993 11.92912871
##  [445]  6.88596566 15.80839697  6.69312889  7.68491453  7.60409965 15.08108638
##  [451]  7.54724785  4.17856285 10.41515920  4.74384313  8.46271787  9.55889571
##  [457] 16.55618278 12.06730265 14.09991179 23.49771839  8.94548699 11.22949785
##  [463] 15.15897431 13.03402652 15.57996109  8.22765773  7.55874013 13.36160518
##  [469] 11.57597649 20.49568888 13.42558675 10.67774925  8.61857661 19.03288126
##  [475] 15.04538577  5.04376241 20.19900124  9.18638032  7.00871123 17.15548351
##  [481]  1.21331929 12.40904726  7.40599332 11.27305349  4.31158731  9.15568804
##  [487] 14.50121985 17.12682308 10.44963651 10.18401080  9.62351756  0.99766081
##  [493] 14.87124299  6.85917390  3.62488020  8.15725179 13.81023292 16.10697345
##  [499] 13.08693368  8.42380490  8.58840217  5.78278069 13.38302585 10.22006047
##  [505] 13.14572348 18.17341487 26.48792364 19.67026593  8.63842711 10.50237296
##  [511] 15.25573549 11.15885522  7.78983836  5.66651278 -1.07398924  7.96496883
##  [517] 15.25846805  9.42382863 14.48222217  5.60516330  6.67913761 10.86408185
##  [523]  7.62109010 10.85597432  5.71495845  2.44885832 10.96475551  2.95134179
##  [529] 22.69193215  8.97781651 16.12909200  3.33476077  8.91973550 15.16633722
##  [535] 12.65434572 15.06752429  3.02949171  6.43867291  4.80110305  6.18839393
##  [541] 11.16050798 16.73676525  9.49326628  8.26630961  9.60297335  9.63975733
##  [547]  0.54496570 17.88341070  9.16900938  8.94025794 13.76396775 11.09844031
##  [553] 14.02181262 13.36069731  8.12736338 16.99063876  5.61494135 10.86426925
##  [559]  7.21374266 16.08636058 12.12329435  7.87764421 13.64749353  7.15524188
##  [565]  3.45436539  4.58490738 19.16535004  8.82726045 13.00574460 15.97119096
##  [571] 11.45197160  6.23206633  6.87813666  7.27660365 16.70259746  0.16723119
##  [577] 11.90022665 18.35789792 -1.30945988 17.96546303 13.22613802  6.14373520
##  [583] 11.40443901 19.60220591 11.31648535  3.74265473 17.69003326 10.08068337
##  [589] 14.70937154  5.61299075 17.44588649 10.15686164  8.66070470  8.46255443
##  [595]  7.30293802 13.08390625 11.27849327  9.93099212  5.11015229 12.93849493
##  [601]  7.86973935 13.63922740  7.62881233 -1.86325171  7.35251334  2.96743520
##  [607] 10.72345468  8.16268603  9.89631507 -2.53328369  7.85151510  8.48513539
##  [613] 15.98834502 11.68244657 15.64514676  3.64315403 10.69265535 12.89163161
##  [619]  6.84929840 11.13002028 10.84177010 11.85536781  6.58246620  8.10115957
##  [625]  7.72179616 20.95214206 13.01024473 10.01664810 17.16277858  3.45592764
##  [631]  8.15898288 13.29549824 13.02065976  6.43504922 10.15021660 16.06451865
##  [637] 11.95752376 20.27520564 19.85740886 12.06363853  6.42193508 14.76427016
##  [643] 15.18882025  8.24175434  6.64646680  4.05527941  8.59925489 17.73747325
##  [649] 13.75176408 14.49959218  2.90913697  7.56144627 12.12270606  3.04930205
##  [655]  9.12598719 11.10793428 10.39734079  9.39730122  9.59237840 12.85903474
##  [661]  7.70296720 16.05696560  6.29969717  0.18586234  8.16324765 10.99463257
##  [667]  6.12557548 20.09990005  9.56769375  3.91032439 10.17522959  8.47262565
##  [673] 12.71968019 12.41192708  8.01181155 15.23213644  9.71472185 14.95220221
##  [679]  4.95922810  8.09103687 -0.72455519  5.58040793 13.09216584 12.05132528
##  [685]  3.60026664 16.92895641 13.97866274  0.03841828 12.99691996  1.45672178
##  [691] 14.77050098  9.29817789  8.73385455  8.79997505  9.56053768 15.86130399
##  [697] 11.64352341  7.35523921 13.41511968  4.43087383  9.07972744  7.67378900
##  [703] 17.89533200  2.21069966 21.72307534 14.41406680 16.45119611  4.20436533
##  [709] 10.21360930 -5.53427481 11.42626466  4.03563512 17.64369526  7.20705015
##  [715] 13.51859423 11.39651909 10.20959877 11.29408022  8.58560614 11.98989813
##  [721]  2.89435949 13.47422193 12.97762437 10.45996183  2.63434378  4.48582090
##  [727]  4.58202745  9.73108223  6.73931211  3.81514104 10.13995704  8.36552858
##  [733]  9.80902586  5.66287272 19.26470686 16.20385906 13.24559670  1.74111370
##  [739]  4.51574205  1.23983730  7.73567064  4.46340689 10.22630045 22.17130730
##  [745]  7.58598715  6.57649231 15.47877977 14.53278269  9.07345486  7.35796034
##  [751] 10.39245755  2.26662341 12.14692665 13.29024684  1.12430890  8.73825683
##  [757] 14.24972549 10.55822071  6.91134747  7.01337980  7.31967289  9.90692639
##  [763]  8.92811429  9.07202221  2.37146177  3.45290341 19.17257724 11.06709635
##  [769]  4.81574768 10.31901018 12.59563462  4.60479249 12.83493903  6.79098913
##  [775] 17.88131878 13.73927158  1.83058897 18.26942346  9.65367211 16.14567988
##  [781]  7.56617488  9.19020232 12.35026396  9.36862663 15.27997024 11.19439836
##  [787] 11.29328026 -0.36737890 10.61239965  9.13176320  6.73058044 11.49584467
##  [793] -1.00818404 14.00763944 15.60976549 12.87940344 11.91433799  8.92446366
##  [799] 13.81617995 11.55787995 12.49133328 10.81101988 12.72062080  9.52979924
##  [805] 18.24300945  3.42259394  5.11954494  1.10437643 12.02248856  0.59366480
##  [811] 16.85650837  2.53619712  8.17110075 14.99745667 16.10543988 14.75531711
##  [817]  7.30779700 14.74696790 10.70078534 11.42274926  5.27379310  4.75689937
##  [823] 10.45394108 13.13709850 13.66159280 12.43625749 12.98628882 21.00006508
##  [829] 10.05806572 10.87295653  8.85128594  9.49366073 13.31274599  8.57077170
##  [835] 10.72069687  8.99031631 11.57421422  8.86743110  1.35607121 18.59419345
##  [841]  6.33484740 10.83052189  7.79220991 17.04439895  2.18599517 12.44126133
##  [847] 11.15217553  9.52214726 14.82400549 21.30290564  8.95174437  6.43982345
##  [853] 12.69853437  5.37569564 17.26960335 13.86412044 15.03683625 14.02879919
##  [859] 13.80887442 12.97709131 19.20362645  8.75928035 11.56950854  8.74298799
##  [865]  5.75603001  0.74749711  7.91028061  3.69946981 -0.75968296 -4.32823520
##  [871] 14.91794240 -1.07622196  2.13070129 10.07977397 12.66712851  7.12088233
##  [877] 11.70444949  7.38991759  7.85983041  8.09235958  5.08101933 -7.88261903
##  [883]  3.19542987  4.19136127 11.55058033  6.78148240 -6.82516341 12.73138280
##  [889]  7.21618885  6.85466385  0.14392060 11.77824575  9.55502323  8.06895470
##  [895] 14.60057107 11.95604999 17.89314197  8.38611371 11.98100515 15.02328153
##  [901] 14.72561449  5.98393955  9.45631728  8.35726549 15.52742458  6.28139218
##  [907] 12.20144561  7.91163771 -1.04685962  8.36768801  6.22504735 11.53385420
##  [913] 18.54061783 15.49390023  5.75653924 16.11602577 13.83203804 10.68941864
##  [919] 10.00468851 10.77493023  7.05369465 11.64368242 12.44039783 12.81684318
##  [925]  8.81502123  8.90437602  9.85250486  6.24772619 -7.16200827  1.89591541
##  [931] -1.95662679  4.44142753  7.24546123 12.21488128  2.85856681 10.67211691
##  [937] 13.49179373  5.34178711  7.91459769 17.18906973  5.99524797  7.31595088
##  [943] 14.82334584 10.13853769  0.08795673 10.83368055 13.72625148 12.52619107
##  [949]  4.61370221 13.56489996 15.92029162  8.55543121  9.69131963 16.60793937
##  [955] 12.49052167  3.62147798  3.44825671  8.61029873 17.29081153  5.64073343
##  [961] 17.44989045  5.56374560 12.19186172 10.48023173  7.05396124  1.48880984
##  [967] 18.84927371  5.85996722 12.52237982 19.16981846 15.07958931 22.26130132
##  [973]  2.32511148  9.19301307 11.46924502 15.10291172 12.66799468 16.07251269
##  [979]  8.99379072 14.24844546 -3.76621936 11.11102718  6.55357152 10.42093984
##  [985] -3.78059602 20.70219013  6.43862775  5.83531407 11.59183229  6.92025437
##  [991]  7.67654213  6.70906563 14.20977949 15.59878440  9.39752965  8.08783637
##  [997]  4.36000363  9.03091390 11.11561033  7.38394324
hist(sdn)

sdn <- as.data.frame(sdn)
ggplot(data = sdn, aes(sdn)) + geom_density(color = 'purple')

Qua đồ thị ta thấy phân phối số ly sữa đậu nành có trung bình là 10 ly và độ lệch chuẩn là 5 với xác suất là hơn 7%

1.1.4 Mô phỏng ngẫu nhiên phân phối nhị thức

Số ly trà tắc xí muội bán được mỗi sáng trong vòng 1000 ngày, chọn ngẫu nhiên 20 người, trong 20 người đó có 20% người gọi trà tắc xí muội

txm <- rbinom(n=1000,size=20,p=0.2)
txm 
##    [1]  2  5  2  2  2  2  4  4  3  2  3  3  8  2  5  3  5  6  5  2  6  6  5  3
##   [25]  3  6  2  3  8  3  4  4  6  4  3  6  4  5  3  5  6  1  4  1  5  3  5  4
##   [49]  3  2  2  3  7  1  2  7  3  4  5  2  2  4  7  4  4  3  2  5  5  4  7  6
##   [73]  4  5  4  4  2  4  5  7  3  5  5  5  2  2  5  5  4  2  5  2  3  6  4  5
##   [97]  3  7  7  7  5  2  7  5  4  3  2  4  2  6  3  4  3  1  5  0  5  4  3  9
##  [121]  2  8  4  4  6  1  6  4  5  1  4  2  3  2  3  2  3  3  4  6  4  3  2  2
##  [145]  4  5  4  8  3  4  2  3  4  6  6  3  4  6  6  6  4  6  5  5  4  7  4  3
##  [169]  4  3  0  4  8  3  3  2  4  4  2  5  7  5  4  2  4  1  3  4  5  4  6  4
##  [193]  5  2  4  5  5  7  3  3  3  2  3  6  4  4  3  4  2  4  6  7  6  3  2  4
##  [217]  1  3  3  5  1  6  6  4  5  5  1  0  3  3  4  2  4  5  2  2  3  2  1  2
##  [241]  4  1  5  4  2  5  4  5  6  5  3  5  3  4  4  5  4  4  4  4  2  4  9  5
##  [265]  5  3  2  3  3  3  9  4  6  8  4  7  2  4  7  5  3  3  2  3  4  4  6  1
##  [289]  3  3  3  5  2  2  0  4  3  3  5  3  6  4  5  6  6  4  7  4 10  3  4  3
##  [313]  3  2  1  6  6  7  3  3  8  2  3  4  3  1  5  4  4  6  4  3  3  4  8  3
##  [337]  5  3  8  3  3  6  3  4  7  4  3  5  3  1  2  6  7  2  4  1  6  6  4  4
##  [361]  1  7  6  4  5  5  3  4  2  1  4  2  3  1  4  3  4  4  5  4  2  2  1  6
##  [385]  8  2  4  3  5  7  0  6  6  2  4  2  2  4  3  5  0  4  5  3  6  6  4  2
##  [409]  2  3  2  6  4  5  7  4  3  6  2  1  2  3  3  3  5  4  3  6  6  3  4  2
##  [433]  8  4  6  5  1  5  3  5  2  6  4  5  3  5  8  3  1  7  1  1  2  6  2  3
##  [457]  5  5  5  3  4  2  3  6  3  4  4  5  4  3  2  3  5  6  3  2  2  6  4  4
##  [481]  6  1  6  2  4  5  7  3  3  3  4  3  3  6  6  6 10  4  3  6  6  5  6  3
##  [505]  4  5  2  5  3  2  2  3  2  4  7  3  6  3  4  1  4  3  4  5  7  4  1  3
##  [529]  4  2  6  4  5  1  6  4  6  1  2  5  7  3  4  2  2  3  7  3  2  2  8  4
##  [553]  2  2  0  6  3  5  2  4  2  1  5  5  5  4  2  2  5  5  3  5  1  6  7  2
##  [577]  7  3  3  9  2  4  5  4  5  4  5  1  2  2  5  6  1  1  4  5  4  8  4  5
##  [601]  4  5  5  4  8  6  5  5  2  4  5  8  6  5  7  4  3 11  1  6  4  2  6  5
##  [625]  4  3  5  5  5  3  3  3  6  2  3  8  6  6  5  4  5  4  3  4  3  4  5  3
##  [649]  6  4  5  2  3  0  3  1  3  3  7  4  5  4  4  4  2  5  3  6  7  4  3  4
##  [673]  4  3  3  3  6  2  2  3  3  7  5  2  5  7  2  4  4  4  5  8  2  1  5  5
##  [697]  2  9  4  3  3  6  4  2  4  2  4  5  3  4  5  7  4  5  0  5  6  4  2  2
##  [721]  4  6  4  5  6  0  9  2  3  3  6  7  2  4  5  3  5  6  5  0  7  5  7  2
##  [745]  4  3  4  4  3  5  1  2  5  5  1  7  4  6  5  5  3  5  5  1  5  3  4  3
##  [769]  4  5  3  5  3  7  3  4  5  3  1  2  5  1  2  5  3  3  3  7  0  4  4  3
##  [793]  3  5  5  4  3  4  2  6  2  3  3  5  2  7  6  3  7  6  3  2  3  2  4  5
##  [817]  1  4  5  1  6  3  4  4  4  5  5  7  1  3  5  3  2  5  5  4  3  4  6  6
##  [841]  6  4  3  6  3  5  4  3  3  5  8  4  1  5  6  2  3  2  3  2  5  4  5  6
##  [865]  7  3  4  7  5  5  4  4  5  3  4  4  6  2  3  4  5  5  4  3  5  3  5  3
##  [889]  2  3  5  4  3  1  3  3  5  5  4  3  2  4  3  4  3  1  2  2  5  1  5  1
##  [913]  5  2  3  1  8  7  3  5  6  3  7  4  6  3  3  2  5  3  4  1  3  1  4  4
##  [937]  7  2  3  5  6  4  4  5  5  5  5  5  5  3  3  1  3  2  4  2  6  4  4  3
##  [961]  4  6  0  3  5  2  2  4  2  3  6  4  2  4  2  5  7  6  3  5  1  4  5  3
##  [985]  3  3  3  4  4  6  7  3  3  7  3  4  1  4  1  3
table(txm)
## txm
##   0   1   2   3   4   5   6   7   8   9  10  11 
##  13  59 140 212 213 176 105  53  20   6   2   1
hist(txm)

20 ngày ngẫu nhiên trong vòng 1000 ngày xác suất có 3-4 ly trà tắc xí muội bán được là cao nhất

1.1.5 Mô phỏng ngẫu nhiên phân phối mũ

Số ly trà đá được mỗi sáng trong vòng 1000 ngày với tỷ lệ là 20%

td <- rexp(n=1000,rate=0.2)
td 
##    [1] 5.117871e+00 8.216068e-01 5.647892e+00 7.510597e+00 1.602531e+00
##    [6] 2.484731e-01 3.019466e+00 7.019193e+00 9.683038e+00 4.906504e+00
##   [11] 4.218387e+00 3.751612e-01 3.857650e+00 1.601318e+01 5.576374e-02
##   [16] 1.505903e+01 5.353562e+00 3.283836e+00 1.229349e+01 5.509594e-01
##   [21] 7.729038e+00 8.390292e+00 1.197824e+01 1.221500e+00 4.003226e+00
##   [26] 6.099912e+00 6.734872e+00 4.732771e-01 2.296419e+00 7.386820e+00
##   [31] 6.834835e-01 4.565683e-01 2.440710e+00 6.066818e+00 2.287539e+00
##   [36] 4.278196e+00 3.364704e+00 3.086344e+00 3.290915e+00 4.540250e+00
##   [41] 1.183058e+01 1.021268e+01 1.461159e+00 1.161320e+01 8.559511e-01
##   [46] 3.578535e+00 1.625433e+01 1.840977e+00 1.161234e-01 6.525448e+00
##   [51] 5.680049e-01 1.203225e+00 9.252250e+00 8.865037e-02 7.896067e-01
##   [56] 4.448610e-01 6.630932e+00 4.446331e+00 2.248072e+00 3.773817e+00
##   [61] 4.513562e-01 6.472037e-01 2.447461e+00 1.527477e+00 8.123112e+00
##   [66] 1.952522e-01 5.750393e+00 3.711768e+00 9.728466e+00 1.394973e+01
##   [71] 2.969651e+00 4.190571e+00 5.802290e-01 6.526373e+00 2.261060e+00
##   [76] 1.082136e+00 1.354242e+01 1.851540e+00 2.400293e+00 2.351386e+00
##   [81] 5.207614e-01 6.929443e+00 1.654632e+01 3.013150e+00 9.008119e-01
##   [86] 1.581818e-01 1.474618e+00 2.101553e+00 4.175041e+00 1.819022e+01
##   [91] 7.402250e-01 1.715246e+00 9.008881e+00 5.021550e+00 4.732803e-01
##   [96] 5.117989e-01 1.867295e+01 9.777649e-01 5.555442e+00 3.432609e+00
##  [101] 1.383767e+00 5.880699e+00 5.998306e-01 1.945362e-01 1.554350e+01
##  [106] 2.564267e+00 2.965449e+00 2.132307e+00 1.940687e+00 7.781648e+00
##  [111] 7.137676e-01 1.025361e+01 6.034306e+00 5.814666e-01 6.879952e+00
##  [116] 3.705562e+00 1.001141e+00 3.775945e+00 1.497909e+01 3.564776e+00
##  [121] 1.411133e+00 1.623782e+00 2.169404e+00 2.069621e+00 6.897728e-01
##  [126] 4.736507e+00 8.443203e-01 4.674572e+00 6.364700e+00 9.244327e-02
##  [131] 3.925875e+00 4.035998e+00 1.782524e+00 9.298516e-02 1.128993e+00
##  [136] 1.981003e+00 7.147797e+00 6.043421e+00 6.405641e+00 7.100634e+00
##  [141] 9.159533e-02 9.988072e+00 5.020863e+00 5.845219e+00 1.452164e+01
##  [146] 3.219900e+00 1.212202e+00 1.870061e+00 6.273617e+00 5.072431e-01
##  [151] 4.952029e+00 3.118456e+00 9.793613e-01 1.946526e+00 3.208580e+01
##  [156] 6.854982e-01 1.267740e+00 4.413256e+00 3.614333e+00 2.839598e+01
##  [161] 2.235533e+00 5.155365e-01 3.790963e+00 6.510748e+00 1.440377e+00
##  [166] 3.738135e+00 5.975959e+00 3.983999e+00 7.175213e-01 6.739613e+00
##  [171] 3.260798e+00 9.292949e-02 2.958487e+00 1.468333e+00 8.098988e+00
##  [176] 3.601898e+00 9.707647e-01 3.320172e+00 1.636140e+00 6.763516e-01
##  [181] 1.150481e+00 6.569495e+00 4.905950e-01 1.136745e-01 1.834778e+00
##  [186] 4.080153e+00 8.006890e+00 6.261840e+00 3.905548e+00 5.246778e+00
##  [191] 3.956704e+00 3.326717e+00 4.392505e+00 1.872223e+01 7.838942e+00
##  [196] 1.227486e+01 4.593279e+00 6.971109e-01 2.105000e+00 1.577427e+00
##  [201] 4.854962e+00 1.662629e+00 7.266287e-02 1.014494e+01 7.762708e+00
##  [206] 9.940788e-01 1.256082e-01 1.560167e+01 4.994298e+00 5.369969e+00
##  [211] 1.644988e+01 3.369936e-02 2.976458e+00 8.905540e+00 4.913320e-01
##  [216] 2.657989e+00 3.474604e-01 1.091181e+00 5.291286e+00 3.002912e+00
##  [221] 1.591642e+01 1.929369e+01 6.293219e+00 2.743177e-01 4.644116e+00
##  [226] 7.665528e+00 1.043552e+00 2.646419e+00 3.436887e+00 5.169812e-01
##  [231] 1.724955e+00 5.822259e+00 5.161890e+00 4.198397e+00 5.220465e-01
##  [236] 1.251110e+00 6.612273e-01 6.315949e+00 1.417889e+00 3.236114e+00
##  [241] 2.240650e+00 1.954354e+01 1.934080e+01 3.089559e+00 2.944218e+00
##  [246] 2.978643e+00 7.294703e+00 2.569678e+00 2.918928e+00 8.002705e-02
##  [251] 1.353080e+01 3.038317e+00 8.535359e-02 1.044166e+01 1.702177e+00
##  [256] 1.911507e+00 4.563348e+00 7.056470e+00 2.150807e+00 4.077894e+00
##  [261] 5.603409e+00 9.509546e+00 4.665703e+00 1.098563e+01 5.164657e+00
##  [266] 1.815504e+01 1.008376e-01 1.936235e+01 2.957151e+00 3.968880e+00
##  [271] 6.979916e+00 4.626124e-01 8.170177e-01 2.121283e+00 3.317126e+00
##  [276] 1.994917e+00 8.918995e+00 1.136658e+01 1.123736e+00 1.688818e+00
##  [281] 5.971003e+00 4.679161e+00 3.007893e+00 1.725184e+01 2.380788e+00
##  [286] 1.011960e+00 1.375801e+01 1.735497e+00 4.095913e+00 1.107053e+00
##  [291] 1.436158e+00 3.402743e+00 2.237677e+00 7.290223e+00 1.451083e+01
##  [296] 6.437167e-01 6.943459e-02 1.705162e+00 1.074130e+01 1.332350e+01
##  [301] 1.471561e+00 4.123166e+00 3.963490e+00 9.771127e-01 1.476521e+00
##  [306] 6.911592e+00 4.598730e-02 2.629574e+00 6.476401e-02 3.835840e+00
##  [311] 1.554988e+00 3.861120e-01 2.958915e-01 2.430737e+00 1.015220e+01
##  [316] 5.635460e+00 1.676994e+01 8.069916e+00 7.002639e-01 1.957160e+01
##  [321] 1.325114e+00 4.952975e+00 2.419179e+00 9.354267e-01 9.372877e+00
##  [326] 5.119343e-01 5.077025e+00 1.196538e+00 1.075488e+00 4.726571e+00
##  [331] 6.991617e+00 2.076928e+00 7.144738e+00 7.175284e+00 1.110945e+01
##  [336] 7.871186e+00 1.666525e+00 1.696260e+00 1.158166e+01 5.331238e+00
##  [341] 1.343082e+01 3.424640e+00 1.972432e+00 2.342627e+00 3.027264e+00
##  [346] 6.670613e+00 1.598339e+01 2.273531e+00 2.623510e-01 1.100216e+00
##  [351] 1.740958e+00 2.469496e+00 4.403364e+00 2.345412e+00 1.912243e+00
##  [356] 6.419186e+00 2.900438e+00 6.949959e+00 1.331492e+01 2.555208e+00
##  [361] 7.412922e+00 1.504562e+00 2.279795e+00 4.382197e+00 1.776581e+00
##  [366] 7.380545e-01 6.833242e+00 9.597230e+00 2.466994e-01 1.775629e+00
##  [371] 2.928042e-01 1.052542e+00 2.315731e+00 9.726772e+00 7.982371e-01
##  [376] 3.033222e+00 7.145110e+00 4.139267e-01 5.981243e+00 1.538311e+01
##  [381] 3.258310e+00 4.252449e-01 2.302155e+00 6.827483e-01 4.257984e+00
##  [386] 1.240087e+00 1.131393e+01 1.817788e+00 7.226862e-01 6.736547e+00
##  [391] 4.951821e+00 6.203528e+00 1.492723e+01 9.142761e-01 6.702922e-01
##  [396] 8.702029e+00 2.062503e+00 1.206980e+01 5.270698e+00 2.033272e+01
##  [401] 1.955120e+00 1.144777e+00 1.447218e+00 5.558150e+00 2.626236e+00
##  [406] 5.556940e+00 1.549958e+00 1.540631e+01 5.834490e+00 1.324908e+00
##  [411] 8.596639e+00 1.836102e-01 6.478332e+00 8.134464e+00 6.565316e+00
##  [416] 4.069798e+00 5.678005e+00 5.384883e+00 3.062819e+00 2.077608e+00
##  [421] 2.627054e+00 7.555109e+00 8.613313e+00 2.529779e+00 1.574234e+01
##  [426] 4.808466e+00 5.025296e+00 9.085828e+00 3.661945e+00 1.206974e+01
##  [431] 1.709855e+01 5.845822e+00 1.320547e+01 2.001907e+00 4.155906e+00
##  [436] 1.137024e+01 7.400850e-01 1.766875e+01 9.205750e+00 9.132726e-01
##  [441] 1.309954e+00 7.025849e+00 2.883302e+00 8.740019e+00 5.973256e+00
##  [446] 5.051054e+00 3.732872e+00 8.602717e+00 1.892507e+00 1.545643e-01
##  [451] 5.620865e+00 1.964911e+00 2.701369e+00 6.962234e+00 5.501010e+00
##  [456] 5.791830e-01 2.702448e+00 4.264209e+00 2.208173e+00 1.445483e-01
##  [461] 7.491592e+00 1.337074e+01 1.802181e+00 7.625532e+00 1.261349e+00
##  [466] 1.369403e+00 2.745607e+01 5.095037e-01 1.781283e+00 3.235884e+00
##  [471] 4.346864e+00 9.745435e+00 5.828225e-01 1.367836e+01 1.647448e+01
##  [476] 6.273148e+00 5.708750e+00 8.055144e+00 2.196377e+00 1.710850e+00
##  [481] 1.341620e+00 3.793025e+00 1.524879e+00 6.135505e+00 3.473438e+00
##  [486] 4.608303e+00 2.501241e-01 4.179127e+00 6.096124e-01 2.890256e+00
##  [491] 1.780621e+00 6.815664e+00 1.978616e+00 4.260548e+00 2.080657e+00
##  [496] 3.108157e+00 1.941810e+00 8.435766e-01 1.515411e+01 7.590954e+00
##  [501] 6.672078e+00 1.906228e+00 2.744055e+00 2.385262e+01 3.083136e+00
##  [506] 5.883023e+00 2.155572e+00 9.773227e+00 2.530080e+00 1.534277e+00
##  [511] 8.859902e+00 5.030758e-01 4.739816e+00 2.103926e+00 4.639526e+00
##  [516] 8.121720e+00 1.690450e+00 8.638704e+00 2.337038e-01 6.927109e+00
##  [521] 8.342956e-01 4.651236e+00 3.681692e+00 2.874260e+00 2.453629e-01
##  [526] 3.325383e+00 3.424183e+00 5.083896e+00 2.153966e+00 1.402334e+00
##  [531] 6.502353e+00 5.800647e+00 1.254611e+00 1.465600e+00 2.067973e-01
##  [536] 2.506646e+00 3.661123e+00 1.439146e+01 4.634627e+00 2.404478e+00
##  [541] 1.874357e+01 5.510278e+00 2.016000e+00 5.406862e-01 6.086289e+00
##  [546] 3.729807e+00 2.065645e-01 2.459390e+00 2.512896e+00 6.408991e+00
##  [551] 1.825747e+00 4.285021e-01 8.980243e-01 1.715327e+00 4.891244e+00
##  [556] 6.046531e+00 4.686131e+00 4.381780e+00 4.312842e+00 1.132962e+01
##  [561] 4.237917e+00 3.650978e+00 9.075387e-01 8.625084e-01 9.471959e-01
##  [566] 6.481998e+00 6.049830e+00 9.481069e-01 5.010196e+00 6.390287e+00
##  [571] 8.897433e-01 2.506873e+00 4.529965e-01 1.167127e+01 1.625017e+00
##  [576] 3.343747e+00 2.552589e-01 7.673762e+00 1.975912e+00 6.758495e+00
##  [581] 3.632075e+00 4.150210e+00 1.667383e+01 1.468451e+00 1.904755e+00
##  [586] 2.467935e+00 7.717059e+00 1.704090e+01 2.482575e+01 1.006764e+01
##  [591] 1.992495e+00 1.537763e+00 1.329949e+01 7.691299e-01 3.466986e+00
##  [596] 1.543693e+00 1.144813e+00 1.069200e+00 8.044451e-01 6.613577e+00
##  [601] 6.156885e-01 4.288635e-01 1.605154e+00 3.880187e-01 3.845394e+00
##  [606] 5.276830e+00 4.492886e+00 5.541346e-04 3.964760e+00 9.987026e+00
##  [611] 7.300243e-02 2.555657e+00 5.860130e+00 1.846179e+00 2.767023e+00
##  [616] 3.643576e+00 4.079984e+00 1.391849e+01 2.779034e+00 7.865375e-01
##  [621] 6.685415e+00 7.747050e+00 3.838030e-01 3.272640e+00 5.693281e+00
##  [626] 4.217106e-01 7.337627e+00 2.567626e-01 8.791320e+00 3.173857e+00
##  [631] 1.112303e+00 2.070696e+00 7.593788e+00 5.091645e-01 4.539492e+00
##  [636] 3.009457e+01 5.100598e-01 4.063627e-01 7.399487e+00 7.205758e+00
##  [641] 1.152428e+00 1.053096e+01 2.086565e+01 5.343410e+00 6.832468e+00
##  [646] 1.576684e+00 1.459909e+01 7.894034e+00 2.316603e+00 9.470620e-02
##  [651] 1.107674e+01 1.271989e+00 6.188508e+00 4.425584e-03 1.650067e+01
##  [656] 6.832727e-01 1.122418e+01 6.699647e+00 1.440894e+01 5.630402e+00
##  [661] 8.443094e+00 5.771815e+00 6.334108e-02 6.828728e-01 1.731687e+01
##  [666] 1.482875e+00 3.593408e+00 4.883334e+00 2.531061e+00 6.677589e+00
##  [671] 1.849414e+00 1.437180e+01 2.802250e+00 5.733098e+00 6.990264e-01
##  [676] 7.199332e+00 4.187359e+00 2.197726e+00 7.733858e-01 8.121895e+00
##  [681] 5.006010e+00 1.810418e+00 4.206492e+00 1.778233e+00 3.456223e+00
##  [686] 8.279933e-01 1.918573e+00 1.063344e+01 2.368824e+00 1.720704e+00
##  [691] 2.680693e+00 1.940388e+00 2.380160e+00 4.336106e+00 1.188183e+01
##  [696] 1.260691e+01 8.757369e+00 1.906120e+01 2.519086e+00 1.055129e+01
##  [701] 1.609344e+01 1.951181e+00 7.619068e+00 2.317741e+00 1.343045e+01
##  [706] 3.046210e+00 4.331526e+00 1.868979e+01 2.656886e-01 3.552035e+00
##  [711] 4.087628e-01 9.994669e-01 3.893064e+00 4.774195e+00 1.878601e+00
##  [716] 1.140147e+01 5.873598e+00 4.462786e-01 1.926146e+00 2.095572e+00
##  [721] 2.760021e+00 2.787598e-01 6.487663e+00 2.789854e+00 9.518470e+00
##  [726] 3.382016e-01 7.458720e-01 3.114142e+00 2.094195e+01 1.459474e+01
##  [731] 1.032403e+01 5.827366e+00 3.115804e-01 6.512793e+00 2.613921e-02
##  [736] 8.156611e-01 8.111530e+00 3.961630e-01 1.675218e+00 7.428847e+00
##  [741] 1.293249e+00 1.358546e+00 5.381836e-01 1.021615e+00 3.475696e+00
##  [746] 4.544417e+00 1.488261e+00 1.773023e+00 5.165378e-01 2.317832e+00
##  [751] 1.400351e+01 1.278061e+01 4.260432e+00 1.187057e+01 1.030281e+00
##  [756] 9.578227e+00 1.252246e+00 3.731014e+00 4.496264e+00 3.561618e+00
##  [761] 1.501407e+00 3.214697e+00 9.083190e-01 4.753686e+00 2.574298e+00
##  [766] 1.042191e+01 8.778685e-01 1.996255e+00 1.222427e+00 3.050329e+00
##  [771] 8.047003e-01 2.596853e-01 2.415196e+00 4.864150e+00 1.162429e+01
##  [776] 3.693116e+00 6.975917e-01 1.892159e+00 6.708019e+00 8.428737e+00
##  [781] 1.101958e+00 3.101633e+00 1.361626e+00 2.802492e+00 2.742986e-02
##  [786] 6.295821e-01 3.278889e+00 6.126113e-01 4.584645e+00 7.731110e+00
##  [791] 6.853334e+00 2.553645e-01 4.833426e+00 9.766289e-01 2.815642e-01
##  [796] 2.933398e+00 6.791694e-01 4.682155e-01 6.782519e+00 1.118533e+00
##  [801] 4.074464e+00 1.327978e+00 6.610970e-01 8.048884e+00 3.734935e+00
##  [806] 2.151780e+00 1.748649e+00 2.454139e+01 2.119262e+00 2.666585e+00
##  [811] 2.603036e+00 8.800692e+00 1.479098e+00 3.693926e-01 1.778240e+00
##  [816] 1.202390e-02 4.253509e-02 8.292750e-01 4.833051e+00 7.965184e+00
##  [821] 1.960985e+00 4.778127e-01 9.586527e-01 9.456378e+00 2.298345e+00
##  [826] 5.720552e+00 4.138837e-01 7.221909e+00 1.982691e+00 3.967332e+00
##  [831] 6.894511e+00 2.309716e-01 6.093293e-01 4.992392e+00 7.510413e+00
##  [836] 4.195158e+00 1.904020e+00 3.154574e+00 3.992692e+00 3.931257e+00
##  [841] 1.563452e+01 7.474641e+00 7.947678e+00 7.363799e+00 1.607916e+00
##  [846] 4.568362e+00 9.132132e+00 1.441885e+00 2.360219e-01 7.706434e-01
##  [851] 5.214533e+00 7.108133e+00 1.203339e+00 9.286197e-01 2.710691e+00
##  [856] 3.535666e+00 4.197076e+00 8.938579e-01 2.353343e+00 7.820631e-01
##  [861] 6.430861e+00 1.343324e+00 1.767359e+00 9.148285e+00 4.082584e-01
##  [866] 1.122748e+01 3.002162e+00 1.368956e-01 5.296866e+00 4.193090e-02
##  [871] 8.331979e+00 6.214408e+00 1.913512e+00 3.722375e+00 1.409978e+01
##  [876] 3.965185e+00 1.335379e+00 3.527461e+00 4.641174e+00 1.729024e+00
##  [881] 1.499119e+00 3.110021e+00 1.931348e+01 8.080471e+00 5.305717e+00
##  [886] 1.601314e+01 1.588003e+00 1.449555e+00 9.712236e+00 1.135469e+00
##  [891] 7.305823e+00 8.367078e+00 1.717389e+01 1.148459e-01 4.129786e+00
##  [896] 1.023632e+00 3.097656e+00 9.576362e+00 5.914842e+00 5.947054e-01
##  [901] 2.028179e+00 5.491674e+00 1.297782e-01 3.829684e+00 9.509972e+00
##  [906] 8.679526e+00 1.170536e+00 1.490571e+00 4.070563e+00 7.973175e+00
##  [911] 3.618252e-01 1.353442e+01 2.128175e+00 3.857670e-01 2.270614e+00
##  [916] 1.266505e+01 2.218804e+00 3.970949e+00 3.812422e+00 7.332045e+00
##  [921] 3.490803e+00 1.480064e+00 1.224875e+00 3.507106e+00 3.191189e-01
##  [926] 9.669214e-01 7.043074e+00 1.971417e+01 2.575322e+00 6.029252e+00
##  [931] 4.247415e+00 8.564271e+00 3.157686e+00 4.398955e+00 7.425632e+00
##  [936] 8.712726e+00 1.090227e+00 6.213259e+00 1.171428e+00 1.211346e+01
##  [941] 4.561731e+00 3.013631e+00 4.635931e+00 6.041445e-03 2.997461e+00
##  [946] 3.940717e-01 9.956861e+00 5.994540e+00 1.287579e+01 1.555637e+00
##  [951] 4.720586e-01 4.487758e+00 9.646879e+00 5.332520e-01 3.645202e+00
##  [956] 4.604202e-01 1.062237e+00 1.114193e+01 1.288362e+00 2.411412e+01
##  [961] 7.371856e+00 1.433002e+00 4.313071e-02 3.374460e+00 3.570283e+00
##  [966] 2.671791e+00 1.001278e+00 7.121358e+00 1.523964e+00 4.067644e+00
##  [971] 5.803494e+00 1.564900e+00 3.203936e+00 7.205730e-01 3.792373e+00
##  [976] 6.421257e+00 1.222427e+00 5.750889e+00 3.209142e+00 4.278668e+00
##  [981] 2.562945e+00 5.288629e+00 4.056093e-01 7.774208e-01 9.862100e+00
##  [986] 5.277152e+00 3.455658e+00 2.731394e+00 4.907407e-01 3.187330e-01
##  [991] 6.591019e+00 3.217431e+00 1.182987e+01 9.125725e+00 1.665867e+00
##  [996] 3.924898e+00 8.075852e+00 4.312704e+00 1.163042e+01 8.380039e-01
hist(td)

td <- as.data.frame(td)
ggplot(data = td,aes(td)) + geom_density(color = 'blue')

Qua đồ thị ta thấy xác suất mà quán không bán được ly trà đá trong vòng 1000 ngày là khoảng 40%

1.2 Thống kê mô tả

1.2.1 Biến bánh căn

library(fBasics)
basicStats(bc)

Thống kê mô tả của biến bánh căn thể hiện trung bình mỗi sáng quán sẽ bán được 47 phần bánh căn trong vòng 1000 ngày, ngày bán được ít nhất sẽ là 45 phần và ngày bán được nhiều nhất là 49 phần.

1.2.2 Biến sữa chua phô mai

scpm <- rpois (n=1000,lambda=10)
scpm
##    [1] 12  8 13 15 12 14  6  8  5 12 13  3 10  9 14 12  9  7 13  4 12  9 11  7
##   [25]  7 13 14 16 12  7 11  9 14 15  7  4  9  9  9 11  6 11 12 15 11  6  6  6
##   [49] 14  6  9 11  5 12 10  6 16 12 11  4 16 16 12  5 10  6  7 14 14  9 14 12
##   [73]  9  9 10 11 11  8 12  7 15  6  7  9  8 10 11  7  9 12  8  9  8 11  7 12
##   [97] 11 11  9 12 14 10 11  6  4  8 13  6 11 10  8 13 10  4 12  8  6  7 10 16
##  [121]  7 15 11  9 12 12  8 12 13 10  8 10  8 16 12 13 12  8 14 12 11 16  8 13
##  [145] 13 14 12 11  6 10  9  6 10  9 11  8 13  8  7  7 10 16  8 10  8  4  6 14
##  [169]  6 20  9 11 12 11 12 15  9 11 11  9 11  9 14  9 10 11 10  5  5 16 12  8
##  [193]  7 13 15 12 11  7 10 14  3  7  9  6  7  8  7  8 18 16 13  9  9  9  9 10
##  [217] 10 10  9  9 11  8 14 12  8 12  7 13 11  7 16  9 12  9  8 11  5 10 11 10
##  [241] 13 15  9 13 18 10 11 11  5  7 11  5 10 12 11 10  8  5  6 12 13 11 11  9
##  [265] 17 14  9  9 12 11 10  8 10  8  8 13 12 14 10 10 13 10 13  6  9  8 12  8
##  [289] 11 12 13  9  8  5  9  5 12  6  8 10  8  7 11 17 12  9 13  7 14  9  4 13
##  [313] 11 12 17  7 12  8  6  7  8 10  8 10 11  6  9 10  5  6  3  8 16 14 13 12
##  [337] 14  9  8 11  9  9 14 12  9 10 12 11 13  8 15 18 17 14 16  9 11 10  8 13
##  [361] 13  8  9  8 12 10 15 11  8 10  8  7 10  9 15  4  9  9  9  9  8  9  4  8
##  [385]  7  8 12  6 12  9  6  8  9  6 13 14 10 11  7  8 13  7 19 10 11 14  6  5
##  [409]  9  7 11 12  5 11 14  8  7  5  6 13 19 12  8 12 12 10  8 15  8 13 17 11
##  [433] 11 11 14 10 12 12 13  7 15 12 14  8 16  5  8 14  9 13 10 10 10 16 14  5
##  [457]  7 15 10  6  8 13 11 14 10 16 20  5 12 12 10  9  7  6  7  4  8  7 11  8
##  [481] 11 12 13 10  9 18  9  6 10  7  6 13  8 18 13  9  7  9 15  8 14  6  6 14
##  [505] 12 10 13 12  9 12 12 10 10 11 13  7  8 11 17  4  9  9  7 13 15 10 10 10
##  [529] 10 13  6 10 14  4  5 12  4 10 14  6  9  8 12  3  9 10 10 15 11 13  7 12
##  [553] 10  6 14  7  7  6 16 14 11 13  7 15 13 12 11 14  7 12  6  8  6  5 16 11
##  [577] 12 12 10  7 12 17  7 10 10  8 10  9 11  8  9 10  9 13 13 16 11 11 11 10
##  [601] 14  6 14 13  9  6 11  6  7  8 12 10 14  8 14 13 13 10  7 10 12  8 16 11
##  [625] 14 10  8 10 17  7 13  7 10 10 15 10 14 10 11 15 10 13 10  6  7  7  7  3
##  [649] 12  8  9 12 11 11 11 10 10 12 15 18 11 11  7  7  8 16 10  2 12 10  7  9
##  [673] 10  8 10  7 10  7  9  4 12  8  7 14  7 11 10  5 12  7 11 13  5 10 13 12
##  [697]  6 11  7 12 10 13  6 14  6  7 10 12  9 14  9 12 11 14  9  7  7 15 11 13
##  [721]  7  4 14 13  8 17 11 13 11 13  7 11 11  6  7  9 15 14  9  6  7  4  5  6
##  [745] 10 16 14  8  7  7 16  7 12  8  9 18 15  6  8 13 14  5  9 10 13  3 11  8
##  [769]  6  9 12  9 13  9 16 16 13 10  9 11 14 12  9 13 16 10 11 10  9 13 12 10
##  [793] 10  7 10  7  4  8 15 10  9  8  9  9  7 11  8  5  8 15  8  7  5 14  9  7
##  [817]  4  4 11 10 15  8 14 12  7  9  9  9 13  3 12 10 12  8 11 14  8 13 10  7
##  [841] 16  9  8  7 12 10  6 11  6 12 11 11 11  6  4  9 11 15  5 13  8 13  5  5
##  [865]  8 16  4 11 10 12  9 11  9 11  8  9 13  7 11 16 10  9 11  8 10 11 12  8
##  [889] 11 13 13 12 13  7 11 14 17  9 15 11 11  8 12 12 14 12  7  5  9 15 10 12
##  [913] 11  8  9  6 15  7  5  8 10  7 11 10 15  7 11  9  9  8 12 10  7 10 10  6
##  [937] 11  9 12  8  8 13  6 10 16 12  9 12  7 10  7 17  9  5 11  9  9  5 15 11
##  [961]  8 11  9 12 11  9  7  7  5  8 13 10 18  7  6 11 10  5 10  5  5  7 10 11
##  [985] 11 12  8  5  8 13 11  7 12  9  6 13 12 13  7 15
basicStats(scpm)

Thống kê mô tả của biến sữa chua phô mai thể hiện trung bình mỗi sáng quán sẽ bán được khoảng 10 hủ trong vòng 1000 ngày, ngày bán được ít nhất sẽ là 2 hủ và ngày bán được nhiều nhất là 22 hủ.

1.2.3 Biến sữa đậu nành

basicStats(sdn)

Thống kê mô tả của biến sữa đậu nành thể hiện trung bình mỗi sáng quán sẽ bán được 10 ly sữa đậu nành trong vòng 1000 ngày, sẽ có những ngày không bán được ly nào và ngày bán được nhiều nhất là khoảng 27-28 ly.

1.2.4 Biến trà tắc xí muội

basicStats(txm)

Thống kê mô tả của biến trà tắc xí muội thể hiện trung bình mỗi sáng quán sẽ bán được 3-4 ly trà tắc xí muội trong vòng 1000 ngày, sẽ có những ngày không bán được ly nào và ngày bán được nhiều nhất là 10 ly.

1.2.5 Biến trà đá

basicStats(td)

Thống kê mô tả của biến trà đá thể hiện trung bình mỗi sáng quán sẽ bán được 4-5 ly trà đá trong vòng 1000 ngày, sẽ có những ngày không bán được ly nào và ngày bán được nhiều nhất là 35 ly.

2 Bài tập về nhà tuần 2

2.1 Chọn và giải thích vấn đề cần mô phỏng

2.1.1 Chọn đối tượng mô phỏng

Mô phỏng doanh thu cho lò chả lụa trong 500 ngày:

2.1.2 Giải thích

  1. Lò chả lụa Châu Mai là 1 lò chả lụa truyền thống lâu đời, được bắt nguồn từ Làng giò chả Ước Lễ và tồn tại được hơn 50 năm

  2. Tại lò thì có các sản phẩm chính: chả lụa, chả bò, giò thủ, chả quế, nem chua,…

  3. Địa chỉ lò: 23 đường số 40, Thị trấn Củ Chi, huyện Củ Chi, TPHCM.

2.1.3 Các sản phẩm được bán tại lò:

  • Chả lụa: 160.000 đồng/kg

  • Chả bò: 260.000 đồng/kg

  • Chả gân: 160.000 đồng/kg

  • Chả quế: 160.000 đồng/kg

  • Giò thủ: 160.000 đồng/kg

  • Chả chiên: 140.000 đồng/kg

  • Bò viên: 180.000 đồng/kg

  • Heo viên: 160.000 đồng/kg

  • Mộc sống: 160.000 đồng/kg

  • Nem chua: 100.000 đồng/kg

2.2 Giải thích các biến

  • clua: Số kg chả lụa bán trong 1 ngày

  • cbo: Số kg chả bò bán trong 1 ngày

  • cgan: Số kg chả gân bán trong 1 ngày

  • cque: Số kg chả quế bán trong 1 ngày

  • gthu: Số kg giò thủ bán trong 1 ngày

  • cchien: Số kg chả chiên bán trong 1 ngày

  • bv: Số kg bò viên bán trong 1 ngày

  • hv: Số kg heo viên bán trong 1 ngày

  • ms: Số kg mộc sống bán trong 1 ngày

  • nc: Số kg nem chua bán trong 1 ngày

Dữ liệu được thu thập từ tháng 2/2022 để thấy được số lượng hàng mà mỗi ngày lò bán được trong vòng 500 ngày

library(readxl)
## Warning: package 'readxl' was built under R version 4.2.3
library(xlsx)
setwd("C:/MPNN2")
data<-read.xlsx("tuần 2.xlsx",2)
data