BT về nhà MPNN
Triệu Nguyễn Mai Thi
2023-06-05
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 errorslibrary(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] 47.51821 45.01771 49.47962 46.47272 47.99421 46.32663 49.52013 49.22804
## [9] 48.72094 48.45171 45.49880 48.47044 45.20580 46.25310 45.24522 45.31056
## [17] 49.68955 49.91778 49.45376 46.01314 46.72926 48.32166 49.92197 48.00984
## [25] 46.47162 48.11096 45.85551 48.37151 46.91832 49.43677 46.52903 47.28280
## [33] 47.63424 49.43467 46.46222 45.08109 46.66270 47.70799 48.66051 47.27553
## [41] 45.31157 46.63154 49.51259 48.35083 48.97065 45.97932 48.22257 46.66538
## [49] 47.24553 45.05417 47.85519 46.32535 47.07479 45.64775 48.91642 48.19184
## [57] 45.72892 45.64195 47.72264 46.95296 48.49970 49.92551 49.52261 48.20863
## [65] 48.12789 47.56262 49.61612 49.35381 49.17334 47.39986 48.34531 45.41298
## [73] 49.83302 45.69214 48.09733 45.18413 46.42604 48.80928 46.17986 49.88610
## [81] 49.76628 47.56551 48.42577 45.10185 46.09065 48.19210 49.46792 49.28025
## [89] 47.82273 47.67024 45.66001 45.42847 47.03326 45.85840 47.39360 45.31728
## [97] 47.01766 45.43590 46.35148 49.07185 49.24087 49.93159 48.81993 49.20594
## [105] 47.74280 45.19944 46.47536 46.77963 46.13782 49.24859 46.57502 45.79122
## [113] 47.39717 46.47794 49.85979 48.43622 48.52603 48.10232 47.93869 45.81736
## [121] 47.40303 47.02230 47.33946 47.20122 46.78277 49.84987 48.68078 48.33505
## [129] 45.25209 48.09405 46.15477 45.95453 45.55034 46.85460 48.55608 45.06729
## [137] 45.38223 48.78771 47.20763 49.17138 47.63203 47.27000 45.00710 45.77331
## [145] 46.14097 47.39628 49.60807 48.73514 46.11747 47.17957 49.17424 46.75923
## [153] 45.60700 46.08858 48.47370 45.91025 49.11798 49.30793 46.95720 45.71185
## [161] 47.67024 47.79416 45.83215 49.80930 48.74357 47.60475 49.19022 49.51393
## [169] 49.53792 48.29501 48.95778 45.61016 45.63083 48.89657 48.12383 45.27075
## [177] 49.01926 46.82510 46.86395 48.41602 48.54539 47.62291 45.72289 46.39746
## [185] 49.70444 49.73314 47.63666 48.65674 46.84632 47.14505 48.64092 46.50329
## [193] 47.51898 47.43234 46.12341 48.23091 45.82957 48.16799 45.20413 49.77022
## [201] 48.81674 47.83116 49.35774 45.14553 47.61098 48.69426 48.51938 49.99149
## [209] 45.54502 48.78702 46.35349 47.08073 48.53574 49.66898 46.90592 47.37905
## [217] 45.85848 45.71933 48.96964 48.24350 45.37073 46.23166 48.48426 47.52692
## [225] 47.83148 45.59091 45.60787 45.26331 49.44656 46.44721 48.99023 48.53624
## [233] 45.05617 47.23616 48.47393 49.40741 49.06371 47.47098 48.91543 45.22863
## [241] 47.39120 45.88954 46.73843 47.85652 47.84909 46.33453 48.65247 47.11726
## [249] 48.04765 47.01704 47.85224 47.40861 47.80949 45.21119 48.88174 45.60320
## [257] 46.52055 48.97402 47.60997 47.14272 46.63688 45.20114 45.66244 48.01693
## [265] 48.39351 46.86406 45.73245 45.16786 47.70370 49.83725 47.27423 48.38762
## [273] 47.81031 45.48569 48.82909 48.70776 48.78025 49.90988 49.10366 46.59461
## [281] 46.53636 45.26215 45.27587 49.14465 45.54115 46.12369 49.20508 48.47070
## [289] 47.05094 49.70412 45.47136 46.35619 46.10397 48.18445 48.62201 45.00478
## [297] 46.88374 45.37214 49.83058 46.38151 48.53615 45.24771 45.97479 46.71351
## [305] 49.77299 47.81656 49.27959 49.94814 47.53435 49.63047 46.76597 46.30273
## [313] 48.18601 49.64051 48.86067 48.16691 46.35022 45.14792 47.81546 45.37236
## [321] 47.90014 46.58707 46.86713 48.69542 49.22553 49.75114 45.38515 45.35844
## [329] 49.91563 48.22580 45.00102 48.80813 46.24530 45.24315 47.09031 49.43119
## [337] 49.33438 46.53388 46.83396 47.22786 47.70923 46.16066 46.37217 47.73702
## [345] 48.65647 48.49712 48.27932 49.49062 45.99990 49.41139 45.61259 45.64296
## [353] 46.57329 47.91940 47.31027 47.46570 49.84373 46.20125 45.62666 48.85277
## [361] 46.30845 48.81072 47.76302 45.41716 45.91030 46.67404 46.94822 49.96426
## [369] 46.57335 45.45518 49.92772 45.87584 49.08285 45.57653 48.57121 47.67633
## [377] 46.97816 46.89930 45.14501 47.65455 45.91328 48.26545 48.20875 49.79800
## [385] 46.79919 46.18576 49.70806 45.78180 45.43430 47.08243 48.81514 49.54304
## [393] 46.13323 46.30775 47.60976 49.16976 48.60821 46.97050 45.13361 45.70767
## [401] 46.74245 48.67195 48.17615 47.48641 45.58508 48.81738 48.25603 49.42452
## [409] 47.17619 48.83190 47.30120 46.46730 47.54181 47.58375 47.02899 45.60900
## [417] 49.70604 46.31658 49.74834 49.38101 47.81783 48.86040 46.33937 48.33997
## [425] 46.34382 47.04374 46.29129 47.76425 48.33810 46.76789 45.27567 47.72701
## [433] 48.59796 46.52849 49.80105 46.25659 47.75775 45.28145 49.24690 47.02388
## [441] 45.88608 46.04840 45.25083 47.60437 46.77736 48.40562 45.85915 48.48626
## [449] 45.00176 46.26552 45.36965 45.85728 46.17536 49.53592 49.02168 49.29648
## [457] 48.31196 49.67008 49.26736 48.14630 47.34175 46.39115 46.69811 49.11504
## [465] 47.32202 46.18818 47.47210 48.61197 49.33763 48.70206 46.87646 48.07476
## [473] 48.60795 48.26305 48.87612 45.01320 48.62662 46.52830 47.16051 45.41753
## [481] 47.73764 48.55295 48.70301 47.01735 47.11684 49.50970 45.74103 47.67230
## [489] 48.61227 47.03036 45.51696 47.36292 45.40222 49.87570 48.33128 48.32537
## [497] 49.79333 45.77740 49.76778 46.80068 45.64351 45.33811 45.77773 48.41385
## [505] 49.87568 48.02012 45.32331 46.45995 47.40336 47.46159 46.79504 47.62237
## [513] 47.86551 48.52166 45.81665 46.35635 45.01661 47.95547 49.75864 49.58914
## [521] 47.96963 47.68756 46.53734 49.53651 46.26048 47.67723 49.63183 47.26946
## [529] 48.92900 46.83124 48.08898 46.06460 48.53576 46.36960 47.84471 46.45039
## [537] 48.34160 49.39160 49.94370 49.33803 45.89486 47.96602 48.05698 49.01111
## [545] 49.06063 46.46825 46.45628 48.54302 48.11532 45.33375 47.39895 49.08965
## [553] 48.78886 46.35391 49.97052 46.61715 46.94536 48.58351 45.44114 47.43014
## [561] 47.05435 46.89132 47.89188 47.19356 46.12311 45.68534 49.14516 47.92878
## [569] 46.15865 48.04173 48.14195 47.19003 47.75150 46.91113 48.09492 47.24538
## [577] 45.49794 45.74751 49.31463 48.97746 46.32230 47.16841 47.75611 47.75404
## [585] 46.69668 46.17105 47.82590 47.45381 49.65609 46.56815 46.56647 47.43643
## [593] 45.35579 48.50847 45.66819 48.00379 48.93974 48.14945 46.79049 48.30426
## [601] 49.18057 47.55331 46.61742 46.44141 48.45897 46.45114 45.32459 47.28523
## [609] 49.15088 47.20815 46.71200 47.03327 48.34480 46.67396 47.88608 49.99124
## [617] 47.48738 45.27256 46.23380 45.48782 45.91920 48.40844 49.76682 45.55386
## [625] 48.07522 46.25359 49.51984 47.14174 49.61799 48.26841 45.95836 45.02796
## [633] 48.91361 48.65174 45.03743 49.13018 46.34079 46.54020 48.63603 47.49980
## [641] 48.78857 47.58033 49.24765 45.42032 45.36887 45.42652 49.68702 48.15855
## [649] 45.27470 45.38647 49.43169 47.67869 46.28532 49.41149 46.87174 45.18256
## [657] 45.88307 49.70635 48.38008 47.29667 48.02780 49.86348 49.45606 47.50898
## [665] 46.59168 48.43591 49.91735 45.28373 48.38105 47.72518 46.86534 46.17776
## [673] 46.63666 48.05690 46.19816 47.37597 49.44296 47.19079 47.91909 49.37863
## [681] 45.02752 45.98389 46.35737 46.87306 48.38361 49.87984 46.66240 46.18816
## [689] 45.98049 46.31734 48.09157 45.68451 47.25918 49.04421 46.87656 46.92334
## [697] 46.42054 47.88486 45.40735 49.84115 47.25119 48.56869 45.35567 46.65853
## [705] 47.72715 47.52228 46.82511 45.33242 46.41452 45.02924 46.33543 45.91743
## [713] 45.13090 47.25335 46.10800 48.86807 45.65066 45.87699 46.91711 46.74108
## [721] 45.84625 47.52411 47.82449 46.59591 47.96468 48.93011 48.19412 46.42654
## [729] 49.94586 49.58342 49.61536 47.53167 46.72022 45.19558 49.57870 45.78125
## [737] 49.97137 45.68217 45.86214 49.72013 47.95278 45.17358 48.16322 47.38525
## [745] 49.52646 45.19898 48.26373 48.99358 47.62311 47.64165 49.31757 46.92453
## [753] 46.23818 48.64915 48.43159 45.29890 47.44847 47.08993 49.21461 47.69488
## [761] 45.05098 49.28605 49.12898 48.78184 45.14150 46.52041 46.60497 46.39634
## [769] 45.09947 49.12620 48.70579 47.30266 48.10763 49.85613 47.65514 49.72706
## [777] 48.83737 46.34797 48.13727 45.39281 47.26160 47.00950 45.20449 46.41380
## [785] 47.84644 48.06366 47.21100 45.72391 49.34612 46.43132 49.11461 47.76210
## [793] 48.37649 45.55040 45.31798 46.13797 47.01169 45.81664 49.11943 49.55642
## [801] 46.95157 45.60278 47.01893 49.92209 46.56754 49.18859 48.66121 47.09165
## [809] 45.47229 48.10706 47.20615 45.16890 45.84314 47.02410 48.60330 45.25622
## [817] 48.26020 48.52536 45.50831 49.98479 46.70931 45.93974 48.23366 46.91220
## [825] 45.42799 47.08849 49.86804 45.41458 48.99557 49.10186 47.15055 49.38199
## [833] 48.13334 47.65554 47.68511 47.63494 46.27026 47.65448 48.19025 47.30240
## [841] 46.04050 49.51833 47.82652 48.29126 45.84669 45.77388 45.94207 45.67225
## [849] 48.42752 49.56310 47.60475 46.67656 49.98682 47.45098 47.24481 48.96954
## [857] 46.58785 47.18123 45.03146 46.44463 45.29525 45.22559 46.98685 48.87854
## [865] 45.46755 45.68075 47.64691 47.60800 48.74810 45.91457 48.22713 49.59096
## [873] 49.63338 49.71747 47.42079 47.00456 48.04308 49.74680 49.62146 49.74601
## [881] 45.39864 45.89871 46.70915 45.47103 47.36024 45.47304 49.94608 47.24824
## [889] 49.13443 47.76794 47.85725 46.21995 46.06774 45.89605 47.37386 45.93816
## [897] 45.72915 48.06582 48.03417 45.15834 49.44406 46.11131 49.20228 47.01809
## [905] 46.49973 49.67893 49.38016 48.16221 49.62128 48.12029 49.09280 46.68339
## [913] 46.20974 47.92621 45.91231 47.09293 48.03697 49.53489 46.00756 45.00624
## [921] 46.73702 47.71537 45.29047 47.91778 47.59828 49.18811 48.96075 49.75301
## [929] 45.44879 46.16831 47.42329 45.34915 47.37845 46.29991 49.50466 48.92949
## [937] 47.30402 46.33704 46.44357 45.00614 48.40114 48.75228 47.04435 49.47290
## [945] 49.54419 49.85467 48.56695 49.75058 46.14213 48.17340 46.90115 47.78868
## [953] 47.78957 45.21379 48.07392 47.05120 45.08080 49.92682 49.60490 49.12992
## [961] 45.78676 49.52123 46.11832 49.24674 46.95840 48.06958 48.98694 46.56850
## [969] 47.09258 48.70579 47.77674 48.61473 45.95562 48.82679 47.39436 45.17962
## [977] 45.68042 48.75038 46.09736 48.73319 46.04334 49.13956 48.28775 48.41796
## [985] 46.52397 49.23152 47.35396 46.86514 46.04827 49.93109 45.02851 49.13221
## [993] 48.68593 47.50012 45.77229 46.27348 45.36878 45.35063 49.11573 45.91349round(runif(n=1000,min=45,max=50),0)## [1] 49 45 48 47 47 49 50 49 45 45 49 46 50 46 50 45 46 46 47 46 46 46 48 49
## [25] 49 49 47 49 48 46 45 48 50 50 48 46 48 47 46 46 45 50 45 46 49 46 49 49
## [49] 48 47 49 48 45 46 50 50 46 49 45 47 46 48 45 48 47 45 47 48 48 47 46 47
## [73] 49 47 47 49 47 46 49 48 45 48 47 49 47 47 49 50 47 49 47 47 48 47 48 46
## [97] 50 48 45 46 46 48 47 50 46 48 47 50 47 50 47 48 45 47 50 45 49 49 46 48
## [121] 48 45 47 47 45 50 46 47 48 45 45 50 47 48 49 49 48 47 48 49 48 47 48 45
## [145] 45 49 46 49 48 46 49 49 49 49 47 47 48 49 49 49 48 50 48 48 49 49 46 46
## [169] 48 46 48 48 49 48 46 49 47 48 47 47 48 48 45 50 46 45 48 49 47 47 45 49
## [193] 49 49 45 48 49 47 50 47 45 46 47 47 47 47 47 49 46 48 48 49 49 49 48 48
## [217] 45 45 48 48 49 46 46 49 45 46 47 50 46 48 50 48 47 46 46 49 47 48 45 50
## [241] 46 47 46 48 46 46 45 49 46 49 46 49 46 46 48 48 47 48 49 46 47 46 48 49
## [265] 48 46 47 49 47 45 47 46 47 49 46 47 49 46 45 48 45 48 46 46 47 49 46 46
## [289] 49 47 46 49 45 47 46 46 47 47 49 49 48 46 50 47 49 46 47 50 48 48 45 47
## [313] 49 48 49 48 50 48 46 49 49 48 50 50 46 46 48 48 49 50 47 46 46 48 49 47
## [337] 47 48 45 49 45 47 49 46 46 50 45 46 49 47 48 48 49 45 48 48 46 49 49 45
## [361] 48 47 48 49 46 46 46 47 46 50 47 49 48 48 48 49 50 45 47 50 47 48 49 45
## [385] 49 49 47 48 47 49 47 47 49 48 46 49 46 48 49 46 45 46 49 48 46 47 49 45
## [409] 50 50 46 47 49 46 48 49 46 49 48 46 45 47 49 47 45 48 49 46 50 48 47 50
## [433] 47 47 46 49 49 48 48 45 47 47 49 47 47 47 45 49 50 48 48 46 49 49 45 48
## [457] 49 46 47 47 49 46 46 49 49 48 48 49 49 49 46 48 45 45 48 46 47 49 48 49
## [481] 47 50 46 45 46 49 46 50 49 48 48 47 49 50 48 48 49 45 46 45 50 49 46 47
## [505] 47 45 49 45 48 45 47 49 49 45 50 50 45 46 50 50 46 48 49 49 48 46 49 46
## [529] 49 49 46 45 47 50 50 49 46 47 48 45 47 47 50 50 47 47 46 48 47 49 49 45
## [553] 49 48 49 46 50 47 48 50 45 49 46 48 48 48 50 50 46 45 45 48 48 47 47 48
## [577] 45 48 49 49 45 50 45 48 47 47 47 47 46 50 48 45 49 48 47 49 49 47 47 47
## [601] 47 48 45 46 47 50 49 49 49 46 49 47 46 45 49 45 49 49 49 50 47 47 49 46
## [625] 47 48 50 50 49 50 50 50 49 46 50 46 46 47 49 45 49 50 46 46 47 50 46 47
## [649] 50 47 46 49 49 46 50 49 47 50 46 48 48 49 46 47 48 45 47 47 48 47 46 49
## [673] 46 47 47 50 47 46 45 50 45 48 46 47 47 49 47 46 49 50 49 49 48 47 48 46
## [697] 45 48 48 49 50 45 48 46 46 45 50 49 50 45 48 46 50 50 48 46 46 49 48 47
## [721] 48 50 47 47 48 45 49 48 49 50 45 47 50 50 47 49 45 49 47 48 50 45 46 46
## [745] 48 50 47 46 48 45 48 48 49 47 47 49 46 46 49 48 47 47 50 46 49 45 47 49
## [769] 48 47 50 47 49 48 47 48 48 45 46 46 47 48 48 46 48 47 46 50 49 50 46 45
## [793] 46 49 48 49 48 48 46 47 46 46 49 49 46 49 49 48 48 49 49 47 47 49 45 46
## [817] 48 49 49 48 47 48 45 50 47 46 47 49 45 47 45 47 48 46 47 49 46 48 50 46
## [841] 48 49 48 48 49 47 49 50 47 46 47 46 47 49 47 48 46 46 48 47 45 49 49 45
## [865] 48 47 49 49 48 50 48 46 49 46 49 46 47 49 45 50 50 46 46 48 47 48 46 49
## [889] 46 45 49 49 48 48 49 48 47 48 49 49 50 49 46 45 45 48 46 48 48 50 50 49
## [913] 48 49 47 49 49 50 46 45 47 46 45 47 49 47 50 48 46 50 50 48 47 47 49 49
## [937] 47 49 48 49 46 46 47 50 49 45 46 50 49 47 48 50 50 50 46 49 46 47 47 46
## [961] 46 47 49 48 45 46 48 45 48 47 49 48 46 46 49 49 49 46 48 45 46 47 47 49
## [985] 47 48 47 50 46 46 50 46 48 46 45 48 48 49 47 46hist(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] 13 7 21 11 7 7 8 10 12 6 12 13 11 10 13 15 14 18 9 13 13 7 10 13
## [25] 10 7 11 8 12 9 12 16 14 10 9 5 18 8 11 8 6 7 6 9 5 11 9 9
## [49] 12 12 10 8 16 11 12 10 10 6 5 5 11 10 5 5 3 11 13 8 5 17 11 9
## [73] 10 8 9 7 6 5 14 6 4 12 10 8 9 8 12 13 12 12 10 10 8 9 8 14
## [97] 11 9 10 10 5 14 15 4 10 11 6 16 9 6 13 8 10 11 11 7 7 4 10 11
## [121] 11 8 5 11 12 13 9 18 12 12 16 11 16 10 13 20 14 14 8 13 12 11 13 15
## [145] 11 7 13 6 10 11 7 12 17 18 4 10 12 15 13 18 18 13 14 10 9 11 14 13
## [169] 13 11 13 10 7 11 8 8 13 6 12 9 13 15 10 10 11 12 15 15 7 12 6 9
## [193] 8 12 13 7 12 11 14 9 11 17 5 13 10 6 15 12 9 7 18 9 9 9 6 8
## [217] 12 6 14 12 16 6 12 9 13 9 16 8 12 5 9 8 11 7 5 12 8 6 10 10
## [241] 14 20 6 10 6 12 11 12 4 12 5 9 7 10 5 9 5 9 4 7 6 13 10 13
## [265] 10 13 11 5 8 8 10 11 13 9 14 11 10 14 10 12 12 10 11 10 12 10 6 8
## [289] 16 9 8 14 16 9 10 9 7 10 9 11 13 7 8 11 6 10 7 10 11 10 10 8
## [313] 8 13 13 13 8 7 13 8 10 8 10 6 7 9 7 8 6 10 7 12 9 11 8 11
## [337] 11 7 7 10 6 12 10 8 7 15 12 9 11 15 15 13 7 9 8 17 8 8 7 11
## [361] 6 8 14 11 15 14 7 8 7 8 9 9 8 4 9 4 11 14 6 12 19 7 14 5
## [385] 14 9 10 9 11 9 13 10 9 14 10 15 9 10 11 16 11 11 4 14 14 6 9 6
## [409] 10 9 5 7 9 9 11 6 11 14 10 12 7 4 11 10 8 3 11 7 13 4 13 11
## [433] 11 19 18 9 7 7 15 10 8 7 9 6 12 15 9 11 13 14 11 9 4 12 10 7
## [457] 14 9 5 7 10 10 5 11 13 5 8 10 7 6 11 13 12 13 7 5 10 13 8 14
## [481] 9 15 6 13 7 8 8 13 9 11 15 8 10 10 2 13 13 20 13 7 10 8 11 10
## [505] 12 13 7 7 9 8 7 9 7 7 11 12 9 7 11 10 9 7 13 5 12 14 11 5
## [529] 7 12 6 5 15 11 7 6 9 13 8 12 15 15 12 11 11 8 8 13 8 9 6 17
## [553] 7 12 14 9 9 8 7 10 8 12 9 7 9 8 9 9 11 15 22 13 13 13 13 11
## [577] 15 6 16 5 12 8 6 12 8 13 13 10 14 7 7 7 6 14 7 8 16 8 17 12
## [601] 9 4 7 11 5 13 4 10 11 9 10 5 11 17 10 7 13 9 7 11 16 13 9 12
## [625] 6 11 8 19 10 9 10 11 7 6 13 8 9 12 5 9 13 11 11 10 13 11 5 8
## [649] 11 11 12 6 7 12 8 6 10 6 16 10 11 19 9 13 6 11 10 8 7 11 8 10
## [673] 11 10 9 9 13 9 7 11 11 3 10 17 6 13 8 7 12 12 12 9 7 11 12 9
## [697] 12 10 6 12 8 9 13 7 10 5 7 6 21 12 9 7 14 12 10 6 11 14 10 5
## [721] 6 10 10 12 9 6 9 9 10 13 7 17 9 11 10 14 13 9 6 8 14 9 9 17
## [745] 8 9 14 12 8 11 19 8 6 10 12 11 5 14 12 16 7 8 5 8 11 14 9 15
## [769] 6 12 16 14 10 9 8 4 14 10 6 9 14 10 10 10 12 8 7 11 8 9 10 14
## [793] 14 14 9 10 14 11 16 10 9 12 9 11 8 8 11 7 8 19 7 8 10 11 9 10
## [817] 10 10 11 7 10 12 8 8 11 7 10 12 16 12 9 10 11 11 9 9 12 9 8 12
## [841] 18 6 17 9 13 17 7 13 10 14 10 11 5 5 12 7 17 9 10 9 10 5 9 12
## [865] 10 12 12 7 4 20 8 8 12 4 12 10 13 10 7 13 10 9 6 8 10 10 9 11
## [889] 8 17 6 11 11 11 8 13 10 10 8 9 11 10 14 7 12 16 6 12 8 14 8 9
## [913] 15 10 12 14 10 9 6 6 6 12 13 12 6 10 11 8 11 12 10 16 7 10 8 9
## [937] 12 15 6 9 5 5 8 13 6 10 8 10 9 10 10 9 14 12 9 8 11 13 10 12
## [961] 7 7 11 15 7 4 9 10 10 6 16 13 8 12 11 11 11 9 7 15 13 12 16 11
## [985] 10 12 10 6 16 6 5 9 20 12 8 19 23 16 14 13scpm <- 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] 0.93045206 15.99900119 9.43956539 9.54942834 11.02762565 5.86755871
## [7] 17.51980386 13.32287102 8.03561125 16.90794497 10.24107874 14.29462551
## [13] -0.79358658 13.29347148 12.39225018 14.22551620 10.14739593 5.08712389
## [19] 11.91020889 16.12827586 14.65192751 7.69468469 12.01939705 3.77521142
## [25] 3.80121180 13.07301242 16.79814282 13.20977354 10.99653097 19.47455013
## [31] 14.32650333 17.89948342 9.62409168 12.83381757 15.74108837 3.39137313
## [37] 9.83082404 8.13643670 11.83643999 11.21027027 17.41533436 9.14745322
## [43] 11.00379649 2.57600065 9.08616330 8.74810858 7.72383386 9.98172472
## [49] 14.43425765 0.66460954 20.37787534 18.74984537 14.00317822 19.96305884
## [55] 16.30402007 3.95111092 14.82690510 5.96798086 15.52293538 7.14329404
## [61] 14.06008860 18.81595255 7.85696256 3.60555447 6.57931404 8.43902610
## [67] 18.60847519 13.21297128 6.15533715 12.88776343 16.01713490 13.97819988
## [73] 11.41935594 14.31202950 14.35849982 4.60787011 13.69172603 18.77761118
## [79] 19.71152501 16.93134035 5.59036869 14.44276888 8.45192104 10.55375473
## [85] -0.90618515 10.48367682 7.00119190 9.66821921 12.61640219 16.50407930
## [91] 16.92775670 8.28859314 1.62909342 14.27134756 20.73716484 11.35454214
## [97] 7.61650679 7.75553071 12.89508156 16.47445659 9.35210981 12.65055637
## [103] 2.71799228 7.63899366 12.78149778 15.38417649 16.76726080 10.62631035
## [109] 22.54384351 9.85384478 11.27310944 6.83618296 10.39777899 1.02409311
## [115] 2.68130084 -0.02239161 26.40048338 8.06846310 5.32536114 8.57327020
## [121] 5.01026154 6.07884843 13.63429540 2.98855688 9.33339586 6.28921459
## [127] 11.62588164 16.17720065 8.00282004 8.65345198 13.30660167 15.70902078
## [133] 15.39831088 13.24236786 13.35755326 2.38652090 4.86054344 10.53252314
## [139] 9.80886239 13.09348895 13.48135849 9.47411359 11.52333410 15.61730073
## [145] 5.08539089 6.05478702 5.10191354 5.55661930 7.64090306 2.33359635
## [151] 11.02011847 4.10489147 13.31529691 12.32927364 6.38109554 12.74898522
## [157] 18.63954360 12.33690916 14.34358815 7.19240899 6.36454227 7.71515631
## [163] 5.20747875 8.28515984 10.75996176 2.39031571 14.93244765 8.64162553
## [169] 17.19226489 8.56186277 15.68472207 12.44189112 4.00112916 14.61489701
## [175] 4.80142737 6.00688126 10.43670709 10.10837936 18.64559138 6.30496669
## [181] 20.31323782 7.46885708 13.43934942 22.13009482 19.03973822 10.21250717
## [187] 20.58751495 14.21137035 -2.59062307 13.81548910 8.11548841 12.63085567
## [193] 11.14896844 5.79599413 20.68691872 2.84966585 3.47917660 8.37843519
## [199] 11.13811612 10.24877756 8.26360673 12.27202184 17.72771913 19.84153464
## [205] 4.20735217 7.27006542 9.61308462 7.58104794 2.83716457 5.57486653
## [211] 3.06165054 9.53124647 6.01976069 9.41484625 13.18077565 25.95392456
## [217] 4.05801767 9.87074481 9.74835540 9.99158243 6.62852066 15.89550552
## [223] 17.92176300 15.79119425 9.29048170 4.93484231 13.39666758 6.82208007
## [229] 3.51236341 5.47505871 12.05982543 5.18291816 20.81500186 12.79780456
## [235] 9.55254235 6.24870208 8.26308487 15.30253390 18.37430057 15.44044106
## [241] 3.29899558 13.33174611 7.73251849 3.23028519 10.78703365 13.69943067
## [247] 11.08827396 9.72065272 11.10058458 13.78389426 8.65402618 16.06251278
## [253] 9.15851678 8.59952997 15.33942206 13.95437794 13.38829259 21.29017643
## [259] 18.16591178 11.91179330 9.98203804 16.05664462 12.98977789 8.42593628
## [265] 15.66983157 10.40504550 14.23996571 17.45542817 11.12242578 9.58014118
## [271] 3.66323028 11.27238849 5.19515796 19.32156874 11.68230391 7.27288878
## [277] 2.77032173 13.14954724 10.99428319 12.50116045 17.31077240 17.52978651
## [283] 13.52525423 8.69793818 8.93441481 13.06889909 8.55431573 14.92699540
## [289] 14.26526853 19.73892135 14.19652564 5.05310856 2.27949257 9.63692994
## [295] 7.34511365 21.22346321 13.08358975 7.89882841 6.21910963 8.70487294
## [301] 8.54486987 7.95000788 1.40190863 10.95818682 22.21893477 7.25697209
## [307] 3.39741658 7.11412676 12.93057773 8.42794354 16.45826072 8.67968579
## [313] 11.95815126 18.98039548 3.72038070 9.39768556 0.65994384 12.43995841
## [319] 19.13711893 14.84661238 8.13410942 10.70239683 6.00952511 12.52611743
## [325] 10.00487669 13.74822384 15.69793373 -1.76716230 1.73679608 18.11381223
## [331] 17.01925719 12.58582804 11.60843733 3.04276032 5.91494649 16.64678425
## [337] 23.44123411 14.06693780 7.02748830 8.93435072 11.20157699 12.35841245
## [343] 14.60093044 11.23959439 3.94974738 8.91054899 13.22463829 1.11610438
## [349] 22.99485431 14.42881038 6.89137265 11.49923978 15.27909013 17.58124561
## [355] 7.95812065 15.37712194 8.60809527 6.90732338 10.11027032 9.84469907
## [361] 7.58345980 10.84335161 2.13531674 10.97018969 3.42187874 4.27942902
## [367] 7.16622285 7.81249916 16.81996107 9.19086483 10.52545623 7.29810872
## [373] 3.84062776 11.22937197 15.51461396 1.91837310 15.30114967 12.60771883
## [379] 1.03178946 12.75515531 5.02600608 15.40743246 4.67448871 6.31598985
## [385] 13.18689454 9.29132125 10.35921458 14.03439159 11.14937344 20.22675047
## [391] 7.26815691 10.33987496 11.81926428 4.14700005 15.52318100 5.44946058
## [397] 15.50578631 1.97385086 1.52198892 6.61500701 13.85281539 14.49103833
## [403] 7.00426522 7.14822205 6.04500527 9.37072202 18.14518636 13.84746227
## [409] 10.17749824 11.04401377 8.33535427 11.68244306 4.49011667 10.79089861
## [415] 14.92866035 11.85339259 9.23102641 18.18691123 7.79159444 14.14491972
## [421] 10.64673570 2.95459980 9.30326184 9.51323795 -3.06030602 15.83413952
## [427] 16.90203630 6.07792439 12.46659362 19.22873887 15.90722886 9.99845813
## [433] 5.31999387 10.29045956 4.21253950 13.42302338 14.18435404 7.44369087
## [439] 11.39906602 18.34283849 3.21194822 3.10680510 15.33371503 11.25919059
## [445] 13.42057060 11.53528309 11.53552063 13.54591576 18.73977450 6.28867564
## [451] 4.30814585 6.89897632 6.13986483 1.38093321 4.26507151 3.65654582
## [457] 8.96350586 12.88778950 14.21696935 15.10699804 12.70854934 9.41041249
## [463] 12.82627733 11.40437239 10.73141409 5.49706112 6.38315876 6.78372527
## [469] 8.45040553 12.30588292 11.64660559 12.75496119 10.05918969 1.29872165
## [475] 17.87880753 6.41778386 15.51705359 6.86979363 22.47576695 10.60837457
## [481] 14.92548492 12.36525172 4.20745346 14.20095461 14.23915444 12.79004526
## [487] 6.01173929 14.51666623 9.73235275 21.77772910 17.04359487 9.15231018
## [493] 8.75859091 7.00606179 15.02798174 13.04746332 9.81317375 12.97945222
## [499] 5.56593298 3.17659376 6.91222635 2.53519817 14.43078851 5.14034948
## [505] 10.49455973 0.51451337 12.12458326 20.17028326 15.57560050 13.57022888
## [511] 13.55337806 4.89820886 12.43127378 11.52610638 10.41619559 17.19311818
## [517] 6.62377779 7.43842199 2.67361129 5.44118410 19.25849211 13.00424338
## [523] 7.32209638 10.86830072 4.15577259 9.42151881 -1.27657547 9.65941686
## [529] 5.47742953 9.74175966 9.88046719 8.29644408 6.54037755 16.15830993
## [535] 17.01104870 12.88246702 10.77978620 20.40341951 10.56707325 15.33908245
## [541] 7.37417367 10.70059672 13.57895789 10.32048796 3.52505872 6.87316182
## [547] 9.43121966 9.72569257 3.44603846 8.51701045 14.65518584 18.12952482
## [553] 8.23062894 4.73657710 9.09965212 7.99196895 16.12015723 5.23635328
## [559] 5.24609677 4.97780678 9.37141704 6.51924318 0.87359267 9.68016271
## [565] 9.69131458 8.34060929 7.97502297 12.85511149 7.16928277 10.13083344
## [571] 13.31077914 6.93413701 7.90649305 8.09642017 7.71101913 9.73322993
## [577] 11.99442886 14.50047960 13.01363176 13.44293406 0.51779867 8.42059244
## [583] 1.02780625 9.26454617 10.56237295 2.49524827 -1.30618387 8.19788072
## [589] 15.19221456 9.73454407 3.95087500 11.89432631 11.00231498 0.28388713
## [595] -0.74975166 9.10329791 16.93750485 13.38877506 3.23276816 15.77528611
## [601] 8.04256985 10.52325692 10.44034980 7.13532813 11.67479820 11.94749622
## [607] 15.18293490 3.18961870 9.05611141 1.77649206 7.45835234 8.65271828
## [613] -0.38156280 7.99316150 9.45166930 11.39355918 7.14800152 17.81941257
## [619] 9.32372730 18.37968833 15.00507768 5.69268247 7.05388986 14.07448083
## [625] 10.97528111 11.32848974 3.20746495 9.10554482 9.29909407 22.79344820
## [631] 14.54958402 7.18804836 11.02753494 13.20245905 13.74162639 13.23096305
## [637] 9.73183576 6.83762680 7.29857081 21.42762302 11.42406469 11.24782762
## [643] 11.38477553 13.11071840 7.77674756 12.61304050 9.82523744 12.57454757
## [649] 6.39195218 11.00054530 19.79611768 16.82416914 5.11771602 7.31853217
## [655] 4.39529084 8.71090325 3.43581217 8.18171781 14.43730550 16.82080345
## [661] 10.38778480 12.52186159 5.68911335 9.31649576 11.27509832 12.23826899
## [667] 18.60314511 15.70065315 11.52463082 18.27285610 5.11338093 0.32176589
## [673] 5.79899798 13.07106577 12.02536710 4.56583378 6.88953208 9.95971054
## [679] 9.86159547 9.76569265 16.37235171 16.65862252 7.93405402 10.91359321
## [685] 8.45847953 10.66386331 14.96048008 11.31656882 12.14837441 2.42519895
## [691] 7.32231238 5.18572245 20.03788461 14.74311816 9.86142736 8.03677113
## [697] 6.63628997 6.64576550 6.78414135 5.01542842 7.83197993 1.65672297
## [703] 9.78494543 15.33407093 7.33159591 17.05146127 16.96247081 9.52718876
## [709] 11.48380988 12.85491192 11.27463747 13.58916102 11.98374755 10.81055826
## [715] 17.44193242 9.50174439 5.14254978 6.80577938 11.83690882 11.20846725
## [721] 12.15053139 10.17320959 3.91792845 10.94356742 18.10649667 5.79868053
## [727] 9.33972249 17.77676220 9.90434649 7.13050192 17.84775478 15.89324985
## [733] 6.80902686 10.23113761 7.94801511 8.17934590 12.04967646 7.52536591
## [739] 15.67149954 14.61094882 7.77536810 8.32322421 11.87336004 11.30861118
## [745] 12.50682909 5.72037604 5.83422406 14.46164308 6.97752945 17.64217521
## [751] 14.24079986 10.06493111 10.78664224 13.07668394 9.29339223 19.40851717
## [757] 6.09290364 1.90311771 9.82049531 7.64494849 11.85852879 9.53060865
## [763] 10.53648478 10.15571003 10.77600272 -1.59580628 12.20469809 13.65319672
## [769] 11.29076868 2.33090095 12.96358396 9.41576925 17.34221429 2.79571214
## [775] 15.76319785 2.30686781 6.61557208 3.23517342 13.09954639 16.71145690
## [781] 8.51641708 6.14258422 22.50707990 12.76643733 9.13048535 13.37928249
## [787] 14.03044579 15.34446544 -4.53368310 14.43538590 13.53979827 11.05924977
## [793] -3.37247200 5.01237334 0.18902447 11.25577680 8.20335753 13.75800884
## [799] 3.74144927 9.78878367 8.09751697 14.36924755 6.93766056 10.35227358
## [805] 9.83991607 18.58475131 18.54716734 12.62250157 5.49586439 9.37504045
## [811] 15.12943417 3.15348895 12.44752114 9.64380114 17.19414626 15.58748324
## [817] 14.70351309 -1.42576091 10.47113596 17.04595315 21.14170550 6.72935654
## [823] 7.61598893 14.97280570 12.32293696 9.71307863 14.15162426 8.64322358
## [829] 15.46672040 18.18108333 17.88895117 8.77973064 1.97992650 1.65995387
## [835] 14.71036704 12.04513296 10.91738452 10.85801781 8.97985109 20.48492545
## [841] 2.36371540 10.69229727 2.21086342 24.90238503 12.70764821 2.26695532
## [847] 9.22881176 9.31244385 5.94850784 9.53648322 5.77162196 0.93729018
## [853] 16.34918309 3.55384604 0.95420580 6.44059787 14.85483791 8.42468849
## [859] 10.72006869 13.46515776 11.51304207 10.71024609 14.45081833 10.35922959
## [865] 5.05850713 8.26506020 5.43972054 14.41512511 9.60427159 15.43292681
## [871] 6.56995677 2.07449335 8.45269813 10.00713677 5.00946624 1.57937650
## [877] 5.76738220 13.65505309 12.16394978 13.37193026 13.69173967 7.36510699
## [883] 10.95612604 7.69206936 9.34992890 8.16912157 12.42153524 2.33304601
## [889] 12.97284591 8.96660094 11.54014318 14.17213626 12.85887550 4.78343468
## [895] 8.56958926 13.31691854 15.20010883 11.22318000 7.87082652 8.23573469
## [901] 17.02951246 2.81817298 11.14879886 11.45119195 13.38828477 14.95868754
## [907] 9.48488419 7.74772997 11.50413628 3.17147534 7.97791629 11.63349346
## [913] 7.55402484 16.93662596 9.32868248 9.51855776 7.65422647 11.39042525
## [919] 2.40628076 9.11389941 10.21331146 8.28926038 13.76182022 21.34495627
## [925] 12.79115190 20.27271702 13.11867890 4.41216662 6.38932556 15.60795770
## [931] 3.52907544 6.00569378 4.18233867 13.38923519 4.33580140 9.81294671
## [937] 13.65404752 13.44891206 10.65507246 8.36933161 18.64999947 8.30674063
## [943] 13.64495199 14.22495434 2.74659443 3.76656285 11.55532371 12.97339731
## [949] 4.70253646 15.74221708 8.60559035 12.13862045 14.04468480 11.77919428
## [955] 6.85172579 7.77541131 9.21421465 7.66790926 18.87261420 12.85590847
## [961] 5.54139960 15.49961693 5.85585260 10.38084233 6.31633509 7.04563631
## [967] 11.44340331 0.44582920 13.61725072 5.93004629 19.78527955 2.44563117
## [973] 5.97764070 9.33734260 9.93908421 4.55702558 13.63670033 9.59308381
## [979] 14.75108345 18.72015090 14.59380851 6.10573903 13.53356953 9.74469965
## [985] 11.11259111 14.35258190 11.85885046 8.41694202 10.03668916 14.65605243
## [991] 17.89971385 5.14502469 16.80939027 -3.50013518 6.64813640 11.29176180
## [997] 16.84149315 8.93485591 12.80635501 14.54905136hist(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] 7 3 5 3 1 2 5 4 3 5 5 4 3 4 7 1 7 5 7 8 6 4 2 2
## [25] 4 3 8 1 1 5 4 6 4 3 2 5 1 3 3 2 3 3 4 3 3 7 4 4
## [49] 7 6 0 4 5 2 8 4 3 1 4 3 3 5 2 5 9 3 5 3 2 3 4 7
## [73] 5 4 3 4 5 5 4 5 5 5 3 3 5 5 5 6 6 4 1 6 2 2 4 2
## [97] 5 8 2 4 3 3 4 5 3 2 3 3 4 5 2 2 5 5 4 4 4 4 4 5
## [121] 3 4 2 3 3 7 6 6 0 6 3 3 2 5 4 5 3 2 5 1 2 6 5 3
## [145] 3 3 5 3 5 5 4 2 4 5 3 5 1 5 2 5 1 5 3 6 5 7 5 5
## [169] 2 3 3 3 5 2 6 6 4 6 3 1 2 5 3 4 7 3 5 5 4 2 5 9
## [193] 4 3 2 4 7 2 4 2 3 2 4 5 6 4 6 5 7 3 7 3 5 5 4 5
## [217] 3 5 3 0 2 3 6 4 7 4 2 5 7 5 4 6 1 4 4 4 4 6 6 7
## [241] 2 5 4 4 3 3 5 1 3 4 6 2 2 1 3 5 3 2 4 2 3 4 8 5
## [265] 6 5 5 5 3 6 2 2 6 6 4 2 4 6 4 6 2 6 4 0 2 2 4 1
## [289] 4 5 5 3 2 5 3 3 4 7 4 4 5 5 6 6 6 1 5 4 3 2 4 6
## [313] 4 5 6 7 4 5 3 4 3 5 4 2 4 2 3 2 5 3 4 5 4 6 5 4
## [337] 2 5 0 3 4 5 4 3 4 4 2 3 2 5 3 5 6 5 4 1 1 4 1 1
## [361] 5 3 8 8 3 2 4 3 3 5 6 3 3 3 3 5 5 3 4 3 3 3 1 2
## [385] 2 2 3 5 4 2 2 2 4 5 2 3 5 2 2 2 0 6 5 3 5 2 3 4
## [409] 1 4 2 7 3 1 2 5 3 1 6 0 4 2 7 3 6 2 4 6 2 6 7 2
## [433] 3 4 4 5 4 4 3 7 6 6 5 6 2 7 3 2 1 7 4 4 2 5 3 6
## [457] 6 3 5 5 4 0 3 4 3 5 3 2 5 6 4 7 3 4 8 5 3 5 2 5
## [481] 5 4 4 6 4 3 3 6 4 5 3 5 3 5 3 3 5 4 3 4 4 4 1 6
## [505] 7 2 6 2 2 2 3 4 7 2 5 3 3 3 3 3 3 5 3 6 4 6 2 6
## [529] 2 3 2 6 5 5 2 3 3 4 4 2 4 7 3 3 4 6 4 2 2 2 4 7
## [553] 5 6 2 2 4 2 10 7 6 7 5 8 5 3 5 2 6 3 4 6 4 6 4 2
## [577] 7 4 6 4 3 1 5 3 4 2 2 2 2 3 5 3 7 3 4 4 6 1 5 2
## [601] 1 5 8 4 7 5 2 5 4 9 6 3 5 4 5 3 5 3 4 4 2 2 5 4
## [625] 3 5 5 3 3 2 4 4 6 2 4 4 3 7 6 3 2 4 5 3 6 6 2 6
## [649] 7 5 7 6 3 3 5 7 2 4 3 4 4 2 6 3 5 4 5 4 3 5 5 5
## [673] 5 6 3 4 3 2 5 4 6 2 3 1 5 5 6 3 1 2 6 5 5 2 4 3
## [697] 1 6 3 5 5 5 3 2 2 2 3 5 6 6 4 5 4 4 7 4 8 4 3 3
## [721] 1 5 1 1 5 2 4 3 0 5 1 3 3 2 3 7 7 6 4 4 4 3 6 4
## [745] 6 4 2 4 7 2 2 3 3 4 2 3 6 1 2 2 5 4 6 3 6 3 3 4
## [769] 4 7 4 2 6 5 8 1 5 4 3 5 3 3 4 9 3 5 0 3 3 3 4 7
## [793] 5 5 6 5 5 4 4 7 5 4 3 5 4 3 4 5 4 3 5 3 3 3 5 7
## [817] 6 4 3 4 6 5 8 6 6 2 1 1 6 3 4 2 5 3 5 4 1 5 6 5
## [841] 2 5 4 6 6 1 3 3 2 6 4 3 5 1 3 3 3 7 3 2 4 4 2 5
## [865] 3 6 4 8 6 4 3 1 3 2 3 4 2 3 3 5 6 3 4 4 2 4 5 1
## [889] 4 3 2 2 3 7 5 4 4 4 3 7 3 2 7 7 3 5 4 3 5 3 6 2
## [913] 3 5 4 4 3 5 6 4 4 2 2 0 2 2 4 4 4 7 8 3 3 1 4 6
## [937] 4 4 5 2 3 3 4 5 5 4 5 2 3 3 3 3 4 7 5 5 3 3 3 1
## [961] 2 7 4 4 6 2 5 7 7 8 8 4 2 1 4 4 6 6 3 5 2 3 4 5
## [985] 3 4 6 4 3 3 5 3 4 4 5 2 4 5 4 6table(txm)## txm
## 0 1 2 3 4 5 6 7 8 9 10
## 11 48 145 217 209 186 106 56 17 4 1hist(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] 1.209045e+01 2.887279e+00 3.283834e+00 5.064195e+00 9.519790e+00
## [6] 7.556840e+00 5.606132e-01 5.828135e+00 4.638581e+00 8.312597e+00
## [11] 4.846968e+00 3.451097e+00 1.518589e+01 1.137002e+01 6.571154e+00
## [16] 5.973225e+00 6.627225e+00 3.794764e+00 8.935944e+00 8.407195e+00
## [21] 7.255037e+00 3.248394e+00 4.164278e+00 1.037753e+00 1.508567e+00
## [26] 8.327723e+00 7.991488e+00 8.575308e+00 5.601014e+00 1.989192e+00
## [31] 7.171077e+00 2.785281e+00 2.295122e+00 3.431325e+00 5.188795e+00
## [36] 5.671114e-01 7.367257e-01 7.537910e-01 5.182480e+00 4.146677e+00
## [41] 4.880429e+00 4.262660e+00 1.392665e+01 5.760774e-01 1.363812e+01
## [46] 5.351079e+00 2.860990e+00 2.074177e-01 4.248946e+00 1.350084e+00
## [51] 8.459998e+00 1.359861e+00 1.378485e+00 1.320697e+00 1.591866e+00
## [56] 5.560313e+00 3.113998e-01 6.265341e-01 5.452726e+00 9.517636e+00
## [61] 6.033258e+00 1.601145e+00 8.390247e+00 9.633896e+00 1.041855e-01
## [66] 1.914626e-01 1.210116e+01 1.088408e+00 2.071157e+00 2.253620e+00
## [71] 9.892759e+00 1.263757e+01 1.278253e-01 1.146958e+01 8.966271e+00
## [76] 4.732025e+00 3.707422e+00 2.458420e+00 5.862235e+00 1.516633e+01
## [81] 1.029456e+01 2.100056e+00 4.196143e-01 4.214570e+00 1.926916e+00
## [86] 1.789923e+00 1.611419e-01 7.509778e-01 1.724970e+00 6.097761e+00
## [91] 1.539500e+01 9.930640e-01 1.609712e+00 2.908681e+00 1.967101e-01
## [96] 2.844368e+00 6.488777e+00 2.098532e+00 1.748906e+00 7.344798e+00
## [101] 3.281183e+00 9.518915e-01 7.284030e+00 3.164376e+00 6.035046e-01
## [106] 1.088331e+01 8.786578e+00 4.132644e+00 1.058022e+00 1.173970e+00
## [111] 1.369016e+00 3.833666e+00 1.208704e+01 1.550773e+01 2.405865e+00
## [116] 7.053691e+00 2.857903e+00 2.496060e+00 3.340373e+00 2.170639e-02
## [121] 1.091230e+00 1.228259e+01 1.003085e+01 1.753787e-01 2.571723e+00
## [126] 2.323936e+00 9.213782e+00 2.641375e+00 7.275063e+00 2.648063e+00
## [131] 9.825265e+00 9.373280e-02 2.834702e+00 1.599131e+00 9.037235e+00
## [136] 3.821583e+00 1.367723e+00 1.861776e+00 2.356109e+00 9.262282e+00
## [141] 1.620220e-01 1.134144e+00 1.353619e+01 4.126539e+00 9.418032e+00
## [146] 3.005103e+00 1.314303e+00 4.904700e+00 2.677183e+00 5.984168e+00
## [151] 2.891304e+00 6.332981e-01 2.186636e+00 8.148966e-01 1.363618e+00
## [156] 3.739931e+00 1.455493e+00 3.494521e+00 1.597150e+00 6.432694e-02
## [161] 1.121620e+01 2.483700e-01 5.221978e+00 8.998788e+00 4.053893e+00
## [166] 4.976266e+00 6.184332e-01 4.429589e+00 1.567284e+01 6.422783e+00
## [171] 6.221430e+00 6.530040e-01 2.753402e+00 2.952350e-01 3.277168e+00
## [176] 9.783350e-01 1.313036e+01 2.378702e+00 2.245377e+00 1.703573e+01
## [181] 4.060166e+00 1.096497e-02 2.505973e+00 4.715971e+00 5.700389e-01
## [186] 3.538480e+00 4.205947e+00 3.382929e+00 1.362711e+00 2.155407e+00
## [191] 4.234293e+00 7.318391e+00 3.286496e+00 3.682335e+00 2.966776e+00
## [196] 2.716292e+00 2.038771e+00 1.081151e+00 4.034506e+00 9.816399e-01
## [201] 1.278048e+00 2.612780e-01 8.572427e-01 5.914979e+00 5.252503e+00
## [206] 1.429037e+00 2.234469e+00 7.416762e+00 4.999900e+00 6.764616e+00
## [211] 5.894733e+00 6.496216e+00 1.157908e+01 3.794609e+00 1.673906e+00
## [216] 3.268187e+00 6.934188e+00 2.850244e+00 2.900543e+00 4.420996e+00
## [221] 7.759288e+00 4.700633e+00 1.542205e+00 2.300521e+00 4.507844e-01
## [226] 6.074939e+00 5.652810e+00 4.788039e+00 4.797502e+00 1.222282e+00
## [231] 8.765195e+00 5.361070e-01 5.588351e-02 4.379871e-01 5.870572e-01
## [236] 3.481670e-01 1.274395e+00 6.324087e+00 1.621068e+00 5.800744e+00
## [241] 1.164336e+01 5.175479e-01 1.570147e+00 1.824330e+00 7.453971e+00
## [246] 7.560765e+00 2.083603e+00 1.628463e+00 4.303256e-01 4.610798e+00
## [251] 4.889705e+00 1.691229e+00 1.426727e+01 3.534973e+00 6.782414e-01
## [256] 1.333008e+00 2.014106e-01 4.532365e+00 4.153392e+00 2.935765e-01
## [261] 9.426842e-01 2.366637e-01 1.749417e+00 2.432037e+00 6.469711e+00
## [266] 1.254404e+00 1.011427e+00 2.185753e+01 1.032580e+01 1.135915e+00
## [271] 6.656219e-01 1.187212e+00 3.119196e+00 1.046045e+01 1.433422e+00
## [276] 1.500768e+01 3.776088e-01 8.528207e+00 5.587850e+00 6.622457e+00
## [281] 9.732417e-01 6.747835e+00 1.025580e+00 5.047045e+00 1.194411e+00
## [286] 4.366148e+00 2.576113e+00 4.241204e+00 1.533909e+00 9.208576e-01
## [291] 5.735115e+00 6.128013e+00 2.187919e+00 7.816686e+00 1.880283e+01
## [296] 1.021583e+01 1.460796e-01 1.182058e+01 4.551695e+00 1.707866e-01
## [301] 1.896669e-01 2.832312e-01 1.042950e-02 3.948673e+00 1.769205e+00
## [306] 1.504978e+01 1.820980e+00 4.464418e+00 4.327768e+00 7.884789e-01
## [311] 3.478601e+00 2.351806e+00 3.437082e+00 4.848083e+00 1.236180e+01
## [316] 2.442819e+00 1.584189e+01 5.160155e+00 2.781710e+00 1.135669e+01
## [321] 8.384133e+00 2.475714e+01 1.554584e+01 8.798445e+00 7.209674e-01
## [326] 6.434695e+00 5.606162e+00 9.491261e+00 8.798697e+00 1.024657e+00
## [331] 5.412627e+00 8.532370e+00 4.775619e+00 6.116358e+00 2.948884e-01
## [336] 9.971926e-01 6.277103e+00 1.185787e+00 1.415892e+00 9.612848e-01
## [341] 6.494541e+00 6.771735e+00 5.517095e+00 1.261461e+01 1.895117e+00
## [346] 5.340991e+00 1.609659e+00 1.051943e+01 2.115445e+00 5.550399e+00
## [351] 3.263614e+00 1.015817e+01 3.078131e+00 2.795636e-01 3.685316e+00
## [356] 1.863913e+00 9.569118e+00 4.977291e-01 7.466305e+00 2.765938e+00
## [361] 6.423743e-01 2.944944e+00 6.111647e+00 9.910057e-01 1.105617e+00
## [366] 4.898191e+00 5.943505e-03 4.674207e+00 1.350822e-01 2.852209e+00
## [371] 1.013314e+00 1.145077e+01 1.366149e-01 3.069325e+00 5.324291e+00
## [376] 9.494341e+00 2.415650e+01 3.970742e+00 8.924190e+00 3.914159e+00
## [381] 7.663674e+00 9.351688e+00 6.152270e-02 3.442187e+00 6.929144e+00
## [386] 1.406735e+01 6.859526e+00 8.748212e-02 2.023991e+00 7.701491e-01
## [391] 3.743507e-01 2.901161e+00 3.618782e+00 1.540886e+00 5.740539e-01
## [396] 7.984007e+00 3.488693e+00 7.647281e+00 1.256205e+00 3.678960e+00
## [401] 3.529348e+00 1.835417e+01 2.538372e+00 3.101785e-01 4.578801e+00
## [406] 2.088627e+00 1.292134e+01 3.156532e+00 3.710800e+00 1.020694e+01
## [411] 8.003796e+00 4.004482e+00 4.479029e+00 1.443903e+00 3.385878e+01
## [416] 3.069517e+00 5.913481e+00 1.452029e+01 1.706077e+00 3.408566e+00
## [421] 1.994550e+00 2.313338e+00 1.921853e+00 1.175143e+00 6.019883e-01
## [426] 2.492187e+00 1.439265e+01 2.517590e+00 2.455340e+00 8.787453e-02
## [431] 4.516243e-01 6.584950e-01 2.696263e+01 2.338326e+00 3.305538e+01
## [436] 2.736191e+00 8.932212e+00 3.505753e+00 1.471435e+01 1.323802e+01
## [441] 8.568992e-01 2.856146e+00 1.097606e+01 2.946308e+00 7.684329e+00
## [446] 8.069878e-02 5.397473e+00 2.105205e+00 2.219285e+00 7.555133e+00
## [451] 4.963363e-01 1.974999e+01 3.022282e+00 1.415259e+01 4.867711e-01
## [456] 4.480718e-01 3.496429e+00 5.277246e+00 5.370307e+00 1.323001e+00
## [461] 1.773239e+00 1.326775e+01 5.589724e-01 6.238010e+00 3.473323e+00
## [466] 5.338667e-01 5.474188e+00 4.748908e+00 4.963286e+00 6.138798e+00
## [471] 4.927841e+00 2.116520e+01 3.279469e+00 1.749481e-01 2.946974e+00
## [476] 2.107529e+00 1.981576e+00 1.044310e+01 2.440335e+00 2.834162e+00
## [481] 7.045806e+00 2.130720e+00 3.902039e+00 3.636904e+00 3.301625e-01
## [486] 2.052236e+00 4.836873e+00 1.482095e+00 7.435063e+00 1.112411e+00
## [491] 1.576222e+00 2.325702e+00 7.711581e+00 3.752731e+00 3.220594e-01
## [496] 7.919391e+00 7.386933e-01 5.896739e-01 3.948605e+00 1.249518e+01
## [501] 1.723311e+00 1.451007e+00 9.911567e+00 1.386885e+01 1.357036e+01
## [506] 5.342721e-01 7.847632e+00 4.900054e+00 1.509243e+01 1.751728e+01
## [511] 1.864906e+00 2.056690e+00 3.290492e-01 2.911365e+00 4.450095e+00
## [516] 5.409640e-01 7.050665e-01 1.017960e+00 4.121382e+00 1.134159e+01
## [521] 2.066830e+00 9.386156e-01 1.452182e-01 2.616098e+00 2.185355e+00
## [526] 2.556265e+00 4.745198e+00 1.330186e+01 1.111908e+01 3.800008e+00
## [531] 7.403039e+00 1.004695e+00 1.096522e+00 2.403394e+00 1.932191e+00
## [536] 1.595093e-01 5.497499e+00 1.148496e+01 1.112612e+01 1.012310e+01
## [541] 1.084506e+00 2.466102e+00 3.898167e+00 4.911849e-01 4.020017e+00
## [546] 3.977698e-01 5.874339e+00 1.022074e+00 8.020909e-01 1.107836e+01
## [551] 3.412656e-01 1.004247e+01 3.947545e+00 1.415953e+00 7.327566e+00
## [556] 1.182058e+01 1.243650e+00 6.496130e+00 3.644176e+00 2.713453e+00
## [561] 1.138988e+01 2.641752e+00 5.724108e+00 1.143680e+00 1.516136e+00
## [566] 1.768897e+00 1.284755e-01 3.847600e+00 8.806444e-01 3.341818e+00
## [571] 1.404687e+00 4.724656e+00 2.517444e-01 2.036228e+00 2.161273e+00
## [576] 1.917352e-01 5.458656e+00 2.252666e+00 9.288981e+00 4.604799e+00
## [581] 6.645571e+00 4.121258e+00 1.459974e+00 2.324146e+00 6.037806e-01
## [586] 1.739272e+00 7.201398e+00 1.101684e+00 1.146510e+01 2.374538e+00
## [591] 1.006663e+01 4.182202e-02 5.033236e+00 2.323220e+00 6.106271e+00
## [596] 1.909715e+00 5.040597e+00 5.175495e-01 4.324171e+00 1.309929e+01
## [601] 7.700039e+00 7.053690e+00 1.174562e+01 1.801986e-01 6.001908e+00
## [606] 4.233629e+00 1.265594e+01 1.413676e+01 1.102303e+01 2.646275e+00
## [611] 3.121819e+00 2.641759e-01 3.639445e+00 1.368255e+01 5.308375e+00
## [616] 2.972787e+00 1.080553e-01 1.964770e+00 3.808266e+00 1.422501e+00
## [621] 1.778131e+01 2.959768e+00 1.199831e+01 1.903884e+00 5.272688e+00
## [626] 6.560371e+00 1.760088e-01 1.616425e-01 4.568118e+00 1.499987e+00
## [631] 4.492085e+00 1.118513e-01 5.512310e-01 2.435768e+00 9.244578e+00
## [636] 1.442330e+01 5.075222e+00 5.509257e+00 9.629978e+00 1.048622e+01
## [641] 4.079512e+00 4.548412e-01 6.543324e+00 1.270292e-01 1.411281e+01
## [646] 1.874900e-01 1.374496e+01 2.322215e+00 3.083432e+00 8.470292e-01
## [651] 3.062293e-01 1.060845e+01 2.023439e+00 2.290938e+00 6.724719e-01
## [656] 9.339749e+00 1.072085e+00 8.386060e+00 7.030664e+00 3.880186e+00
## [661] 6.236615e+00 1.257549e+01 1.377336e+00 9.197957e+00 1.516884e+01
## [666] 5.087179e+00 3.477759e+00 9.987503e-04 3.036375e+00 4.649821e+00
## [671] 1.539637e+00 1.033551e+01 7.157342e-01 5.357436e+00 6.344234e-01
## [676] 6.029358e+00 1.304699e+00 2.040165e+00 1.549179e+01 6.800022e-01
## [681] 1.151333e+01 3.397661e-01 2.175210e+00 1.517916e+00 1.355409e+00
## [686] 1.729056e+01 8.009257e-02 3.309273e+00 4.109368e+00 3.533324e+00
## [691] 5.947942e+00 1.383206e+00 1.531898e+00 1.709765e+00 7.197188e-01
## [696] 1.641333e-01 2.725120e-01 1.563084e+00 6.964879e+00 3.666693e-01
## [701] 1.097781e+00 3.707195e+00 2.079459e+00 1.820039e+01 4.066777e+00
## [706] 5.368224e-01 3.574931e+00 4.267038e+00 1.467452e+01 4.454903e+00
## [711] 8.263556e+00 6.487596e-02 1.248377e+00 2.786471e+00 1.525400e-01
## [716] 2.689146e+00 8.722394e-01 1.198632e+01 2.072612e+00 1.073609e+01
## [721] 1.289008e-01 1.230564e+00 1.201390e+00 1.324078e+01 4.574541e+00
## [726] 1.960379e+00 2.934724e+00 4.122670e-01 3.312270e-01 5.906771e+00
## [731] 1.624834e+00 1.705073e+00 4.670704e+00 7.307900e+00 1.327413e+00
## [736] 1.163539e+01 3.650410e+00 2.309957e+00 8.075534e-01 7.728998e+00
## [741] 8.836885e+00 1.151950e+00 2.826832e+00 1.145744e+00 1.759403e-01
## [746] 5.083061e+00 8.812800e+00 9.046384e-01 1.720896e+00 1.233698e+01
## [751] 2.726631e+00 3.511076e-01 7.352960e-01 7.270840e+00 1.066472e+00
## [756] 2.138902e+00 3.596160e+00 1.874934e+01 1.633628e+01 1.983560e+00
## [761] 1.034642e+00 1.571866e+01 6.212164e-01 6.525511e-01 2.684571e+00
## [766] 5.032142e+00 7.193154e+00 1.221994e+01 2.606151e+00 8.698162e+00
## [771] 4.138655e+00 2.177433e-01 1.699096e-01 2.637942e+00 3.996575e+00
## [776] 2.576363e-01 1.293615e+01 3.961267e+00 6.442534e+00 3.136779e+00
## [781] 4.625878e+00 1.612214e+01 2.827507e+00 7.546258e+00 1.324868e+00
## [786] 5.497874e+00 5.356038e+00 5.018085e-01 1.161078e+00 3.557874e+00
## [791] 2.291590e+00 1.890916e+01 1.012506e+00 1.118586e+01 8.433856e-01
## [796] 9.345220e+00 9.103620e+00 9.785739e+00 3.012665e+00 6.004400e+00
## [801] 2.290751e+00 4.427753e+00 1.924676e-01 1.322392e+01 9.419687e+00
## [806] 3.692837e-01 7.333999e-01 3.002844e+00 3.242780e+00 8.994349e+00
## [811] 2.800482e+00 1.324244e+00 1.922233e+00 8.473239e+00 1.339189e+01
## [816] 1.586272e+01 2.405246e+00 3.951317e+00 4.287100e+00 5.456432e+00
## [821] 6.021533e+00 2.272192e+00 7.939023e+00 2.080752e+00 5.739882e+00
## [826] 5.652120e-02 6.282246e+00 6.983165e-01 1.133910e+01 4.577802e+00
## [831] 6.735720e+00 6.233689e-01 9.496350e+00 4.155144e+00 8.478687e+00
## [836] 2.674597e+00 3.636841e+00 1.150123e+01 2.999950e+00 2.222633e+00
## [841] 1.758132e-01 5.364332e+00 4.124898e+00 6.040473e-01 2.204933e+01
## [846] 1.856004e+00 4.785101e+00 7.944605e+00 8.721568e+00 3.323160e-01
## [851] 6.091907e+00 8.188491e+00 3.849123e+00 1.517250e+01 3.656623e+00
## [856] 1.797588e+00 2.099397e+00 8.322303e-02 1.356651e+01 4.244259e-01
## [861] 5.139465e+00 2.253072e+00 9.602753e+00 3.925551e-01 1.478621e+01
## [866] 4.374226e+00 6.589435e+00 1.893243e-01 6.396374e+00 4.965167e+00
## [871] 3.829190e+00 2.169338e+00 3.778812e+00 8.350495e+00 7.759941e-01
## [876] 8.699205e+00 8.411386e+00 5.267861e-01 1.427211e-01 6.324882e+00
## [881] 8.340405e+00 6.337912e-01 2.500185e+00 6.790296e+00 6.338030e+00
## [886] 4.099193e+00 1.613521e+00 3.407031e+00 8.075337e+00 4.742691e+00
## [891] 7.224629e+00 3.886203e+00 5.296701e+00 1.015344e+01 1.059909e+01
## [896] 5.444125e+00 1.597745e+01 1.374979e-01 1.294460e+00 9.582643e+00
## [901] 2.341173e+00 2.492027e+00 1.334643e+00 2.481718e+00 2.106338e+00
## [906] 1.950891e+01 2.173427e-02 1.044216e+01 2.637878e+00 3.584693e+00
## [911] 1.451673e+00 2.540619e+00 7.994664e+00 8.487300e+00 1.940283e+01
## [916] 1.381949e+00 8.142628e+00 1.577053e+01 8.638294e-01 1.155211e-02
## [921] 2.641838e+00 3.774135e+00 1.602646e+00 5.510088e+00 3.982847e-01
## [926] 6.232327e+00 3.176376e+00 1.887990e+00 5.094888e+00 1.323778e+00
## [931] 1.931854e+00 3.790298e+00 1.820360e+00 8.819860e-01 2.668319e+00
## [936] 4.983555e-01 2.108572e+00 5.887675e+00 2.272290e+00 2.752882e-01
## [941] 7.114408e-01 2.054769e+00 4.897992e+00 1.013986e+01 3.317128e+00
## [946] 3.662682e+00 1.031073e+01 1.089865e+01 3.530912e+00 7.158289e-01
## [951] 3.526520e+00 1.203318e+01 1.724604e+00 7.420409e+00 4.952338e+00
## [956] 1.788944e+01 9.862901e-02 6.321753e+00 7.783178e-01 2.825642e-03
## [961] 2.485140e+00 2.770784e+00 2.735054e+00 1.317688e+01 7.164443e-01
## [966] 3.565049e+00 2.948863e+00 1.789172e+00 2.290598e+00 1.391760e+00
## [971] 1.001492e+01 6.493654e-01 2.557373e+00 5.756254e+00 5.616737e-01
## [976] 2.626073e+00 1.242184e+01 8.217208e-01 1.188393e+00 2.369981e+00
## [981] 2.348703e+00 1.632064e+01 1.662494e+00 2.483408e-01 3.067879e+00
## [986] 5.072165e-01 2.010588e-01 1.735217e+01 9.607583e-01 5.618917e+00
## [991] 4.804153e+00 4.424414e+00 5.508945e+00 1.310977e+00 4.492195e-01
## [996] 2.182515e+00 5.930871e+00 1.484182e+01 1.511136e+01 1.986054e-01hist(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] 13 10 11 7 7 10 8 12 13 7 9 5 5 13 12 12 14 9 9 13 9 17 4 9
## [25] 7 10 10 11 8 11 8 6 4 7 8 13 7 9 7 11 8 11 10 10 4 8 11 8
## [49] 13 11 9 5 6 11 9 10 14 6 7 8 3 10 11 9 7 10 7 16 14 8 13 12
## [73] 13 7 8 16 15 11 12 12 13 6 7 11 7 12 11 11 11 13 8 12 16 15 5 7
## [97] 15 5 12 11 12 11 13 10 8 8 10 12 17 11 10 16 8 8 13 13 5 10 15 16
## [121] 9 9 8 12 7 13 7 16 11 11 9 9 9 16 10 10 14 9 6 10 8 12 10 5
## [145] 14 10 12 16 11 11 9 11 8 11 8 8 6 5 8 6 11 8 6 8 7 12 10 11
## [169] 6 15 8 9 13 7 8 13 11 14 8 5 11 9 13 17 9 7 10 5 8 14 4 9
## [193] 6 13 8 10 6 16 6 8 11 9 15 6 9 11 9 14 12 15 8 13 14 14 9 10
## [217] 9 7 9 15 20 8 8 6 8 13 16 8 9 12 8 12 15 11 12 9 10 9 3 2
## [241] 8 10 13 15 9 11 17 8 11 6 11 11 10 3 16 10 10 14 11 10 8 7 11 10
## [265] 12 9 12 10 11 13 10 9 12 10 5 4 10 11 7 11 9 9 8 7 14 4 10 11
## [289] 7 12 12 7 12 9 8 6 13 11 9 7 11 19 8 15 5 7 8 7 7 7 10 12
## [313] 17 11 8 8 8 7 13 6 16 6 12 13 10 13 10 9 10 13 8 12 19 7 8 8
## [337] 11 9 10 10 12 14 5 10 15 6 13 5 14 18 10 5 7 9 12 12 5 11 4 10
## [361] 10 11 6 9 10 6 11 8 8 6 9 9 11 7 11 8 5 11 9 12 8 8 12 10
## [385] 10 8 15 11 14 7 9 11 7 8 12 7 9 5 13 10 9 9 6 14 10 8 11 6
## [409] 11 12 4 13 9 8 11 14 7 11 6 9 16 19 9 9 9 14 11 8 10 6 12 9
## [433] 8 9 9 9 12 13 12 13 10 8 15 9 8 13 11 8 13 12 10 11 11 11 13 13
## [457] 10 13 7 11 10 12 4 7 15 5 9 13 11 5 11 9 15 7 12 13 11 12 6 6
## [481] 8 12 7 8 4 11 10 13 15 8 9 14 11 9 11 12 13 7 10 8 19 10 6 6
## [505] 7 9 14 16 11 14 10 13 9 5 8 7 16 8 8 12 9 9 5 8 4 12 15 10
## [529] 8 8 11 10 9 7 10 11 12 9 15 13 14 10 10 12 13 12 10 11 8 19 14 8
## [553] 13 6 10 13 11 9 7 9 13 17 16 7 13 10 16 9 13 9 13 6 7 6 10 9
## [577] 15 7 10 10 10 8 12 11 15 9 13 8 10 14 10 11 10 6 13 5 12 10 11 16
## [601] 5 11 11 5 9 6 6 8 13 8 8 15 4 12 9 12 10 9 9 9 11 6 15 13
## [625] 13 8 10 7 10 12 8 13 9 9 7 16 11 5 7 8 12 10 4 8 15 14 11 6
## [649] 9 3 10 8 7 15 11 8 6 9 12 10 12 6 9 7 10 13 14 12 7 11 10 8
## [673] 7 12 8 7 13 12 10 10 13 12 11 8 13 9 14 6 5 9 9 10 7 9 6 10
## [697] 10 9 8 6 14 4 10 7 9 13 21 7 11 7 5 8 12 8 12 10 13 13 7 9
## [721] 11 19 8 15 13 10 13 13 9 13 13 9 9 15 7 5 7 10 12 11 12 6 7 10
## [745] 10 12 6 8 9 8 4 5 14 6 9 10 14 11 8 6 7 7 7 9 10 13 10 13
## [769] 11 6 10 7 14 5 3 9 15 12 12 19 12 6 10 7 12 4 10 6 7 13 8 10
## [793] 14 7 8 8 11 6 11 6 15 12 13 8 11 11 10 11 14 10 11 9 14 3 10 14
## [817] 5 11 9 11 6 13 8 13 11 16 10 11 15 13 8 9 12 11 9 8 10 9 11 9
## [841] 7 10 5 8 11 15 6 7 13 14 15 15 12 11 12 11 13 10 6 11 7 12 11 10
## [865] 9 17 12 11 11 4 12 12 6 5 11 10 4 15 7 3 6 5 10 17 12 12 5 7
## [889] 7 7 8 9 8 15 13 12 8 9 13 12 12 10 11 8 13 7 12 8 9 18 9 14
## [913] 10 8 10 9 13 16 11 11 10 7 10 9 11 12 8 7 7 9 9 11 13 11 10 13
## [937] 9 7 9 16 7 8 13 10 10 14 14 7 10 9 12 4 15 9 11 15 8 8 9 8
## [961] 9 13 14 7 11 9 11 8 10 12 10 12 7 12 5 11 12 8 8 13 14 12 14 10
## [985] 5 11 11 8 9 9 11 10 7 20 8 14 8 9 14 6basicStats(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:
Doanh thu = Giá bán * Số lượng
2.1.2 Giải thích
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
Hình thức bán: chỉ bán cho khách hàng mang đi, bỏ sỉ cho các cửa hàng, tiểu thương bán tại chợ và các quán ăn, bệnh viện, nhà thờ,…Tại lò thì có các sản phẩm chính: chả lụa, chả bò, giò thủ, chả quế, nem chua,…
Địa chỉ: 23 đường số 40, Thị trấn Củ Chi, huyện Củ Chi, TPHCM.
2.2 Dữ liệu nghiên cứu
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
2.2.1 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.3 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
setwd("C:/MPNN2")
library(xlsx)
library(readxl)## Warning: package 'readxl' was built under R version 4.2.3data <-read.xlsx("tuần 2.xlsx",2)
data