Mơ hình nghiên cứu trong mục này sử dụng số liệu mảng không gian, gồm 63 tỉnh /thành phố theo thời gian từ năm 2004 đến năm 2014. Số liệu của các biến số trong mơ hình được tổng hợp từ hai nguồn: số liệu điều tra Doanh nghiệp hàng năm của Tổng cục Thống kê và số liệu điều tra PCI (năng lực cạnh tranh cấp tỉnh).
Các kiểm định cho thấy tồn tại hiện tượng sai số không gian và trễ không gian (Phụ lục PL10). Do đó, nếu sử dụng các mơ hình số liệu gộp hoặc số liệu mảng sẽ bỏ qua sự tương quan theo không gian địa lý (kết quả ước lượng xem trong phần phụ lục PL9). Vì vậy chúng tơi tiến hành ước lượng các mơ hình số liệu mảng khơng gian, bao gồm 4 mơ hình: Mơ hình sai số khơng gian (SEM), mơ hình có trễ khơng gian (SAR), mơ hình Dublin khơng gian (SDM), mơ hình tự hồi quy khơng gian (SAC).
Các biến số được giải thích trong bảng 4.23.
Bảng 4.23 Giải thích các biến số trong mơ hình nghiên cứu
Tên biến Giải thích các biến Xu hướng ảnh hưởng
LI Là biến phụ thuộc, đại diện cho chuyển
dịch cơ cấu lao động giữa các ngành
trong mỗi tỉnh và được đo bằng chỉ số
Lilien theo cơng thức 2.2
5¢/3 Là biến độc lập thể hiện mức khác biệt
trong thu nhập bình quân của người lao
động giữa các ngành trong mỗi tỉnh. Biến
này được đo bởi độ lệch chuẩn của thu nhập trong các ngành cấp 1 chia cho mức thu nhập trung bình của lao động trong tỉnh.
Khác biệt thu nhập là một động
lực quan trọng thúc đẩy chuyển
dịch lao động. Mức khác biệt
càng cao thì người lao động
càng có động lực thuyên.
5¢/3 2 Là bình phương của KBTN. Biến này đưa vào nhằm kiểm
soát quan hệ phi tuyến giữa sự khác biệt này và CDCCLĐ.
Tên biến Giải thích các biến Xu hướng ảnh hưởng
tự nhiên của vốn bình quân theo lao động của các doanh nghiệp trên địa bàn tỉnh.
thể tạo thêm nhiều việc làm, thúc
đẩy lao động dịch chuyển.
ln-!ỳ Là biến độc lập đo bằng logarit cơ số tự
nhiên của tổng số lao động đang làm việc tại các doanh nghiệp trên địa bàn tỉnh.
Quy mô lao động càng lớn thì
thơng tin có xu hướng đến với
người lao động chậm hơn, do vậy
quy mơ lao động có thể cản trở
quá trình chuyển dịch lao động.
(/ ( Là biến độc lập thể hiện cho chất lượng
của thị trường lao động, được đo bằng
chỉ số thị trường lao động - là thành
phần con của chỉ số PCỊ
Biến này có xu hướng tác động cùng chiều với chuyển dịch lao
động và có ảnh hưởng lan tỏạ
G_LABOR Đo tốc độ tăng lao động hàng năm trong
khối doanh nghiệp tại mỗi địa phương.
Nó đại diện cho xu hướng của cầu lao
động từ phía các doanh nghiệp trên địa
bàn tỉnh.
Tốc độ tăng lao động hàng năm
có xu hướng cùng chiều với CDCCLĐ giữa các ngành.
Nguồn: Nghiên cứu của tác giả
Sau khi tính được giá trị của các biến số theo từng tỉnh từ số liệu Điều tra doanh nghiệp hàng năm, chúng tôi nối dữ liệu theo các năm để được một file số liệu mảng theo 63 tỉnh/ thành phố và theo thời gian 11 năm, từ năm 2004 đến năm 2014.
Sau đây là thống kê mơ tả một số biến trong mơ hình nghiên cứụ
Bảng 4.24 Thống kê mô tả các biến số trong mơ hình nghiên cứu
Variable Số quan sát Mean Std. Dev. Min Max
LI 693 1,879906 1,853291 0 12,27548 KBTN 693 0,3055882 0,2884572 0,0878145 4,879329 lnCAP 693 5,918101 0,7332821 4,270606 8,267736 lnSIZE 693 10,96801 1,104359 8,343316 14,70607 G_LABOR 693 0,0842096 0,1314927 -0,6115508 1,888314 DTLD 693 4,085245 2,164907 0 9,596658
Nguồn: Tính tốn của tác giả từ số liệu nghiên cứu với phần mềm STATA Thiết lập ma trận trọng số khơng gian
Các mơ hình hồi quy theo số liệu mảng không gian đều yêu cầu xây dựng một ma trận trọng số theo không gian. Một số phương pháp xây dựng ma trận này đã được
giới thiệu ở chương 2. Trong mơ hình nghiên cứu này, chúng tơi xây dựng ma trận trọng số không gian = ZgM[<1×<1 cho 63 tỉnh, đây là một ma trận vng cỡ 63x63, được chuẩn hóa theo hàng (tức là tổng các phần tử theo hàng của ma trận bằng 1).
Các phần tử của ma trận trọng số được xác định như sau:
gM = £
!(0 < 9M < 9) ∑ !(0 < 9M<1 < 9)
M 8ℎ 0 < 9M < 9, ≠ ¥
0 á 6ườ&$ ℎợ| ị& ạ
Trong đó, hàm chỉ báo: !Z0 < 9M < 9[ nhận giá trị bằng 1 nếu khoảng cách từ tỉnh i đến tỉnh j là 9M nằm trong khoảng (0,d), và nhận giá trị bằng 0 trong các trường hợp khác.
Khoảng cách giữa hai tỉnh i và j được tính bằng khoảng cách giữa hai điểm đại diện, đặt tại trung tâm hành chính của mỗi tỉnh theo cơng thức:
9M = G(?M− ?)+ (aM − a)
ở đó, (?, a) là tọa độ địa lý của điểm i, xác định bởi vĩ độ và kinh độ đo bằng số thập phân.
Ngưỡng khoảng cách d được lựa chọn theo đặc điểm riêng của từng nghiên cứu, trong nghiên cứu này ngưỡng khoảng cách d được chọn bằng 5 đơn vị để phù hợp với đặc điểm địa lý của Việt Nam (Nguyen Khac Minh, Pham Anh Tuan, Nguyen Viet Hung, 2015).