Tổng hợp Xác suất thống kê - Mai Huy Sơn Bài 1: Gọi mu1 mức tiêu thụ xăng trung bình động cũ Gọi mu2 mức tiêu thụ xăng trung bình động H0: mu1 mu2 (dự đốn có sở) > t.test(Cu,Moi,mu=0,alt="g") p-value = 0.9467 > 0.05 => chấp nhận H0, bác bỏ H1 Vậy mức ý nghĩa 5%, dự đoán nhà sản xuất khơng có sở Bài 3: a) > t.test(Thang1,conf.level=0.96) [77.73251 , 139.76749] b) Gọi p tỷ lệ người dùng điện thoại lớn 300 tháng H0: p >= 10% H1: p < 10% > sum(Thang1>300) [1] > length(Thang1) [1] 20 > prop.test(0,20,p=0.1,correct=T,alt="l") p-value = 0.1318 > 0.05 => chấp nhận H0, bác bỏ H1 Vậy mức ý nghĩa 5%, tỷ lệ người dùng điện thoại lớn 300 tháng chiếm 10% c) Gọi mu tiền điện thoại trung bình tháng H0: mu = 200 H1: mu # 200 > t.test(Thang1, mu=200, alt="t") p-value = 3.249e-06 < 0.05 => bác bỏ H0, chấp nhận H1 Vậy mức ý nghĩa 5%, tiền điện thoại trung bình tháng 200 d) Gọi mu1 tiền điện thoại trung bình tháng nhóm ly Gọi mu2 tiền điện thoại trung bình tháng nhóm có gia đình H0: mu1 mu2 > n1=subset(Thang1, TinhTrangHonNhan=="Da ly hon") > n2=subset(Thang2, TinhTrangHonNhan=="Dang co gia dinh") > t.test(n1, n2, mu=0, alt="g") p-value = 0.5617 > 0.05 => chấp nhận H0, bác bỏ H1 Tổng hợp Xác suất thống kê - Mai Huy Sơn Vậy mức ý nghĩa 10% trung bình người dùng tiền điện thoại tháng nhóm ly khơng cao nhóm có gia đình e) Gọi mu1, mu2, mu3 tiền sử dụng điện thoại tháng trung bình nhóm chưa lập gia đình, ly hơn, có gia đình H0: mu1 = mu2 = mu3 H1: tồn mui # muj với i,j = 1:3 > anova(lm(Thang1~TinhTrangHonNhan)) p-value = 0.538 > 0.05 => chấp nhận H0, bác bỏ H1 Vậy mức ý nghĩa 5%, tiền điện thoại trung bình loại tình trạng nhân Nguồn biến Tổng bình phương thiên Giữa nhóm 5291 Nội nhóm 69923 Tổng cộng 75214 Bậc tự Phương sai Tỷ số F 17 19 2645.3 4113.1 0.6431 Fth=qf(1-0.05,2,17) [1] 3.591531 F-value < Fth => chấp nhận H0, bác bỏ H1 Bài 4: a) > table(GioiTinh,KhuVuc) KhuVuc GioiTinh HaiDao MienNui NongThon ThanhPho Nam Nu 6 13 13 10 10 26 16 b) H0: khu vực sống giới tính khơng có mối liên hệ với H1: khu vực sống giới tính có mối liên hệ với > a=table(GioiTinh,KhuVuc) > chisq.test(a) p-value = 0.7067 > 0.05 => chấp nhận H0, bác bỏ H1 Vậy mức ý nghĩa 5%, khu vực sống giới tính khơng có mối liên hệ với Bài 5: Gọi mu trọng lượng trung bình gói đường kính xuất H0: mu = H1: mu < Tổng hợp Xác suất thống kê - Mai Huy Sơn Cách 1: > x=rep(c(0.9,0.95,0.98,1,1.01,1.02),c(5,10,15,10,5,5)) > t.test(x,mu=1,alt="l") p-value = 8.578e-06 < 0.05 => bác bỏ H0, chấp nhận H1 Vậy mức ý nghĩa 5%, trọng lượng trung bình gói đường xuất nhỏ kg Cách 2: > tth=qt(1-0.05,49) > tth [1] 1.676551 tkd = -4.7663 tkd < -tth => bác bỏ H0, chấp nhận H1 Bài 6: H0: loại xe ô tô ưa chuộng H1: loại xe ô tô không ưa chuộng > x=c(32,35,45,40) > p0=c(0.25,0.25,0.25,0.25) > chisq.test(x,p=p0) p-value = 0.4612 > 0.05 => chấp nhận H0, bác bỏ H1 Vậy mức ý nghĩa 5%, loại xe ô tô ưa chuộng Bài 7: Tìm ước lượng điểm cho tỷ lệ cho số IQ>120 > sum(ChiSoIQ>120) [1] 18 > length(ChiSoIQ) [1] 200 > 18/200 [1] 0.09 Tìm ước lượng điểm cho trung bình IQ > mean(ChiSoIQ) [1] 99.97 Tìm khoảng tin cậy 95% cho số IQ trung bình Tổng hợp Xác suất thống kê - Mai Huy Sơn > t.test(ChiSoIQ,conf.level=0.95) [ 97.71837 , 102.22163] Tìm khoảng tin cậy 96% cho tỷ lệ IQ>130 > sum(ChiSoIQ>130) [1] > length(ChiSoIQ) [1] 200 > prop.test(4,200,correct=T,conf.level=0.96) [0.006148167 , 0.055846647] Tìm ước lượng điểm cho phương sai số IQ (tính phương sai bình thường var) > var(ChiSoIQ) [1] 260.7529 Bài 8: Gọi md1, md2, md3 trung vị trọng lượng bao gạo đo máy H0: md1 = md2 = md3 H1: tồn mdi # mdj với i,j = 1:3 > m1=c(20.2,20.3,19.8,19.7,20.3,20.4) > m2=c(20.5,20.1,20.7,19.9,19.2,19.9) > m3=c(19.5,19.7,19.8,19.4,19.7,19.6) > kruskal.test(list(m1,m2,m3)) p-value = 0.03477 < 0.05 => bác bỏ H0, chấp nhận H1 Vậy mức ý nghĩa 5%, trọng lượng bao gạo đo máy khác Bài 9: a) > table(VungGui) VungGui KhuVuc1 KhuVuc2 KhuVuc3 NoiTinh 75 67 62 96 b) H0: tỷ lệ bưu phẩm gửi vùng không thay đổi H1: tỷ lệ bưu phẩm gửi vùng thay đổi Tổng hợp Xác suất thống kê - Mai Huy Sơn > x=c(75,67,62,96) > p0=c(0.2,0.25,0.3,0.25) > chisq.test(x,p=p0) p-value = 0.0002492 => bác bỏ H0, chấp nhận H1 Vậy mức ý nghĩa 5%, tỷ lệ bưu phẩm gửi vùng thay đổi Bài 10: Gọi mu1 huyết áp trung bình bệnh nhân trước dùng Timolol Gọi mu2 huyết áp trung bình bệnh nhân sau dùng Timolol H0: mu1 - mu2 30 > t.test(truoc,sau,mu=30,alt="g",paired=T,var.equal=T) p-value = 0.9943 => chấp nhận H0, bác bỏ H1 Vậy mức ý nghĩa 5%, kết luận hàm lượng timolol trước điều trị lớn sau điều trị 30 Bài 11: Gọi p1 tỷ lệ nam yêu thích thể thao mức Gọi p2 tỷ lệ nữ yêu thích thể thao mức H0: p1 p2 > sum(GioiTinh=="Nam" & TheThao=="3") [1] 31 > sum(GioiTinh=="Nu" & TheThao=="4") [1] 37 > table(GioiTinh) GioiTinh Nam Nu 227 248 > x=c(31,37) > n=c(227,248) > prop.test(x,n,correct=F,alt="g") p-value = 0.6527 => chấp nhận H0, bác bỏ H1 Tổng hợp Xác suất thống kê - Mai Huy Sơn Vậy mức ý nghĩa 5%, tỷ lệ nam thích thể thao mức khơng cao tỉ lệ nữ thích thể thao mức Bài 12: Gọi p1 tỷ lệ sản phẩm bj lỗi từ dây chuyền thứ Gọi p2 tỷ lệ sản phẩm bị lỗi từ dây chuyền thứ hai H0: p1 = p2 H1: p1 # p2 > x=c(50,60) > n=c(600,650) > prop.test(x,n,correct=F,alt="t") p-value = 0.5758 => chấp nhận H0, bác bỏ H1 Vậy mức ý nghĩa 5%, chất lượng dây chuyền Bài 12’: H0: phân phối số lần bán hàng thực tuân theo pp nhị thức H1: phân phối số lần bán hàng không tuân theo pp nhị thức > x=c(10,41,60,20,6,3) > p0=dbinom(0:5,5,0.4) > chisq.test(x,p=p0) p-value = 0.0356 => bác bỏ H0, chấp nhận H1 Vậy mức ý nghĩa 5%, phân phối số lần bán hàng không tuân theo pp nhị thức Bài 13: H0: chiều cao người dân việt nam có tuân theo phân phối chuẩn H1: chiều cao người dân việt nam không tuân theo phân phối chuẩn Bài 15: Gọi md1 trung vị hàm lượng cholesterol đo trước tham gia chương trình Gọi md2 trung vị hàm lượng cholesterol đo sau tham gia chương trình Tổng hợp Xác suất thống kê - Mai Huy Sơn H0: md1 >= md2 H1: md1 < md2 > truoc=c(265,240,258,295,251,245,287,314,260,279) > sau=c(229,231,227,240,238,241,234,256,247,239) > wilcox.test(truoc,sau,mu=0,alt="l",paired=T) p-value = 0.9979 => chấp nhận H0, bác bỏ H1 Vậy mức ý nghĩa 5%, không đủ sở cho tham gia chương trình hàm lượng cholretrol máu giảm Bài 16: a) Kiểmđịnh xem mức độ đồng số IQ nam có cao mức độ đồng số IQ nữ hay không mức ý nghĩa 5% Gọi S12 phương sai số IQ nam Gọi S22 phương sai số IQ nữ H0: S12 >= S22 H1: S12 < S22 > var.test(nam,nu,mu=0,alt="l") p-value = 0.6484 => chấp nhận H0, bác bỏ H1 Vậy mức ý nghĩa 5%, mức độ đồng số IQ nam không cao mức độ đồng số IQ nữ b) Xem độ biến động Trí Tuệ IQ nam có thấp độ biến động Trí Tuệ IQ nữ hay không? Gọi S12 phương sai số IQ nam Gọi S22 phương sai số IQ nữ H0: S12 >= S22 H1: S12 < S22 ... Vậy mức ý nghĩa 5%, không đủ sở cho tham gia chương trình hàm lượng cholretrol máu giảm Bài 16: a) Kiểm định xem mức độ đồng số IQ nam có cao mức độ đồng số IQ nữ hay không mức ý nghĩa 5% Gọi... không tuân theo pp nhị thức Bài 13: H0: chiều cao người dân việt nam có tuân theo phân phối chuẩn H1: chiều cao người dân việt nam không tuân theo phân phối chuẩn Bài 15: Gọi md1 trung vị hàm... => chấp nhận H0, bác bỏ H1 Vậy mức ý nghĩa 5%, khu vực sống giới tính khơng có mối liên hệ với Bài 5: Gọi mu trọng lượng trung bình gói đường kính xuất H0: mu = H1: mu < Tổng hợp Xác suất thống