Xác suất thống kê trong y học giới thiệu phần mềm r

XÁC SUẤT THỐNG KÊ TRONG Y HỌC THÀNH VIÊN NHÓM 6 Nguyễn Thị Lâm Oanh 7 Trịnh Thị Nam Phương 1 Bùi Minh Quân 8 Đỗ Thị Mai Phương 2 Đặng Thu Nhung 9 Cao Thị Yến Nhi 3 Đào Thị Phương 4 Phùng Minh Quân 5 Nguyễn Thị Thu Quyên Bài đọc thêm số 2: Giới thiệu phần mềm R Sau khi nghiên cứu toàn bộ lý thuyết xác suất thống kê, chúng ta cần một công cụ để tính toán nhanh và chính xác các công thức Phần mềm R thỏa mãn được các yêu cầu đó Đây là phần mềm được các nhà nghiên cứu thống kê trên thế giới hợp tác với nhau để cho tất cả thành viên trong ngành thống kê học và toán học trên thế giới có thể sử dụng một cách thống nhất và hoàn toàn miễn phí (Sau đây nhóm em xin trình bày về cách sử dụng phần mềm R) Phần mềm R được ứng dụng cho • Phần ước lượng • Phần kiểm định • Trong tương quan hồi quy 1 Hướng dẫn cài đặt phần mềm Truy cập đường link: https://cran.r-project.org/ Click: Download for Windows => install R for the first time => Download R-4.3.1 for Windows Chọn Save • Chọn Run as administator • Chọn Finish 2 Cách nhập dữ liệu trong R Để phần mềm có thể tính toán, trước hết chúng ta cần nhập dữ liệu Ví dụ: Trọng lượng của một loại sản phẩm A là đại lượng ngẫu nhiên phân phối theo quy luật chuẩn với độ lệch tiêu chuẩn là 1 gam Cân thử 25 sản phẩm loại này ta có kết quả: Hãy nhập dữ liệu trên vào R Chúng ta có thể sử dụng function có tên c như sau > Trongluong Sanpham Dulieu Dulieu 3 Sử dụng phần mềm R cho phần ước lượng Dữ liệu giống như ví dụ trên, với độ tin cậy 1-α=95%, hãy ước lượng trọng lượng trung bình, trọng lượng trung bình tối đa, trọng lượng trung bình tối thiểu cho sản phẩm nói trên Đây là trường hợp ước lượng cho tổng thể có phân phối chuẩn, chưa biết phương sai nên ta dùng hàm t.text Trong cửa sổ R ta thực hiện các lệnh sau: Khai báo các biến trongluong-c(18,19,20,21) sanpham-c(3,5,15,2) > DuLieu=rep(trongluong,sanpham) (thể hiện mối quan hệ giữa 2 đại lượng trọng lượng và sản phẩm ở trên) Ước lượng chi phi trung bình thì chúng ta dùng hàm t.test với công thức như sau: >t.test(DuLieu,alt="" conf level=095) Trong đó ước lượng đối xứng nên ta dùng đối thiết là t và nhập độ tin cậy tương ứng 0,95 Và ấn enter ta thu được kết quả như sau: Ước lượng trọng lượng tối thiểu Ước lượng trọng lượng tối đa 4 Sử dụng phần mềm R cho phần kiểm định Ví dụ: Hàm lượng vitamin B12 trong thuốc tiêm của một cơ sở sản xuất là đại lượng ngẫu nhiên phân phối chuẩn với hàm lượng trung bình theo quy định là 200 g /ml Nghi ngờ thuốc không có đủ hàm lượng vitamin B12 như quy định, người ta kiểm tra ngẫu nhiên 25 ống thuốc thu được kết quả về hàm lượng( tính theo g /ml) Với mức ý nghĩa α=0,01 hãy kết luận về điều nghi ngờ nói trên Trong cửa sổ R ta thực hiện các lệnh sau: >Hamluong=seq( 185,210, by=5) >Tanso=c( 2,4,6,9,3,1) >Dulieu=rep( Hamluong, Tanso) >t.test( Dulieu, mu=200, alt= “t”, conf.level=0.99) Ta được kết quả như sau: Đọc kết quả: Số bậc tự do là 24; Vì p-value= 0.02537>0.01 nên chưa có cơ sở bác bỏ H0, chấp nhận H1, nghĩa là với mức ý nghĩa α=0,01 ta nói rằng ý kiến nghi ngờ thuốc không có đủ hàm lượng vitamin B12 như quy định là không đúng 5 Sử dụng phần mềm R cho trong tương quan hồi quy Ví dụ: Theo dõi sự phụ thuộc giữa mức suy giảm hàm lượng đường trong mía X (%) và thời gian chế biến (Y) ta thu được kết quả sau đây : Nhập số liệu vào R và ta vẽ một biểu đồ tán xạ như sau: >hamluongduongthoigianchebiendataplot(hamluongduong ~ thoigianchebien)