STATA GIAI PHAP MOI CHO PHAN TICH SO LIEU NGHIEN CUU LAM NGHIEP

en.t-1 : $:i3 nH:H Tgp chf F{}RESTRt' ${II*iH(:ffi AIUl} "il"lit'HN (}[,t ]f ;Y J{XjRlq,-{fi., {}H \'IET\.I \I F{ TRT{S'TRV I I}iI}'IIRSITY MTJC LUC T1p chi: ilHoA HtlS & c0ilG un ffiHf /.s.t-\i 1859 Trang Nguy6n Thi Minh H[ng, fguy6n Minh Thu Ph6n lQp vd tuy6n - 10 chon m6t s6 chung vi khu6n Lactic c6 kha ndng sinh t6ng hqp Amvlase vd Bacteriocm mnrEp - 3828 ' Nguy6n Thi Hdng G6m, D6 Quang.Trung, Nguy6n Thi Minh 11- 15 Hdng, Hd Hf,i Ninh Buoc d6u nh6n gi6ng hoa Tu1'lip vitng (Tulipa gesneriana) bdng k! thu{t nu6i cFry In Vitro Nguy6n Quj'nh Trang, Vfr Thi HuQ, Khudt Thl Hii Ninh, 16 - 2l Nguy€n Thi Tho Nhdn gi6ng In Vitro lan Phi Di6p tim (Dendrobium anosmun) Cong ngh€srnh hgc & Giatng elytdng Duong Tryng Hi6u, Nguy6n Thanh Hir Thuc trang nuong 16'rd m6 hinh phuc h6i nmg tr6n ddt nuong rdy dreo hucrng n6ng l6m k6t hqp tai vung h6 YCn LAp, tinh Quang Ninh 22- 29 vfl Quang Nam, LC XuAn Thfng o6.q.nn Tufln M6t so d* die- 30 - 37 sinh hoc vd sinh thai hoc loai Md sa pa t4i Vuon quoc gia Hodng Lien- hnh LAo Cai Ltun sinh Quan lf nhg & M6i Tdi nguydn trttutg L€ Dinh Phuong, O6lntr Tudn M6t s6 aa" Aeto noi moc vd tar srnh 38 - {5 cua loai Cioian hatlMichelia tonkinensis A.Cher) tar Vuon quoc gia Ben En- Thanh H6a Bii Dinh Dric, Bni V[n Bfc Nghi6n cuu dd mat bien phap r-at ly, ccr Cong nghiQp rimg Kinh ti' & chinh s6ch ttng dgng COng nghQ thdng tin ' so J \-{\I rl.i 2013 Ir 16 - 52 gicri phong tru Sdu r6m lum l6ng (Dasychra crcuthq Coliennette) hai Th6ng md rii't4i LQc Binh - Lang Son D0ng Thanh Hei, Vfl Ti0n Thinh, Ph4m Anh Cudng, TrAn Vin Dffng, Giang Trgng Toin, Nguy€n Chf Thhnh, Nguy6n Xuin Trud'ng, Biri Hda Binh D6nh gi6 muc d0 uu tien ket noi cua c6c khu nrng dac dung tai miAn B6c b6i cAnh bien dOi hflu Ph4m Hfru Khrlnh, Vfr Tii5n Thinh Du b6o bi6n d6ng qudn Bd t6t (Bos gaurtts) tai Vuon qu6c gia C6t Ti0n 53 - 61 th6 62 - 66 s6 67 - 72 Da 73 - 81 Nguy'6n Hftu Quin, Nguy6n ThC NhA, Ph4m Hfi'u Hing M6t dac di6m sinh hoc cua Bo que hai Lu6ng (Baculum apicalis Chen et He) tai Thanh H6a Nguy6n Vinh Thanh, Trdn Trung Thhnh,.Cn Thi Khuy€n dang thdnh ph6n lodi thu vd hi6n trang m6t s6 loai thu o khu rung thu6c xd Y6n D6ng, huy6n Y6n M6, tinh Ninh Binh Vfi Ti0n Thinh Nghidn cuu t0p tinh hoat d6ng cia lodi.Don 32-85 (Alherurus mqcrouru,\ Linnaeus, 1758) dieu kien nuoi nh6t Gi5y pf€p so: hop 96 nhua -li_ JJ = -,lr-rl l- -^-4-.rC.L C'gJ '* Jg^ tnQng ' r^ -: Quich Vin Thi0m, TrAn Vin Chri Nghi6n cuu anh huong nhi6t 86 - 91 dQ ep ddu voi phun toi d0 b6n keo vd d0 b6n u6n cua v6t li6u phirc 7='2 Trinh Hidn Mai, Nguyel Phan Thi6t, Phan Tirng Hung.Anh \alo huong 92- l0l cna thoi gian ngdm t6m d6n ch6t luqng g5 Keo Lai (Acacia mangium x lccrcia auncuffirmis) vd g5 M6 (Mangtietia conifera) bi6n tinh bang hat SiO: Dao Lan Phuong Hoirn thi6n phuong ph6p xac dinh gi6 tr tar san ngan 102 - I r a tiu san dar han cua doanh nghidp nhdnuoc lrong lirh luc n6ng l6m nSluip l0 hui :^^ ^^ Irh6n Lt uu VUrry h, -':r$ig rnai Van Xudn +DonG-l-fAt'CI ' Chu Thi Thu, Ph4m Thanh Qu6 Quan 11 nhir nuoc doi voi :ganh tai ngu1.6n thi6n nhi6n vd m6i truong b[ng cric cong cu kinh :- EIS r Kinh nghiOm th6 gioi r a Viet Nam ,/ \am Thfng Ifng dung phan mdm Moodle xdy t:.ong bai giang tryc tuy6n tai truong Dai hoc Ldm nghiOp Phung :: Biri \I4nh Hung, Nguy6n Thi Bfch Phuqng Stata rn r ;ho phin tich so lieu nghien cuu l6m nghicp * gidi lll - l2l dung 122 - 127 ph6p 128 - 132 Ứng dụng công nghệ thông tin STATA – GIẢI PHÁP MỚI CHO PHÂN TÍCH SỐ LIỆU NGHIÊN CỨU LÂM NGHIỆP Bùi Mạnh Hưng, Nguyễn Thị Bích Phượng ThS Trường Đại học Lâm nghiệp TĨM TẮT Stata phần mềm thống kê phát triển từ năm 1985 Stata phần mềm thể sức mạnh vượt trội phân tích số liệu nghiên cứu nói chung Lâm nghiệp nói riêng Stata ứng dụng nhiều lĩnh vực như: Điều tra rừng, Lâm học, khoa học đất, Quản lý bảo vệ tài nguyên rừng, Chế biến lâm sản… Nó thực phân tích từ đơn biến đến đa biến, từ định tính đến định lượng Stata cung cấp nhanh chóng thông tin mẫu phân bố thực nghiệm thông qua lệnh tính tốn đặc trưng mẫu, lệnh vẽ biểu đồ xây dựng bảng biểu Stata cịn có khả phân tích việc so sánh mẫu quan sát tìm điều kiện tốt cho đại lượng lâm nghiệp Chương trình so sánh tiêu chuẩn tham số phi tham số, từ hai mẫu đến nhiều mẫu, từ mẫu độc lập đến mẫu liên hệ Hơn nữa, phân tích phương sai hai nhân tố ưu việt Stata Một chức ưu việt Stata phân tích tương quan hồi quy Stata phân tích tuyến tính biến đến đa biến Chương trình phân tích tương quan phi tuyến với loại hàm khó thực như: Gompertz, Schumacher, Koller, Verhulst-Robertso Vì vậy, việc khai thác sử dụng Stata phân tích số liệu nghiên cứu lâm nghiệp thực cần thiết vấn đề đáng quan tâm Từ khóa: Đặc trưng mẫu, phân tích phương sai, phân tích số liệu, thống kê, tương quan, so sánh, Stata I ĐẶT VẤN ĐỀ Stata phần mềm phân tích thống kê bắt đầu phát triển vào năm 1985 tập đoàn Stata (Wiki, 2012) Stata chữ viết tắt cụm từ “Statistics and data” tạm dịch “Thống kê số liệu” Phần mềm sử dụng rộng rãi để phân tích số liệu nghiên cứu nhiều lĩnh vực khác như: kinh tế, xã hội học, y tế… (Wiki, 2012) Stata liên tục phát triển, hồn thiện nâng cấp chức trình lệnh từ phiên phiên 12 Stata phần mềm mạnh phân tích số liệu nói chung, cịn ứng dụng để phân tích số liệu nghiên cứu lâm nghiệp Tuy nhiên, qua trình khai thác, ứng dụng Stata phân tích số liệu, chứng minh sức mạnh lĩnh vực Stata đem lại nhiều lời giải cho vấn đề khó khăn đặt phân tích số liệu lâm nghiệp Những vấn đề số lượng dung lượng mẫu thường lớn Ngoài ra, lĩnh vực lâm nghiệp có nhiều chun mơn sâu khác như: Điều tra rừng, Lâm sinh, Quản lý tài nguyên rừng, Chế biến lâm sản… Mỗi 128 chun mơn sâu lại có u cầu riêng phân tích Tuy nhiên, với hệ thống trình lệnh phong phú, Stata ứng dụng lĩnh vực khoa học Nó sử dụng để phân tích đơn biến đến đa biến, từ phân tích định tính đến phân tích định lượng, từ tính tốn đặc trưng mẫu đến lập phân bố thực nghiệm, so sánh mẫu quan sát, phân tích phương sai tương quan hồi quy Trong khuôn khổ báo tác giả trình bày khả ứng dụng Stata đáp ứng yêu cầu người sử dụng cơng tác phân tích số liệu nghiên cứu lâm nghiệp II PHƯƠNG PHÁP NGHIÊN CỨU Kế thừa số liệu lĩnh vực lâm nghiệp như: Điều tra rừng, Lâm học, Chế biến lâm sản, Quản lý tài nguyên rừng… để phục vụ cho việc chạy mô hình phân tích khác sau Tham khảo tài liệu, nghiên cứu phần mềm ứng dụng, phân tích số liệu khác có liên quan, tham khảo ý kiến chuyên gia lĩnh vực nghiên cứu khoa học lâm nghiệp, phân tích số liệu nghiên cứu lâm nghiệp nói riêng lĩnh vực khác nói chung TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ LÂM NGHIỆP SỐ (KỲ I) - 2013 Ứng dụng công nghệ thông tin Xây dựng qui trình kỹ thuật phần mềm Stata để thực nội dung nghiên cứu chạy thử quy trình lệnh với số liệu mơ hình thu thập cơng tác ngoại nghiệp III KẾT QUẢ NGHIÊN CỨU 3.1 Cung cấp thông tin mẫu lâm nghiệp phân bố thực nghiệm Yêu cầu phân tích số liệu lâm nghiệp thông tin thống kê mẫu, biến như: đường kính, chiều cao, độ ẩm, độ dày ván gỗ… thông tin phân bố thực nghiệm Những thông tin thực cần thiết để nhà nghiên cứu có nhìn khái quát, nhận định sơ mẫu nghiên cứu Đây sở ban đầu để đề xuất biện pháp kỹ thuật nói chung kỹ thuật lâm sinh nói riêng, nhằm nâng cao chất lượng rừng, cải thiện hiệu quản lý tài nguyên rừng, hay điều chế rừng theo mục đích kinh doanh phù hợp (Nguyễn Hải Tuất Nguyễn Trọng Bình, 2005) Với lệnh “Summary and descriptive statistics” (Tổng hợp thống kê mơ tả), Stata nhanh chóng tính tốn xuất kết đặc trưng mẫu như: dung lượng mẫu, trung bình mẫu, phương sai, độ lệch, độ nhọn… Stata cung cấp thông tin vị trí, biến động hình dạng phân bố Stata tính tốn cho mẫu nhỏ mẫu qua lập phân bố thực nghiệm Lệnh “Tables of summary Statistics (tables)” (Bảng thống kê) dễ dàng lập phân bố thực nghiệm cho biến rời rạc số tái sinh/ô, số sâu/cành hay số sinh viên/lớp Riêng với biến liên tục, Stata phần mềm đánh giá ưu việt để lập phân bố thực nghiệm Khi sử dụng Stata dễ dàng điều chỉnh số lượng tổ chia, cự ly tổ giá trị cận dưới, cận Đây ưu điểm trội Stata, phù hợp với phân tích số liệu nghiên cứu lâm nghiệp Hình 01 Hộp thoại histogram để điều chỉnh số tổ, cự ly tổ để lập phân bố thực nghiệm Lệnh “Distributional graphs” “Graphics” (biểu đồ) dễ dàng kiểm chứng mức độ phù hợp phân bố Hình 02 Các lệnh vẽ biểu đồ Stata thực nghiệm theo phân bố lý thuyết khác như: phân bố chuẩn, phân bố Chisquared, phân bố đối xứng… dựa vào phương TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ (KỲ I) - 2013 129 Ứng dụng cơng nghệ thơng tin pháp biểu đồ Ngồi ra, Stata trang bị hệ thống lệnh vẽ biểu đồ hình 02 Hệ thống lệnh giúp xây dựng nhiều loại biểu đồ khác nhau, từ biểu đồ dạng 2D, đến biểu đồ không gian 3D, từ biểu đồ dạng đường, cột đến biểu đồ tần số chuỗi thời gian Những biểu đồ xây dựng từ thơng tin biến Đó sở trực quan để đề xuất kết luận đối tượng nghiên cứu 3.2 Khả so sánh tìm điều kiện tốt cho đại lượng lâm nghiệp Trong phân tích số liệu lâm nghiệp, có nhiều trường hợp phải so sánh mẫu thu so sánh đường kính, chiều cao vị trí địa hình khác nhau, độ xác loại dụng cụ, loại thiết bị quan trắc khác hay tìm hiểu ảnh hưởng mật độ đến sinh trưởng địa… Những so sánh phân tích nhằm cho thấy mẫu có khác biệt hay khơng, nhân tố có thực ảnh hưởng tới kết thí nghiệm hay khơng? Và cơng thức nào, điều kiện tốt cho sinh trưởng phát triển rừng đại lượng lâm nghiệp khác Để giải vấn đề Stata có hệ thống lệnh so sánh Các lệnh giúp nhà khoa học lâm nghiệp thực tiêu chuẩn thống kê cách nhanh chóng, đơn giản Stata cung cấp tiêu chuẩn tham số phi tham số tiêu chuẩn U Mann Whitney, tiêu chuẩn Kruskal and Wallis, tiêu chuẩn Wilcoxon Từ việc so sánh mẫu đến nhiều mẫu, từ mẫu độc lập đến mẫu liên hệ Có thể tóm tắt lệnh so sánh từ mẫu trở lên Stata sau: So sánh đại lượng lâm nghiệp Các mẫu độc lập mẫu k mẫu (k >2) Tiêu chuẩn tham số Các mẫu liên hệ mẫu k mẫu (k >2) Tiêu chuẩn phi tham số Sơ đồ 01 Sơ đồ lệnh so sánh mẫu trở lên Stata Bên cạnh lệnh so sánh, Stata cịn mạnh phân tích phương sai để tìm hiểu ảnh hưởng hay nhiều nhân tố như: địa hình, độ dốc, độ ẩm, ánh sáng, phân bón… tới sinh trưởng phát triển đối tượng nghiên cứu lâm nghiệp rừng, nấm, độ cứng dẻo ván ép… Stata cung cấp lệnh để phân tích phương sai nhân tố, hai nhân tố không lặp nhân tố có lặp lại Những lệnh giúp người phân tích trả lời câu hỏi: nhân tố có thực ảnh 130 hưởng tới số liệu đối tượng nghiên cứu hay khơng? Ví tuổi có ảnh hưởng tới hệ số biến đổi bon hay khơng? Ánh sáng có ảnh hưởng tới sinh trưởng vườn ươm hay khơng? Địa hình có ảnh hưởng tới sinh trưởng rừng hay khơng? Nồng độ thuốc có ảnh hưởng tới sinh trưởng nấm hay khơng? Hay nhiệt độ ép có ảnh hưởng tới độ bền ván ép hay không? Tuy nhiên, hạn chế nhỏ Stata phân tích phương sai khơng tích hợp tiêu TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ LÂM NGHIỆP SỐ (KỲ I) - 2013 Ứng dụng công nghệ thông tin chuẩn Duncan Spss Duncan tiêu chuẩn sử dụng để so sánh kiểm tra khác biệt nhóm (David, 2008) Dựa vào việc phân nhóm tiêu chuẩn Duncan, nhà nghiên cứu tìm công thức tốt cho đối tượng nghiên cứu Ví độ tàn che tốt cho sinh trưởng địa, nồng độ thuốc tốt để kích thích rễ, nhiệt độ tốt cho ván ép… Dưới ví dụ quy trình sử dụng Stata để phân tích phương sai nhân tố Quy trình phân tích phương sai nhân tố Statistics\ Linear models and related \ ANOVA/ MANOVA\ Analysis of variance and convariance Trong hộp thoại anova - Analysis of variance and covariance khai báo d0 vào Dependent variable, khai báo cong_thuc vào Model Ok Dưới ví dụ kết phân tích phương sai để tìm hiểu ảnh hưởng độ tàn che tới sinh trưởng đường kính lồi Gội nếp giai đoạn vườn ươm giai đoạn 3-4 tháng tuổi (Bùi Mạnh Hưng, 2011) Bảng 01 Kết phân tích phương sai nhân tố Source Between groups Within groups Total Analysis of Variance SS df MS F 2.07758343 4.57257639 494 692527809 009256228 6.65015982 497 013380603 Bartlett's test for equal variances: chi2(3) = Prob > F 74.82 22.3961 0.0000 Prob>chi2 = 0.000 Nguồn: Bùi Mạnh Hưng (2011), ĐH quốc gia Úc 3.3 Khả phân tích cung cấp thơng tin mối quan hệ đại lượng lâm nghiệp Trong phân tích số liệu lâm nghiệp, thường xuyên phải phân tích mối quan hệ đại lượng lâm nghiệp mối tương quan hệ số chuyển đổi bon (BEF) với đường kính, chiều cao cây, phân tích mối quan hệ chiều cao lửa với độ ẩm vật rơi rụng… Dựa vào kết việc phân tích giúp nhà nghiên cứu dễ dàng xác định đại lượng khó đo đếm, khó tính tốn thể tích cây, hệ số chuyển đổi sinh khối, hình số… dựa vào đại lượng dễ đo đếm, dễ tính tốn đường kính ngang ngực, chiều cao….Khi phân tích mối quan hệ, thường có số câu hỏi cần phải trả lời là: mối quan hệ tuân theo dạng hàm nào? mức độ quan hệ đại lượng chặt hay không chặt? Để giải vấn đề Stata có trang bị lệnh phân tích tương quan hồi quy; phân tích tuyến tính phi tuyến Phân tích tuyến tính sử dụng lệnh “Linear regression” sau đưa biến độc lập (X) phụ thuộc (Y) vào Các thông tin tương quan như: hệ số tương quan, tham số phương trình, ước lượng kiểm tra tồn tham số tổng thể bảng đây: Lệnh “Nonlinear regression least square” tạm dịch “Tương quan phi tuyến phương pháp bình phương nhỏ nhất” giúp nhà phân tích thăm dị nhanh mối quan hệ đại lượng theo số loại hàm cơ khác nhau: Linear, logistic, Gompertz hàm Exponential Và cho biết TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ LÂM NGHIỆP SỐ (KỲ I) - 2013 131 Ứng dụng công nghệ thông tin mức độ liên hệ đại lượng theo dạng hàm Một chức khác lệnh giúp phân tích theo dạng hàm phi tuyến Thơng tin cần có biến, tham số giá trị bắt đầu tham số (David, 2011) Stata tự động tính tốn dừng lại tổng bình phương sai số dư bé hệ số xác định lớn (David, 2011) IV KẾT LUẬN Stata phần mềm mạnh, ứng dụng để phân tích số liệu nói chung số liệu nghiên cứu lâm nghiệp nói riêng Stata giúp giải đáp nhiều vấn đề khó khăn đặt phân tích số liệu lâm nghiệp Stata giúp nhà nghiên cứu tính toán nhận giá trị đặc trưng mẫu cách nhanh chóng, từ đặc trưng vị trí tới đặc trưng hình dạng phân bố Stata cung cấp cho người sử dụng lệnh ước lượng số trung bình tổng thể, lệnh so sánh mẫu Từ so sánh mẫu đến nhiều mẫu, từ so sánh mẫu độc lập đến mẫu liên hệ Từ nhà nghiên cứu kết luận mẫu hay khơng Hơn nữa, Stata cịn có khả ứng dụng tốt phân tích phương sai hai nhân tố, phân tích tương quan tuyến tính phi tuyến Vì vậy, chương trình đem lại sở khoa học cần thiết cho kết luận, đề xuất để góp phần vào việc quản lý phát triển tài nguyên rừng bền vững tương lai TÀI LIỆU THAM KHẢO Bùi Mạnh Hưng (2011), Phân tích ảnh hưởng nhân tố sinh thái tới sinh trưởng hai loài địa miền Bắc Việt Nam, Luận văn thạc sỹ, ĐH Quốc gia Úc, Canberra, Úc Bùi Mạnh Hưng (2012), Nghiên cứu khai thác ứng dụng phần mềm Stata phân tích số liệu nghiên cứu Lâm nghiệp, Đề tài cấp sở năm 2012, ĐH Lâm nghiệp, Hà nội, Việt Nam David, L.(2008), Online Statistics: An Interactive Multimedia Course of Study, University of HoustonDowntown, Texas, USA Available from: (Accessed 14 November, 2011) David G., 2011, Nonlinear, (Accessed 15 November, 2011) Nguyễn Hải Tuất Nguyễn Trọng Bình (2005), Khai thác sử dụng Spss để xử lý số liệu nghiên cứu lâm nghiệp, Nhà xuất Nông nghiệp, Hà nội Wiki, 2012 , xem ngày 20/11/2012 STATA – A NEW SOLUTION FOR FORESTRY DATA ANALYSIS Bùi Manh Hung, Nguyen Thi Bich Phuong SUMMARY Stata is a statistical software which has developed from 1985 Stata has shown superior strength in data analysis in many fields, especially in forestry Stata has been applied in many fields of forest science, such as forest inventory, silviculture, soil science, natural resources management, forest product processing and so on It can analyze from the univariate to multi-variable and from qualitative to quantitative analysis Stata can quickly provide information about the sample and the frequency distribution through descriptive statistics menu, graph and table menus Stata also shows the strength of analysis to compare observed samples and find the best conditions for the forestry variables The program can analyze by using parametric and non parametric tests, from two samples to multiple samples, from independent samples to paired samples Moreover, Stata also has very powerful ability about analysis of variance with one or two factors One of functional advantages of Stata is to analyze regression It can analyze linear regression from one variable to multi-variable Stata also can carry out non-linear correlation analysis for many difficult function types, such as Gompertz, Schumacher, Koller, Verhulst-Robertso and so on Therefore, the use of Stata to analyze research data in forestry is really necessary and to be a concerned issue Key words: ANOVA, comparison, data analysis, descriptive statistics, regression Stata, statistics Người phản biện: TS Vũ Khắc Bảy Ngày nhận bài: 10/6/2013 Ngày phản biện: 28/6/2013 Ngày định đăng: 20/9/2013 132 View publication stats TẠP CHÍ KHOA HỌC VÀ CÔNG NGHỆ LÂM NGHIỆP SỐ (KỲ I) - 2013 ... Thfng Ifng dung phan mdm Moodle xdy t:.ong bai giang tryc tuy6n tai truong Dai hoc Ldm nghiOp Phung :: Biri I4nh Hung, Nguy6n Thi Bfch Phuqng Stata rn r ;ho phin tich so lieu nghien cuu l6m nghicp... , xem ngày 20/11/2012 STATA – A NEW SOLUTION FOR FORESTRY DATA ANALYSIS Bùi Manh Hung, Nguyen Thi Bich Phuong SUMMARY Stata is a statistical software which has developed from 1985 Stata. .. trí tới đặc trưng hình dạng phân bố Stata cung cấp cho người sử dụng lệnh ước lượng số trung bình tổng thể, lệnh so sánh mẫu Từ so sánh mẫu đến nhiều mẫu, từ so sánh mẫu độc lập đến mẫu liên hệ