Giới thiệu dự án

Ngành công nghiệp chế biến, chế tạo (CBCT) giữ vai trò là động lực then chốt thúc đẩy tăng trưởng kinh tế Việt Nam. Theo số liệu của Tổng cục Thống kê, ngành này chiếm tới 49,9% tổng số lao động trong khu vực doanh nghiệp và đóng góp hơn 85% tổng kim ngạch xuất khẩu quốc gia (đạt 281,5 tỷ USD năm 2020). Tuy nhiên, giai đoạn 2019–2021 đánh dấu bước ngoặt phức tạp khi đại dịch Covid-19 bùng phát, gây ra cú sốc kép: đứt gãy chuỗi cung ứng nguyên vật liệu đầu vào và suy giảm tổng cầu tiêu dùng toàn cầu.

Vấn đề cốt lõi đặt ra là các doanh nghiệp CBCT niêm yết tại Việt Nam chưa có công cụ lượng hóa chính xác mức độ lãng phí nguồn lực đầu vào và khoảng cách so với biên giới hạn sản xuất tối ưu (Production Possibility Frontier - PPF). Việc thiếu các phân tích định lượng đa chiều khiến các quyết định phân bổ chi phí quản lý, lao động và cấu trúc tài chính trong thời kỳ khủng hoảng thiếu cơ sở khoa học vững chắc.

Mục tiêu cụ thể của dự án nghiên cứu gồm:

  1. Hệ thống hóa cơ sở lý thuyết về hiệu quả kỹ thuật (Technical Efficiency - TE) và các phương pháp tiếp cận biên sản xuất hiện đại.
  2. Xây dựng tập dữ liệu bảng từ 277 doanh nghiệp CBCT niêm yết trên các sàn HOSE, HNX và UPCoM thuộc hệ thống Vietstock Finance trong giai đoạn 2019–2021.
  3. Ứng dụng đồng thời mô hình phi tham số DEA (Data Envelopment Analysis) và mô hình tham số SFA (Stochastic Frontier Analysis) để tính toán điểm số TE.
  4. Ước lượng mô hình kinh tế lượng hồi quy tác động cố định (FEM) nhằm xác định mức độ ảnh hưởng của 13 biến vi mô và vĩ mô đến hiệu quả kỹ thuật, từ đó đề xuất hệ thống giải pháp nâng cao hiệu quả vận hành.

Phạm vi nghiên cứu tập trung vào 277 công ty cổ phần thuộc 19 nhóm ngành cấp II của ngành CBCT tại 3 miền Bắc, Trung, Nam trong giai đoạn 2019–2021. Giới hạn nghiên cứu là dữ liệu dựa trên báo cáo tài chính đã kiểm toán, chưa phản ánh đầy đủ nhóm doanh nghiệp vừa và nhỏ (SMEs) chưa niêm yết.


Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Các phương pháp đánh giá hiệu quả doanh nghiệp truyền thống tại Việt Nam chủ yếu dựa trên các chỉ số tài chính đơn lẻ hoặc hồi quy OLS thuần túy. Những cách tiếp cận này bộc lộ nhiều hạn chế khi phân tích hiệu quả kỹ thuật đa biến.

Tiêu chí so sánh Phân tích tỷ số tài chính (Financial Ratios) Hồi quy OLS truyền thống Phân tích bao dữ liệu (DEA) Phân tích biên ngẫu nhiên (SFA)
Bản chất mô hình Đơn biến, tính toán tĩnh Tham số, hồi quy trung bình Phi tham số, quy hoạch tuyến tính Tham số, kinh tế lượng nâng cao
Xử lý đa đầu vào / đầu ra Kém (chỉ so sánh từng cặp) Kém (chỉ 1 biến phụ thuộc) Rất mạnh (nhiều input/output) Tốt (mô hình hóa hàm sản xuất)
Tách biệt sai số ngẫu nhiên Không có Không tách được phi hiệu quả Không (gộp sai số vào phi hiệu quả) Tách biệt nhiễu $v_i$ và phi hiệu quả $u_i$
Yêu cầu giả định hàm Không Bắt buộc giả định tuyến tính Không cần giả định dạng hàm Giả định hàm Cobb-Douglas/Translog

Yêu cầu phân tích hệ thống được chuẩn hóa theo khung MoSCoW:

  • Must have: Thu thập đầy đủ 5 biến đầu vào ($\ln LD$, $\ln NV$, $\ln QL$, $\ln TS$, $\ln BH$) và 1 biến đầu ra ($\ln DT$); kiểm định đa cộng tuyến (VIF) và phương sai sai số thay đổi (White test).
  • Should have: Ước lượng đồng thời DEA hiệu suất thay đổi theo quy mô (VRS) và SFA ước lượng hợp lý cực đại (MLE).
  • Could have: So sánh mức độ khác biệt TE theo vùng miền (Bắc, Trung, Nam) và cấu trúc sở hữu.
  • Won't have: Mô hình hóa các biến động phi tài chính theo thời gian thực (real-time stream processing).

Thiết kế hệ thống và mô hình toán học

Hệ thống xử lý và phân tích dữ liệu được xây dựng trên nền tảng công nghệ chuyên dụng:

  • Công cụ định lượng: Stata MP phiên bản 17.0, DEAP Version 2.1 (Coelli), Python 3.10 (thư viện SciPy, Pandas, PyDEA).
  • Hệ cơ sở dữ liệu: Bộ dữ liệu bảng dạng cấu trúc (Panel Data) gồm 831 quan sát ($277 \text{ DN} \times 3 \text{ năm}$).

Mô hình hàm sản xuất biên ngẫu nhiên (SFA)

Mô hình SFA dạng Cobb-Douglas được thiết lập:

$$\ln(DT_{it}) = \beta_0 + \beta_1 \ln(LD_{it}) + \beta_2 \ln(NV_{it}) + \beta_3 \ln(QL_{it}) + \beta_4 \ln(TS_{it}) + \beta_5 \ln(BH_{it}) + v_{it} - u_{it}$$

Trong đó:

  • $DT_{it}$: Doanh thu thuần của doanh nghiệp $i$ tại năm $t$.
  • $LD_{it}, NV_{it}, QL_{it}, TS_{it}, BH_{it}$: Lần lượt là Chi phí lao động, Vốn chủ sở hữu, Chi phí quản lý, Tổng tài sản, Chi phí bán hàng.
  • $v_{it} \sim \text{iid } N(0, \sigma_v^2)$: Nhiễu thống kê ngẫu nhiên.
  • $u_{it} \sim \text{iid } N^+(0, \sigma_u^2)$: Thành phần phi hiệu quả kỹ thuật không âm.

Mô hình quy hoạch tuyến tính DEA (Dạng VRS - Input-oriented)

$$\min_{\theta, \lambda} \theta$$ Ràng buộc: $$\sum_{j=1}^{n} \lambda_j x_{ij} \le \theta x_{io}, \quad \forall i=1,\dots,m$$ $$\sum_{j=1}^{n} \lambda_j y_{rj} \ge y_{ro}, \quad \forall r=1,\dots,s$$ $$\sum_{j=1}^{n} \lambda_j = 1; \quad \lambda_j \ge 0, \quad \forall j=1,\dots,n$$

Mô hình hồi quy các yếu tố ảnh hưởng (Second-Stage Econometric Model)

$$TE_{it} = \gamma_0 + \gamma_1 \text{tuoi}{it} + \gamma_2 \text{vungmien}i + \gamma_3 \text{gioitinh}i + \gamma_4 \text{trinhdo}i + \gamma_5 \text{cautrucvon}{it} + \gamma_6 \ln(ROS{it}) + \gamma_7 \ln(ROE{it}) + \gamma_8 \ln(ROA{it}) + \gamma_9 \text{tysuatQL}{it} + \gamma{10} \text{tysuatBH}{it} + \gamma{11} \text{sanCK}i + \gamma{12} \text{TdCovid}t + \gamma{13} \text{mucphobien}i + \varepsilon{it}$$

Methodology

Quy trình nghiên cứu áp dụng mô hình thác nước kết hợp kiểm định vòng lặp:

  1. Giai đoạn 1 (Thu thập & Tiền xử lý): Trích xuất số liệu từ 277 báo cáo tài chính kiểm toán, xử lý giá trị khuyết và logarit hóa các biến quy mô.
  2. Giai đoạn 2 (Xác định kích thước mẫu): Áp dụng công thức Green (1991): $N \ge 50 + 8p$. Với $p = 13$ biến độc lập, kích thước mẫu tối thiểu là $50 + 8 \times 13 = 154$. Mẫu thực tế $N = 277$ doanh nghiệp đạt chuẩn tin cậy thống kê cao.
  3. Giai đoạn 3 (Ước lượng tham số & Phi tham số): Chạy giải thuật Simplex quy hoạch tuyến tính trên DEAP 2.1 và ước lượng MLE hàm SFA trên Stata 17.0.
  4. Giai đoạn 4 (Kiểm định khuyết tật mô hình): Thực hiện kiểm định Breusch-Pagan, White Test, Hausman Test để chọn giữa Pooled OLS, FEM và REM.

Implementation và kết quả

Development process

Quá trình phân tích dữ liệu được mô đun hóa thông qua các đoạn mã thực thi kinh tế lượng chuẩn mực.

* ====================================================================
* STATA SCRIPT: UOC LUONG HAM BIEN NGAU NHIEN SFA VA HOI QUY FEM
* Author: Nguyen Minh Tuan (2022) - Thesis Data Analytics
* ====================================================================

* 1. Khai bao du lieu bang (Panel data setup)
iis id_dn
tis nam

* 2. Uoc luong mo hinh bien san xuat SFA bang Maximum Likelihood
frontier lndt lnld lnnv lnql lnts lnbh, distribution(half) emean

* 3. Trich xuat chi so hieu qua ky thuat TE (Technical Efficiency)
predict te_sfa, te

* 4. Kiem dinh lua chon mo hinh: Fixed Effects vs Random Effects
xtreg te_sfa tuoi vungmien gioitinh trinhdo cautrucvon lnros lnroe lnroa tysuatql tysuatbh sanck tdcovid mucphobien, fe
estimates store fixed_effects

xtreg te_sfa tuoi vungmien gioitinh trinhdo cautrucvon lnros lnroe lnroa tysuatql tysuatbh sanck tdcovid mucphobien, re
estimates store random_effects

hausman fixed_effects random_effects

* 5. Kiem dinh khuyet tat va hieu chinh sai so chuan vung (Robust standard errors)
xtreg te_sfa tuoi cautrucvon lnros lnroe lnroa tysuatql tysuatbh tdcovid, fe vce(robust)

Testing và validation

Kết quả chẩn đoán mô hình kinh tế lượng đảm bảo tính vững của các ước lượng:

  • Kiểm định đa cộng tuyến: Giá trị Hệ số Phóng đại Phương sai (VIF) trung bình đạt $2,14$ (tất cả các biến đều có $\text{VIF} < 5,0$), loại trừ nguy cơ đa cộng tuyến nghiêm trọng.
  • Kiểm định Breusch-Pagan Lagrangian Multiplier: Cho giá trị $p\text{-value} < 0,001$, bác bỏ mô hình Pooled OLS thông thường.
  • Kiểm định White & Modified Wald: Xác nhận sự hiện diện của hiện tượng phương sai sai số thay đổi ($p < 0,01$). Nghiên cứu đã xử lý triệt để bằng phương pháp ước lượng FEM hiệu chỉnh sai số chuẩn vững (Robust Standard Errors).
                            PHÂN BỐ DOANH NGHIỆP THEO LĨNH VỰC

Kết quả đạt được

Đặc điểm mẫu nghiên cứu phản ánh sát cấu trúc ngành công nghiệp Việt Nam:

  • Độ tuổi hoạt động: Trung bình đạt 14,26 năm ($\text{Min} = 3$, $\text{Max} = 33$, độ lệch chuẩn $5,192$), cho thấy các doanh nghiệp có bề dày kinh nghiệm và tính ổn định nhất định.
  • Phân bố không gian: Khu vực miền Nam chiếm ưu thế với 133 doanh nghiệp (47,50%), miền Bắc và miền Trung chiếm 147 doanh nghiệp (52,50%).
  • Ước lượng hàm SFA: Các biến Chi phí lao động ($\beta_1 = 0,284$, $p < 0,01$), Chi phí quản lý ($\beta_3 = 0,192$, $p < 0,05$) và Tổng tài sản ($\beta_4 = 0,315$, $p < 0,01$) có tác động tích cực đáng kể nhất tới doanh thu thuần.
  • Tác động của Covid-19: Biến giả đại dịch $\text{TdCovid}$ mang hệ số âm có ý nghĩa thống kê ($\gamma_{12} = -0,042$, $p < 0,01$), phản ánh hiệu quả kỹ thuật trung bình của toàn ngành sụt giảm $4,2%$ trong giai đoạn phong tỏa so với mức trước dịch năm 2019.
  • Cấu trúc vốn và Tỷ suất chi phí: Cấu trúc đòn bẩy tài chính cao và tỷ suất bán hàng ($\text{tysuatBH}$) tương quan nghịch với TE ($\gamma_{10} = -0,118$, $p < 0,05$), chứng minh việc thâm dụng chi phí phân phối trung gian trong dịch làm xói mòn hiệu quả vận hành.

Đổi mới và đóng góp

Nghiên cứu mang lại những đóng góp học thuật và thực tiễn mang tính hệ thống so với các công trình tiền nhiệm:

  1. Phương pháp tiếp cận đối sánh (Methodological Triangulation): Khác với công trình của Lê Ngọc Danh & Ngô Thị Thanh Trúc (2021) chỉ sử dụng đơn lẻ DEA-Tobit, hoặc Hoàng Thanh Hiền & Huỳnh Thị Diệu Linh (2022) chỉ dùng SFA cho khối doanh nghiệp nhà nước, nghiên cứu này kết hợp đồng thời DEA và SFA trên cùng một mẫu doanh nghiệp niêm yết, tạo cơ chế đối chứng chéo và loại bỏ sai số đo lường.
  2. Cập nhật dữ liệu thời gian thực trong khủng hoảng: Đóng đính khoảng trống nghiên cứu của giai đoạn biến động 2019–2021, cung cấp bằng chứng thực nghiệm về tính dễ tổn thương của chuỗi cung ứng doanh nghiệp sản xuất trước cú sốc ngoại sinh.
  3. Bóc tách chi tiết cấu trúc chi phí quản trị: Chứng minh định lượng rằng tối ưu hóa tỷ suất chi phí quản lý ($\text{tysuatQL}$) và tái cấu trúc đòn bẩy nợ có thể cải thiện hiệu quả kỹ thuật thêm từ $8,5%$ đến $12,3%$ trong điều kiện phục hồi sản xuất.

Ứng dụng thực tế và triển khai

Nghiên cứu mở ra khung hành động cụ thể cho việc triển khai tối ưu hóa tại các nhà máy và cơ quan quản lý:

                LỘ TRÌNH 4 BƯỚC TỐI ƯU HÓA HIỆU QUẢ KỸ THUẬT (TE)
  • Ứng dụng cho doanh nghiệp sản xuất: Sử dụng hệ số biên từ hàm SFA làm trọng số để phân bổ ngân sách. Doanh nghiệp ngành thực phẩm và hóa chất cần tập trung tự động hóa để cắt giảm tỷ suất chi phí quản lý vận hành từ mức $15%$ doanh thu xuống dưới $9%$.
  • Ứng dụng công nghệ: Triển khai hệ thống hoạch định tài nguyên doanh nghiệp (ERP SAP S/4HANA hoặc Oracle NetSuite) tích hợp mô-đun phân tích hiệu quả sản xuất theo thời gian thực.
  • Chính sách vĩ mô: Cung cấp cơ sở định lượng để Bộ Công Thương và các sở ban ngành địa phương thiết kế các gói tín dụng hỗ trợ lãi suất cho nhóm doanh nghiệp sản xuất có chỉ số TE dưới mức trung bình ngành nhưng có năng lực phục hồi tốt.

Hạn chế và hướng phát triển

  • Hạn chế về dữ liệu: Nghiên cứu mới khai thác số liệu thứ cấp của 277 doanh nghiệp niêm yết trên thị trường chứng khoán, chưa bao quát được toàn diện nhóm doanh nghiệp vừa và nhỏ (SMEs) chiếm tỷ trọng lớn trong nền kinh tế.
  • Giới hạn biến số: Mô hình chưa đưa vào các biến phi tài chính như tỷ lệ số hóa dây chuyền, mức độ phát thải carbon (ESG) hay tỷ lệ đầu tư R&D do hạn chế công bố thông tin.
  • Hướng mở rộng:
    1. Ứng dụng mô hình DEA-Malmquist động để phân tích sự thay đổi năng suất các yếu tố tổng hợp (TFP) theo chuỗi thời gian dài 5–10 năm.
    2. Kết hợp thuật toán học máy (Random Forest Regression, XGBoost) với SFA để dự báo sớm nguy cơ suy giảm hiệu quả kỹ thuật trong điều kiện đứt gãy chuỗi cung ứng.

Đối tượng hưởng lợi

                 LỢI ÍCH ĐỊNH LƯỢNG CHO CÁC BÊN LIÊN QUAN

Câu hỏi thường gặp

1. Yêu cầu kỹ thuật phần mềm để tái lập mô hình nghiên cứu này là gì?

Để tái lập nghiên cứu, hệ thống cần cài đặt:

  • Stata MP 16.0 trở lên (khuyến nghị Stata 17.0) hỗ trợ gói lệnh frontierxtreg.
  • Phần mềm DEAP 2.1 do Tim Coelli phát triển chạy trên môi trường Windows/DOS.
  • R (phiên bản 4.2+) với các package frontier, Benchmarking, plm để kiểm chứng độc lập.

2. Sự khác biệt căn bản giữa DEA và SFA trong tính toán điểm số TE là gì?

DEA là phương pháp phi tham số, tạo ra một đường biên bao dữ liệu bằng quy hoạch tuyến tính và quy toàn bộ khoảng cách từ quan sát tới đường biên là phi hiệu quả. SFA là phương pháp tham số kinh tế lượng, phân tách sai số thành hai phần: nhiễu ngẫu nhiên $v_i$ (do thời tiết, biến động vĩ mô) và thành phần phi hiệu quả kỹ thuật $u_i$. SFA tránh được hiện tượng thổi phồng mức phi hiệu quả của DEA khi dữ liệu có nhiều nhiễu.

3. Làm thế nào để xử lý hiện tượng phương sai sai số thay đổi trong mô hình?

Trong nghiên cứu này, kiểm định White và Modified Wald phát hiện phương sai sai số thay đổi ở mức ý nghĩa $1%$. Giải pháp xử lý là áp dụng phương pháp ước lượng tác động cố định với sai số chuẩn vững:

xtreg te_sfa `bien_doc_lap`, fe vce(robust)

Tùy chọn vce(robust) điều chỉnh ma trận hiệp phương sai của các ước lượng, đảm bảo kiểm định $t$ và khoảng tin cậy không bị chệch.

4. Chi phí quản lý và chi phí bán hàng tác động như thế nào đến điểm số TE?

Kết quả hồi quy thực nghiệm chỉ ra rằng chi phí quản lý đóng vai trò đầu vào tạo ra doanh thu ($\beta_3 = 0,192$), nhưng tỷ suất chi phí bán hàng trên doanh thu lại tương quan nghịch với hiệu quả kỹ thuật ($\gamma_{10} = -0,118$). Điều này chứng minh việc gia tăng chi phí quản trị đúng mức giúp chuẩn hóa quy trình, nhưng lạm dụng chi phí trung gian bán hàng trong khủng hoảng sẽ làm sụt giảm TE.

5. Khung thời gian và lộ trình tối ưu hóa TE cho một nhà máy chế biến mất bao lâu?

Lộ trình tối ưu hóa điển hình kéo dài 9 đến 12 tháng, chia làm 3 giai đoạn:

  1. Giai đoạn 1 (Tháng 1-3): Thu thập dữ liệu định mức, phân tích DEA/SFA nội bộ để định vị điểm lãng phí.
  2. Giai đoạn 2 (Tháng 4-8): Tái cấu trúc quy trình, số hóa hệ thống quản lý chi phí qua ERP/MES.
  3. Giai đoạn 3 (Tháng 9-12): Đánh giá lại chỉ số TE định kỳ và cân bằng lại cấu trúc vốn.

Kết luận

Công trình nghiên cứu đã giải quyết thành công bài toán lượng hóa hiệu quả kỹ thuật của 277 doanh nghiệp ngành công nghiệp chế biến, chế tạo niêm yết tại Việt Nam trong giai đoạn biến động Covid-19 (2019–2021). Bằng việc kết hợp hài hòa giữa phương pháp phi tham số DEA và mô hình kinh tế lượng tham số SFA, nghiên cứu cung cấp bức tranh toàn cảnh về thực trạng phân bổ nguồn lực, chỉ rõ các nhân tố trọng yếu tác động đến năng suất kỹ thuật như chi phí quản lý, cấu trúc vốn, vị trí địa lý và tác động ngoại sinh của dịch bệnh.

Kết quả nghiên cứu khẳng định rằng việc nâng cao hiệu quả kỹ thuật không chỉ phụ thuộc vào quy mô tổng tài sản hay vốn chủ sở hữu, mà cốt lõi nằm ở năng lực tinh gọn bộ máy quản trị, tối ưu hóa tỷ suất chi phí phân phối và chủ động ứng dụng công nghệ số trong vận hành. Đây là tài liệu tham khảo giá trị cho các nhà hoạch định chính sách, các nhà nghiên cứu kinh tế lượng ứng dụng và các nhà quản trị doanh nghiệp trên hành trình nâng cao năng lực cạnh tranh trong kỷ nguyên số.