(Tiểu luận) báo cáo bài tập lớn môn xác xuất thống kê chủ đề 3

41 7 0
(Tiểu luận) báo cáo bài tập lớn môn xác xuất thống kê chủ đề 3

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

Thông tin tài liệu

  ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC BÁCH KHOA *** BÁO CÁO BÀI TẬP LỚN MÔN: XÁC XUẤT THỐNG KÊ CHỦ ĐỀ Giảng viê n hướ ng dản: Hoả ng Vản Hả Nho m sinh viê n thư c hiê  n: Nho m 15 – L10 Họ tên MSSV Nhiệm vụ Trả n Nguyê n Xuả n Tiê n 2012189 Soả n bả o cả o , ly thuyêt   Trướng Bả o Khả nh 2011396 Soả n bả o cả o, ly thuyêt   Nguyê n Cướ ng Thinh 2014596 Soả n Codê Rstudio hoảt đo  ng Phu ng Minh Đư c 2013012 Soản Codê Rstudio hoảt đo  ng Hồ Chí Minh, ngày 08 tháng năm 2022 h   LỜI CẢM ƠN Trong thớ i giản ho c tả p bo  mo n Xảc Xuảt Thong Kê ớ lớ p, chu ng êm đả co ho  i tiêp xu c vả lảm quên với nhiê u kiên thư c, lả cớ sớ đê chu ng êm co thê hoả n thảnh bả i tả  p lớ n nảy Đả y cu ng lả nhưng kiên thư c quy bảu phu c vu cho quả trình hoc tả p, lảm viê  c sảu nả y cu ả chu ng êm Ngoải rả, chu ng êm cảm thảy bản thản co sư tiên bo  viê  c chu đo  ng ho c tả p, tìm kiêm thong tin, trảu doi kì nảng lảm viê  c nhom, tảo moi quản hê  gản kêt với cả c bả n nho m lớp Đê co đước kêt quả nảy lả nhớ sư tả n tả m quả trình giảng dảy, truyên đảt  kiên thư c ớ lớ p vả hướng dả n chu ng êm quả trình thư c hiê  n bả i tả p lớ n cu ả Thả y Hoảng Vả n Hả Chu ng êm xin gư i lớ i cả m ớn sảu sả c vả chản thảnh đên thảy h   Mục lục Hoạt động 1: Cơ sở lý thuyết .3 3.1 3.1.1 Một số khái niệm 3.1.2 Mô tả liệu đồ thị 3.1.3 Mô tả liệu định lượng .5 3.1.4 Các phân phối thường gặp thống kê 3.2 Thống kê mô tả Phân ch phương sai – Kiểm định ANOVA (Analysic of Variance) .12 3.2.1 Phân ch phương sai yếu tố (One way ANOVA) 12 3.2.2 Phân ch phương sai hai yếu tố (Two way ANOVA) 16 THỰC HIỆN: 20 4.1 Đọc liệu, thực thống kê mô tả thực kiểm định: 20 4.2 Phân ch phương sai nhân tố (one way anova): .22 4.3 Phân ch phương sai nhân tố (two way anova): .27 Hoạt động 2: 30 5.1 Đọc liệu, thực thống kê mô tả thực kiểm định: 31 5.2 Phân ch phương sai nhân tố (one way anova): .33 5.3 Phân ch phương sai nhân tố (two way anova): .35 TÀI LIỆU THAM KHẢO 37 h   Hoạt động 1: Tả  p tin Diêt.csv (cung cảp bớ i Đả i ho c Shêfiêld, Anh) chư ả tho ng tin vê mo  t thư nghiê  m vê hiê  u quả cu ả cả c chê đo  ản kiê ng viê  c giả m cả n nả ng đoi vớ i nhưng ngưới trướng thả nh Mo  t ngướ i thảm giả sê đước ả p du ng mo  t bả chê đo  ản kiê ng khả c nhảu vong tuả n lê Cả n nả ng cu ả ngướ i thảm giả sê đước ghi nhả  n trước vả sảu kêt thu c thư nghiê  m đê đả nh giả hiê  u quả cu ả tưng chê đo  ả n kiê ng Cảc biên chình bo  dư liê  u: • Pêrson = so thư tư cu ả ngướ i thảm giả thư nghiê  m • gêndêr = giớ i tình cu ả ngướ i thảm giả (1 = nảm, = nư ) • Agê = tuoi (nảm) • Hêight = chiê u cảo (cm) • prê.wêight = cản nả ng trước ảp du ng chê đo  ả n kiê ng (kg) • Diêt = chê đo  ản kiê ng (3 chê đo  khảc nhảu) • wêight6wêêks = cản nả ng sảu tuả n ản kiê ng Cảc bướ c thư c hiê  n: Đo c dư liê  u (Import dảtả): Diêt.csv Lả m sảch dư liê  u (Dảtả clêảning): NA (dư liê  u khuyêt) Lả m ro dư liê  u: (Dảtả visuảlizảtion) (ả) Chuyên đoi biên (nêu cản thiêt) (b) Thong kê mo tả: dung thong kê mảu vả dung đo thi t.têst: Du ng mo  t kiêm đinh phu hớ p cho hải biên prê.uêight vả uêight6wêêks ANOVA mo  t nhản to: Chê đo  ản kiê ng Diêt nảo hiê  u quả nhảt viê  c giả m cản ANOVA hải nhản to: Chê đo  ản kiê ng Diêt vả giới tình gêndêr ả nh hướ ng thê nảo đên viê  c giả m cản wêightLOST? h   Cơ sở lý thuyết  3.1 Thống kê mô tả 3.1.1 Một số khái niệm -Tong thê (populảtion): tả p hớp tảt cả nhưng phả n tư mảng đả c trưng quản tả m hảy cản nghiê n cư u -Mảu (sảmplê): lả mo  t tả p đước chon rả tư tong thê -Thảm so (pảrảmêtêr): lả mo  t đả  c trưng cu thê cu ả mo  t tong thê -Thong kê (stảtistic): lả mo  t đả c trưng cu thê cu ả mo  t mả u -Chon mo  t mảu ngả u nhiê n: Mo  t mảu ngả u nhiê n (rảndom sảmplê) go m n phả n tư đướ c cho n rả tư mo  t tong thê phải tho ả cảc điê u kiê  n sảu: Moi phả n tư ton g thê phải đướ c cho n ngảu nhiên vả đo  c lả  p Moi phản tư tong thê co khả nả ng đướ c cho n nhảu (xảc suảt  đước chon bảng nhảu) Mo i mảu cớ n cu ng co cu ng khả nảng đướ c cho n tư tong thê -Phướng phả p cho n mảu ngảu nhiê n đớn giản (simplê rảndom sảmpling): Đả nh so cảc phản tư cu ả tong thê tư đên N Lả  p cảc phiêu cu ng đả nh so vả y Tro n đê u cả c phiêu, sảu đo cho n co hoản lả i n phiêu Cảc phả n tư cu ả tong thê co so thư tư phiêu lảy rả sê đướ c cho n lả m mảu      -Thong kê mo tả (dêscriptivê stảtistics): lả quả trình thu thả p, tong hớ p vả xư ly dư liê  u đê biên đoi dư liê  u thảnh tho ng tin    Thu thả  p dư liê  u: khả o sảt, đo đảc, Biêu diê n dư liê  u: du ng bảng vả đo thi Tong hớp dư liê  u: tình cảc thảm so mảu trung bình mảu (sảmplê mêản), phướng sải mảu (sảmplê vảriảncê), trung vi (mêdiản), 3.1.2 Mô tả liệu đồ thị  Giới thiệu: -Viê  c mo tả dư liê  u bả ng đo thi sê cho tả mo  t cải nhìn tong quản vê dư liê  u trướ c vảo phản tìch cu thê -Cả c loả i đo thi đướ c sư du ng sê phu thuo  c vảo dảng biên cản phản tìch h   -Trong phản nảy, tả sê khảo sảt chu yêu vê đo thi to chưc tản so (histogrảm) -Phản bo tản so (frêquêncy distribution): lả mo  t bảng chưả cảc khoảng đước phản nho m thêo dư liê  u quản trả c vả cả c tả n so tướng ư ng cu ả dư liê  u nảm bê n tưng khoả ng Phản bo tản so cho phêp mo tả phả n phoi cu ả dư liê  u, xêm xê t tình đoi xưng/bảt đoi xư ng, tả p trung/phản tản cu ả dư liê  u nhả n dả ng phả n phoi chuản (bêll-shảpêd), xảc đinh modê (unimodảl, bimodảl)       -Phản bo tản so đướ c mo tả bảng đo thi to chư c tản so (histogrảm) Xây dựng phân bố tần số: -Trong mo  t bảng phản bo tản so:   Moi nhom co bê ro  ng bảng nhảu Bê ro  ng cu ả moi nho m đước xả c đinh bớ i: Giátrịlớnnhất −Giá trịbé nhất  Số khoảngcầnchia Cả c khoảng kho ng tru ng nhảu Nên chon so khoảng toi thiêu ≥ -Đo thi biêu diên bảng phản bo tản so goi lả đo thi  to chưc tản so (histogrảm)   - Chiả dư liê  u thảnh bảo nhiê u khoảng lả tot? lả quả trình "thư" vả "sải", đo thi to chư c tản so kho ng đướ c quả "lớ m chớ m", co nhiê u đình vả khong co dảng "khoi", chì rả đướ c sư biên thiên dư liê  u -Dả ng điê  u cu ả phả n phoi (shảpê of thê distribution) go i lả đoi xư ng (symmêtric) nêu cảc giả tri quản trảc cả n bả ng xung quảnh trung tả m    h   -Dảng điê  u cu ả phản phoi go i lả bảt đoi xưng (skêwêd) nêu dư liê  u quản trả c kho ng phản bo đoi xư ng xung quảnh trung tảm -Sư du ng đo thi histogrảm đê nhả n biêt phản phoi xảc suảt cu ả mo  t đả i lướng ngảu nhiên 3.1.3 Mô tả liệu định lượng Các độ đo trung tâm: -Trung bình (mêản) lả đả i lướ ng thướ ng đướ c sư du ng nhảt đê đo giả tri trung tảm cu ả dư liê  u Vớ i mo  t tong thê co N phả n tư, trung bình tong thê tì nh bớ i:   N   x ∑ = i  μ=   N   x + x 2+ … + x N  i  =  N  Vớ i mo  t mảu cớ n, trung bình mảu n  x =   x ∑ =  x n i  x + x 2+ … + x n  = n Trung bình bi ả nh hướ ng bớ i cảc giả tri ngoải lải (outliêrs) -Trong mo  t tả p dư liê  u đước sả p xêp thêo thư tư tả ng dản, trung vi (mêdiản) lả giả tri "chình giư ả" cu ả dư liê  u (50% bê n trê n, 50% bê n dưới)   Trung vi kho ng bi ảnh hướ ng bớ i cảc điêm ngoả i lải (outliêrs) Vi trì cuả trung vi: sảp xêp dư liê  u thêo thư tư tả ng dản, go i i lả vi trì cu ả trung vi i= n +1 + Nêu i chản, trung vi = Xi + Nêu i lê , trung vi¿  X [ i ] + X [ i ]+ , vớ i [i] lả phản nguyê n cu ả i h   -Modê (yêu vi) lả mo  t đả i lướ ng đê đo xu hướ ng trung tảm cu ả dư liê  u, lả giả tri thướng xảy rả nhảt, kho ng bi ả nh hướ ng bớ i cả c điêm ngoả i lải, co thê sư du ng cho cả dư liê  u đinh tình vả dư liê  u đinh lướ ng, co thê co nhiê u modê hoả c khong to n tả i modê Các độ đo biến thiên: -Khoảng biên thiên (rảngê) lả đo  đo sư biên thiê n đớn giản nhảt, lả đo  chê nh lê  ch giưả giả tri lớ n nhảt vả bê nhảt cu ả dư liê  u quản trả c  Khoảng biến thiên= X max − X min   Bo phả n bo cu ả dư liê  u Bi ả nh hướ ng bớ i cả c điêm ngoải lải -Tả co thê loải bo cả c điêm ngoải bả ng cảch sư du ng khoả ng tư phản vi (IntêrQuảrtilê Rảngê-IQR) (hảy co n đướ c go i lả đo  trải giưả )  Co ng thư c tì nh khoả ng tư phả n vi: IQR = Q3 − Q1 vớ i Q1 lả phản vi thư (mưc 25%) vả Q3 lả phản vi thư (mư c 75%) cu ả dư liê  u   Cảc điêm Q1, Q2, vả Q3 đướ c go i lả cảc điêm tư phản vi Co ng thư c tì m phản vi: Sảp xêp dư liê u thêo thư tư tả ng dả n, go i Q 1, Q2 (trung vi), Q3 lản lướ t lả phả n vi thư 1, vả cu ả dư liê  u Vi trì cuả Q1, Q2 vả Q3 đướ c xảc đinh sảu: Vi trì phản vi thư nhảt = 0.25(n + 1) Vi trì phản vi thư hải = 0.5(n + 1) Vi trì phản vi thư bả = 0.75(n + 1) vớ i n lả so giả tri quản trảc -Đo thi ho  p:  Đê biêu diê n khoảng tư phản vi  vả cảc điêm ngoải lải (outliêrs) tả dung đo thi ho  p (boxplot) h    Khi vê nhiê u đo thi boxplot cu ả nhiê u tả  p dư liê  u khả c nhảu bê n cảnh nhảu, tả co n co thê so sả nh đướ c đo  phả n tản vả so sả nh giả tri trung tảm (trung bình/trung vi) cu ả cả c tả p dư liê  u nả y -Phướng sải (Vảriảncê) lả trung bình cuả bình phướng đo  lê  ch cả c giả tri so với trung bình   Phướng sải phản ảnh đo  phả n tản hảy sư biên thiê n cu ả dư liê  u Phướng sải ton g thê  N  ( x − μ ) ∑ = i σ  = i  N  vớ i N lả so phả n tư cu ả tong thê, µ lả trung bì nh tong thê, x i lả giả tri thư i cu ả biên x  Phướng sải mả u n ( X  − X ) ∑ = i S= i n −1 vớ i  X  lả trung bình mả u, n lả cớ mả u, Xi lả giả tri thư i cu ả biên X h   -Đo  lê  ch tiê u chuản (Stảndảrd dêviảtion) đước du ng đê đo sư biên thiê n, biêu diê n sư biên thiê n xung quảnh trung bình,  Co cu ng đớn vi đo vớ i dư liê  u goc Đo  lê  ch chuản cu ả tong thê, ky hiê  u lả σ:  Đo  lê  ch chuản mảu:  -Hê  so biên thiê n (Coêficiênt of Vảriảtion) đước sư du ng đê so sả nh sư biên thiê n cuả hải hảy nhiêu tả  p dư liê  u, co thê đo ớ cả c đớn vi khảc nhảu    Đo moi liên hê  giưả sư biên thiê n vả trung bình Đớn vi tình bảng % Co ng thưc:  S CV =  100%  X  3.1.4 Các phân phối thường gặp thống kê Phân phối chuẩn tắc: -Biên ngảu nhiên liên tuc X nhả  n giả tri khoảng (−∞, +∞) đướ c go i lả co phản phoi chuản thảm so µ, σ nêu hảm mả  t đo  xả c suảt co dảng -Biên ngả u nhiên X đướ c go i lả co phả n phoi chuản tảc nêu no co phả n phoichuản vớ i thảm so µ = vả σ 2 = 1, ky hiê  u X ∼   (0, 1) -Thêo quy ướ c, hả m phản nhiê n chuản ho ả đước ky phoi cu ả biên ngả u hiê  u lả Φ(x), tưc lả h   Nhận xét: Biêu đo QQ-plot  cho tả thảy như ng giả tri quản sảt đả phả n nảm trê n đướ ng thả ng kì vo ng cu ả phả n phoi chuản đo biên weightLOST ớ chê đo  ản kiê ng tuản thêo phản phoi chuản Ngoải rả, tả co thê dung hảm shapiro.test  đê kiêm trả: shảpiro.têst(diêt2$wêightLOST) Nhận xét: Giả thiêt H_0: Cản nả ng giả m ớ chêt đo  ản kiê ng tuản thêo phản phoi chuản Giả thiêt H_1: Cản nả ng giả m ớ chê đo  ản kiê ng kho ng tuản thêo phản phoi chuản Vì Pr(>F) = 0,8722 > mưc y nghìả α=0,05α=0,05 nên tả chảp nhả n giả thiêt H_0 Vả  y biên weightLOST ớ chê đo  ản kiê ng tuả n thêo phản phoi chuản Đoi vớ i chê đo  ản kiê ng 3: diêt3 = subsêt(Diêt1, Diêt == 3) qqnorm(diêt3$wêightLOST) qqlinê(diêt3$wêightLOST) 25 h   Nhận xét: Biêu đo QQ-plot  cho tả thảy như ng giả tri quản sảt đả phả n nảm trê n đướ ng thả ng kì vo ng cu ả phả n phoi chuản đo biên weightLOST ớ chê đo  ản kiê ng tuản thêo phản phoi chuản Ngoải rả, tả co thê dung hảm shapiro.test  đê kiêm trả: shảpiro.têst(diêt3$wêightLOST) Nhận xét: Giả thiêt H_0: Cản nả ng giả m ớ chê đo  ản kiê ng tuản thêo phản phoi chuản Giả thiêt H_1: Cản nả ng giả m chê đo  ản kiê ng khong tuản thêo phản phoi chuản Vì Pr(>F) = 0,372 > mư c y nghìả α=0,05α=0,05 nên tả chảp nhả  n giả thiêt H_0 Vả  y biên weightLOST ớ chê đo  ản kiê ng tuả n thêo phản phoi chuản Tả du ng hảm leveneTest  đê kiêm trả giả đinh tình đong nhảt cuả cảc phướng sải: library(cảr) #Yêu cầu gói lệnh 'car' để sử dụng hàm 'leveneTest'  lêvênêTêst(wêightLOST ~ ảs.fảctor(Diêt),dảtả = Diêt1) #Thực kiểm định tính đồng phương sai Nhận xét: Giả thiêt H_0: Phướng sải cản nả ng giả m ớ cảc phướng phả p ả n kiê ng bả ng nhảu Giả thiêt H_1: Co ì t nhảt phướng phảp ản kiêng co phướng sải cản nả ng giả m khảc nhảu Vì Pr(>F) = 0.6313 > mưc y nghìả α=0,05α=0,05 nên tả chảp nhả n giả thiêt H_0 Vả  y phướng sải cả n nả ng giảm ớ cả c phướng phảp ản kiê ng lả nhảu 26 h   ả Du ng lê  nh aov đê thưc hiê  n phả n tìch Anovả mo  t nhản to: ảnovả1 = ảov(wêightLOST ~ ảs.fảctor(Diêt),dảtả = Diêt1) #Thực phân tích  phương sai summảry(ảnovả1) #Tóm tắt kết mơ hình phân tích phương sai Nhận xét: Vì Pr(>F) < mư c y nghìả α=0,05α=0,05 nên bảc bo giả thiêt H_0, chảp nhả  n giả thiêt H_1 Vả y co ìt nhảt phướng phảp ản kiêng co cản nả  ng trung bình giảm khảc nhảu, tưc mư c đo  hiê  u quả đoi vớ i viê  c giảm cả n giưả cảc phướng phả p ả n kiê ng lả khảc nhảu ả Du ng hả m TukeyHSD  vả plot  đê thưc hiê  n cả c so sả nh bo  i (multiplê compảrisons) sảu phản tìch phướng sải: TukêyHSD(ảnovả1) #Thực so sánh bội plot(TukêyHSD(ảnovả1)) #Vẽ đồ thị so sánh bội Nhận xét: Sự khác cân nặng giảm trung bình nhóm chế độ ăn kiêng 1: Giả thiêt H_0: Trung bình cản nả ng giảm ớ nho m chê đo  ản kiê ng vả chê đo  ả n kiê ng bảng nhảu Giả thiêt H_1: Trung bình cản nả ng giảm ớ nho m chê đo  ản kiê ng vả chê đo  ả n kiê ng khả c nhảu Tả nhả n thảy p ảdj = 0,9987711 > mư c y nghìả α=0,05α=0,05, nên tả chảp nhả  n giả thiêt H_0 Vả  y trung bình cản nả  ng giảm ớ nho m chê đo  ản kiê ng vả chê đo  ả n kiê ng bảng nhảu 27 h   Sự khác cân nặng giảm trung bình nhóm chế độ ăn kiêng 1: Giả thiêt H_0: Trung bình cản nả ng giảm ớ nho m chê đo  ản kiê ng vả chê đo  ả n kiê ng bảng nhảu Giả thiêt H_1: Trung bình cản nả ng giảm ớ nho m chê đo  ản kiê ng vả chê đo  ả n kiê ng khả c nhảu Tả nhả n thảy p ảdj = 0,0188047 < mư c y nghìả α=0,05α=0,05, nê n tả bảc bo giả thiêt H_0, chảp nhả n giả thiêt H_1 Vả  y trung bình cản nả ng giả m ớ nho m chê đo  ả n kiê ng vả chê đo  ả n kiê ng khảc nhảu Mả  t khảc tả dưả vả o giả tri diff = 1,848148 > nên tả co thê kêt luả n trung bình cản nả ng giả m ớ nho m chê đo  ản kiê ng cảo hớn so vớ i nho m chê đo  ản kiê ng Sự khác cân nặng giảm trung bình nhóm chế độ ăn kiêng 2: Giả thiêt H_0: Trung bình cản nả ng giảm ớ nho m chê đo  ản kiê ng vả chê đo  ả n kiê ng bảng nhảu Giả thiêt H_1: Trung bình cản nả ng giảm ớ nho m chê đo  ản kiê ng vả chê đo  ả n kiê ng khả c nhảu Tả nhả n thảy p ảdj = 0.0152020 < mư c y nghìả α=0,05α=0,05, nê n tả bảc bo giả thiêt H_0, chảp nhả n giả thiêt H_1 Vả  y trung bình cản nả ng giả m ớ nho m chê đo  ả n kiê ng vả chê đo  ả n kiê ng khảc nhảu Mả  t khảc tả dưả vả o giả tri diff = 1.880148 > nên tả co thê kêt luả n trung bình cản nả ng giả m ớ nho m chê đo  ản kiê ng cảo hớn so vớ i nho m chê đo  ản kiê ng Kết luận: phướng phảp co cản nả ng giảm trung bình cảo hớn phướng phảp vả Vả  y phướng phảp lả phướng phả p ả n kiê ng co hiê  u quả nhảt viê  c giảm cản 4.3 Phân tích phương sai nhân tố (two way anova): ả Thưc hiê  n phản tìch phướng sải hải nhả n to đê đảnh giả sư ảnh hướ ng chế độ ăn kiêng vả giới tính đên sư giả m cản: ảnovả2 = ảov(wêightLOST ~ ảs.fảctor(gêndêr)*ảs.fảctor(Diêt),dảtả = Diêt1) #Thực phân tích phương sai summảry(ảnovả2) #Tóm tắt kết mơ hình phân tích phương sai Nhận xét:: Giả thiêt H_0: Cản nả ng trung bình giảm giư ả cả c giớ i tình bả ng nhảu Giả thiêt H_1: Cản nả ng trung bình giảm giu ả cả c giớ i tình lả khảc nhảu 28 h   Vì Pr(>F)= 0,82062 > mư c y nghìả α=0,05α=0,05 nên tả chảp nhả n giả thiêt H_0 Vả  y giới tính kho ng ảnh hướ ng đên viê  c giảm cả n Giả thiêt H_0: Cản nả ng trung bình giảm giư ả cả c phướng phả p ả n kiê ng Giả thiêt H_1: Co ì t nhảt phướng phảp ản kiêng co cản nả ng trung bình giảm khảc nhảu Vì Pr(>F) = 0,00546 < mưc y nghìả α=0,05α=0,05 nên bảc bo giả thiêt H_0, chảp nhả n giả thiêt H_1 Vả  y chế độ ăn kiêng co ảnh hướ ng đên viê  c giảm cả n ả Phả n tìch sư tướng tả c giưả chế độ ăn kiêng  vả giới tính đên sư giảm cả n: Giả thiêt H_0: Kho ng co sư tướng tảc giưả giớ i tì nh vả phướng phả p ả n kiê ng Giả thiêt H_1: Co sư tướng tảc giưả giớ i tì nh vả phướng phả p ả n kiê ng Tư kêt quả cu ả phản tìch phướng sải, tả co Pr(>F) = 0,04884 < mư c y nghìả α=0,05α=0,05 nên bảc bo giả thiêt H_0, chảp nhả n giả thiêt H_1 Vả y co sư tướng tả c giưả giớ i tình vả chê đo  ản kiê ng Ngoải rả, tả co thê vê đo thi  đê đả nh giả sư tướng tả c: intêrảction.plot(Diêt1$Diêt,Diêt1$gêndêr,Diêt1$wêightLOST,typê="b",col=c(2:3), lêg.bty="o",lêg.bg="bêigê",lwd=2,pch=c(18,24)) Nhận xét:: Tả nhả n thảy đướ ng thảng trê n đo thi cảt nhảu, tư c co sư tướng tảc giưả giớ i tình vả phướng phả p ả n kiê ng 29 h   30 h   Hoạt động 2: Tả  p tin Folds5x2_pp.ods (cung cảp bớ i Hêysêm Kảyả, Khoả Ky thuả t Mả y tình, Đải ho c Bog ziỗi, TR-34342, Bờikt, Istnbul, Tho Nhỡ Ky vả Pìnảr Tu fêkci, Khoả Ky thuả  t Çorlu, Đả i ho c Nảmìk Kêmảl, TR-59860 Çorlu, Têkirdảg , Tho Nhì Ky ) chưả 9568 điêm dư liê  u đướ c thu thả p tư Nhả mảy điê  n chu trì nh hon hớ p vong nảm (2006 - 2011), nhả mả y điê  n nả y đướ c thiêt lả p đê lả m viê  c vớ i đảy đu phu tả i.  Cảc biên so mo i trướng xung quảnh trung bình hảng giớ Giả tri trung bình đướ c lảy tư cả c cảm biên khảc nhảu đả t xung quảnh nhả mảy ghi lả i cả c biên xung quảnh moi giảy Cả c biên đướ c đưả rả mả kho ng cả n chuản ho ả Cảc biên chình bo  dư liê  u: • AT = Nhiê  t đo  mo i trướng xung quảnh (oC) • V = chản kho ng xả (cmHg) • AP = ảp suảt mo i trướ ng xung quảnh (mbảr) • RH = Đo  ảm tướng đoi (%) • PE = sả n lướ ng điê  n nả ng thưc thêo giớ (MW) Cảc bướ c thư c hiê  n: Đo c dư liê  u (Import dảtả): dư liê  u đước chuyên sảng dả ng csv => Folds5x2_pp csv Lả m sảch dư liê  u (Dảtả clêảning): NA (dư liê  u khuyêt) Lả m ro dư liê  u: (Dảtả visuảlizảtion) Thong kê mo tả: du ng thong kê mả u vả du ng đo thi t.têst: Du ng mo  t kiêm đinh phu hớp ANOVA mo  t nhản to: ANOVA hải nhản to: nhiê  t đo  xung quảnh vả chản kho ng xả ả nh hướ ng đên sản lướ ng điê  n 31 h   5.1 Đọc liệu, thực thống kê mô tả thực kiểm định: ả Đo c dư liê  u vả o R: Đo c dư liê  u tư ilê Foldsx2_pp.csv vả hiên thi 6 do ng dư liê  u đảu tiê n Sêtwd(“C:/Usêrs/Admin/Dêsktop/Nêw foldêr/Activê2/CCPP") dảả sêtwd("C:/Usêrs/Admin/Dêsktop/Nêw foldêr/Activê2/CCPP") > dảả hêảd(dảả)   ì AT V AP RH PE 14.96 41.76 1024.07 73.17 463.26 25.18 62.96 1020.04 59.08 444.37 5.11 39.40 1012.16 92.14 488.56 20.86 57.32 1010.24 76.64 446.48 10.82 37.50 1009.23 96.62 473.90 26.27 59.44 1012.23 58.77 443.67 b Lảm sảch dư liê  u: Loả i bo giả tri khuyêt (NA) plảcê_NA=ảpply(is.nả(dảả),2,which) nêw_dảả=nả.omit(dảả) ảnyNA(dảả) c Lảm ro dư liê  u Gả n cả c biên cho dư liê  u: ATF)< 2ê-16 bê hớn mư c y nghìả α=0,05, nên tả bảc bo H_0, chảp nhả  n giả thiêt H_1 Vả  y Nhiệt độ ả nh hướ ng đên sản lướ ng điê n Giả thiêt H_0: Sản lướ ng điê  n kho ng thảy đoi giư ả cả c điê u kiê  n chản khong xả Giả thiêt H_1: Sản lướ ng điê  n co sư thảy đoi giư ả cả c điê u kiê  n chản kho ng xả Vì Pr(>F) = 6.85ê-11 bê hớn mưc y nghìả α=0,05α=0,05 nên bảc bo giả thiêt H_0, chảp nhả n giả thiêt H_1 Vả  y Chân không xả  co ảnh hướ ng đên sản lướ ng điê n 37 h   38 h   TÀI LIỆU THAM KHẢO Slidê bả i giảng thảy Hoả ng Vản Hả Hoảng Tro ng & Chu Nguyê n Mo  ng (2008) Thống kê ứng dụng kinh tế xã hội, Nxb Thong Kê Tong Đình Quy Giáo trình Xác suất thống kê , Nxb Bảch Khoả – Hả No  i Stêphảniê Glên , ANOVA Test: Definition, Types, Examples, SPSS, truy cả  p tư: https://www.stảtisticshowto.com/probảbility-ảnd-stảtistics/hypothêsistêsting/ảnovả/ Bài 4: Phân tích phương sai (ANOVA) , 2020, truy cả  p tư: http://ibf.iuh.êdu.vn/wp-contênt/uploảds/2020/10/BAI-4.-PHAN-TICHPHUONG-SAI.pdf  Tải liê  u hướng dả n sư du ng Rstudio:   Nguyên Vản Tuản, Phân tích số liệu biểu đồ Rstudio, truy cả p tư: https://crản.r-projêct.org/doc/contrib/Intro_to_R_Viêtnảmêsê.pdf  Kê nh Youtubê: Lêảrn to SCIENCE 39 h

Ngày đăng: 04/04/2023, 09:35

Tài liệu cùng người dùng

Tài liệu liên quan