MPNN
Nguyễn Phương Nhung
2023-07-09
1 .TUẦN 4: XÂY DỰNG MÔ HÌNH MÔ PHỎNG
1.1
2 .TUẦN 3: XÁC ĐỊNH PHÂN PHỐI CÁC BIẾN
2.1 Dữ liệu nghiên cứu
setwd("C:/PTDLDT")
library(readxl)## Warning: package 'readxl' was built under R version 4.2.3library(xlsx)
data <- read.xlsx("dulieucp.xlsx", sheetIndex = 1)
str(data)## 'data.frame': 121 obs. of 7 variables:
## $ Date: chr "01/03/2023" "01/04/2023" "01/05/2023" "01/06/2023" ...
## $ ACB : num 22700 22650 22850 23150 23100 ...
## $ FPT : num 80000 80000 80600 80300 80300 80400 80900 80700 80000 80800 ...
## $ GVR : num 14600 14850 14900 14400 14400 ...
## $ NVL : num 14650 14000 13700 13600 13850 ...
## $ SSI : num 18900 19000 19300 18900 19000 18900 19300 19300 19200 19400 ...
## $ VNM : num 78500 79300 81200 79800 80700 81000 79400 80000 80600 79400 ...summary(data)## Date ACB FPT GVR
## Length:121 Min. :21400 Min. :77500 Min. :13950
## Class :character 1st Qu.:24000 1st Qu.:79400 1st Qu.:14800
## Mode :character Median :24550 Median :80700 Median :15500
## Mean :24157 Mean :81295 Mean :15974
## 3rd Qu.:25000 3rd Qu.:83100 3rd Qu.:16700
## Max. :26350 Max. :87300 Max. :19600
## NVL SSI VNM
## Min. :10250 Min. :18250 Min. :65500
## 1st Qu.:12700 1st Qu.:19650 1st Qu.:69200
## Median :13600 Median :21200 Median :74400
## Mean :13307 Mean :21584 Mean :73100
## 3rd Qu.:14300 3rd Qu.:22850 3rd Qu.:76400
## Max. :15600 Max. :26600 Max. :81300sd(data$ACB)## [1] 1269.64Biến ACB có giá trị trung bình là 24157 tức giá đóng cửa trung bình của cổ phiếu ACB trong giai đoạn 01/03/2023 đến 30/06/2023 là 24.157 VNĐ. Bên cạnh đó giá trị cao nhất của cổ phiếu ACB là 26.350 VNĐ và thấp nhất là 21400 VNĐ, độ lệch chuẩn của mã ACB là xấp xỉ 1270 VNĐ.
sd(data$FPT)## [1] 2358.737Biến FPT có giá trị trung bình là 81295 tức giá đóng cửa trung bình của cổ phiếu FPT trong giai đoạn 01/03/2023 đến 30/06/2023 là 81.295 VNĐ. Bên cạnh đó giá trị cao nhất của cổ phiếu FPT là 87.300 VNĐ và thấp nhất là 77.500 VNĐ, độ lệch chuẩn của mã ACB là xấp xỉ 2359 VNĐ.
sd(data$GVR)## [1] 1561.824Biến GVR có giá trị trung bình là 15974 tức giá đóng cửa trung bình của cổ phiếu GVR trong giai đoạn 01/03/2023 đến 30/06/2023 là 15.974 VNĐ. Bên cạnh đó giá trị cao nhất của cổ phiếu GVR là 19.600 VNĐ và thấp nhất là 13.950 VNĐ, độ lệch chuẩn của mã GVR là xấp xỉ 1562 VNĐ.
sd(data$NVL)## [1] 1373.665Biến NVL có giá trị trung bình là 13307 tức giá đóng cửa trung bình của cổ phiếu NVL trong giai đoạn 01/03/2023 đến 30/06/2023 là 13.307 VNĐ. Bên cạnh đó giá trị cao nhất của cổ phiếu GVR là 15.600 VNĐ và thấp nhất là 10.250 VNĐ, độ lệch chuẩn của mã GVR là xấp xỉ 1374 VNĐ.
sd(data$SSI)## [1] 2275.744Biến SSI có giá trị trung bình là 21584 tức giá đóng cửa trung bình của cổ phiếu SSI trong giai đoạn 01/03/2023 đến 30/06/2023 là 21.584 VNĐ. Bên cạnh đó giá trị cao nhất của cổ phiếu SSI là 26.600 VNĐ và thấp nhất là 18.250 VNĐ, độ lệch chuẩn của mã SSI là xấp xỉ 2276 VNĐ.
sd(data$VNM)## [1] 4505.238Biến VNM có giá trị trung bình là 73100 tức giá đóng cửa trung bình của cổ phiếu VNM trong giai đoạn 01/03/2023 đến 30/06/2023 là 73.100 VNĐ. Bên cạnh đó giá trị cao nhất của cổ phiếu VNM là 81.300 VNĐ và thấp nhất là 65.500 VNĐ, độ lệch chuẩn của mã VNM là xấp xỉ 4505 VNĐ. Chứng tỏ rằng trong thời gian này biến động giá cổ phiếu VNM là khá đáng kể.
2.2 Xác định phân phối các biến đầu vào
2.2.1 Mã chứng khoán ACB
Đồ thị Histogram:
hist(data$ACB)
Đồ thị QQ plot:
library(stats)
qqnorm(data$ACB)
qqline(data$ACB)
Kiểm định Shapiro-Wilk Test:
Với giả thuyết H0: Dãy số liệu có phân phối Chuẩn
shapiro.test(data$ACB)##
## Shapiro-Wilk normality test
##
## data: data$ACB
## W = 0.83704, p-value = 3.103e-10Kết quả cho thấy giá trị tới hạn P_value = 3.103e-10 < 5% nên chưa có cơ sở để chấp nhận giả thuyết Ho tức giá cổ phiếu ACB không có tuân theo phân phối chuẩn.
2.2.2 Mã chứng khoán FPT
Đồ thị Histogram:
hist(data$FPT)
Đồ thị QQ plot:
library(stats)
qqnorm(data$FPT)
qqline(data$FPT)
Kiểm định Shapiro-Wilk Test:
Với giả thuyết H0: Dãy số liệu có phân phối Chuẩn
shapiro.test(data$FPT)##
## Shapiro-Wilk normality test
##
## data: data$FPT
## W = 0.94099, p-value = 4.646e-05Kết quả cho thấy giá trị tới hạn P_value = 4.646e-05 < 5% nên chưa có cơ sở để chấp nhận giả thuyết Ho tức giá cổ phiếu FPT không có tuân theo phân phối chuẩn.
2.2.3 Mã chứng khoán GVR
Đồ thị Histogram:
hist(data$GVR)
Đồ thị QQ plot:
library(stats)
qqnorm(data$GVR)
qqline(data$GVR)
Kiểm định Shapiro-Wilk Test:
Với giả thuyết H0: Dãy số liệu có phân phối Chuẩn
shapiro.test(data$GVR)##
## Shapiro-Wilk normality test
##
## data: data$GVR
## W = 0.88784, p-value = 4.452e-08Kết quả cho thấy giá trị tới hạn P_value = 4.452e-08 < 5% nên chưa có cơ sở để chấp nhận giả thuyết Ho tức giá cổ phiếu GVR không có tuân theo phân phối chuẩn.
2.2.4 Mã chứng khoán NVL
Đồ thị Histogram:
hist(data$NVL)
Đồ thị QQ plot:library(stats)
qqnorm(data$NVL)
qqline(data$NVL)
Kiểm định Shapiro-Wilk Test:
Với giả thuyết H0: Dãy số liệu có phân phối Chuẩn
shapiro.test(data$NVL)##
## Shapiro-Wilk normality test
##
## data: data$NVL
## W = 0.93264, p-value = 1.305e-05Kết quả cho thấy giá trị tới hạn P_value = 1.305e-05 < 5% nên chưa có cơ sở để chấp nhận giả thuyết Ho tức giá cổ phiếu NVL không có tuân theo phân phối chuẩn.
2.2.5 Mã chứng khoán SSI
Đồ thị Histogram:
hist(data$SSI)
Đồ thị QQ plot:
library(stats)
qqnorm(data$SSI)
qqline(data$SSI)
Kiểm định Shapiro-Wilk Test:
Với giả thuyết H0: Dãy số liệu có phân phối Chuẩn
shapiro.test(data$SSI)##
## Shapiro-Wilk normality test
##
## data: data$SSI
## W = 0.92973, p-value = 8.555e-06Kết quả cho thấy giá trị tới hạn P_value = 8.555e-06 < 5% nên chưa có cơ sở để chấp nhận giả thuyết Ho tức giá cổ phiếu SSI không có tuân theo phân phối chuẩn.
2.2.6 Mã chứng khoán VNM
Đồ thị Histogram:
hist(data$VNM)
Đồ thị QQ plot:
library(stats)
qqnorm(data$VNM)
qqline(data$VNM)
Kiểm định Shapiro-Wilk Test:
Với giả thuyết H0: Dãy số liệu có phân phối Chuẩn
shapiro.test(data$VNM)##
## Shapiro-Wilk normality test
##
## data: data$VNM
## W = 0.94751, p-value = 0.0001333Kết quả cho thấy giá trị tới hạn P_value = 0.000133 < 5% nên chưa có cơ sở để chấp nhận giả thuyết Ho tức giá cổ phiếu VNM không có tuân theo phân phối chuẩn.
2.3 Chuyển đổi dữ liệu tuân quy luật phân phối chuẩn
2.3.1 Chuyển đổi biến ACB
ACB <- sample(data$ACB, 3)
hist(ACB)
Kiêm định phân phối chuẩn:
shapiro.test(ACB)##
## Shapiro-Wilk normality test
##
## data: ACB
## W = 0.75, p-value < 2.2e-16Kết quảt trên cho thấy giá trị tới hạn P_value = 0.6369 > 5% tức có cơ sở để thừa nhận giả thuyết Ho nên giá cổ phiếu ACB tuân theo quy luật phân phối chuẩn
2.3.2 Chuyển đổi biến FPT
FPT <- sample(data$FPT, 5)
hist(FPT)
Kiểm định phân phối chuẩn:
shapiro.test(FPT)##
## Shapiro-Wilk normality test
##
## data: FPT
## W = 0.95646, p-value = 0.7832Kết quả trên cho thấy giá trị tới hạn P_value = 0.9605 > 5% tức có cơ sở để thừa nhận giả thuyết Ho nên giá cổ phiếu FPT tuân theo quy luật phân phối chuẩn.
2.3.3 Chuyển đổi biến GVR
GVR <- sample(data$GVR, 5)
hist(GVR)
Kiểm định phân phối chuẩn:
shapiro.test(GVR)##
## Shapiro-Wilk normality test
##
## data: GVR
## W = 0.80738, p-value = 0.09295Kết quả trên cho thấy giá trị tới hạn P_value = 0.159 > 5% tức có cơ sở để thừa nhận giả thuyết Ho nên giá cổ phiếu GVR tuân theo quy luật phân phối chuẩn.
2.3.4 Chuyển đổi biến NVL
NVL <- sample(data$NVL, 5)
hist(NVL)
Kiểm định phân phối chuẩn:
shapiro.test(NVL)##
## Shapiro-Wilk normality test
##
## data: NVL
## W = 0.94239, p-value = 0.6829Kết quả trên cho thấy giá trị tới hạn P_value = 0.5092 > 5% tức có cơ sở để thừa nhận giả thuyết Ho nên giá cổ phiếu NVL tuân theo quy luật phân phối chuẩn.
2.3.5 Chuyển đổi biến SSI
SSI <- sample(data$SSI, 5)
hist(SSI)
Kiểm định phân phối chuẩn:
shapiro.test(SSI)##
## Shapiro-Wilk normality test
##
## data: SSI
## W = 0.85853, p-value = 0.223Kết quả trên cho thấy giá trị tới hạn P_value = 0.3811 > 5% tức có cơ sở để thừa nhận giả thuyết Ho nên giá cổ phiếu SSI tuân theo quy luật phân phối chuẩn.
2.3.6 Chuyển đổi biến VNM
VNM <- sample(data$VNM, 5)
hist(VNM)
Kiểm định phân phối chuẩn:
shapiro.test(VNM)##
## Shapiro-Wilk normality test
##
## data: VNM
## W = 0.98345, p-value = 0.9522Kết quả trên cho thấy giá trị tới hạn P_value = 0.07119 > 5% tức có cơ sở để thừa nhận giả thuyết Ho nên giá cổ phiếu VNM tuân theo quy luật phân phối chuẩn.
3 .TUẦN 2: GIẢI THÍCH VẤN ĐỀ CẦN MÔ PHỎNG
3.1 Vấn đề cần mô phỏng
Thị trường chứng khoán Việt Nam hiện đang là kênh đầu tư khá hấp dẫn và nhiều tiềm năng. Trong những năm gần đây, thị trường hoạt động cực kì sôi động và đã có những bước tiến vượt bậc. Hàng loạt các công ty chứng khoán ra đời, các nhà đầu tư trong và ngoài nước rót vốn ồ ạt, các chỉ số thị trường liên tục phá kỷ lục. Việc lựa chọn kênh đầu tư để sinh lời bằng việc đầu tư vào cổ phiếu đã được rất nhiều người lựa chọn hiện nay. Việc xây dựng một danh mục đầu tư cổ phiếu an toàn và mang lại lợi nhuận vẫn là vấn đề mà giới đầu tư quan tâm đến hàng đầu.
Với nhiều nhà đầu tư nhóm cổ phiếu VN30 là thuật ngữ quen thuộc, mô tả tập hợp các cổ phiếu dẫn đầu ngành, lĩnh vực và thị trường. Nhóm cổ phiếu VN30 là tập hợp 30 mã cổ phiếu đạt các tiêu chuẩn giá trị vốn hóa lớn, tính thanh khoản cao, tỷ lệ free-float thấp hơn 5%.Rổ VN30 được coi là nhóm cổ phiếu blue chip, có tính thanh khoản cao, giá trị vốn hóa lớn, thường được nhà đầu tư lựa chọn đầu tư vì có tính an toàn cao, biến động cổ phiếu thấp. Các cổ phiếu bên trong rổ VN30 cho thấy xu hướng thị trường đang nghiêng về lĩnh vực nào … Tạo nên sức hấp dẫn lớn với nhiều nhà đầu tư bởi khả năng tăng trưởng ổn định. Các cổ phiếu bên trong rổ VN30 cho thấy xu hướng thị trường đang nghiêng về lĩnh vực nào. Từ đó giúp nhà đầu tư đánh giá được tổng quan thị trường, lựa chọn cổ phiếu đầu tư trong lĩnh vực phù hợp.
Với những đặc điểm về giá trị vốn hóa và tính thanh khoản cao, có nên đầu tư cổ phiếu VN30 không là mối quan tâm của nhiều người. Thị trường chứng khoán có hàng nghìn mã cổ phiếu, việc lựa chọn mã cổ phiếu phù hợp với mục đích sinh lời không hề đơn giản. Vì vậy trong bài mô phỏng này, tác giả sẽ thực hiện “Mô phỏng danh mục đầu tư và tỷ suất sinh lợi của 6 cổ phiếu thuộc nhóm VN30”, cụ thể hơn là mô phỏng ngẫu nhiên giá đóng cửa của 6 cổ phiếu trong rổ VN30 là ACB, FPT, GVR, NVL, SSI, VNM vào ngày 3/7/2023 từ đó có thể mô phỏng được giá trị và tỷ suất sinh lợi của danh mục đầu tư gồm 6 cổ phiếu trên trong ngày 03/07/2023 với số lượng đầu tư cổ phiếu là 1000 cp mỗi mã.
3.2 Mô hình nghiên cứu
Để mô phỏng được giá trị, tỷ suất sinh lời của danh mục đầu tư 6 cổ phiếu thuộc nhóm VN30, tác giả xin đề xuất mô hình như sau:
TGTDMĐT = 1000.ACB+1000.FPT+1000.GVR+1000.NVL+1000.SSI+1000.VNMTrong đó:
Biến phụ thuộc:
TGTDMĐT: Tổng giá trị danh mục đầu tư 6 cổ phiếu thuộc nhóm VN30 (đơn vị: VNĐ)
Biến độc lập:
ACB: Giá đóng cửa cổ phiếu ACB ngày 04/07/2023 (đơn vị: VNĐ)
FPT: Giá đóng cửa cổ phiếu FPT ngày 04/07/2023 (đơn vị: VNĐ)
GVR: Giá đóng cửa cổ phiếu GVR ngày 04/07/2023 (đơn vị: VNĐ)
NVL: Giá đóng cửa cổ phiếu NVL ngày 04/07/2023 (đơn vị: VNĐ)
SSI: Giá đóng cửa cổ phiếu SSI ngày 04/07/2023 (đơn vị: VNĐ)
VNM: Giá đóng cửa cổ phiếu VNM ngày 04/07/2023 (đơn vị: VNĐ)
3.3 Mô tả các biến trong mô hình nghiên cứu
ACB: là mã cổ phiếu của ngân hàng TMCP Á Châu, ACB là lã cổ phiếu thứ 10 của ngành ngân hàng có mặt trong rổ chỉ số VN30, mã cổ phiếu ACB chính thức được niêm yết trên Sở giao dịch chứng khoán Hồ Chí Minh (Ho Chi Minh Stock Exchange - Hose) từ ngày 09/12/2020 với giá tham chiếu 26.400 đồng/cp.
FPT: là mã cổ phiếu Công ty Cổ phần FPT.Ngày 13/12/2006 mã cổ phiếu FPT được niêm trên trên sàn Hose. Mã cổ phiếu FPT là một trong những cổ phiếu đầu ngành công nghệ thông tin trên thị trường chứng khoán Việt Nam.
GVR: là mã cổ phiếu của Tập đoàn công nghiệp cao su Việt Nam. Lĩnh vực kinh doanh chính là sản xuất plastic và cao su tổng hợp. Tháng 2/2018, GVR tiến hành phát hành cổ phiếu lần đầu ra công chúng (IPO), cổ phiếu GVR nằm trong rổ VN30, tuy có tuổi đời gia nhập sàn chứng khoán thấp hơn những cổ phiếu khác trong VN30 nhưng vẫn là cổ phiếu thu hút được nhiều nhà đầu tư chứng khoán.
NVL: là mã cổ phiếu của CTCP Tập đoàn Đầu tư Địa ốc No Va (NovaLand), Tập đoàn Nova hoạt động trong lĩnh vực sản xuất kinh doanh thuốc thú y, thuốc thủy sản, xây biệt thự cho thuê nhưng thế mạnh, mũi nhọn là kinh doanh bất động sản. Ngày 28/12/2016 mã cổ phiếu NVL của Tập đoàn NovaLand chính thức được niêm yết trên sàn Hose và là một trong những mã cổ phiếu thuộc ngành bất động sản lớn trong rổ VN30.
SSI: là mã cổ phiếu của CTCP Chứng khoán SSI. SSI hoạt động trên các lĩnh vực dịch vụ tài chính lớn bao gồm Dịch vụ Chứng khoán Khách hàng Cá nhân, Quản lý quỹ đầu tư, Dịch vụ Ngân hàng đầu tư,…và là một trong những công ty chứng khoán lớn nhất Việt Nam được niêm yết trên sàn Hose và lọt vào danh sách VN30.
VNM: là mã cổ phiếu nổi tiếng được Công ty cổ phần Sữa Việt Nam (Vinamilk) phát hành. Nhiều năm nay, Vinamilk luôn giữ vị trí số 1 trong ngành sữa Việt Nam đặc biệt trong nhóm ngành tiêu dùng. VNM là một trong những mã cổ phiếu lớn được niêm yết từ khá lâu trên sàn Hose (ngày 19/01/2006) và là mã cổ phiếu có giá trị lớn trong nhóm VN30.
Tác giả sẽ mô phỏng ngẫu nhiên giá đóng cửa của các mã cổ phiếu trên ở ngày 04/07/2023 thông qua các giá trị về tham số (trung bình, phương sai…) được tính toán từ giá lịch sử của từng mã cổ phiếu từ ngày 01/03/2023 đến ngày 30/06/2023. Dữ liệu lịch sử của từng mã Cổ phiếu được tổng hợp từ nguồn investing.com (https://www.investing.com/)
4 .TUẦN 1: Mô phỏng ngẫu nhiên các biến ngẫu nhiên
4.1 Mô phỏng biến Mì Quảng:
Số tô Mì Quảng bán được trong tháng có phân phối đều trong khoảng 1050-1500 tô/tháng
mq <- runif(1000, min = 1050, max = 1500)
mq <- round(runif(1000, min = 1050, max = 1500),0)
mq## [1] 1326 1258 1123 1348 1480 1299 1440 1099 1280 1459 1263 1290 1271 1429
## [15] 1325 1178 1194 1418 1154 1224 1343 1362 1179 1244 1093 1353 1219 1449
## [29] 1363 1425 1079 1059 1453 1311 1491 1443 1196 1098 1101 1202 1283 1204
## [43] 1214 1387 1395 1280 1361 1273 1384 1336 1468 1187 1309 1367 1306 1144
## [57] 1227 1298 1341 1381 1357 1152 1154 1240 1087 1189 1078 1305 1452 1409
## [71] 1431 1158 1094 1237 1105 1250 1337 1337 1155 1158 1276 1319 1261 1413
## [85] 1472 1247 1071 1138 1338 1117 1375 1296 1399 1229 1106 1377 1288 1342
## [99] 1269 1205 1382 1462 1206 1402 1473 1059 1169 1371 1103 1452 1158 1473
## [113] 1235 1283 1076 1109 1458 1354 1078 1338 1384 1278 1230 1259 1060 1304
## [127] 1436 1297 1301 1490 1375 1414 1492 1066 1395 1244 1287 1165 1472 1409
## [141] 1372 1237 1380 1149 1194 1207 1482 1488 1210 1075 1293 1460 1376 1479
## [155] 1107 1440 1241 1428 1415 1457 1145 1134 1151 1316 1206 1413 1405 1314
## [169] 1380 1380 1416 1100 1166 1215 1333 1415 1487 1224 1377 1146 1417 1189
## [183] 1179 1298 1399 1461 1216 1336 1254 1436 1450 1112 1371 1133 1233 1353
## [197] 1311 1315 1329 1153 1438 1474 1254 1383 1109 1471 1265 1375 1291 1204
## [211] 1153 1144 1402 1331 1476 1471 1129 1446 1239 1128 1374 1416 1368 1242
## [225] 1328 1460 1139 1457 1434 1052 1244 1071 1444 1264 1254 1437 1071 1263
## [239] 1193 1124 1271 1474 1216 1092 1109 1100 1094 1127 1354 1222 1132 1266
## [253] 1471 1194 1361 1438 1174 1172 1406 1139 1389 1156 1294 1178 1368 1179
## [267] 1393 1335 1430 1445 1254 1137 1066 1473 1223 1161 1247 1173 1058 1497
## [281] 1401 1426 1470 1103 1338 1419 1099 1133 1097 1060 1078 1110 1431 1282
## [295] 1274 1252 1344 1232 1437 1250 1139 1313 1171 1206 1469 1233 1464 1140
## [309] 1137 1347 1247 1310 1262 1180 1221 1161 1292 1258 1126 1115 1492 1291
## [323] 1288 1146 1181 1336 1349 1425 1126 1093 1445 1181 1127 1178 1206 1232
## [337] 1439 1212 1309 1322 1173 1088 1352 1495 1177 1498 1305 1326 1416 1473
## [351] 1071 1400 1325 1345 1366 1437 1278 1311 1071 1129 1443 1234 1345 1222
## [365] 1204 1352 1218 1401 1494 1362 1252 1221 1267 1397 1378 1465 1341 1445
## [379] 1140 1425 1352 1085 1419 1282 1227 1213 1176 1351 1347 1161 1461 1236
## [393] 1171 1154 1235 1387 1245 1198 1402 1159 1336 1061 1272 1369 1468 1424
## [407] 1257 1318 1200 1408 1491 1183 1275 1091 1326 1143 1460 1208 1225 1286
## [421] 1083 1100 1484 1216 1360 1282 1284 1268 1379 1100 1319 1265 1053 1408
## [435] 1197 1085 1112 1257 1492 1166 1174 1089 1400 1203 1141 1094 1168 1233
## [449] 1097 1488 1413 1069 1244 1340 1075 1285 1109 1297 1098 1379 1143 1160
## [463] 1261 1125 1134 1194 1071 1133 1114 1452 1368 1370 1257 1280 1333 1087
## [477] 1163 1052 1290 1418 1194 1054 1353 1361 1420 1159 1388 1174 1317 1167
## [491] 1177 1288 1153 1056 1179 1463 1227 1350 1072 1101 1259 1172 1370 1423
## [505] 1257 1089 1430 1165 1412 1228 1254 1299 1453 1188 1185 1170 1204 1154
## [519] 1174 1382 1069 1396 1497 1080 1331 1447 1500 1168 1352 1445 1410 1290
## [533] 1100 1399 1128 1351 1401 1460 1487 1301 1055 1378 1264 1251 1245 1204
## [547] 1259 1462 1373 1443 1354 1402 1107 1372 1391 1209 1253 1382 1168 1499
## [561] 1427 1430 1151 1139 1393 1301 1086 1345 1279 1313 1222 1311 1310 1180
## [575] 1093 1414 1146 1219 1362 1451 1311 1310 1404 1455 1419 1287 1074 1366
## [589] 1320 1273 1390 1185 1300 1397 1420 1377 1062 1225 1160 1454 1367 1427
## [603] 1137 1097 1067 1127 1382 1171 1138 1167 1311 1213 1468 1343 1180 1223
## [617] 1053 1063 1227 1428 1483 1417 1471 1351 1080 1055 1294 1299 1238 1435
## [631] 1073 1463 1441 1300 1392 1286 1157 1421 1230 1317 1370 1080 1427 1216
## [645] 1130 1434 1454 1119 1242 1497 1112 1282 1153 1130 1318 1361 1333 1068
## [659] 1443 1400 1119 1069 1297 1063 1257 1359 1387 1349 1421 1439 1272 1220
## [673] 1457 1478 1053 1401 1246 1287 1251 1122 1483 1214 1059 1154 1445 1385
## [687] 1457 1388 1490 1429 1130 1405 1258 1116 1206 1282 1207 1288 1323 1373
## [701] 1204 1232 1125 1072 1439 1208 1391 1118 1168 1202 1202 1282 1173 1159
## [715] 1219 1403 1174 1456 1107 1486 1473 1296 1163 1251 1438 1232 1464 1082
## [729] 1059 1111 1065 1396 1067 1482 1263 1293 1151 1303 1176 1085 1500 1358
## [743] 1418 1213 1112 1244 1489 1196 1135 1119 1208 1462 1431 1084 1368 1295
## [757] 1304 1224 1366 1310 1389 1400 1481 1294 1331 1089 1416 1470 1288 1478
## [771] 1103 1262 1050 1411 1404 1400 1200 1226 1462 1470 1139 1382 1217 1398
## [785] 1499 1233 1351 1051 1155 1132 1133 1264 1218 1342 1374 1405 1296 1403
## [799] 1396 1166 1237 1381 1124 1293 1472 1216 1142 1417 1072 1411 1085 1277
## [813] 1209 1183 1279 1084 1196 1331 1238 1309 1388 1396 1279 1086 1085 1294
## [827] 1498 1436 1262 1324 1439 1379 1097 1066 1219 1321 1430 1067 1056 1229
## [841] 1065 1085 1234 1463 1113 1363 1167 1159 1346 1341 1309 1378 1450 1407
## [855] 1421 1123 1252 1296 1476 1499 1333 1486 1403 1105 1421 1209 1110 1215
## [869] 1246 1100 1156 1458 1283 1227 1346 1438 1246 1124 1359 1137 1287 1463
## [883] 1071 1407 1484 1331 1310 1336 1275 1284 1488 1325 1151 1087 1091 1146
## [897] 1250 1341 1054 1094 1456 1315 1312 1338 1199 1168 1287 1328 1440 1114
## [911] 1212 1094 1230 1294 1389 1161 1127 1372 1100 1115 1393 1084 1313 1112
## [925] 1233 1314 1079 1161 1306 1153 1391 1260 1386 1425 1345 1471 1217 1204
## [939] 1217 1466 1182 1474 1452 1183 1141 1312 1298 1137 1203 1360 1414 1480
## [953] 1251 1367 1183 1432 1437 1165 1397 1209 1226 1149 1360 1387 1488 1354
## [967] 1272 1181 1176 1361 1172 1169 1227 1261 1200 1087 1272 1357 1050 1116
## [981] 1086 1496 1183 1437 1180 1296 1401 1186 1450 1075 1298 1196 1468 1074
## [995] 1181 1083 1369 1156 1369 1367hist(mq)
library(fBasics)
basicStats(mq)## mq
## nobs 1.000000e+03
## NAs 0.000000e+00
## Minimum 1.050000e+03
## Maximum 1.500000e+03
## 1. Quartile 1.167750e+03
## 3. Quartile 1.388000e+03
## Mean 1.276606e+03
## Median 1.279500e+03
## Sum 1.276606e+06
## SE Mean 4.047416e+00
## LCL Mean 1.268664e+03
## UCL Mean 1.284548e+03
## Variance 1.638157e+04
## Stdev 1.279905e+02
## Skewness -2.237200e-02
## Kurtosis -1.194310e+00Số phần Mì quảng mà quán bán được trung bình trong 1 tháng khoảng 1270 phần với độ lệch chuẩn là 133 phần. Trong đó, quán bán được nhiều nhất là khoảng 1498 phần và bán được ít nhất khoảng 1050 phần.
4.2 Môp phỏng biến Cao Lầu:
Số phần Cao Lầu bán được trong tháng có phân phối đều trong khoảng 600-900 tô/tháng:
cl <- runif(1000, min=600, max=900)
cl <- round(runif(1000, min = 600, max = 900),0)
cl## [1] 645 663 622 883 737 796 705 611 826 889 855 883 893 808 768 633 781 899
## [19] 836 751 821 609 772 820 659 752 725 730 722 870 814 855 730 706 639 666
## [37] 771 606 713 816 668 651 832 667 811 890 685 761 742 835 618 861 768 884
## [55] 687 723 754 860 741 636 658 817 609 887 836 865 760 645 752 678 714 768
## [73] 741 623 756 850 721 840 675 669 625 632 604 634 646 871 862 831 800 644
## [91] 753 654 842 892 801 634 797 796 864 617 755 812 716 707 641 724 614 819
## [109] 775 602 872 809 616 765 843 606 658 748 881 603 791 782 681 785 718 739
## [127] 665 641 612 789 611 792 660 863 784 616 651 863 647 729 778 768 809 804
## [145] 635 619 621 649 745 778 886 657 871 616 684 676 652 681 787 723 695 705
## [163] 648 863 782 687 791 668 758 898 790 815 869 825 637 617 603 892 762 647
## [181] 605 823 811 618 643 772 841 855 643 705 785 608 809 632 680 870 609 668
## [199] 750 787 634 711 836 642 727 841 739 662 644 802 707 867 832 772 708 741
## [217] 614 607 755 866 787 767 789 633 681 835 853 833 870 724 823 872 869 613
## [235] 731 690 895 731 870 806 795 730 842 882 809 804 701 841 820 641 684 895
## [253] 641 650 765 699 674 729 754 895 827 634 750 694 861 768 619 673 722 624
## [271] 737 690 875 664 750 775 819 683 607 855 632 796 734 726 635 871 771 665
## [289] 758 737 760 783 659 876 794 728 817 669 673 807 702 736 855 845 785 701
## [307] 823 778 868 715 766 694 855 693 708 694 862 863 621 629 651 850 856 627
## [325] 681 610 779 843 803 789 834 814 812 633 724 629 672 773 635 868 636 893
## [343] 884 618 783 756 674 879 880 638 874 786 684 618 847 717 733 633 711 647
## [361] 647 804 651 859 703 705 775 888 809 882 701 809 731 666 787 736 824 632
## [379] 617 720 763 794 784 814 616 655 818 890 675 814 607 866 881 761 819 633
## [397] 757 624 677 635 780 680 806 815 794 654 647 652 601 663 669 846 731 688
## [415] 673 894 605 727 718 849 614 829 863 717 874 600 805 875 658 753 615 809
## [433] 784 887 633 807 848 829 789 705 848 739 769 850 615 794 668 683 657 664
## [451] 773 776 892 874 850 747 714 603 716 836 741 850 768 871 656 807 794 737
## [469] 720 628 787 877 798 630 859 603 649 741 725 795 867 884 642 832 609 864
## [487] 874 741 686 745 629 803 608 816 865 875 829 602 604 754 703 648 670 890
## [505] 685 826 702 649 800 869 639 702 705 630 700 731 630 862 784 892 747 659
## [523] 773 899 752 619 807 765 831 723 709 799 722 734 628 645 759 688 768 774
## [541] 692 690 855 893 690 662 700 660 706 893 669 859 665 856 784 764 614 896
## [559] 679 676 612 633 807 674 760 839 688 718 863 721 865 810 681 886 741 865
## [577] 690 810 654 872 846 872 812 688 779 730 850 728 726 765 877 870 805 751
## [595] 761 695 726 780 829 662 626 879 802 869 771 838 607 616 747 632 887 858
## [613] 670 737 695 869 805 863 748 714 777 772 814 873 735 697 673 735 849 754
## [631] 821 797 789 612 653 714 897 696 678 856 672 656 894 894 703 856 692 804
## [649] 717 763 731 661 890 687 633 795 734 814 846 820 732 745 808 622 823 672
## [667] 855 874 637 764 601 888 663 628 745 843 647 657 665 714 867 702 892 655
## [685] 790 889 800 830 871 830 657 726 766 851 867 649 802 639 698 823 723 793
## [703] 690 791 634 775 681 691 809 711 732 897 661 855 702 708 724 895 703 660
## [721] 832 692 864 843 791 740 746 662 618 846 736 603 736 643 746 832 803 864
## [739] 716 680 675 700 666 794 650 686 713 821 739 855 875 876 678 767 601 677
## [757] 606 662 749 877 657 717 881 801 894 701 845 624 745 686 821 736 679 808
## [775] 749 619 688 843 602 813 700 682 751 882 803 854 653 607 702 610 675 747
## [793] 642 827 655 674 606 676 622 757 825 786 770 616 656 865 668 691 809 661
## [811] 889 726 878 873 873 815 665 653 868 744 832 852 863 843 682 701 765 842
## [829] 716 812 889 656 859 746 619 656 783 751 624 602 706 725 671 829 628 835
## [847] 743 809 842 800 805 704 889 784 701 605 745 704 741 797 671 675 875 770
## [865] 802 749 603 750 717 647 657 723 610 698 755 810 827 702 890 623 883 781
## [883] 861 831 680 842 647 763 631 804 689 603 647 831 738 805 886 716 836 703
## [901] 628 734 745 719 668 813 610 686 756 770 805 767 812 790 870 679 812 883
## [919] 790 605 802 615 695 890 799 808 800 739 608 684 749 844 765 767 739 820
## [937] 699 768 836 609 636 833 807 730 794 813 894 899 894 691 744 800 613 758
## [955] 785 893 881 609 613 646 735 696 607 601 822 896 785 802 879 611 807 734
## [973] 759 728 612 704 864 705 876 710 612 835 781 793 645 681 882 714 773 710
## [991] 635 744 656 881 620 871 616 711 773 736hist(cl)
basicStats(cl)## cl
## nobs 1000.000000
## NAs 0.000000
## Minimum 600.000000
## Maximum 899.000000
## 1. Quartile 672.000000
## 3. Quartile 820.000000
## Mean 747.566000
## Median 746.000000
## Sum 747566.000000
## SE Mean 2.763147
## LCL Mean 742.143762
## UCL Mean 752.988238
## Variance 7634.982627
## Stdev 87.378388
## Skewness 0.019067
## Kurtosis -1.207907Trung bình quán bán được khoảng 751 phần Cao Lầu trong 1 tháng với độ lệch chuẩn là 87 phần. Trong đó, số phần Cao lầu nhiều nhất quán bán được trong 1 tháng là khoảng 900 phần và ít nhất khoảng 600 phần.
4.3 Mô phỏng biến Bánh mỳ:
Số bánh mỳ Hội An bán được trong tháng có phân phối Poisson với tham số lambda=150
library(ggplot2)
bmy <- rpois(1000, lambda = 150)
bmy <- table(bmy)
bmy <- as.data.frame(bmy)
ggplot(bmy,aes(bmy, Freq)) + geom_col(fill= 'red') + geom_point() + geom_line(aes(as.integer(bmy), Freq), color ='green')
bmy <- rpois(1000, lambda = 150)
basicStats(bmy)## bmy
## nobs 1.000000e+03
## NAs 0.000000e+00
## Minimum 1.080000e+02
## Maximum 1.960000e+02
## 1. Quartile 1.410000e+02
## 3. Quartile 1.580000e+02
## Mean 1.499820e+02
## Median 1.500000e+02
## Sum 1.499820e+05
## SE Mean 3.927850e-01
## LCL Mean 1.492112e+02
## UCL Mean 1.507528e+02
## Variance 1.542800e+02
## Stdev 1.242095e+01
## Skewness 1.159300e-01
## Kurtosis 9.068600e-02Số bánh mỳ mà quán bán được trung bình trong 1 tháng khoảng 150 phần với độ lệch chuẩn là 12 chiếc. Trong đó, quán bán được nhiều nhất là khoảng 204 chiếc và bán được ít nhất khoảng 114 chiếc.
4.4 Mô phỏng biến Cơm hến:
Số phần Cơm Hến bán được trong một tháng có phân phối chuẩn với kỳ vọng Mu=150 và sigma^2=15
ch <- rnorm(1000, mean = 150, sd = 15)
ch <- round(rnorm(1000, mean = 150, sd = 15),0)
ch## [1] 174 155 133 145 131 158 147 165 158 167 134 145 159 144 175 149 145 144
## [19] 162 180 140 132 165 150 135 179 167 140 174 162 158 170 178 163 172 153
## [37] 161 166 133 146 140 155 163 148 153 120 162 147 136 136 142 143 143 168
## [55] 146 157 181 143 155 139 191 149 118 177 135 162 147 158 165 145 131 169
## [73] 152 118 154 121 142 145 166 152 173 158 160 138 167 158 148 129 156 143
## [91] 164 140 140 164 147 154 107 151 130 164 163 134 166 127 148 164 135 155
## [109] 140 113 175 151 156 146 170 158 128 159 157 142 159 155 144 160 150 158
## [127] 146 136 157 122 168 161 165 173 137 166 147 145 161 147 131 140 115 150
## [145] 134 153 133 150 149 171 144 142 138 173 178 159 144 158 143 168 130 150
## [163] 140 150 139 158 139 168 168 137 157 125 172 162 164 168 147 138 151 154
## [181] 156 134 165 150 151 141 123 157 134 165 156 144 147 129 125 147 133 132
## [199] 158 155 129 154 148 142 128 166 156 151 171 157 178 162 161 158 114 154
## [217] 150 166 152 164 163 147 174 163 138 149 187 161 183 159 152 184 143 148
## [235] 155 140 166 118 147 144 124 143 137 158 156 132 148 131 154 154 129 132
## [253] 168 145 153 151 159 128 147 135 169 192 131 163 115 157 161 134 162 189
## [271] 133 154 146 182 152 153 157 138 146 142 149 148 130 179 164 142 145 129
## [289] 157 101 157 165 147 140 193 169 132 167 163 173 146 153 169 147 162 158
## [307] 148 131 164 126 139 151 153 149 159 121 153 156 160 181 162 176 123 142
## [325] 127 145 146 159 168 139 162 156 142 174 152 149 144 127 152 124 152 154
## [343] 148 143 159 158 158 153 173 133 150 147 161 168 145 149 117 155 147 137
## [361] 137 130 154 130 162 159 148 151 130 166 141 154 141 173 147 162 137 167
## [379] 161 139 155 140 167 129 166 164 161 151 141 158 143 175 167 146 134 169
## [397] 139 145 134 164 183 150 163 158 127 161 170 164 181 154 163 165 138 145
## [415] 155 127 147 150 166 149 150 149 162 145 156 147 176 141 140 150 142 123
## [433] 130 131 162 177 154 147 150 137 132 152 156 130 156 160 158 129 148 133
## [451] 128 159 134 171 138 135 147 163 148 150 154 162 166 137 138 145 141 153
## [469] 136 155 131 166 156 149 140 144 156 151 131 148 161 144 163 157 153 136
## [487] 149 150 139 131 158 166 145 168 138 148 164 146 147 147 158 163 149 141
## [505] 144 132 144 151 165 152 132 115 174 147 155 165 157 140 149 136 127 165
## [523] 146 144 133 126 154 163 159 145 162 136 159 163 169 159 145 125 181 148
## [541] 159 163 162 175 131 149 166 149 116 150 146 126 156 144 144 132 179 154
## [559] 131 146 131 147 151 162 133 166 169 171 154 153 144 146 142 154 156 141
## [577] 131 161 161 146 154 132 131 178 144 143 149 144 154 159 153 154 145 143
## [595] 162 163 160 139 152 114 168 161 176 160 137 167 152 162 157 127 180 147
## [613] 142 173 146 144 165 170 159 151 154 142 115 116 159 165 161 146 123 162
## [631] 157 139 173 138 146 117 162 176 141 127 161 155 167 156 148 170 161 148
## [649] 159 159 146 162 134 162 157 148 158 155 151 146 156 151 159 158 166 162
## [667] 151 158 143 152 153 153 152 133 138 133 111 146 170 170 149 181 138 155
## [685] 144 148 152 153 126 167 151 149 148 152 157 150 139 153 165 144 172 141
## [703] 143 154 122 163 180 144 160 150 114 117 159 154 176 183 166 144 142 122
## [721] 152 161 154 124 172 130 129 148 146 157 196 175 139 151 139 140 145 138
## [739] 134 162 151 153 160 113 157 121 122 125 153 150 154 143 139 155 158 146
## [757] 162 172 125 182 159 143 138 154 157 119 132 158 160 131 114 139 152 146
## [775] 150 148 102 141 129 165 160 160 171 132 156 168 158 150 163 107 165 182
## [793] 145 152 170 138 123 135 163 127 145 131 154 161 158 160 183 100 136 132
## [811] 110 138 149 132 150 128 182 162 159 153 141 136 137 119 164 140 183 144
## [829] 133 173 131 162 144 166 153 167 157 159 153 147 162 137 137 170 164 146
## [847] 131 163 164 140 151 153 127 155 127 144 143 166 145 182 142 149 137 125
## [865] 153 121 147 155 126 114 128 147 157 153 143 110 141 143 142 136 159 146
## [883] 123 167 119 173 169 159 146 179 159 157 151 148 151 166 166 184 143 148
## [901] 169 135 144 108 138 134 151 162 136 137 129 145 124 134 119 166 137 158
## [919] 143 147 153 135 148 129 157 126 157 145 150 165 162 159 175 180 164 149
## [937] 153 147 142 168 153 160 140 140 171 144 126 157 170 163 152 138 134 188
## [955] 112 155 144 137 133 167 145 149 184 137 128 149 149 147 153 148 168 162
## [973] 194 165 174 164 131 182 146 138 160 163 149 126 149 152 147 139 170 179
## [991] 149 136 163 140 148 143 131 146 134 140hist(ch)
basicStats(ch)## ch
## nobs 1000.000000
## NAs 0.000000
## Minimum 100.000000
## Maximum 196.000000
## 1. Quartile 140.000000
## 3. Quartile 161.000000
## Mean 150.040000
## Median 150.000000
## Sum 150040.000000
## SE Mean 0.490888
## LCL Mean 149.076709
## UCL Mean 151.003291
## Variance 240.971371
## Stdev 15.523253
## Skewness -0.149040
## Kurtosis 0.076862Trung bình quán bán được khoảng 149 phần Cơm hến trong 1 tháng với độ lệch chuẩn là 14 phần. Trong đó, số phần Cơm hến nhiều nhất quán bán được trong 1 tháng là khoảng 197 phần và ít nhất khoảng 108 phần.
4.5 Mô phỏn biến Nước Mót:
Số ly Nước Mót bán được trong tháng có phân phối nhị thức với tham số n=6000 và xác suất bán được p=0.5 (trong 1 tháng chọn ngẫu nhiên 6000 người vô quán và trong 6000 người đó sẽ có khoảng 3000 người gọi Nước mót).
nuoc <- rbinom(1000, size = 6000, prob = 0.5)
nuoc## [1] 3008 3000 2946 2980 2988 2953 2992 3098 3069 3070 3018 2926 3006 3019
## [15] 3043 3002 2981 3029 3042 3012 2968 2991 3055 3027 3014 2952 3064 2977
## [29] 3014 2981 3024 3005 3002 3033 2915 3018 2998 3002 3033 2958 3051 3047
## [43] 2989 2990 2977 3017 3011 2964 3016 2999 2960 3043 3055 2992 3030 2983
## [57] 3055 2944 2993 3013 2972 3010 3027 2974 3019 3040 3005 2970 2943 2970
## [71] 2951 2992 2976 3044 3002 2969 2963 2970 3026 3033 2956 3029 2985 3037
## [85] 3052 2948 3006 3042 3002 2971 3028 2943 2972 2985 3030 2907 3029 2972
## [99] 2992 2969 2963 2931 2920 2987 3033 2990 2984 2983 3007 3042 3084 3013
## [113] 2993 3027 2981 3016 3020 3040 2976 2995 2943 2927 3057 2954 3044 3014
## [127] 3073 2995 3050 3028 2964 3041 3067 3065 2992 3006 3034 3003 2960 2973
## [141] 3034 3026 3008 2993 2994 2939 2963 3009 2970 2984 2908 3053 3039 3066
## [155] 3020 3012 3002 3047 2907 3002 2986 3035 3010 3005 2963 3000 3004 3000
## [169] 3062 2965 3018 2989 2934 3043 3004 2950 2961 2971 3045 2979 2959 2907
## [183] 3084 2923 2963 3001 3019 2959 3050 3038 2946 3002 2957 3014 2987 3002
## [197] 2960 2998 3020 2969 2948 3006 2948 3025 2987 3022 3102 3020 2960 2983
## [211] 3009 3024 3004 2972 2961 2956 2977 2927 2956 2963 2956 3017 3057 3008
## [225] 2995 3046 2998 3003 2976 2987 3005 3023 2980 3046 3016 3089 2959 3063
## [239] 2985 2981 3018 3015 2977 2959 3039 3021 2973 2940 3005 2908 3099 2945
## [253] 2969 3027 3066 2905 2994 2955 3068 2988 3069 3027 2991 3076 2986 2983
## [267] 2978 2974 3029 3032 2938 3048 3048 3037 2986 3040 2939 3040 2941 3003
## [281] 3067 2966 2998 3030 2986 2997 3057 3000 3005 3041 2955 3020 3038 2997
## [295] 2899 3031 2961 2916 2989 3066 2974 3073 2959 3026 3034 2957 2999 2959
## [309] 2987 3018 2989 2934 2993 3025 3032 3009 2950 2967 2987 2955 3006 3061
## [323] 3003 2984 3019 3064 2972 3103 3064 2914 3084 3022 2993 3025 2985 3019
## [337] 2978 3006 2982 3032 2997 3010 3050 2953 3061 3051 2942 3018 2926 2943
## [351] 2994 2966 2934 3041 3013 2937 2974 3022 2995 2985 3027 2987 3012 3030
## [365] 3000 3036 3056 3000 3022 2920 3055 3060 3022 2968 2923 3026 3032 2988
## [379] 3007 2999 2974 3011 2981 3035 2927 2975 3024 2966 3052 3034 2981 2866
## [393] 3058 3004 2997 3061 2980 2975 3080 2953 3034 2998 2970 3036 2993 2965
## [407] 3010 3013 3042 2982 3031 2992 3062 3017 2938 2977 3027 3030 3018 3026
## [421] 3024 3013 3012 3001 3034 2986 2992 2965 2993 3066 2998 2989 2976 2990
## [435] 2900 3007 3004 3039 2961 3005 3049 3022 2990 3012 3079 2932 2988 2976
## [449] 2977 2944 3016 3022 3032 3037 3046 3013 3017 2983 3053 2951 2965 2948
## [463] 2987 3054 2973 3019 2963 3007 2976 3058 3051 3022 3008 2968 2978 3117
## [477] 2916 2992 3044 2936 2974 2984 3012 2978 3059 2977 2982 3020 2941 2997
## [491] 2952 2966 2978 2972 2995 3013 2968 2981 2948 3010 3021 3027 2934 3039
## [505] 3063 2955 3018 3069 2964 2987 2988 3032 3024 3035 3012 3045 2985 3086
## [519] 3046 3047 3044 2995 3000 2965 2959 2964 3011 2981 2989 2974 3012 2888
## [533] 2986 2984 3024 2984 3010 3021 3002 2979 2968 3018 2945 2994 3055 3050
## [547] 2997 3008 2970 3028 2990 2979 3094 2948 2967 2988 2958 2988 3115 3024
## [561] 3061 2969 3022 3001 2960 2973 3024 2981 3026 3065 2979 2936 2963 2966
## [575] 3008 2926 2987 2989 3010 2998 2991 2974 2947 3010 3013 2990 3057 3055
## [589] 2896 2995 2996 2987 3025 3008 3033 3006 2978 3030 2987 3025 2976 2950
## [603] 2986 2993 2986 2987 2978 3020 3002 2958 3053 2946 3017 3013 3040 2935
## [617] 3001 2975 2956 3092 2995 3035 2895 2990 2960 3011 2992 2980 2999 3001
## [631] 2967 2957 3038 3063 3022 3024 2986 2954 3021 2960 2931 3024 3040 3010
## [645] 3031 3018 3015 3042 2985 3001 3051 2974 3029 3027 2999 2923 3013 3046
## [659] 3012 2987 3004 2962 2957 3022 3008 3071 3023 2953 3047 2981 2977 3030
## [673] 3044 2944 3002 2979 3051 2998 3023 3072 3093 3057 3035 3004 3042 2972
## [687] 3029 3007 2938 3052 3000 2989 3091 2934 3041 2975 2898 2962 2982 3013
## [701] 3008 3049 2978 2968 2985 2999 2920 2955 3073 3035 2935 2977 2983 3031
## [715] 2963 2996 2979 2993 2919 3007 3018 3027 2913 3031 3014 2984 3030 3031
## [729] 2969 3002 3017 2969 3016 2994 3038 3039 3028 3001 3047 3023 2962 2989
## [743] 2990 2979 2941 3040 2999 2924 3016 2951 2986 3014 2965 2953 3032 3005
## [757] 3040 2981 2990 3044 3033 2867 3062 2986 2959 3007 3005 3003 2938 3044
## [771] 3050 3071 3037 2996 3023 2996 2945 2997 2967 3026 3000 3020 3064 2979
## [785] 2969 2991 2940 3013 3059 3032 3005 3000 3013 2993 3027 3003 2980 2919
## [799] 3005 3025 2972 3041 3057 2998 2965 2991 2942 3065 3033 2947 3016 3010
## [813] 3005 2946 3000 2995 2973 2987 3023 3021 3008 2973 3006 3010 2988 3015
## [827] 2967 3096 2953 3057 3020 3021 2967 3022 2993 3008 2987 3016 3060 2937
## [841] 3012 2948 3045 2946 3010 3007 3062 3062 2972 2980 2971 2992 2931 2952
## [855] 3045 3047 2972 2997 2982 3088 2955 3035 2936 3010 2965 3018 3000 2986
## [869] 2950 2998 2951 3000 3019 3000 3008 2961 3014 2992 3043 3012 3022 2973
## [883] 2939 3027 3003 3002 2921 2941 2972 2966 3024 3015 2980 2962 2995 3025
## [897] 2977 3003 3052 2986 3033 2951 2912 2959 3022 2978 2967 2997 2937 3010
## [911] 2985 3000 3003 2969 2963 3011 3073 2946 2923 3029 2978 3040 3011 2991
## [925] 3047 3055 2988 3021 3011 2965 2996 3029 2979 3033 2980 3069 3064 3038
## [939] 2968 2976 3019 3031 3048 2993 2999 3092 2951 2984 3035 2986 3003 3000
## [953] 2955 2966 2994 3054 2960 3006 3081 3001 2988 2902 3063 2951 2967 3003
## [967] 3055 2995 2930 3009 2973 3022 2990 2974 3021 2908 2989 3041 2910 3065
## [981] 3011 3012 2981 3036 2979 2978 3005 3010 2983 3068 3041 2952 3036 2979
## [995] 3014 2997 3019 2990 3045 3034hist(nuoc)
basicStats(nuoc)## nuoc
## nobs 1.000000e+03
## NAs 0.000000e+00
## Minimum 2.866000e+03
## Maximum 3.117000e+03
## 1. Quartile 2.974000e+03
## 3. Quartile 3.026000e+03
## Mean 2.999633e+03
## Median 3.000000e+03
## Sum 2.999633e+06
## SE Mean 1.245602e+00
## LCL Mean 2.997189e+03
## UCL Mean 3.002077e+03
## Variance 1.551524e+03
## Stdev 3.938939e+01
## Skewness -8.649700e-02
## Kurtosis 2.412300e-02Số ly Nước mót mà quán bán được trung bình trong 1 tháng cho 6000 khách khoảng 2999 ly với độ lệch chuẩn là 38 ly. Trong đó, quán bán được nhiều nhất là khoảng 3124 ly và bán được ít nhất khoảng 2873 ly.
4.6 Mô phỏng biến tuổi thọ:
Tuổi thọ của bếp điện có phân phối mũ với tuổi thọ trung bình 6 năm nên phân phối mũ với tham số lambda=1/6.
tt <- rexp(1000, rate = 1/6)
tt## [1] 0.059188124 1.276885861 0.814081002 0.088095662 1.318525457
## [6] 2.075574657 4.726301774 9.899359524 0.179508241 2.507723118
## [11] 1.616340643 13.202891047 0.078330464 2.484191987 5.972946323
## [16] 0.617555690 24.751010451 2.390356222 0.346432857 11.839794471
## [21] 11.300307737 11.082223978 2.905172057 0.851322736 7.110432721
## [26] 12.687492358 0.700267065 9.897309076 1.035633944 11.191897787
## [31] 13.802746117 0.686769051 9.817249171 4.179114300 44.658842926
## [36] 6.738653037 8.104654055 17.519932239 6.648334489 1.505021006
## [41] 1.649714283 0.582358570 1.330925823 7.262840099 4.392433971
## [46] 1.794306477 5.501703794 11.678254867 3.724098662 3.509477030
## [51] 0.461958045 12.777810867 0.977525440 0.612910573 0.895499120
## [56] 0.125399596 8.154362507 0.244767553 0.977092900 8.731462266
## [61] 0.448850404 17.887845172 1.573284800 8.525020288 2.421167300
## [66] 3.674140119 9.071581483 6.085429363 4.726703477 3.533916545
## [71] 1.304598690 8.590035952 2.757432463 1.100185177 8.620425091
## [76] 10.165096588 5.607251672 0.119552216 1.146925016 12.494715045
## [81] 1.698723922 9.648543566 1.515159811 4.560280069 6.965082317
## [86] 1.970083471 4.970717042 1.146542475 5.351346027 5.853161441
## [91] 1.159352981 1.045222512 0.479813551 1.016417485 12.718244465
## [96] 3.553762397 6.445198305 19.931903287 6.954641946 5.117002762
## [101] 5.690795677 10.299510039 3.164617379 0.409911565 3.044880985
## [106] 4.672386383 6.956889524 4.477813912 0.052703842 6.520222412
## [111] 0.383577352 17.694071039 9.080348624 2.776411769 3.480771032
## [116] 20.119371950 9.463369511 16.113051239 9.666015020 7.095194969
## [121] 7.593703547 5.045218606 14.506515685 0.201726926 1.301468773
## [126] 0.132608159 9.898729264 15.770184036 5.291613566 8.776967358
## [131] 2.274096083 9.232874593 0.481573032 22.401134497 3.065827386
## [136] 11.921237848 4.931369450 1.188101870 4.391110722 4.101583131
## [141] 1.908785686 3.524631548 1.489392861 1.410504182 5.229778591
## [146] 3.815691946 10.989573709 13.551849540 1.375061357 22.938835623
## [151] 4.664682157 4.760436522 11.132412560 3.991100189 1.184645947
## [156] 2.456704454 0.712137806 9.390174195 7.629091831 14.120504099
## [161] 4.416479624 1.809150998 9.505608942 0.983700145 11.872988864
## [166] 2.266189529 8.719015874 9.049333959 10.511841036 4.140998734
## [171] 12.986650015 0.600726611 1.575960321 5.269384080 0.064947923
## [176] 13.655176620 1.938488164 2.837475311 2.280510337 1.788359707
## [181] 3.961331078 4.823850212 1.842877012 3.531954582 2.416195247
## [186] 14.139336425 0.552579382 1.929880335 3.514709112 15.173669244
## [191] 7.240587802 5.590305250 0.668421533 1.171880281 4.530187396
## [196] 17.651828768 0.905862367 2.112415287 1.773149251 0.647951173
## [201] 6.013026016 5.866412049 5.700713542 3.743787252 13.217820430
## [206] 2.421863040 20.602854340 4.788052766 0.987120626 11.262450087
## [211] 0.007012352 0.102886112 3.485818855 4.339697999 7.135659160
## [216] 9.746213425 13.960128077 11.000421289 4.094216905 1.613746924
## [221] 0.098340567 8.716020757 10.419145066 7.752273140 1.247553536
## [226] 2.212987167 5.413059158 5.821053732 3.075059499 23.548435571
## [231] 2.482183693 6.713863397 9.391148306 1.222003232 3.925438634
## [236] 4.143419534 0.584770206 8.298524922 1.184365520 16.917440771
## [241] 0.956183868 2.136774500 2.017040108 1.538854807 5.179569614
## [246] 1.598689579 13.753080468 0.114831957 5.715346634 4.672868851
## [251] 0.749494188 24.127662184 2.145320035 8.759152539 3.500267431
## [256] 10.731065969 14.125604901 1.169628440 0.193150487 2.665944401
## [261] 1.098738465 7.821352020 0.184181807 5.289130699 4.361805445
## [266] 1.656213017 0.910851568 0.030387886 4.973323725 6.184811571
## [271] 0.771995204 0.704840414 17.954262881 4.647839528 3.262513309
## [276] 2.529240555 1.518190353 22.970162788 0.206740189 17.914489686
## [281] 0.235030896 0.156245970 10.803317029 1.812833245 4.934318940
## [286] 4.574441274 10.797635413 0.464210384 2.950735291 12.067975904
## [291] 1.823022248 5.505690040 4.289832840 10.905394493 7.189713696
## [296] 0.986066492 3.836171717 1.199146358 1.561841662 20.576964675
## [301] 4.843094336 5.525826009 1.444778388 8.802520362 7.923826882
## [306] 8.716494843 2.664628641 1.137897382 2.244192473 1.552650928
## [311] 22.365665617 6.914058307 0.093976540 4.714474305 5.153768057
## [316] 6.144880883 1.025406805 0.889078651 11.421383164 0.694086176
## [321] 1.701728848 5.649270857 4.180057210 2.188088876 7.029756842
## [326] 4.149065992 1.923332336 1.882187102 0.445080433 1.797593093
## [331] 1.304918689 2.612953492 4.786080686 9.889710154 1.444862889
## [336] 7.721183242 2.340613070 1.446814364 5.393523857 1.195926148
## [341] 0.076673694 2.822532612 0.785856807 9.854423571 7.749034127
## [346] 9.460338480 8.343441566 2.579157827 9.770375348 6.586230671
## [351] 1.958019605 11.094209910 4.158369170 9.780471190 0.842479873
## [356] 1.264649210 1.444673189 12.862764510 5.282030336 3.448507355
## [361] 0.143686168 3.956185665 9.423839729 0.512383782 1.196373763
## [366] 0.087680778 4.036847790 8.380150102 7.877149025 1.229685254
## [371] 5.625059343 5.332934350 1.416964743 3.393404461 36.360384373
## [376] 4.389852214 4.700341572 18.608139629 1.392712299 3.283748229
## [381] 0.730119003 2.935817565 4.681450719 5.616302813 6.237276071
## [386] 17.040666549 5.306508520 4.923522306 13.848914556 0.277700209
## [391] 3.191060350 4.433253550 0.720617584 2.521425454 4.186150884
## [396] 24.278228524 13.631615292 4.700789513 5.000129705 3.440778701
## [401] 0.797948839 2.307353636 11.321963800 1.589581423 12.341619830
## [406] 3.688622757 6.047316793 12.719249430 11.917294375 5.811098959
## [411] 2.341472802 3.318454456 4.457451867 3.680543401 2.799547473
## [416] 6.136295775 2.264944760 4.225930839 0.548648684 0.056713851
## [421] 2.399979481 0.038775885 2.889768161 0.978332293 8.164905471
## [426] 12.792643348 0.410799145 0.158898739 4.071124521 0.306712796
## [431] 1.527494523 8.505000124 15.100122363 4.659021194 3.406645143
## [436] 1.525659089 1.073484477 2.909507892 5.296999551 12.174278195
## [441] 17.715170457 6.818308130 2.730505753 0.473767018 3.672281315
## [446] 17.195057987 14.324827690 2.536160409 0.198820858 9.329064330
## [451] 5.088552929 23.251676886 13.535029459 0.221543437 4.877380640
## [456] 12.933413218 1.766249727 9.972786589 8.890761204 1.009493226
## [461] 2.700506797 19.681464984 10.367256369 11.722376324 18.957420587
## [466] 6.201711712 7.342217568 9.298475465 2.037968625 1.433251163
## [471] 6.301149588 7.555788941 1.302101612 0.172565439 4.635106601
## [476] 10.273224659 0.414772213 0.463528428 2.161413928 17.879026740
## [481] 0.032475645 18.436599567 8.225208621 1.148468326 13.594011794
## [486] 1.438297170 5.234765818 7.729795022 0.276880403 8.477546207
## [491] 2.037184808 0.374745710 1.159191362 9.837698657 36.102803616
## [496] 3.501846822 3.793846037 2.352101211 2.047284077 9.069828630
## [501] 3.778954254 3.144275514 1.216010942 3.856439053 1.663217051
## [506] 2.236679521 2.658083710 0.487459375 12.503522727 5.756163164
## [511] 2.900407009 11.420954883 10.955254610 7.276730960 1.667197991
## [516] 5.892856566 1.680088956 0.041197739 2.761792663 6.688806463
## [521] 0.028772905 0.020576529 0.158974771 9.979289755 3.900050274
## [526] 7.212652785 2.154599009 2.099706447 4.086980662 4.249232869
## [531] 26.450398317 4.743653710 5.270304584 3.849638835 14.371182648
## [536] 7.321381174 12.369800291 10.936687488 12.641407003 1.197290607
## [541] 1.743499838 4.640348532 13.896154056 3.635691558 2.830874576
## [546] 8.920566316 3.687045090 9.363333418 4.206597699 2.561255233
## [551] 27.899654677 9.203148532 7.834753381 0.681627000 18.866792231
## [556] 3.960380095 26.503381780 4.166915289 1.843661869 5.662480810
## [561] 3.826929140 2.004050612 0.920434758 13.593027714 0.089619965
## [566] 10.385722875 0.352684805 1.990382810 1.846367746 7.140478184
## [571] 0.735043106 6.725901134 7.971577579 4.126347587 9.031522929
## [576] 1.737764096 2.564895028 12.759455200 3.827339150 1.355941750
## [581] 1.431127575 11.008194621 8.878536622 0.797525352 3.227859939
## [586] 0.002699763 4.555375139 0.445976833 11.407729766 1.811344961
## [591] 0.758229548 4.458432711 13.868345825 0.906349931 7.020689763
## [596] 6.411726566 11.097025814 6.304529937 0.494539848 16.120819877
## [601] 15.371826130 2.409833883 4.813720071 4.919087408 2.296212103
## [606] 0.901984385 25.185380381 8.503871552 1.693566131 5.540078088
## [611] 2.849890728 10.656536210 9.918568240 11.793700363 8.435625649
## [616] 0.748666564 31.561619917 5.929339790 0.172377602 1.452675164
## [621] 4.331877228 7.788014011 5.006042852 2.557952221 1.922979833
## [626] 3.233272791 10.708022017 1.315040049 2.366934380 5.840781837
## [631] 5.441668692 0.498381353 15.865597050 24.333454850 26.103226057
## [636] 1.102622294 1.350689773 8.593952534 4.289423506 2.259795574
## [641] 2.157048298 0.724701044 2.409302607 18.956598063 3.513278668
## [646] 1.694537216 16.232642244 1.868766856 0.597930578 4.663989893
## [651] 11.095974650 4.925002468 4.791203845 0.232594647 9.652319669
## [656] 6.244354740 1.396698827 13.547636748 5.059657255 2.143668552
## [661] 4.174497211 23.018990221 6.740705659 3.676183261 2.757581815
## [666] 0.561882617 5.628064818 4.302116342 1.965311721 15.336221374
## [671] 29.390500706 1.366016813 1.335224298 1.362060986 5.091507649
## [676] 2.980757480 8.746604919 6.406695562 8.709511835 1.248220565
## [681] 10.546927770 18.648180528 1.001114540 8.840652476 1.788227848
## [686] 1.152580275 8.400848478 5.104459109 6.179508651 0.910963682
## [691] 7.050722010 17.362070094 5.756917080 3.958166426 24.929082724
## [696] 9.982994567 4.700856451 4.273723114 1.911559794 3.957706285
## [701] 3.093024957 8.575330942 10.045246191 1.836221505 0.673668412
## [706] 3.464624508 0.407305936 0.968922168 13.346132520 0.126983656
## [711] 4.189079645 6.916493685 10.055479348 10.435262672 4.207314474
## [716] 0.368018369 11.060458224 0.416772597 8.362635041 22.950660946
## [721] 5.347856858 8.632297664 8.105046624 8.274489378 0.214435553
## [726] 9.775726061 3.799008133 3.031640268 21.540860437 13.568820561
## [731] 1.036049097 3.441617260 0.707637744 8.686892229 7.731614370
## [736] 20.260577648 9.663625788 26.114289990 6.492333133 6.948866905
## [741] 1.379725011 5.902420727 1.875106573 1.752707690 6.736938257
## [746] 1.453873417 7.009533645 1.645664373 3.460031895 5.111137831
## [751] 1.275827764 4.444740336 1.718953673 2.592125114 5.382398936
## [756] 5.638209343 5.988908300 0.096032537 6.771776380 2.713631163
## [761] 5.572304243 8.271048987 8.676280876 0.683640595 5.060515237
## [766] 4.449756313 17.689771895 3.942453755 1.792317308 2.521678007
## [771] 23.236826452 0.612517683 4.698493469 1.943196964 0.884502746
## [776] 22.321846637 4.572889813 3.247322895 0.426175947 7.239744236
## [781] 14.367467141 0.817062261 3.070411664 9.198292433 2.173551206
## [786] 0.383341779 10.855950694 4.901293343 0.131391923 7.800433203
## [791] 4.006247291 0.686207176 3.142127456 4.074352671 2.557440944
## [796] 5.989758677 5.835750650 2.093480865 2.045567807 3.375628690
## [801] 4.704027688 5.336730765 10.469919528 1.039117905 0.641923808
## [806] 0.237354226 11.198657792 2.359785983 3.916786886 1.106460910
## [811] 2.366513335 3.337468095 1.572505590 6.648217575 8.510348360
## [816] 2.653733594 8.410886082 9.813425939 11.861501198 11.419007141
## [821] 2.577966638 1.123386287 11.585895359 5.671281248 10.923023555
## [826] 0.881149258 0.637081652 2.355020086 6.258436988 1.687417117
## [831] 0.852028797 0.468555171 0.859730346 0.641727568 2.773776127
## [836] 9.589287467 9.391282137 0.032561188 4.701354174 5.024244734
## [841] 6.741305932 3.487232354 9.366228211 1.754620800 0.079089641
## [846] 10.302276193 5.201662364 5.628036061 0.045085665 2.585101616
## [851] 3.732201445 0.897353729 3.351133880 0.605371197 5.655612275
## [856] 1.675349968 0.997889208 1.437113363 12.456103742 0.341092214
## [861] 0.830115940 1.907836438 2.347819046 2.189467938 9.921723301
## [866] 1.016997873 11.973485950 0.723168796 3.301580211 2.511607497
## [871] 2.799043749 16.573562632 4.490928214 8.840190924 12.462574895
## [876] 4.822419617 12.398191809 7.319729905 5.980389224 1.676451095
## [881] 1.470189581 3.767944235 2.118442219 0.919111823 4.659161602
## [886] 6.589110452 20.760340153 3.181147699 3.654086220 6.571702134
## [891] 5.333781363 3.288725453 13.785837728 7.682390805 6.399810265
## [896] 8.256409630 10.250471223 4.775479686 15.981487121 19.525960653
## [901] 7.753949420 8.040055884 15.877106275 3.062048477 20.859377601
## [906] 0.749685475 6.596791094 5.222531794 0.605015676 12.629194461
## [911] 7.159839006 9.875292231 3.907461746 1.931900684 1.952228124
## [916] 4.697495827 15.545473571 4.452816406 11.868543852 2.599219632
## [921] 3.656105840 3.762759190 1.814459845 4.007913929 5.803763028
## [926] 23.319908592 3.909801018 4.058644786 5.626962840 0.869496680
## [931] 0.152748571 6.748860083 0.396147280 2.924481492 2.866109854
## [936] 2.633709960 0.336121512 1.611519932 0.424603913 8.119228843
## [941] 9.561159230 11.282106958 3.650325503 13.925587322 14.418275210
## [946] 0.877560382 1.492096933 2.231399476 0.082675556 9.702911127
## [951] 10.767396613 8.627270489 3.148281505 2.042088364 1.062793256
## [956] 0.882660189 12.184299085 1.742379192 1.876901784 3.258271731
## [961] 0.634603037 1.714790753 4.046915900 0.850876537 3.243506150
## [966] 0.926267543 6.618553634 5.783113232 11.200218121 11.612374886
## [971] 1.977577466 7.910665121 0.401476960 15.345137086 1.242709234
## [976] 2.385494982 3.238568567 0.306476926 3.834000338 6.815701848
## [981] 8.055066973 2.642994649 2.163547098 1.684827603 9.599042247
## [986] 1.994135083 7.998732123 0.386667193 1.023257369 10.657971571
## [991] 1.132731554 2.102140239 7.926151156 4.831608573 0.529396883
## [996] 5.097353421 2.618758854 22.309294435 0.848296216 4.686542189hist(tt)
basicStats(tt)## tt
## nobs 1000.000000
## NAs 0.000000
## Minimum 0.002700
## Maximum 44.658843
## 1. Quartile 1.613190
## 3. Quartile 8.417071
## Mean 5.787611
## Median 4.097900
## Sum 5787.611185
## SE Mean 0.182211
## LCL Mean 5.430050
## UCL Mean 6.145172
## Variance 33.200991
## Stdev 5.762030
## Skewness 1.894417
## Kurtosis 5.088628Tuổi thọ trung bình bếp điện nấu cuả quán là khoảng 6 năm với độ lệch chuẩn là 6,4 năm . Trong đó tuổi thọ ngắn nhất của bếp điện là khoảng 0 năm và dài nhất khoảng 43 năm.