Методы моделирования

Министерство образования и  науки, молодежи и спорта Украины 

Сумской государственный университет

Кафедра моделирования сложных  систем

 

 

 

 

 

 

 

 

 

 

 

 

 

Курсовая работа

по дисциплине «Моделирование экономических, экологических и социальных процессов»

по теме: «Исследование заданного социально-экономического процесса с помощью практического применения эконометрического моделирования»

 

 

 

 

 

 

 

 

 

Выполнила:                  студентка

группы ИН-83

Чуйко И.Ю.

 

Проверил:         Назаренко А. М.

 

 

 

 

 

 

Сумы 2012 

Содержание

 

Постановка задачи 3

Введение 5

1 Линейные регрессионные уравнения 6

2 Корреляционный анализ системы 13

3 Многофакторные производственные функции 23

Список используемой литературы 30

Приложение. Код программы 31

 

 

Постановка задачи

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

По ходу выполнения работы должны быть раскрыты следующие вопросы.

  1. Теория применения линейных производственных функций (ПФ) вида

 

при моделировании социально-экономических  процессов.

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

а) построение однофакторной, двухфакторной и трехфакторной ПФ и оценивание их качества приближения с помощью коэффициента детерминации ;

б) формулировка выводов, объясняющих полученные результаты.

  1. Теоретические основы проведения корреляционного анализа систем «показатель — факторы» и «факторы».
  2. Проведение корреляционного анализа систем «показатель-факторы» и «факторы» на данных своего варианта.
  3. Проведение регрессионного анализа линейной трехфакторной ПФ.
  4. Теория применения ПФ типа Кобба-Дугласа вида

 

при моделировании  социально-экономических процессов.

  1. Построение ПФ типа Кобба-Дугласа на примере данных своего варианта и оценивание ее качества приближения с помощью коэффициента детерминации .
  2. Проведение корреляционного и регрессионного анализа линеаризированной ПФ типа Кобба-Дугласа.
  3. Сравнение линейной модели и функции типа Кобба-Дугласа, выявление наиболее адекватной модели и описание ее с экономической точки зрения.

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

Для заданного варианта необходимо выполнить задания курсовой работы (приведенные выше), придерживаясь  требований по содержанию и правил оформления, которые даны ниже.

Вариант № 99

 

454.72

496.59

588.72

597.19

664.02

749.32

790.78

862.93

928.91

 

84.02

93.70

110.91

116.34

123.48

144.62

155.95

160.84

176.90

 

68.59

89.88

97.99

115.53

135.02

145.68

171.22

179.79

211.15

 

26.38

32.19

37.72

41.63

49.21

53.51

59.70

68.07

73.92


 

 

Введение

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

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

 

 

1 Линейные регрессионные уравнения

Будем исходить из того, что  между объясняемой и объясняющими переменными выбрана линейная связь. Имеем

 

Здесь – фиктивная переменная, введённая для удобства; слагаемое отражает влияние на других факторов, ошибки измерений, ошибки выбора модели.

Пусть с целью исследования линейной связи проведена выборка  объёма . Тогда для наблюдаемых величин можно записать

 

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

Если  – возможные оценки (приближенные значения) параметров , то функция регрессии, соответствующая модели, имеет вид

 

Отклонения выборочных данных от неё определяются величинами

 

Отметим, что несмещенной  оценкой неизвестной дисперсии  возмущений является выборочная дисперсия

 

где – число связей, накладываемых функцией регрессии на выборку. Следуя положениям § 3.3 [1], заключаем, что , т.е. общее число связей равно числу оценок, от которых зависит функция регрессии.

Критерием выбора оценок в математической статистике является условие минимума дисперсии, которое при фиксированном значении эквивалентно условию минимума функции ошибок

 

Имеем

 

В результате для определения  МНК-оценок приходим к системе линейных уравнений с неизвестными

 

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

Таким образом, функция регрессии, соответствующая МНК-оценкам , имеет вид

 

Несмещенной оценкой неизвестной  дисперсии  является МНК-оценка

 

 

Матричный способ оценки

Процесс оценивания регрессионной  модели при ³ является довольно громоздким, поскольку приходится вычислять большое число сумм и решать системы уравнений с тремя и более неизвестными, что без использования ЭВМ весьма затруднительно. Если же в распоряжении пользователя ЭВМ имеются стандартные программы, позволяющие осуществлять действия над матрицами, то регрессионный анализ значительно упрощается.

Запишем регрессионное соотношение для наблюдаемых величин в развёрнутом виде

 

Здесь для всех . Матричная запись системы уравнений такова:

 

где

 

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

 

имеет ранг, равный , и, следовательно, существует обратная матрица .

Нетрудно заметить, что  система линейных уравнений, из которой  определяются МНК-оценки , может быть записана в виде

 

откуда находим вектор-столбец  искомых МНК-оценок. Имеем

 

Таким образом, вектор-оценку можно определять двумя способами: либо решая систему линейных уравнений, либо пользуясь формулой.

Далее получаем

 

 

 

 

 

 

 

Коэффициент детерминации

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

 

Можно показать, что 

££.

Значение  (в процентах) означает, что линейная модель объясняет всей дисперсии показателя, остальные не обусловлены линейной моделью.

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

 

 

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

Для однофакторного регрессионного уравнение вида

 

имеем следующее:

.

Проведем вычисления по формуле

 

       

 


тогда:

 

Таким образом, оцененное  однофакторное уравнение регрессии  имеет вид

.

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

 

Полученное значение указывает на то, что модель объясняет исходных значений, а остальные носят случайный характер.

Рассмотрим теперь двухфакторное  регрессионное уравнение вида

.

Проведем вычисления по формуле

 

      


тогда:

 

Таким образом, оцененное  двухфакторное уравнение регрессии  имеет вид

.

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

 

Полученное значение указывает  на то, что данное уравнение объясняет  исходных значений, а остальные 0,74 носят случайный характер.

Рассмотрим теперь трехфакторное  регрессионное уравнение вида

.

Проведем вычисления по формуле

 

       


тогда:

 

само уравнение будет  иметь вид

.

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

 

Полученное значение указывает  на то, что данное уравнение объясняет  исходных значений, а остальные 0,13 носят случайный характер.

Анализируя полученные результаты можно заметить, что при добавлении факторов в модель ее адекватность повышается. Например, добавление к двуфакторной модели позволило улучшить качество приближения на 0 но в то же время добавление второго фактора улучшило приближение на

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

 

  2 Корреляционный  анализ системы

Изучение проблемы спецификации переменных (выбора факторов) следует  начинать с корреляционного анализа  связи между переменными. Будем  различать корреляционный анализ системы  «показатель-факторы» и корреляционный анализ системы «факторы».

 

Корреляционный анализ системы  «показатель-факторы»

Запишем матрицу исходных данных в виде

 

Матрица имеет размерность . Тогда транспонированная матрица имеет размерность , и для матрицы размерности получаем

 

Элементы полной корреляционной матрицы 

 

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

 

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

 

Частный коэффициент корреляции характеризует истинную корреляцию показателя и фактора при исключении влияния других факторов.

Значимость выборочных частных  коэффициентов корреляции проверяется по двустороннему критерию

 

где критическое значение вычислено по таблице Фишера-Иейтса при уровне значимости и числе степеней свободы . Выполнение критерия для коэффициента означает наличие стохастической связи между показателем и фактором , и переменную следует включить в число существенных переменных.

Может случиться так, что  для одного или нескольких факторов неравенство не выполняется. В подобных ситуациях многие исследователи  исключают соответствующие факторы  из модели. Однако к исключению из модели переменных следует относиться осторожно. Ведь таким образом из модели может  быть исключена переменная, оказывающая  существенное влияние на показатель . Может случиться ситуация, когда несколько частных коэффициентов корреляции окажутся незначимыми, и исключенная вместе с другими переменная с начальными незначимым коэффициентом могла бы иметь значимый коэффициент в сокращенной модели, если бы она была там оставлена.  При изменении состава системы «показатель-факторы» будут изменяться матрицы и , а значит, другие значения будут иметь и частные коэффициенты корреляции.

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

С помощью корреляционного  анализа системы «показатель-факторы» мы только определяем кандидатов на исключение из числа объясняющих переменных. Какие из них должны быть исключены  из модели может подсказать корреляционный анализ системы «факторы».

 

Корреляционный анализ системы  «факторы»

Рассмотрим матрицу

 

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

 

через элементы, которой выражаются выборочные частные коэффициенты корреляции системы «факторы». Имеем

 

Значимость коэффициентов проверяется по двустороннему критерию при уровне значимости a и числе степеней свободы . Если для одного или нескольких коэффициентов выполняется , это будет означать наличие в модели проблемы мультиколлинеарности (зависимости между факторами). В этом случае спецификация переменных должна быть пересмотрена заново. Количество объясняющих переменных необходимо уменьшить.

Каким образом осуществлять новую спецификацию переменных? Представляется разумным последовательно уменьшать  количество переменных на единицу и на каждом шаге заново проводить корреляционный анализ системы “показатель-факторы” и системы “факторы”. Здесь на каждом этапе могут возникать следующие ситуации.

  1. Только один из коэффициентов оказался значимым. В этом случае из двух переменных и из модели исключается та, которой соответствует незначимый выборочный коэффициент корреляции. Случай, когда обе переменные ранее оказались существенными маловероятен. Если же обе переменные ранее попали в список кандидатов на исключение из модели, то должна сработать интуиция исследователя. Можно отработать также оба варианта.
  2. Сразу несколько коэффициентов значимы. Здесь может возникнуть большая свобода в принятии решения о том, какую из переменных удалить. Если есть сомнения в правильности действий, необходимо просчитывать разные варианты.

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

В решении проблемы спецификации модели может помочь регрессионный  анализ. Иногда, вместо того, чтобы исключать  из модели одну из двух или нескольких переменных, следует подумать о сохранении в модели сразу всех переменных, например, путем линейного ограничения  на параметры. Количество объясняющих  переменных будет уменьшено, однако данные будут использованы по всем переменным. Такая процедура, как  правило, приводит к повышению точности модели.

Регрессионный анализ модели

Будем исходить из линейной модели

 

которая на базисных данных в векторной форме принимает вид

 

Оцененная регрессия 

 

где вычислены ее основные характеристики.

Вектор МНК-оценок удовлетворяет системе линейных уравнений

 

и может быть получен по формуле

 

Оцененная ковариационная матрица  оценок коэффициентов регрессии имеет вид

 

где –оценки дисперсии остатков, равная

 

На диагонали матрицы  находятся оценки дисперсий оценок , корни из которых дают стандартные ошибки коэффициентов регрессии

 

Оцененная ковариационная матрица  значений на базисных  данных имеет вид

 

На ее диагонали находятся  оценки ,, корни из которых представляют собой стандартные ошибки значений регрессии на базисных данных.

Регрессионный анализ  позволяет  проверить значимость объясняющих  переменных в уравнении регрессии. Для этого можно использовать двусторонний критерий значимости коэффициентов  регрессии

 

где – критическое значение, вычисленное по таблице Стьюдента при уровне значимости α и числе степеней свободы . Этот критерий эквивалентен критерию Фишера, и может быть осуществлена проверка результатов корреляционного анализа. В случае необходимости можно использовать односторонние критерии значимости коэффициентов регрессии, если из экономических соображений или из опытных данных имеются большие аргументы в пользу положительного или отрицательного направления действия объясняющей переменной на показатель .

Предположим, что нами оценена  модель (2.10)

 

Значимость регрессионной  модели в целом определяется согласно критерию

 

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

 

Пусть корреляционный и регрессионный  анализы модели, а также сокращенных  версий, получающихся после последовательного  удаления из исходной модели на наш  взгляд несущественных переменных привели  нас к модели

 

Значимость суммарного вклада исключенных переменных можно проверить с помощью -статистики

 

которую в данном случае можно  переписать в эквивалентной форме

 

Если выполняется неравенство

 

где – критическое значение, равное соответствующему квантилю распределения Фишера  при уровне значимости и числах степеней свободы и , то мы делаем вывод о значимости совместного вклада исключенных переменных в “длинную” регрессию. В этом случае “короткая” регрессия должна быть скорректирована, спецификация переменных должна быть пересмотрена заново.

Проведение корреляционного анализа

Проведем корреляционный анализ системы «показатель-факторы» для заданного варианта задания.

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

Матрица из примет вид:

 

Зная  можно рассчитать матрицу :

 

Используя которую можно  рассчитать полную корреляционную матрицу :

 

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

Систему «показатель-факторы» будем характеризовать выборочными  частными коэффициентами корреляции . Для их вычисления рассчитаем значения матрицы :

 

Зная ее рассчитаем искомые  частные коэффициенты корреляции:

 

Для того, чтобы определить, какие факторы значимо влияют на изучаемый показатель найдем . При уровне значимости и числе степеней свободы , . Тогда из критерия , видим, что из трех факторов только материальные затраты на выпуск продукции не значимо влияют на стоимость выпущенной продукции (при исключении влияния других факторов).

Проанализируем теперь систему  «факторы».

Матрица будет иметь вид:

 

Тогда, рассчитав матрицу , через элементы которой выражаются выборочные частные коэффициенты корреляции системы «факторы», получим:

 

и найдем матрицу

 

элементы которой указывают  на наличие в линейной модели проблемы мультиколлинеарности, поскольку значение превосходит

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

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

.

В этом случае система «показатель-фактор»  будет иметь вид .

Проводя расчеты, аналогично приведенным выше, получим:

    1. Частные коэффициенты корреляции  будут равны: 

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

    1. Матрица выборочных частных коэффициентов корреляции системы «факторы» будет равна:

 

элементы которой указывают  на наличие в линейной модели проблемы мультиколлинеарности, поскольку значение превосходит

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

Убедимся теперь, что регрессионный  анализ линейной модели приведет к  аналогичным результатам. Модель имеет  вид:

.

Рассчитав отклонение функции  регрессии от заданных значений показателя , с помощью формулы (2.16) найдем

 

 

Воспользовавшись результатами предыдущих расчетов из главы один, найдем

 

 

 

Откуда     

Тогда уравнение регрессии  перепишется в виде

 

 

 

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

 

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

Проверим значимость регрессионной  модели в целом. Получим:

 

что намного превосходит  и значит уравнение регрессии в целом значимо. 

3 Многофакторные производственные  функции

Функция

 

называется производственной функцией Кобба-Дугласа. Для нее  эластичность замещения ресурсов , т.е. при увеличении капиталовооруженности производства на 1% предельная норма их замещения также увеличивается на 1% независимо от объемов используемых ресурсов.

Для функции Кобба-Дугласа  предельные производительности ресурсов

 

 

Коэффициенты эластичности выпуска по затратам ресурсов

 

 

Как и следовало ожидать, сумма коэффициентов эластичности

 

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

 

В процессе производства, описываемом  функцией Кобба-Дугласа, ресурсы взаимозаменяемы, причем с ростом затрат одного ресурса  высвобождается все большее количество другого. Линии уровня функции имеют вид, представленный на рисунке.

Формула Кобба-Дугласа  является частным случаем более общей формулы:

 

где показатели эластичности выпуска по затратам капитала и труда 

 

не связаны между собой. Можно предположить, что обе величины a и b находятся между и . Показатели a и b больше нуля, так как увеличение затрат ресурсов и должно вызвать рост выпуска . С другой стороны, a и b меньше единицы, так как разумно предположить, что увеличение затрат и приводит к более медленному росту выпуска продукции (затраты других ресурсов предполагаются постоянными).

Производственная функция  Кобба–Дугласа (3.5) является однородной, причем

 

Если по отдельности показатели эластичности a и b указывают на процентное увеличение (или уменьшение) выпуска при однопроцентных колебаниях величин капитала и труда , то сумма ab отразит уже общую реакцию производства на указанные изменения затрат. Степень однородности nab характеризует эффект от масштаба производства.

Пусть,  например, затраты  и возросли в пропорции l, т.е. расход и капитала, и труда увеличился на . Тогда новый уровень выпуска

.

При ab имеем постоянный эффект от масштаба производства ( увеличивается в той же пропорции, что и , и ), наблюдается постоянная отдача факторов. Если ab, то рост производства превышает отметку. В подобных случаях (при возрастании и в некоторой пропорции выпуск растет в большей пропорции) говорят о возрастающем эффекте от масштаба производства или положительном эффекте масштаба. Его экономический смысл заключается в повышении эффективности производства под воздействием концентрации ресурсов, кооперации труда. Если же ab, то рост выпуска не достигает отметки. Налицо убывающий эффект от масштаба производства ( растет в меньшей пропорции, чем и ) или отрицательный эффект масштаба: наращивание затрат ресурсов оборачивается снижением их продуктивности. Подобное бывает, например, при переконцентрации производства, когда его размеры выходят за оптимальные границы.