СТАТИСТИЧЕСКОЕ ИССЛЕДОВАНИЕ ДИНАМИКИ ПРОИЗВОДСТВА И РЕАЛИЗАЦИИ САХАРНОЙ СВЕКЛЫ, МЕТОДЫ АНАЛИЗА ВЛИЯНИЯ РАЗЛИЧНЫХ ФАКТОРОВ НА ЕЕ УРОЖАЙНОСТЬ - Студенческий научный форум

VIII Международная студенческая научная конференция Студенческий научный форум - 2016

СТАТИСТИЧЕСКОЕ ИССЛЕДОВАНИЕ ДИНАМИКИ ПРОИЗВОДСТВА И РЕАЛИЗАЦИИ САХАРНОЙ СВЕКЛЫ, МЕТОДЫ АНАЛИЗА ВЛИЯНИЯ РАЗЛИЧНЫХ ФАКТОРОВ НА ЕЕ УРОЖАЙНОСТЬ

Гайнетдинова Г.Р. 1
1Башкирский государственный аграрный университет
 Комментарии
Текст работы размещён без изображений и формул.
Полная версия работы доступна во вкладке "Файлы работы" в формате PDF
Производство, потребление, качество продуктов питания в совокупности составляют одну из фундаментальных основ жизнедеятельности человека, определяющих производственную независимость государства. Именно сахарная свекла определяет продовольственную независимость нашей страны [4].

Сахарная свекла - важнейший экспортный продукт. Чем больше зерновые ресурсы страны, тем крепче внешнеполитические позиции государства на мировой арене.

От того, насколько интенсивно в стране будет наращиваться производство сахарной свеклы, зависит реальность преодоления кризиса в АПК, улучшение продовольственного обеспечения населения [1].

Проведем анализ состава и структуры посевных площадей сахарной свеклы в динамике. Данные представлены в Приложение А.

По данным Приложение А видно, что в 2015 году по сравнению с 2007 годом произошло увеличение посевных площадей на 13,2 %, в том числе «Доминика» на 47,8% и «Льговская» на 17,2% , а площадь «Рамонская» снизилась на 14,9 %.

В продолжительной динамике изучено изменение посевной площади, валового сбора, объема реализации и урожайности продукции растениеводства по группам однородных культур (таблица 1).

Таблица 1 Динамика производства продукции растениеводства

Показатели

Годы

Темп роста,

базисный, %

 

2007 г.

2008 г.

2009 г.

2010 г.

2011 г.

2012 г.

2013 г.

2014 г.

2015 г.

 

Посевная площадь, га-всего в т.ч.

566,8

573,7

581,3

589,5

598,4

608,0

618,4

629,5

641,5

113,2

Доминика

196,4

206,2

216,6

227,4

238,8

250,7

263,2

276,4

290,2

147,8

Рамонская

257,6

252,5

247,4

242,5

237,6

232,9

228,2

223,7

219,2

85,1

Льговская

112,7

115,0

117,3

119,6

122,0

124,5

126,9

129,5

132,1

117,2

Валовой сбор, тыс. ц-всего, в т.ч.

73432,1

73776,9

66629,5

66817,4

68090,0

75888,7

76604,1

77684,3

81568,5

111,1

Доминика

22730,5

23639,7

17020,6

17701,4

19432,2

27655,1

28761,3

30199,4

31709,4

139,5

Рамонская

31567,8

30620,8

29702,2

28811,1

27946,8

27108,4

26295,1

25506,3

26781,6

84,8

Льговская

19133,7

19516,4

19906,7

20304,9

20710,9

21125,2

21547,7

21978,6

23077,6

120,6

Урожайность с 1 га, ц-всего, в т.ч.

129,6

128,6

114,6

113,3

113,8

124,8

123,9

123,4

127,2

98,1

Доминика

115,7

114,6

78,6

77,8

81,4

110,3

109,3

109,3

109,3

94,4

Рамонская

122,5

121,3

120,1

118,8

117,6

116,4

115,2

114,0

122,2

99,7

Льговская

169,8

169,7

169,7

169,8

169,8

169,7

169,8

169,7

174,7

102,9

Реализовано, тыс. ц,

60459,8

61668,9

62902,3

64160,4

65443,6

66752,4

68087,5

69449,2

70838,2

117,2

Уровень товарности, %

82,3

83,6

94,4

96,0

96,1

88,0

88,9

89,4

86,8

105,5

По данным таблицы 1 можно сделать вывод, что в 2015 году относительно 2007 года посевная площадь возросла на 13,2 %, урожайность с 1 га уменьшилась на 1,9 %.

Применив способ разложения составных относительных показателей на простые, изучено изменение уровня производства сахарной свеклы на 100 га пашни. Данные представлены в Приложение Б.

Ур= У × d1 × d2 (1)

Данная взаимосвязь прослеживается и при расчетах индексов данных показателей.

Jyp = Jy ×Jdl ×Jd2, (2)

где Jyp - индекс уровня производства сахарной свеклы на 100 га пашни;

Jyp = Ур последнего периода (3)

Ур первого периода

Проведен анализ изменения показателей уровней ряда динамики урожайности, рассчитав для этого абсолютные, относительные и средние показатели ряда динамики по формулам, приведенным в таблице 2.

Таблица 2 Способы расчета показателей ряда динамики

Показатели

Способы расчета

базисный

цепной

Абсолютный прирост, ц

б = Уi – У1

 

Средний абсолютный прирост, ц

б=

ц =

Темп роста, %

%

%

Средний темп роста, %

%

 

Темп прироста, %

= %

 

Средний темп прироста, %

   

Абсолютное содержание 1% прироста, ц

-

 

где У1 – начальный уровень ряда динамики,

Уп – конечный уровень ряда динамики,

n - число лет.

Показатели рассчитываются двумя способами: базисным и цепным. Расчеты показателей ряда динамики представим в виде таблицы 3.

Таблица 3 Показатели ряда динамики урожайности

Годы

Урож-ть,ц с 1 га

Абсолютный прирост

Темп роста, %

Темп прироста, %

Абсолютное содержание 1% прироста

Базис.

Цепн.

Базис.

Цепн.

Базис.

Цепн.

2007

129,6

             

2008

128,6

-2,4

-1,0

99,2

99,2

-0,8

-0,8

1,3

2009

114,6

-1,4

-14,0

88,4

89,1

-11,6

-10,9

1,1

2010

113,3

12,6

-1,3

87,4

98,9

-12,6

-1,1

1,1

2011

113,8

13,9

0,5

87,8

100,4

-12,2

0,4

1,1

2012

124,8

13,4

11

96,3

109,7

-3,7

9,7

1,2

2013

123,9

2,4

-0,9

95,6

99,3

-4,4

-0,7

1,2

2014

123,4

3,3

-0,5

102,7

99,6

2,7

-0,4

1,2

2015

127,2

3,8

3,8

103,1

103,1

3,1

3,1

1,3

В среднем

122,1

-0,3

-0,3

99,8

99,8

-0,2

-0,2

 

1. Средний абсолютный прирост:

ц,

ц,

К=n-1, где К – количество цепных абсолютных приростов,

n - количество уровней ряда.

2. Средний коэффициент роста:

=;

К=n-1, где К – количество цепных коэффициентов роста

  1. Средний темп роста:

  1. Средний темп прироста:

  1. Средний уровень ряда динамики:

Для выявления тенденции в рядах динамики используем следующие методы выравнивания:

1. Механическое выравнивание: а) укрупнения периодов; б) расчета скользящих средних.

2. Плавного уровня: а) среднего абсолютного прироста; б) среднего темпа роста.

3. Аналитическое выравнивание уровней ряда методом наименьших квадратов (МНК) по уравнению прямой.

Для исследования интенсивности изменения явления во времени рассчитываются показатели ряда динамики урожайности сахарной свеклы.

Проведен анализ тенденции ряда динамики.

Таблица 4 Выравнивание ряда динамики методом укрупнения периодов

Годы

Урожайность, ц с 1 га, У

По 3-х летиям

Периоды

Сумма

Средние

2007

129,6

-

-

-

2008

128,6

2007 -2009

372,8

124,3

2009

114,6

-

-

-

2010

113,3

-

-

-

2011

113,8

2010 -2012

351,9

117,3

2012

124,8

-

-

-

2013

123,9

-

-

-

2014

123,4

2013 - 2015

374,5

124,8

2015

127,2

-

-

-

Рисунок 1 Выравнивание ряда динамики методом укрупнения периодов

Методом укрупнения периодов за 2007 – 2015 годы выявлена нестабильная тенденция увеличения урожайности с 1 га.

Таблица 5 Выравнивание ряда динамики методом скользящей средней

Годы

Урожайность, ц с 1 га, У

По скользящим 3-х летиям

Периоды

Сумма

Средние

2007

129,6

-

-

 

2008

128,6

2007 -2009

372,8

124,3

2009

114,6

2008 -2010

356,5

118,8

2010

113,3

2009 -2011

341,7

113,9

2011

113,8

2010-2012

351,9

117,3

2012

124,8

2011 -2013

362,5

120,8

2013

123,9

2012 -2014

372,1

124,0

2014

123,4

2013 - 2015

374,5

124,8

2015

127,2

-

-

-

Рисунок 2 Выравнивание ряда динамики методом скользящей средней

Методом скользящей средней выявлена тенденция уменьшения урожайности с 1 га за 2007-2010 годы, с 2010 года наблюдается увеличение урожайности с 1 га в ООО «Нива»

Таблица 6 Выравнивание ряда динамики по среднему абсолютному приросту и среднему коэффициенту роста

Годы

Урожайность, ц с 1 га, Yi

Порядковый номер года t

Yt=129,6-0,3*(t-1)

Yt = 129,6*1,0511(t-1)

2007

129,6

1

129,6

129,6

2008

128,6

2

129,3

129,3

2009

114,6

3

129,0

129,1

2010

113,3

4

128,7

128,8

2011

113,8

5

128,4

128,6

2012

124,8

6

128,1

128,3

2013

123,9

7

127,8

128,1

2014

123,4

8

127,5

127,8

2015

127,2

9

127,2

127,5

Рисунок 3 Выравнивание ряда динамики по среднему абсолютному приросту и среднему коэффициенту роста

Таблица 6 и рисунок 3 показывают, что методом выравнивания ряда динамики по среднему абсолютному приросту с 2007 по 2015 годы выявлено незначительное изменение, причем за 2008-2011 годы наблюдается увеличение урожайности с 1 га.

Исходя из этого, для проявления тенденции динамики используем уравнение прямой

, (4)

где - - выравненное значение урожайности,

a, b – неизвестные параметры,

a – значение выравненной урожайности для центрального в динамическом ряду года, содержательной интерпретации не имеет,

b – ежегодный прирост (снижение) урожайности,

t – значения дат.

Для определения неизвестных параметров a и b в соответствии с требованием способа наименьших квадратов необходимо решить систему нормальных уравнений:

(5)

Система упрощается, если воспользоваться способом отсчета времени от условного начала (таблица 7).

Таблица 7 Аналитическое выравнивание урожайности по уравнению прямой методом наименьших квадратов

Годы

Урожайность, ц с 1 га,У

t

t^2

Уx t

Выравненный уровень урож-ти, ц с 1 га,

(У- )

(У- )^2

2007

129,6

-4

16

-518,4

122,2

7,4

55,1

2008

128,6

-3

9

-385,8

122,3

6,3

40,2

2009

114,6

-2

4

-229,2

122,3

-7,7

59,9

2010

113,3

-1

1

-113,3

122,4

-9,1

83,2

2011

113,8

0

0

0

122,5

-8,7

75,7

2012

124,8

1

1

124,8

122,6

2,2

4,9

2013

123,9

2

4

247,8

122,7

1,2

1,5

2014

123,4

3

9

370,2

122,7

0,7

0,4

2015

127,2

4

16

508,8

122,8

4,4

19,2

Итого

1099,2

0

60

4,9

1102,5

-3,3

340,1

Поскольку , то система уравнений принимает вид:

, (6)

тогда

, (7)

. (8)

Подставляя в формулы (7) и (8) рассчитанные величины, получим:

,

Подставив полученные коэффициенты а и b в формулу (4), получим уравнение прямой

.

При правильном выборе уравнения сумма фактических значений урожайности должна максимально приближаться к сумме расчетных значений урожайности , что и получилось 1099,2=1099,2

Таблица 8 Выравнивание ряда динамики аналитическим методом по линейной функции

Годы

Урожайность, ц с 1 га, Yi

Порядковый номер года t

Линейная функция

Yt

Yi - Yt

(Yi - Yt)2

2007

129,6

1

121,8

7,8

60,7

2008

128,6

2

121,9

6,7

45,0

2009

114,6

3

122,0

-7,4

54,4

2010

113,3

4

122,1

-8,8

76,7

2011

113,8

5

122,1

-8,3

69,5

2012

124,8

6

122,2

2,6

6,7

2013

123,9

7

122,3

1,6

2,6

2014

123,4

8

122,4

1,0

1,0

2015

127,2

9

122,5

4,7

22,4

2016

x

10

122,5

x

x

2017

x

11

122,6

x

x

Итого

1099,2

x

х

x

338,9

Рисунок 4 Выравнивание ряда динамики аналитическим методом по линейной функции

Таблица 9 Выравнивание ряда динамики аналитическим методом по логарифмической функции

Годы

Урожайность, ц с 1 га,

Yi

Порядковый номер года, t

Логарифмическая функция

Yt

Yi - Yt

(Yi - Yt)2

2007

129,6

1

124,7

4,9

24,0

2008

128,6

2

123,4

5,2

26,5

2009

114,6

3

122,7

-8,1

65,9

2010

113,3

4

122,2

-8,9

79,2

2011

113,8

5

121,8

-8,0

63,9

2012

124,8

6

121,5

3,3

11,1

2013

123,9

7

121,2

2,7

7,4

2014

123,4

8

120,9

2,5

6,0

2015

127,2

9

120,7

6,5

41,8

2016

x

10

120,5

x

x

2017

x

11

120,4

x

x

Итого

1099,2

x

x

0,1

325,8

Рисунок 5 Выравнивание ряда динамики аналитическим методом по логарифмической функции

Таблица 10 Выравнивание ряда динамики аналитическим методом по полиномиальной функции

Годы

Урожайность, ц с 1 га, Yi

Порядковый номер года, t

Полиномиальная функция

Yt

Yi - Yt

(Yi - Yt)2

2003

129,6

1

129,1

0,5

0,2

2004

128,6

2

123,7

4,9

23,9

2005

114,6

3

119,9

-5,3

27,9

2006

113,3

4

117,6

-4,3

18,6

2007

113,8

5

116,9

-3,1

9,7

2008

124,8

6

117,8

7,0

49,3

2009

123,9

7

120,2

3,7

13,7

2010

123,4

8

124,2

-0,8

0,6

2011

127,2

9

129,8

-2,6

6,5

2012

x

10

136,9

x

x

2013

x

11

145,6

x

x

Итого

1099,2

x

x

0

150,5

Рисунок 6 Выравнивание ряда динамики аналитическим методом по полиномиальной функции

Таблица 11 Выравнивание ряда динамики аналитическим методом по степенной функции

Годы

Урожайность, ц с 1 га, Yi

Порядковый номер года, t

Cтепенная функция

Yt

Yi - Yt

(Yi - Yt)2

2007

129,6

1

124,4

5,2

26,6

2008

128,6

2

123,2

5,4

28,7

2009

114,6

3

122,5

-7,9

63,1

2010

113,3

4

122,0

-8,7

76,5

2011

113,8

5

121,7

-7,9

61,9

2012

124,8

6

121,4

3,4

11,9

2013

123,9

7

121,1

2,8

7,9

2014

123,4

8

120,9

2,5

6,4

2015

127,2

9

120,7

6,5

42,6

2016

x

10

120,5

x

x

2017

x

11

120,3

x

x

Итого

1099,2

x

x

1,3

325,6

Рисунок 7 Выравнивание ряда динамики аналитическим методом по степенной функции

Таблица 12 Выравнивание ряда динамики аналитическим методом по экспоненциальной функции

Годы

Урожайность, ц с 1 га, Yi

Порядковый номер года, t

Экспоненциальная функция

Yt

Yi - Yt

(Yi - Yt)2

2007

129,6

1

121,5

8,1

65,0

2008

128,6

2

243,1

-114,5

13105,4

2009

114,6

3

364,6

-250,0

62509,0

2010

113,3

4

486,2

-372,9

139022,6

2011

113,8

5

607,7

-493,9

243933,9

2012

124,8

6

729,2

-604,4

365342,9

2013

123,9

7

850,8

-726,9

528347,8

2014

123,4

8

972,3

-848,9

720656,1

2015

127,2

9

1093,9

-966,7

934420,0

2016

x

10

1215,4

x

x

2017

x

11

1336,9

x

x

Итого

1099,2

x

x

-4370,1

3007402,7

Рисунок 8 Выравнивание ряда динамики аналитическим методом по экспоненциальной функции

Отбор наиболее адекватной функции проведем с помощью среднеквадратического отклонения:

  1.  
    1. Линейная функция:

  1.  
    1. Логарифмическая функция:

  1.  
    1. Полиномиальная функция:

  1.  
    1. Степенная функция:

  1.  
    1. Экспоненциальная функция:

Наиболее адекватной функцией будет – полиномиальная функция, так как у нее среднеквадратическое отклонение наименьшее.

= 0,7821t2 – 7,7398t + 136,06

Полученное уравнение тренда позволяет установить устойчивость выявленной тенденции и составить прогноз. Для этого рассчитаем показатели колеблемости:

1) Размах колеблемости:

R=(Уi-)max - (Уi-)min , (9)

2) Среднее квадратичное отклонение от тренда:

(10)

3) Коэффициент колеблемости:

, (11)

4) Коэффициент устойчивости:

Куст= 100% - Vy, (12)

где - фактическое значение ряда динамики;

-средний уровень ряда динамик;

У - средний уровень ряда динамики;

По отобранной функции в качестве тренда определим показатели колеблемости и сделаем вывод о возможности прогнозирования.

1. Размах колеблемости:

2. Среднее абсолютное отклонение:

3. Дисперсия колеблемости

=

4. Среднеквадратическое отклонение тренда

ц

5. Относительный размах колеблемости

6. Относительное линейное отклонение

7. Коэффициент колеблемости

8.Коэффициент устойчивости уровня ряда динамики

Так как коэффициент устойчивости больше 50%, то уровни ряда динамики устойчивы и данное уравнение тренда подходит для расчета прогноза на перспективу.

Выполним интервальный прогноз на 2 года:

,

где

- интервальный прогноз,

- табличное значение Стьюдента,

при ,

Интервальный прогноз на 2016 год:

Интервальный прогноз на 2017 год:

Таким образом, если выявленная тенденция по полиномиальной функции сохранится, то в следующие два года с вероятностью 95% можно ожидать увеличение урожайности, причем в 2016 году урожайность на 1 га будет составлять от 666,87 до 681,13 ц, а в 2017 году – от 665,74 до 682,26 ц.

Основным фактором изменения уровня производства продукции растениеводства является урожайность, поэтому следует изучить влияние данного фактора на увеличение валового сбора.

Для этого проведем индексный анализ валового сбора и средней урожайности сахарной свеклы по двум последним годам (таблица 7). Индексный анализ проведем тремя системами индексов.

Таблица 13 Индексный анализ валового сбора и средней урожайности культур

Виды культур

Базисный период

Отчетный период

Валовой сбор, ц

урож-ть, ц с 1 га

пос.площ, га

урож-ть, ц с 1 га

пос.площ., га

базис

отчет

условный

 

У0

П0

У1

П1

У0 П0

У1 П1

У0 П1

Доминика

115,7

196,4

77,8

290,2

22730,5

22589,9

33586,5

Рамонская

122,5

257,6

118,8

219,2

31567,8

26042,9

26862,0

Льговская

169,8

112,7

169,8

132,1

19133,7

22427,1

22427,3

Итого

408,1

566,7

366,4

641,5

73432,0

71059,8

82875,9

В среднем

136,0

188,9

122,1

213,8

24477,3

23686,6

27625,3

Индексный анализ проведем по следующей схеме:

Изменение: Абсолютное Относительное

Валового сбора

(17)

Размера посевных площадей

(18)

Средней урожайности

(19)

(20)

Урожайности отдельных культур

(21)

Посевных площадей

(22)

Структуры посевных площадей

(23)

Подставив в формулы (17) - (23) рассчитанные величины, получим:

Изменение: Абсолютное Относительное

Валового сбора

Размера посевных площадей

Средней урожайности

Урожайности отдельных культур

Посевных площадей

Структуры посевных площадей

Проверим правильность расчетов. Взаимосвязь между рассчитанными показателями:

, (24)

(25)

(26)

(27)

Подставив числовые значения в формулы (24), получим:

, ,

,

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

  1. Афанасьев, В.Н., Маркова А.И. Статистика сельского хозяйства. - М.: Финансы и статистика, 2002 .-270с.

  2. Ефимова, М.Р., Петрова Е.В. Общая теория статистики. - М.: Инфра-М, 2002.-416с.

  3. Зинченко, А.П. Сельскохозяйственная статистика с основами социально-экономической статистики. - М.: Издательство МСХА, 1998. - 427с.

  4. Рафикова, Н.Т. Основы статистики. – М.: Финансы и статистика, 2005. – 512с.

  5. Рафикова, Н.Т. Основы статистики. Уфа, 1995. - 136с.

  6. Система ведения агропромышленного производства в Республике Башкортостан. - Уфа: Гилем ,1997. - 416 с.

Просмотров работы: 739