Bài viết trình bày tình hình nuôi cá tra và dịch bệnh gan thận mủ trên cá tra năm 2014 tại An Giang và Đồng Tháp, thu thập thông tin về dịch bệnh gan thận mủ trên cá tra được các cơ sở nuôi ghi chép trong năm 2014 và thu thập thông tin về các yếu tố nguy cơ theo biểu mẫu để tổng hợp phân tích mối liên quan đến bệnh gan thận mủ tại các cơ sở nuôi cá tra.
KHOA HỌC KỸ THUẬT THÚ Y TẬP XXIII SỐ - 2016 MỘT SỐ ĐẶC ĐIỂM DỊCH TỄ BỆNH GAN THẬN MỦ Ở CÁ TRA TẠI CÁC TỈNH AN GIANG VÀ ĐỒNG THÁP NĂM 2014 Bùi Thị Việt Hằng1,2, Nguyễn Thế Hiền1, Nguyễn Thị Lan Hương1, Võ Đình Chương1, Nguyễn Thị Việt Nga1, Kim Văn Vạn2, Nguyễn Văn Long1 TÓM TẮT Nghiên cứu số đặc điểm dịch tễ bệnh gan thận mủ cá tra tiến hành 600 sở nuôi cá 23 huyện thuộc tỉnh An Giang Đồng Tháp Kết nghiên cứu cho thấy: - Dịch bệnh xuất hầu hết huyện điều tra; thường từ tháng – tháng 11 hàng năm, với tỷ lệ trung bình 45,8% (95% CI 41,78 – 49,86) hộ điều tra; tỷ lệ ao điều tra có bệnh 37,30% (95% CI 34,57 – 40,08); tỷ lệ diện tích thả ni bị bệnh 34,63% (95% CI 31,43 – 37,95); - Kết phân tích yếu tố nguy cho thấy: Tỷ số chênh nguy xảy dịch bệnh 5,01 (95% CI 2,18 – 12,95) lần mật độ 200 con/m2 so với mật độ từ 10-50 con/m2; tỷ số chênh 2,21 (95% CI 1,40 – 3,55) lần sở thả nuôi không vét bùn đáy ao so với sở ni có vét bùn; tỷ số chênh 3,58 (95% CI 2,34 – 5,48) lần sở nuôi không khử trùng dụng cụ nuôi so với sở ni có khử trùng dụng cụ Từ khóa: Cá tra, Bệnh gan thận mủ, Tỷ lệ nhiễm bệnh, Yếu tố nguy cơ, Tỉnh An Giang, Đồng Tháp Some epidemiological characteristics of Bacillary necrosis of Pangasius in An Giang, Dong Thap Provinces in 2014 Bui Thi Viet Hang, Nguyen The Hien, Nguyen Thi Lan Huong, Võo Dinh Chuong, Nguyen Thi Viet Nga, Kim Van Van, Nguyen Van Long SUMMARY Study on some epidemiological characteristics of Bacillary necrosis of Pangasius was conducted at 600 Pangasius culture farms in 23 districts of An Giang, Dong Thap provinces The studied result showed that: Epidemic occurred in most of the surveyed districts, normally from May to November yearly, with the average disease infection rate of the investigated farms was 45.8% (95% CI 41.78 – 49.86) Of which, the infection rate of pond was 37.30% (95% CI 34.57 – 40.08), the infection rate of culture areas was 34.63% (95% CI 31.43 – 37.95) The result of risky analysis indicated that the odds ratio of disease outbreak between the stocking rate of 200 fish/m2 and 10-50 fish/m2 was 2.21 (95% CI 1.40 – 3.55); the odds ratio between the culture ponds without removing mud and with removing mud was 2.21 (95% CI 1.40 – 3.55) This ratio between the farms disinfected farming tools and the farms did not disinfected farming tools was 3.58 (95% CI 2.34 – 5.48) Keywords: Pangasius, Bacillary necrosis, Infection rate, Risky factors, An Giang, Dong Thap provinces Cục Thú y Khoa Thủy sản, Học viện Nông nghiệp Việt Nam 66 KHOA HỌC KỸ THUẬT THÚ Y TẬP XXIII SỐ - 2016 ĐẶT VẤN ĐỀ Cá tra sản phẩm nông nghiệp đề nghị đưa vào chương trình sản phẩm quốc gia Chính phủ ban hành nghị định riêng quy định nuôi, chế biến xuất sản phẩm cá tra (Nghị định số 36/2014/ NĐ-CP ngày 29/4/2014 Chính phủ) Điều khẳng định vị trí ngành nuôi cá tra phát triển kinh tế nước ta Hiện nay, tổng diện tích ni cá tra ước đạt 4.900 ha, tập trung chủ yếu tỉnh An Giang, Đồng Tháp, Vĩnh Long, Cần Thơ Hậu Giang Theo báo cáo địa phương năm gần cho thấy, tình hình dịch bệnh cá tra có chiều hướng gia tăng mạnh, tác động tiêu cực đến q trình ni, xuất cá tra nước ta Đã có số nước yêu cầu Việt Nam phải có kế hoạch phịng, chống dịch bệnh cá tra, bệnh gan thận mủ, xuất huyết Mặc dù nghề nuôi cá tra nước ta phát triển mạnh, đem lại nguồn ngoại tệ lớn nhờ việc xuất khẩu, việc phòng, chống, nghiên cứu điều tra dịch tễ dịch bệnh cá tra chưa thực quan tâm triển khai nhiều Cụ thể, việc phòng, chống, nước ta chưa có quy định cụ thể cho dịch bệnh cá tra; trước thực theo quy định Thông tư số 17/2014/TT-BNNPTNT ngày 20/6/2014 Bộ Nông nghiệp Phát triển nơng thơn quy định phịng chống dịch bệnh cho động vật thủy sản Tuy nhiên, hiệu phòng, chống hạn chế quy định chưa cụ thể, chưa thực phù hợp dịch bệnh cá tra Mặt khác, thiếu điều tra, nghiên cứu dịch tễ bệnh cá tra, nên biện pháp chưa có tính thuyết phục, chưa dựa sở khoa học nên tính khả thi cịn hạn chế Vì vậy, chúng tơi tiến hành đề tài “Nghiên cứu số đặc điểm dịch tễ bệnh gan thận mủ cá tra tỉnh An Giang Đồng Tháp năm 2014” II NỘI DUNG, NGUYÊN LIỆU VÀ PHƯƠNG PHÁP NGHIÊN CỨU 2.1 Nội dung nghiên cứu - Tình hình ni cá tra dịch bệnh gan thận mủ cá tra năm 2014 An Giang Đồng Tháp - Thu thập thông tin dịch bệnh gan thận mủ cá tra sở nuôi ghi chép năm 2014 - Thu thập thông tin yếu tố nguy theo biểu mẫu để tổng hợp phân tích mối liên quan đến bệnh gan thận mủ sở nuôi cá tra 2.2 Nguyên liệu - Số liệu nuôi cá tra ổ dịch bệnh gan thận mủ Chi cục Thú y nêu thu thập báo cáo Cục Thú y năm 2013 2014 - Bộ phiếu điều tra Cục Thú y thiết kế hướng dẫn việc thu thập thông tin, liệu - Số liệu địa lý chi tiết đến cấp xã năm 2011 Nhà xuất Tài nguyên Môi trường Bản đồ Việt Nam xây dựng cung cấp (2012) - Phần mềm vẽ đồ ArcGIS 9.3 phần mềm phân tích thống kê R (R Development Core Team, 2012) gói phân tích tương ứng epiR (Stevenson, 2012), lme4 (Bates and Sarkar, 2007) 2.3 Phương pháp Nghiên cứu cắt ngang (cross-sectional study) ứng dụng để thu thập số liệu sở nuôi ghi chép năm 2014 2.3.1 Lựa chọn địa điểm nghiên cứu Tại tỉnh điều tra (An Giang Đồng Tháp), Chi cục Thú y lập danh sách tất sở nuôi cá tra địa bàn tỉnh, sau chọn ngẫu nhiên 300 sở nuôi cá tra để vấn, thu thập thông tin 67 KHOA HỌC KỸ THUẬT THÚ Y TẬP XXIII SỐ - 2016 2.3.2 Thu thập thông tin yếu tố nguy Thông tin yếu tố nguy thu thập câu hỏi thông tin chung cho 607 sở nuôi cá tra (chuẩn bị ao nuôi, chọn giống thả giống ao ni, chăm sóc quản lý ao, dịch bệnh ao nuôi) 2.3.3 Quản lý phân tích số liệu - Thơng tin liệu từ thực địa nhập vào file MS Excel - Phân tích mơ tả tình hình ni cá tra, dịch tễ học mô tả (theo không gian, thời gian theo đối tượng) theo phương pháp dịch tễ, thống kê thường quy - Tỷ lệ sở nuôi cá tra bị bệnh tính theo phương pháp Fleiss (Fleiss, 1981) So sánh tỷ lệ sở bị bệnh phân tích, dựa vào số Chi-square - Phân tích đa tầng nhiều biến (multilevel analysis) áp dụng để định lượng yếu tố nguy theo phân tầng khác (Dohoo et al., 2001; Dohoo, 2003; Long, 2013) , cụ thể theo bước sau: + Bước 1: Xác định “Đầu ra” sở có cá tra bị bệnh gan thận mủ + Bước 2: Phân tích sàng lọc (bivariate analyses) để xác định mối liên hệ “Đầu ra” yếu tố nguy Phương pháp kiểm tra Wald (Agresti, 2007) sử dụng để xác định mối liên hệ biến nguy với “Đầu ra” Tất biến nguy có mối liên hệ với “Đầu ra” mặt thống kê sinh học có giá trị P < 0.2 giữ lại để đưa vào phân tích đa biến (Bước 3) Những biến nguy có mối liên hệ với “Đầu ra” giá trị P > 0.2 loại bỏ, theo nguyên tắc, loại bỏ biến có giá trị P cao cịn biến có P < 0.2 + Bước 3: Phân tích đa biến (multivariate analyses): Tất biến có P < 0.2 xác định bước đưa vào mơ hình đa biến chạy mơ hình phân tích Chỉ biến nguy có mối liên hệ mặt thống kê sinh học với “Đầu ra” có giá trị P < 0.05 giữ lại để tính ảnh hưởng biến đa cấp bước + Bước 4: Yếu tố ảnh hưởng ngẫu nhiên đưa vào mơ hình hồi quy logic hỗn hợp (mixed-effects logistic regression model) để phân tích đánh giá mức độ ảnh hưởng đo yếu tố cấp sở nuôi cá tra, cấp xã cấp tỉnh lên khả sở nuôi cá tra bị bệnh gan thận mủ Dựa cách tiếp cận này, chuyển dạng logic (logit transform) sở bị bệnh i xã j tỉnh k, pijk, mơ hình hóa theo chức tuyến tính loạt yếu tố ảnh hưởng m β1… βm ảnh hưởng yếu tố nguy cấp độ sở nuôi, cấp xã Xj tỉnh Pk tính theo cơng thức sau: m pijk log = β + ∑ β m x mijk + X j + Pk + ε ijk i =1 1 − pijk Kết mơ hình cuối thể tỷ số chênh điều chỉnh cho biến nguy Tỷ số chênh (odds ratio, OR) lớn cho thấy sở nuôi cá tra phơi nhiễm với yếu tố nguy có khả bị bệnh tăng lên Ngược lại, OR < cho thấy nguy cơ sở nuôi cá tra bị bệnh giảm OR = cho thấy khơng có ảnh hưởng yếu tố lên nguy cơ sở nuôi bị bệnh 68 Đường cong đặc trưng thể hoạt động thu nhận (Receiver Operating Characteristic, ROC) vẽ để biểu thị khả dự đoán mơ hình logic Vùng đường cong ROC có giá trị dao động từ 0-1 để mô tả khả dự đốn mơ hình sở bị bệnh không bị bệnh Giá trị đường cong ROC lớn, cho thấy khả mô hình đốn xác cao (Hosmer and KHOA HỌC KỸ THUẬT THÚ Y TẬP XXIII SỐ - 2016 Lemeshow, 2000) Phân tích thống kê tiến hành cách sử dụng phân tích thống kê lme4 package (Bates and Sarkar, 2007) chương trình R 2.15.2 (R Development Core Team, 2012) III KẾT QUẢ NGHIÊN CỨU 3.1 Tình hình ni cá tra tỉnh An Giang Đồng Tháp 3.1.1 Nuôi cá tra tỉnh An Giang - Năm 2014, diện tích ni cá tra 1.269 ha, 95,35% so kỳ năm 2013 Cá tra nuôi khoảng 1.793 sở, hộ dân Trong đó, ni tập trung chủ yếu huyện, gồm: Chợ Mật độ nuôi sở điều tra Mới, Châu Phú, Châu Thành, Thoại Sơn Phú Tân, với tổng diện tích từ 100 đến 200 - Diện tích vùng ni doanh nghiệp 538 ha, chiếm 42,53% (tăng 70,79%), cịn lại nơng hộ chiếm khoảng 57,47% so với kỳ năm 2013; doanh nghiệp có vùng ni địa bàn tỉnh xây dựng mơ hình ni cá tra thương phẩm xuất đạt tiêu chuẩn an tồn chất lượng theo tiêu chuẩn quốc tế Theo tổng diện tích cá tra chứng nhận đạt tiêu chuẩn quốc tế 240 ha, với sản lượng 105.000 tấn/năm, chiếm tỷ lệ khoảng 40% tổng sản lượng nguyên liệu xuất (hình 1) Mức độ bệnh sở điều tra Hình Bản đồ thể mật độ nuôi cá tra, dịch bệnh gan thận mủ hộ điều tra tỉnh An Giang năm 2014 3.1.2 Nuôi cá tra tỉnh Đồng Tháp - Năm 2014, cá tra nuôi 700 sở 81 xã thuộc 12 huyện/ thị/ thành phố thuộc tỉnh Đồng Tháp, với tổng diện tích khoảng 2.070 (Bảng 2) Trong đó, diện tích ni cá tra thương phẩm chiếm khoảng gần 50% (1.000 ha) - Số liệu tổng hợp cho thấy Đồng Tháp tỉnh có tổng diện tích mặt nước ni cá tra cao nước (chiếm khoảng 28% tổng diện tích ni nước), với tổng sản lượng cá tra đạt khoảng 386.610 (chiếm khoảng 27% tổng sản lượng cá tra nước) So với tỉnh An Giang, diện tích ni sản lượng cá tra tỉnh Đồng Tháp cao gần gấp đơi; có nhiều sở nuôi cá tra tỉnh Đồng Tháp theo hướng thâm canh (hình 2) 69 KHOA HỌC KỸ THUẬT THÚ Y TẬP XXIII SỐ - 2016 Mật độ nuôi sở điều tra Mức độ bệnh sở điều tra Hình Bản đồ thể mật độ nuôi cá tra, dịch bệnh gan thận mủ hộ điều tra tỉnh Đồng Tháp năm 2014 3.2 Tình hình dịch bệnh gan thận mủ cá tra 3.2.1 Dựa số liệu Chi cục Thú y tỉnh báo cáo - Tại tỉnh An Giang: Kết điều tra Cục Thú y cho thấy dịch bệnh gan thận mủ cá tra xuất tỉnh vào năm 2014, Chi cục Thú y không báo cáo số liệu theo quy định hành - Tại tỉnh Đồng Tháp: Dịch bệnh gan thận mủ xuất 22 xã thuộc huyện (Châu Thành, Hồng Ngự, Lai Vung, Thanh Bình, Tháp Mười thị xã Hồng Ngự), với tổng số 349/4.530 ha, chiếm khoảng 7,96% diện tích ni huyện + Đặc điểm dịch tễ không gian: Dịch bệnh xuất hai huyện Thanh Bình (245 ha) huyện Lai Vung (88 ha) cao so với huyện khác + Đặc điểm dịch tễ thời gian: Dịch bệnh gan thận mủ xuất nhiều giai đoạn từ tháng 1-3/2014 giai đoạn từ tháng 6-9/2014 70 + Đặc điểm dịch tễ đối tượng cá tra: Dịch bệnh xuất chủ yếu cá tra nuôi thương phẩm, theo hướng công nghiệp, nuôi tập trung 3.2.2 Dựa số liệu Cục Thú y tổ chức điều tra thu thập - Đặc điểm dịch tễ không gian: Dịch bệnh xuất hầu hết huyện điều tra hai tỉnh (Bảng 2) với tỷ lệ 45,8% (95% CI 41,78 – 49,86) hộ điều tra có dịch bệnh; tỷ lệ ao điều tra có bệnh 37,30% (95% CI 34,57 – 40,08); tỷ lệ diện tích thả ni bị bệnh 34,63% (95% CI 31,43 – 37,95) Các tỷ lệ tỉnh Đồng Tháp cao so với tỉnh An Giang Tại tỉnh Đồng Tháp, kết điều tra cho thấy huyện Lai Vung địa bàn có diện tích bị bệnh nhiều (phù hợp với số liệu Chi cục Thú y), tiếp huyện Thanh Bình - Đặc điểm dịch tễ thời gian: Dịch bệnh xuất từ đầu tháng kéo dài đến hết tháng 11/2014, từ đầu tháng có dịch bệnh xuất tỉnh Đồng Tháp KHOA HỌC KỸ THUẬT THÚ Y TẬP XXIII SỐ - 2016 - Đặc điểm dịch tễ đối tượng cá tra: Tỷ lệ cá bị bệnh 9,52% (95% CI 9,52 – 9,52) số cá giống thả ni; tỷ lệ cá chết bệnh 49,16% (95% CI 49,15 – 49,17) Như vậy, thấy mức độ lây lan bệnh gan thận mủ không cao (9,52% cá thả nuôi bị bệnh), cá bị bệnh khả chết lại cao (49,16%), xem bảng bảng Bảng Tình hình dịch bệnh xuất huyết bệnh gan thận mủ cá tra nuôi tỉnh điều tra Chỉ số so sánh Số huyện có bệnh Số xã có bệnh Số hộ ni Số ao ni Diện tích ni (ha) Số giống thả ban đầu (con) Đồng Tháp 30 306 772 682 536.728.700 An Giang 32 301 448 164 59.146.000 - ao 39 583 1018 595 528.808.700 - 10 ao 19 18 118 170 27.218.000 > 10 ao 84 81 39.848.000 I Tỉnh II Số ao nuôi III Số ao lắng 39 560 1029 519 507.709.700 20 42 155 290 76.220.000 22 19 8.912.000 1 14 18 3.033.000 IV Thời gian xuất bệnh sau thả ngày 12 24 20 19.981.000 - 10 ngày 36 223 504 343 372.678.900 11 - 20 ngày 14 22 70 69 28.000.000 > 20 ngày 21 76 67 18.039.000 329 546 347 157175800 607 1,220 846 595.874.700 Không bệnh Tổng cộng 62 3.3 Phân tích xác định yếu tố nguy liên quan đến dịch bệnh “cơ sở có cá tra bị bệnh gan thận mủ” Một số mối liên quan điển sau: 3.3.1 Phân tích nhị biến xác định mối liên quan yếu tố nguy với “cơ sở có dịch bệnh hay khơng có dịch bệnh” - Tỷ số chênh (OR) bị bệnh gan thận mủ sở thả nuôi với mật độ 200 con/m2 5,01 (95% CI 2,18 - 12,95) lần so với sở nuôi với mật độ từ 10-50 con/m2 Kết trình bày bảng Kết phân tích nhị biến (bảng 3) cho thấy, tổng số 21 yếu tố nguy hỏi thu thập thơng tin, có 18 yếu tố nguy có mối liên hệ sinh học, thống kê (P < 0.2) với - Tỷ số chênh bị bệnh gan thận mủ sở thả nuôi không vét bùn đáy ao 2,21 (95% CI 1,40 - 3,55) lần so với sở nuôi có vét bùn đáy ao 71 72 123 An Giang 11 6 - 10 ao > 10 ao 24 1 69 380 24 32 399 150 305 Số ao bệnh 22 21 11 - 20 ngày > 20 ngày Tổng cộng 278 223 - 10 ngày Không bệnh 12 ngày 455 66 44 323 22 Khoảng thời gian bị bệnh 251 Số ao lắng 261 - ao Số ao nuôi 155 Số hộ bệnh Đồng Tháp Tỉnh Chỉ số so sánh 293 62 42 170 19 107 181 19 32 242 56 237 Diện tích bị bệnh (ha) 56,739,130 4610562 4119598 47067970 941000 1000 620000 10648860 45469270 3502973 3225850 50010307 11982070 44757060 Số cá bệnh (con) 27,892,344 1050528 1209896 25384360 247560 100 38000 3302479 24551765 714446 664015 26513883 1300910 26591434 Số cá chết (con) 45,8 (41,78 - 49,86) 100,00 (77,16 - 100,00) 100,00 (78,05 - 100,00) 100,00 (97,54 - 100,00) 100,00 (63,97 - 100,00) 100,00 (1,26 - 100,00) 50,00 (6,76 - 93,24) 57,14 (40,96 - 72,28) 44,82 (40,65 - 49,05) 100,00 (42,13 - 100,00) 61,11 (35,75 - 82,70) 44,77(40,68 - 48,91) 40,86 (35,26 - 46,65) 50,65 (44,91 - 56,39) Tỷ lệ hộ bệnh (95% CI) 37,30 (34,57 - 40,08) 86,84 (77,13 - 93,51) 62,86 (50,48 - 74,11) 64,09 (59,73 - 68,28) 91,67 (73,00 - 98,97) 7,14 (0,18 - 33,87) 22,73 (7,82 - 45,37) 44,52 (36,54 - 52,70) 36,93 (33,97 - 39,96) 28,57 (19,24 - 39,47) 27,12 (19,35 - 36,08) 39,19 (36,18 - 42,27) 33,48 (29,12 - 38,06) 39,51 (36,04 - 43,06) Tỷ lệ ao bệnh (95% CI) 34,63 (31,43 - 37,95) 92,42 (83,32 - 97,45) 60,53 (47,99 - 72,12) 49,53 (44,12 - 54,95) 95,93 (76,22 - 99,95) 7,17 (0,36 - 29,68 ) 18,76 (4,72 - 43,36) 36,87 (31,30 - 42,71) 34,88 (30,78 - 39,15) 23,54 (14,81 - 34,30) 18,74 (13,17 - 25,43) 40,66 (36,69 - 44,73) 34,06 (26,84 - 41,86) 34,76 (31,18 - 38,47) Tỷ lệ diện tích bệnh (95% CI) 9,52 (9,52 - 9,52) 25,56 (25,54 - 25,58) 14,71 (14,70 - 14,72) 12,63 (12,63 - 12,63) 4,71 (4,70 - 4,72) 0,03 (0,03 - 0,04) 6,96 (6,94 - 6,97) 13,97 (13,96 - 13,98) 8,96 (8,95 - 8,96) 8,79 (8,78 - 8,80) 11,85 (11,84 - 11,86) 9,46 (9,45 - 9,46) 20,26 (20,25 - 20,27) 8,34 (8,34 - 8,34) Tỷ lệ cá bệnh (95% CI) Tỷ lệ cá chết bệnh (95% CI) 49,16 (49,15 - 49,17) 22,79 (22,75 - 22,82) 29,37 (29,33 - 29,4) 53,93 (53,92 - 53,94) 26,31 (26,23 - 26,38) 10,00 (8,21 - 12,03) 6,13 (6,07 - 6,19) 31,01 (0,99 - 31,04) 54,00 (53,99 - 54,010 20,40 (20,36 - 20,43) 20,58 (20,54 - 20,62) 53,02 (53,01 - 53,03) 10,86 (10,84 - 10,87) 59,41 (59,40 - 59,42) Bảng Tình hình dịch bệnh xuất huyết bệnh gan thận mủ cá tra nuôi tỉnh điều tra KHOA HỌC KỸ THUẬT THÚ Y TẬP XXIII SỐ - 2016 KHOA HỌC KỸ THUẬT THÚ Y TẬP XXIII SỐ - 2016 Bảng Kết phân tích sàng lọc nhị biến (bivariate analysis) yếu tố nguy đầu (cơ sở ni cá tra có dịch bệnh gan thận mủ) TT Yếu tố nguy Tỉnh Diện tích thả nuôi (ha) Số lượng ao nuôi Ao lắng Mật độ thả (con/m2) Độ sâu ao Tuổi ao nuôi Thả xen kẽ Khử trùng dụng cụ 10 Cấp thoát nước riêng 11 Vét bùn đáy ao 12 Phơi đáy 13 14 15 16 17 18 19 20 21 Chlorine khử trùng đáy ao Lọc nước ao Diệt tạp chlorine Sử dụng BKC Sử dụng vôi bột Khử trùng trước thả Khử trùng sau thả Thay nước ao Hút bùn đáy ao Tổng Phân loại Cơ sở điều tra Cơ sở có dịch bệnh OR (95% CI) P An Giang 301 123 1,00 Tham chiếu Đồng Tháp 306 155 1,49 (1,08 - 2,05) 0,01568 10 17 13 4,07 (1,42 - 14,60) 0,01515 ao 344 126 1,00 Tham chiếu ao 116 59 1,79 (1,17- 2,74) 0,0072 > ao 147 93 2,98 (2,00 - 4,47)