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

Đâ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

XÁC SUẤT

THỐNG KÊ

TRONG Y HỌC

5 Nguyễn Thị Thu Quyên

6 Nguyễn Thị Lâm Oanh

7 Trịnh Thị Nam Phương

8 Đỗ Thị Mai Phương

9 Cao Thị Yến Nhi

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).

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 <- c (18,19,20,21)

> Sanpham <- c (3,5,15,2)

> Dulieu <-data.frame(trongluong,sanpham)

> 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

Hình vẽ: Minh họa mối liên hệ giữa thời gian chế biến và hàm lượng đường

Để tính hệ số tương quan mẫu,

ta dùng lệnh cor(x,y) và thu được kết quả sau:

>cor(thoigianchebien,

hamluongduong)

[1] 0.8236616

Trong khi đó tính bằng tay theo công thức ta được :

Để xác định hệ số a và b của phương trình đường thẳng hồi quy

mô tả tương quan giữa thời gian chế biến(Y) và hàm lượng đường (X) :y=ax+b ta dùng lệnh

>lm(thoigianchebien ~ hamluongduong)

Vậy ta có phương trình hồi quy tuyến tính thu được là :

y = 0,3028x – 5,9044

CREDITS: This presentation template was created by Slidesgo , including icons by Flaticon , and infographics & images by Freepik and illustrations