Bài giảng phương pháp thí nghiệm trong chăn nuôi và thú y tập 1 part 9

Chia sẻ: Ahfjh Kasjfhka | Ngày: | Loại File: PDF | Số trang:7

Thêm vào BST

Báo xấu

215
lượt xem 45
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Tin sinh học (bioinformatics) là một lĩnh vực khoa học sử dụng các công nghệ của các ngành toán học ứng dụng, tin học, thống kê, khoa học máy tính, trí tuệ nhân tạo, hóa học và hóa sinh (biochemistry) để giải quyết các vấn đề sinh học. Một thuật ngữ thường được dùng thay thế cho tin sinh học là sinh học tính toán (computational biology). Tuy nhiên, tin sinh học thiên về việc phát triển các giải thuật, lý thuyết và các kĩ thuật thống kê và tính toán để giải quyết các bài toán bắt nguồn từ...

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Bài giảng phương pháp thí nghiệm trong chăn nuôi và thú y tập 1 part 9

3.11.5. Ví d • Ví d v b ng tương liên 2×2. Hawkins và c ng s (1993) ñã ti n hành nghiêm c u nh hư ng c a vi c thi n ñ n s xu t hi n b nh ti u ñư ng chu t. Bi t r ng t l m c m c b nh này chu t là 24% ñ i v i con ñ c và 73 % con cái. T 100 chu t thí nghi m, chia m t cách hoàn toàn ng u nhiên v 1 trong 2 cách s lý thi n và không thi n. S chu t 2 lô thí nghi m ñư c theo dõi cho ñ n 140 ngày tu i và ti n hành l y m u nghiên c u t 42 ngày tu i. B nh ti u ñư ng ñư c xác ñ nh ñ i v i chu t có hàm lư ng ñư ng trong máu l n hơn 200 mg/ dl. K t qu thí nghi m ñư c ghi l i b ng sau: Cách x lý Thi n Không thi n T ng M c b nh 26 12 38 Không m c b nh 24 38 62 T ng s 50 50 100 Câu h i ñ t ra: T l m c b nh 2 cách x lý có khác nhau không? L i gi i H0: π1 = π2, Gi thi t H1: π1 ≠ π2, Nh p t n su t quan sát (O) vào b ng tương liên 2×2 Cách x lý Thi n Không thi n T ng M c b nh 26 12 38 Không m c b nh 24 38 62 T ng s 50 50 100 Ư c tính t n su t theo lý thuy t (E): Cách x lý Vacxin Thi n Không thi n T ng 50 × 38 50 × 38 M c b nh = = 19 = = 19 38 100 100 50 × 62 50 × 62 Không m c b nh = = 31 = = 31 62 100 100 T ng s 50 50 100 Tính giá tr χ2 th c nghi m theo công th c sau ( O − E − 0,5)2 χ =∑ 2 E 57
( 26 − 19 − 0,5) 2 ( 38 − 31 − 0,5) 2 2 χ= +L+ = 7,17 19 31 df = (S hàng - 1) ×(S c t – 1) = (2 - 1) × (2 - 1) = 1 Xác ñ nh giá tr P b ng cách so sánh giá tr χ2 = 7,17 v i df = 1 ta có 0,005 < P < 0,01 Rút ra k t lu n P< 0,05 bác b gi thi t H0; ch ng t r ng t l chu t sau khi thi n m c b nh ñái ñư ng cao hơn so v i chu t không b thi n (π1 = 26/50 = 0,52; π2 = 12/50 = 0,24). • Ví d v b ng tương liên 4×3. ðây là các lo i vacxin phòng b nh ñư c so sánh v i ñ i ch ng. V i các m c ñ : không, trung bình và nhi m b nh n ng ñư c ghi lai sau 24 tháng. S li u ñư c trình bày b ng sau: M c ñ b nh Vacxin Không Trung bình N ng T ng ð i c h ng 100 (137.3) 71 (42.6) 29 (20.1) 200 A 146 (133.9) 32 (41.6) 17 (19.6) 195 B 149 (132.5) 28 (41.2) 16 (19.3) 193 C 146 (137.3) 37 (42.6) 17 (20.1) 200 T ng s 541 168 79 788 B ng cho ta th y t n su t O (E) - có nghĩa là t n su t quan sát và ư c tính. Chúng ta ki m tra H0 r ng không có m i liên quan nào gi a m c ñ b nh và tiêm phòng vacxin, nghĩa là t t c các ñ u có cùng m c ñ nh hư ng. • N u H0 ñúng thì t n su t ư c tính ñư c tính như sau: Nhóm ñ i ch ng, không nhi m b nh (ô ñ u tiên trong b ng): E1 = [(T ng s hàng th nh t) × (T ng s c t th nh t)] /T ng s toàn b quan sát = (200 × 541) /788 = 137,3 Chúng ñư c bi u hi n trong d u (....) c a b ng. Tương t như v y ta có th tính ñư c các giá tr còn l i. • Tính χ2 th c nghi m 58
(O − E ) 2 χ2 = ∑ E = 10.14 + 18.86 + 3.99 + 1.10 + 2.21 + 0.33 + 2.05 + 4.20 + 0.58 + 0.55 + 0.75 + 0.46 = 45.2 • B c t do (df): df = (S hàng - 1) ×(S c t - 1) = (4 - 1) × (3-1) = 6 Xác ñ nh giá tr P P = P(χ 2 > 45.2) = 0.000 6 K t l u n: Có m i liên h có ý nghĩa gi a các lo i vacxin và m c ñ nhi m b nh - có nghĩa là các m c ñ khác nhau c a nghi m th c cho m c ñ mi n d ch khác nhau (P=0,000). Ki m tra l i t n su t O so v i E ta th y nhóm ñ i ch ng có s trư ng h p m c b nh cao hơn (không có vacxin). χ2 45,20 Áp d ng Minitab Chúng ta s nh p s li u vào WorkSheet c a Minitab như sau: Data Display Row Khong TB Nang 1 100 71 29 2 146 32 17 3 149 28 16 4 146 37 17 Lưu ý r ng ta ph i ng m hi u r ng row1 = ñ i ch ng, row2, row3, row4 tương ng v i các vacxin A, B, C. 59
Th c hi n phép th b ng Minitab Chi-Square Test: Khong, TB, Nang Stat > Tables > Chi-Square Test… Expected counts are printed below observed counts Khong TB Nang Total 1 100 71 29 200 137.31 42.64 20.05 2 146 32 17 195 133.88 41.57 19.55 3 149 28 16 193 132.50 41.15 19.35 4 146 37 17 200 137.31 42.64 20.05 Total 541 168 79 788 Chi-Sq = 10.138 + 18.863 + 3.994 + 1.098 + 2.205 + 0.332 + 2.054 + 4.201 + 0.580 + 0.550 + 0.746 + 0.464 = 45.224 DF = 6, P-Value = 0.000 Ta cũng có k t qu tương t như ph n tính tay. ( O − E − 0,5)2 trang 56 χ = ∑ 2 Lưu ý: Công th c ch s d ng ñ i v i nh ng thí E nghi m có s lư ng quan sát h n ch . Các ph n m m th ng kê không th c hi n theo công th c này; vì v y trư c khi x lý s li u nên xem xét và quy t ñ nh có s d ng ñư c ph n m m th ng kê hay không!!! 60
3.12. Bài ki m tra s 4 Khi ti n hành th nghi m 2 lo i vacxin (A và B); ñ i v i vacxin A, sau khi ti n hành thí nghi m trên 1350 gà thí nghi m quan sát th y 250 con ch t; tương t v i vacxin B, t 2535 gà th y 500 con ch t. M t bác s thú y quy t ñ nh ch n vacxin A ñưa vào phòng b nh cho gà tr i chăn nuôi c a mình. Anh (ch ) cho bi t quy t ñ nh ñó ñúng hay sai, vì sao? 61
4. Ph l c B NG S NG U NHIÊN 81 37 66 40 77 65 29 99 77 42 92 78 15 25 07 76 79 24 21 84 48 03 48 91 03 57 56 56 42 76 57 27 60 60 16 30 76 96 94 49 86 49 52 63 66 70 80 71 09 64 84 36 03 54 53 39 36 30 69 27 73 59 16 61 43 18 86 80 19 42 23 78 86 08 44 08 55 51 12 97 10 46 82 01 40 55 50 91 24 12 34 43 20 37 71 52 13 25 67 31 63 34 98 49 54 23 60 36 10 40 08 12 34 46 59 82 91 74 60 92 18 40 40 07 42 21 10 22 39 57 86 80 03 29 64 96 73 84 72 47 59 86 66 45 91 17 29 15 92 05 97 60 76 48 44 58 89 64 01 26 30 99 69 70 16 08 76 29 74 90 18 42 43 71 47 22 10 21 08 69 14 49 02 64 25 44 27 12 36 82 67 84 58 21 61 72 45 23 63 43 99 76 35 87 72 35 14 61 70 33 94 30 18 23 70 30 80 72 72 04 50 42 77 64 94 44 17 80 67 98 72 15 00 52 41 76 16 85 33 23 10 38 18 55 57 31 38 12 97 80 91 47 94 45 67 92 31 55 16 91 46 52 61 13 33 04 30 47 97 11 30 03 87 98 33 06 29 77 56 41 29 21 02 78 61 84 33 50 43 75 42 28 40 16 12 42 03 44 10 28 83 59 26 14 81 77 04 94 98 12 33 71 07 29 35 25 86 82 52 43 87 22 31 54 76 04 80 79 92 37 97 31 53 34 10 57 19 48 32 86 73 53 23 83 40 45 57 33 18 29 13 61 64 03 38 09 01 88 13 14 29 32 83 46 27 05 18 31 46 93 59 83 90 79 53 91 47 02 26 90 70 71 37 04 12 71 30 23 31 51 92 96 09 93 08 52 94 79 45 34 87 29 28 54 53 54 33 39 22 61 46 98 84 24 28 71 42 75 98 07 83 78 88 92 75 35 07 41 70 05 83 13 45 06 24 89 75 66 06 27 69 26 97 35 72 95 58 30 84 12 70 41 36 92 05 62 89 01 62 31 07 82 88 94 99 80 07 37 94 52 15 26 90 39 39 51 53 40 98 78 55 80 29 81 32 27 28 59 29 74 27 46 15 47 00 47 94 04 03 43 80 73 03 69 35 68 22 77 82 26 83 58 62 71 77 88 00 70 45 58 45 69 97 79 98 33 45 64 83 62 20 36 34 64 67 29 08 47 56 72 25 15 57 13 07 95 01 02 02 70 86 74 56 14 94 33 49 73 62 71 82 87 56 32 99 86 35 13 22 12 25 90 89 20 82 87 46 23 14 27 00 98 13 94 00 85 09 30 97 98 72 40 81 87 33 96 58 28 08 64 61 99 16 38 11 08 28 65 70 71 79 51 31 38 27 99 64 57 99 98 79 93 50 34 41 50 21 49 74 52 03 52 53 24 89 53 96 19 31 06 36 19 99 62 65 08 46 68 44 96 73 98 65 41 72 37 46 27 11 41 88 27 35 22 39 59 19 39 65 55 59 20 25 48 23 61 78 35 48 89 24 20 27 94 31 17 47 50 37 11 15 19 46 34 23 80 37 60 30 50 54 55 44 08 73 05 63 52 47 43 82 40 98 97 92 13 46 31 02 67 83 93 99 35 06 85 63 39 04 12 93 91 86 88 63 68 62 75 91 38 64 64 87 77 53 05 29 76 06 23 88 81 10 33 02 86 86 93 12 00 74 72 31 23 20 17 06 56 26 91 86 60 48 28 08 93 56 03 26 44 81 76 68 15 22 70 38 56 71 59 69 38 45 64 79 98 69 02 11 90 62
B NG XÁC SU T C A PHÂN B TIÊU CHU N HOÁ Các giá tr trong b ng là c a phân b chu n v i trung bình b ng 0 và ñ l ch chu n là 1. ng v i m i giá tr z trong là giá tr P, P(Z < z). z P z P z P z P -4,00 0,00003 -1,50 0,0668 0,00 0,5000 1,55 0,9394 -3,50 0,00023 -1,45 0,0735 0,05 0,5199 1,60 0,9452 -3,00 0,0013 -1,40 0,0808 0,10 0,5398 1,65 0,9505 -2,95 0,0016 -1,35 0,0885 0,15 0,5596 1,70 0,9554 -2,90 0,0019 -1,30 0,0968 0,20 0,5793 1,75 0,9599 -2,85 0,0022 -1,25 0,1056 0,25 0,5987 1,80 0,9641 -2,80 0,0026 -1,20 0,1151 0,30 0,6179 1,85 0,9678 -2,75 0,0030 -1,15 0,1251 0,35 0,6368 1,90 0,9713 -2,70 0,0035 -1,10 0,1357 0,40 0,6554 1,95 0,9744 -2,65 0,0040 -1,05 0,1469 0,45 0,6736 2,00 0,9772 -2,60 0,0047 -1,00 0,1587 0,50 0,6915 2,05 0,9798 -2,55 0,0054 -0,95 0,1711 0,55 0,7088 2,10 0,9821 -2,50 0,0062 -0,90 0,1841 0,60 0,7257 2,15 0,9842 -2,45 0,0071 -0,85 0,1977 0,65 0,7422 2,20 0,9861 -2,40 0,0082 -0,80 0,2119 0,70 0,7580 2,25 0,9878 -2,35 0,0094 -0,75 0,2266 0,75 0,7734 2,30 0,9893 -2,30 0,0107 -0,70 0,2420 0,80 0,7881 2,35 0,9906 -2,25 0,0122 -0,65 0,2578 0,85 0,8023 2,40 0,9918 -2,20 0,0139 -0,60 0,2743 0,90 0,8159 2,45 0,9929 -2,15 0,0158 -0,55 0,2912 0,95 0,8289 2,50 0,9938 -2,10 0,0179 -0,50 0,3085 1,00 0,8413 2,55 0,9946 -2,05 0,0202 -0,45 0,3264 1,05 0,8531 2,60 0,9953 -2,00 0,0228 -0,40 0,3446 1,10 0,8643 2,65 0,9960 -1,95 0,0256 -0,35 0,3632 1,15 0,8749 2,70 0,9965 -1,90 0,0287 -0,30 0,3821 1,20 0,8849 2,75 0,9970 -1,85 0,0322 -0,25 0,4013 1,25 0,8944 2,80 0,9974 -1,80 0,0359 -0,20 0,4207 1,30 0,9032 2,85 0,9978 -1,75 0,0401 -0,15 0,4404 1,35 0,9115 2,90 0,9981 -1,70 0,0446 -0,10 0,4602 1,40 0,9192 2,95 0,9984 -1,65 0,0495 -0,05 0,4801 1,45 0,9265 3,00 0,9987 -1,60 0,0548 0,00 0,5000 1,50 0,9332 3,50 0,99977 -1,55 0,0606 4,00 0,99997 M t vài giá tr t i h n c a z: 0,80 0,90 0,95 0,975 0,99 0,995 0,999 P 0,842 1,282 1,645 1,960 2,326 2,576 3,090 z 63