经济学之粮食生产模型

一、粮食生产模型

1、选择变量和模型的关系形式

研究目的：研究我国近十年粮食产量与主要影响因素之间的定量关系。

（1）确定模型应包含的变量

被解释变量：粮食总产量（万吨），y

解释变量：X1农用化肥施用量（万公斤）

X2粮食播种面积（千公顷）

X3成灾面积（千公顷）

X4农业机械劳动力（万千瓦）

X5农业劳力（万人）

变量选择依据：国家粮食政策对粮食生产的积极扶持作用，在样本区间内具有一致性；

我国粮食统一平价收购，，农民面临的是“无限的需求”。

因此，影响粮食产量的主要因素是：资本和劳动力

农业生产的特点决定了资本主要是：土地和化肥

考虑到我国的传统农业耕作方式有一定的变革，一方面从事非农产业的农民增加；另一方面，农业机械化水平有所提高，表示劳动力的变量为：农机动力和农业劳动力。

（2）确定模型的数学形式和参数范围

被解释变量和解释变量之间的关系是投入产出关系，因此选用C-D生产函数，设定模型为：

lnY=a0+a1lnX1+a2lnX2++++u

解释变量的个数从2~5个，可根据实际情况选择最有结果；

偏回归系数为相应变量的产出弹性，理论值应介于0~1之间；

截距a0为常熟虚拟变量，包含政策、气候等难以量化的影响因素；

u为随即扰动项，描述模型之外的因素对模型的干扰。

2、采集样本数据

本模型使用时间序列数据，数据来源于《中国农村统计年鉴》（1984——1996）（见表2.10.1）

其中：由于统计数据没有粮食生产的专用数据，农用化肥施用量（万公斤）、成灾面积（千公顷）、农业机械总动力（万千瓦）、农业劳力（万人）均为农业总体数据，尤其是农业劳力（万人）为农林牧副渔也总劳动力，只有几个年份有种植业、农业劳力数据，占农林牧副渔

也劳力的比重为0.8~0.85之间。表2.10.1

obs Y

粮食总产

量（万吨）

X1

农用化肥施用

量（万公斤）

X2

粮食播种面积

（千公顷）

X3

成灾面积（千

公顷）

X4

农业机械劳动力

（万千瓦）

X5

农业劳力（万

人）

1983387281659.811404716209.31802231645.1 1984407311739.8112884152641949731685 1985379111775.810884522705.32091330351.5 1986391511930.6110933236562295030467 1987402981999.311126820393.72483630870 1988394082141.511012323944.72657531455.7 1989407552357.111220524448.72806732440.5 1990446242590.311346617819.32870833330.4 1991435292805.1112314278142938934186.3 1992442662930.211056025894.73030834037 1993456493151.9110509231333181733258.2 1994445103317.9109544313833380232690.3 1995466623593.7110060222673611832334.5 1996504543827.91125482123438546.932260.4将原始数据格式转换为文本文件格式，并保存为：f:\sq\lzn2101.txt（记住样本区间、各列对

应的变量名称）

38728 1659.8 114047 16209.3 18022 31645.1

40731 1739.8 112884 15264 19497 31685

37911 1775.8 108845 22705.3 20913 30351.5

39151 1930.6 110933 23656 22950 30467

40298 1999.3 111268 20393.7 24836 30870

39408 2141.5 110123 23944.7 26575 31455.7

40755 2357.1 112205 24448.7 28067 32440.5

44624 2590.3 113466 17819.3 28708 33330.4

43529 2805.1 112314 27814 29389 34186.3

44266 2930.2 110560 25894.7 30308 34037

45649 3151.9 110509 23133 31817 33258.2

44510 3317.9 109544 31383 33802 32690.3

46662 3593.7 110060 22267 36118 32334.5

50454 3827.9 112548 21234 38546.9 32260.4

进入TSP窗口，进行如下操作：

Create a 1983 1996

Read(o) F:\sq\ lzn2101.txt y x1 x2 x3 x4 x5

write(o) F:\sq\ lzn2101.txt y x1 x2 x3 x4 x5

save F:\sq\ lzn2101.txt

load F:\sq\ lzn2101.txt

3、参数估计结果及统计检验

（1）ls lny c lnx23 lnx5 lnx123 lnx423

把化肥施用量改为单位面积化肥施用量、农机劳力改为单位面积农机劳力，能够在一定程度上避免可能的多重共线性。用OLS法估计，得到如下回归方程：

Create a 1983 1996

Read(o) F:\sq\ lzn2101.txt y x1 x2 x3 x4 x5

write(o) F:\sq\ lzn2101.txt y x1 x2 x3 x4 x5

save F:\sq\ lzn2101.txt

load F:\sq\ lzn2101.txt

genr lny=log(y)

genr lnx23=log(x2-x3)（不受灾的面积）

genr lnx5=log(x5)

genr lnx123=log(x1/(x2-x3)) （单位面积化肥施用量）

genr lnx423=log(x4/(x2-x3)) （单位面积农机劳力）

smpl 1983 1995

K=4,n=13,查t\F\DW等统计量的临界值，进行统计检验；

根据经济意义和t检验结果，依次分别剔除变量LNX5、LNX423，用OLS法估计，得到如下回归方程：

（2）ls lny c lnx23 lnx5 lnx123