Исследование зависимости случайных составляющих стохастической производственной функции при оценке технической эффективности Текст научной статьи по специальности «Математика»

С. А. Айвазян, М. Ю. Афанасьев, В. А. Руденко

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

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

ключевые слова: эконометрическая модель; производственный потенциал; производственные факторы; интеллектуальный капитал; копула-функции; нормальная копула-функция; зависимость компонентов ошибки. JEL classification: C50; C51; C52; D24.

мпирическое оценивание производственной функции или функции затрат являет-

ся стандартной эконометрической задачей. Понятие граничной производственной

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

Понятие «стохастической границы» впервые было введено в работах (Meeusen, van den Broek, 1977; Aigner et al., 1976). Стохастическая граница, наряду с детерминированными составляющими модели, описывающими поведение основных факторов производства, включает в себя и случайную составляющую, моделирующую сопутствующие факторы, оказывающие воздействие на производственный процесс.

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

1. Введение

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

В работе (Aigner et al., 1976) впервые рассматривалась модель вида

У, = f(X, ß) + e,, (1)

где y t — выпуск продукции 7-й компании, xt — вектор входных данных, ß — вектор параметров. При этом предполагалось, что ошибка e7 представлена в виде суммы двух случайных величин e7 = n 7 + u , где n 7 — независимые одинаково распределенные случайные величины, имеющие нормальное распределение с параметрами (0; ouv ), а случайные величины u 7 независимы от n 7 и удовлетворяют условию u7 < 0. Основное внимание уделялось случаю, когда в качестве распределения величин —u 7 используется усеченное в нуле нормальное распределение N + (0; о и). Кроме того, предлагалось рассматривать и другие варианты, например, когда величины — u 7 имеют показательное распределение. Важно отметить, что в данном случае выполняется неравенство y7 < f (x7, ß) + v7 и, таким образом, граница в этой модели является действительно стохастической.

Обоснование приведенной спецификации модели заключается в том, что процесс производства подвержен влиянию двух экономически различимых факторов. Компонента ut < 0 обеспечивает нахождение показателей выпуска каждой из компаний ниже стохастической границы y7 < f (xt, ß) + v i. При этом предполагается, что эти отклонения от границы являются результатом влияния факторов, находящихся под контролем компаний, таких как экономическая неэффективность распределения ресурсов, усилий производителя и его рабочих и т. п. Но сама граница может быть различной как для разных компаний, так и в течение времени для одной и той же компании. В данной интерпретации компонента n 7 отражает внешние благоприятные или, напротив, неблагоприятные факторы, на которые компании повлиять не могут (например, климатические условия, экономическое состояние региона и т. п.). Ошибки в измерениях при наблюдениях также составляют одну из частей компоненты n . Интерес представляет сравнение относительного вклада компонентов n 7 и u 7 при оценивании дисперсий о n и о2 соответственно. Кроме того, эта спецификация предполагает, что

оценка эффективности производства должна осуществляться на основе отношения

y

f (x,, ß) + V,'

а не с помощью отношения

y,

f (х,, Ь)'

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

В данной работе будем исследовать класс моделей производственной функции компании, в общем случае имеющий вид

=Ь0 (х] Г1 Г... к Г (2)

где Я — объем производства, зафиксированный для ,-й компании, , =1,..., п; х],..., х, — объемы факторов производства; п — число компаний. Тогда yi = 1пЯ,, Xi = (х , . . . , хр) ,

ь = (Ьо, Ь],..., Ьр), е,= V-и,,

* = f(x,b) = lnbo + b)x) +... + b,xp +e,. (3)

I

<u ¡5

Случайная величина Vi подчиняется (0; oV )-нормальному распределению; случай- ^

ч оа

ная величина Ц распределена в соответствии с усеченным в нуле нормальным законом, имеющим среднее значение т и дисперсию оЦ. Тогда, в соответствии с концепцией сто- £ хастической границы, мерой технической эффективности производства является величина TEi = Е(ехр{—Ц.}| е.). Предпосылкой, существенно упрощающей оценку параметров функции (3) и используемой в стандартных методах максимизации функции правдоподобия, яв- ^ ляется предположение о том, что случайные величины У{ и и. статистически независимы. ®

Вопрос о том, является ли это предположение существенным для рассматриваемого клас- ^

са моделей, остается открытым. Напомним, что случайная составляющая У объясняет внеш- ^

00

ний «статистический шум» в модели, а Ui — внутреннее систематическое снижение эффективности компании, называемое технической неэффективностью. Однако может случиться, что данные величины являются зависимыми. Так, например, при некоторых внешних экономических изменениях может оказаться, что крупные компании, сумевшие «удержаться на плаву» и воспользоваться сложившейся ситуацией, смогли не только остаться в лидерах по оценкам технической эффективности, но и повысить свою достижимую стохастическую границу, в то время как маленькие компании будут вынуждены снижать свою техническую эффективность. При этом величины Уi и и. могут оказаться зависимыми.

Также вопросы о соответствии реальности возникают в тех случаях, когда при нахождении оценок параметров оказывается, что У i близко к нулю (т. е. дисперсия оУ близка к 0). Возможно, такой эффект возникает из-за искажения оценок вследствие неверной предпосылки о независимости величин У и и .

о

2. Проверка гипотезы о независимости случайных величин V, и U, с помощью копула-функций

Целью данного исследования является проверка статистической гипотезы о независимости случайных величин V и U . Будем проверять эту гипотезу с помощью аппарата копула-функций, позволяющего описывать законы многомерного распределения вероятностей. В ряде работ уже были произведены расчеты для моделей стохастической граничной функции с помощью построения копул. Так, в (Amsler et al., 2014; Shi, Zhang, 2011) исследуется зависимость компонент Uit между собой в разные моменты времени для случая панельных данных. В (Lai, Huang, 2013; Carta, Steel, 2012) рассматриваются модели с множеством оценок выпусков компаний в предположении зависимости компонент Uik между выпусками1. В работе (Smith, 2008) рассмотрена задача оценки зависимости V и U, однако трудно судить о качестве построенных моделей, т. к. ни принятого критерия качества, ни истинных значений параметров, описывающих эту зависимость, не приводится.

В соответствии с (Айвазян, Фантаццини, 2014; Благовещенский, 2012) приведем некоторые определения и теорему, необходимые для дальнейшего изложения.

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

Определение 1. Функция C(ul, uu, ..., uk), определенная на единичном k-мерном кубе Ik = [0;l]k (т.е. ut £ [0;1], i = 1, U, ..., k), называется копула-функцией, если она обладает следующими свойствами:

1) область значений функции — единичный интервал [0;1];

2) если ut = 0 хотя бы для одного значения i £ {l,U,...,k}, то C(u1, uU, ., uk) = 0;

3) C(1, .,1, u7,1,...,1) = ut для всех ut £ [0;1];

4) C(u1, uU, ..., uk) является k-возрастающей функцией в том смысле, что для всех (a1, au, ..., ak ) £ [0;1]k и (b1, bu, ..., bk ) £ [0;1]k с ai < bt справедливо равенство:

UU

1)^+..+7k C(u1. ,uU^,...,ukit) > 0, где uj1 = üj и uj2 = bj для всех j £ {1, 2, ..., k}.

i1=1 7k =1

Другими словами, копула — это k-мерное распределение на k-мерном кубе I = [0; 1] , все маржинальные распределения которого являются равномерными на отрезке I = [0; 1].

Теорема 1 (Теорема Склара). Пусть H(•) — k-мерная функция распределения с частными распределениями F1,...,Fk. Тогда существует k-мерная копула-функция C(•) такая, что для всех действительных x1,xu,...,xk выполняется:

H (x, xu,..., xk) = C (F1( x), Fu (xu),..., Fk (xk)). (4)

Если все частные функции распределения непрерывны, то копула-функция определена единственным образом. Обратно, если C (•) — копула-функция, а F1,..., Fk — функции распределения, то функция H(•), определяемая выражением (4), является совместной функцией распределения с частными распределениями F1,..., Fk .

Доказательство можно найти в (Sklar, 1996; Joe, 1997; Благовещенский, 2012).

Следуя (Айвазян, Фантаццини, 2014), рассмотрим два типа копула-функций: эллиптические и архимедовы.

Определение 2. Пусть X — k-мерный случайный вектор и S£ Rkxk — неотрицательно определенная симметричная матрица. Если существуют m£ Rk и R ^ R+ такие, что характеристическая функция вектора X — m имеет вид2

у x—m(t) = exP(it T m) u

для любого t £ Rk, то X называют случайным вектором, имеющим распределение эллиптического типа с параметрами m и S.

Копула-функция C(•) называется эллиптической, если она соответствует распределению эллиптического типа.

Дадим определение архимедовых копул в двумерном случае.

Определение 3. Рассмотрим непрерывную, строго убывающую и выпуклую функцию f(u), определенную при u £ [0;1] и удовлетворяющую условию f(1) =0. Определим псевдообратную функцию f[-1](t) соотношением:

Г— 11/ч if—1 (t) для 0 < t <f(0) ,

fr 1 (t) = \ , где f (t) — обычная обратная функция к f (u).

10 для t > f (0)

2

2 i — мнимая единица.

Двумерная функция С: [0;1]2 ^ [0;1], определенная как С(и1,и2) = ф[ 1] (ф(«) + ф(и2)), называется архимедовой копула-функцией с генератором ф(-).

I

ф

s

Архимедовы копулы могут быть представлены в явном аналитическом виде и, как прави- Q-ло, допускают простую вычислительную реализацию, вследствие чего применяются в раз- „^ личных областях, например, в банковской и сельскохозяйственной сферах (Пеникас, Сима- ® кова, 2009; Салмин, 2013; Кучмент, Демидов, 2013). <5

Эллиптические копулы, как правило, применяются в актуарной математике, страховании j|

vT

и финансовой сфере (Frees, Valdez, 1998; Rosenberg, Schuermann, 2006, и др.). Некоторые из них имеют ряд общих свойств с многомерным нормальным распределением. Кроме то- ^ го, преимуществом копул этого типа является тот факт, что они позволяют учесть всю ко- ^ вариационную матрицу, описывающую зависимость нескольких переменных, в отличие от S3

QQ

архимедовых копул, оценивающих лишь один параметр «связи».

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

Определение 4. Многомерная нормальная (гауссова) копула-функция определяется соотношением:

C(u,u2,...,uk;S) = Фк (Ф-1 (u),Ф-1 (u2),...,Ф-1 (uk);S),

где Фк (•,...,•; S) — k-мерная нормальная функция распределения с нулевым вектором средних значений и ковариационной матрицей S, Ф-1 (•) — обратная функция одномерного стандартного нормального распределения.

В данной работе будем строить модель (3) с зависимыми величинами U и V с помощью двумерной нормальной копула-функции (Айвазян, Фантаццини, 2014).

л Ф (Mj ) Ф (U

C (Ui, u2, r) = 2 J J

^) ^) х2 - 2rxj + y2

I-2 J J exP

2ж\11- r -» -»

dx dy = , ч

(5)

о

2(1- г2)

= Ф2 (Ф-1 («1), Ф-1 («2); г),

где Ф2 (•) — функция двумерного нормального распределения с нулевым средним значением, единичными дисперсиями компонентов и корреляцией г, Ф- (и1) — обратная функция одномерного стандартного нормального распределения.

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

Согласно теореме Склара СШгт^ (у), Еи (и)) = Н (V, и), где СЫогт (•) — нормальная копула-функция, FV (у) — распределение случайной величины V, т. е. N (0; ст^), Ри (и) — распределение случайной величины и, т. е. ^(т;о2и), Н(V,и) — совместная функция распределения и и V с частными распределениями FV (у) и Fи (и).

Плотность двумерной нормальной копула-функции в общем случае имеет вид:

Ф1, «2) = ехр(-\? (2-1 - 1)С1 (6)

где £ = (Ф ] (и]), Ф ] (и2)) — вектор, компонентами которого являются обратные функции

/] И

одномерного стандартного нормального распределения, 2 = рица, I — единичная матрица.

,г 1, V V

ковариационная мат-

Совместная плотность распределения величин и и V при этом выражается как

f (х, у) = с(и, и2 )■ ^ (х)- /и (у), (7)

где /у(х) и/и(х) — одномерные плотности распределения величин V и и соответственно. В рассматриваемой модели получаем следующую совместную плотность распределения:

f ( x, У) =

Xexp

1

1

1

Vl- r2 2ж°и °V 1-Ф(-/М / Su )

•exp

1 x2 1 (y-m)2

2 ct2 2

о

1

2(1- r2 )

x2 r2 2rx

V

о

Ф-

о

ф(( y -m)/ s и )-Ф(-т/^и) 1-Ф(-т/сти )

X

+

+r2

/ /

Ф-1 \ v

ф(( y -m)/ Ou )-Ф(-т/ои )

1-ф(-т/ои )

(8)

Подробные расчеты приведены в Приложении 1. С помощью (8) найдем плотность е = V — и.

о у+г (Г (г)

Г(г) = Р(е < г) = P(V — и < г) = / / Дх,у)dхdy, /е(г) = ^^ = //(у + г,у)(у.

о -о (2 о

Обозначим в модели (3) zi = 1пЬ0 + Ь]х] +... + Ьрхр. Тогда yi = г, + е,. Пусть Н7 (0 — плотность распределения У (случайной величины, реализациями которой являются у,). Тогда Ну (уг) = / (уг — г,). »

Функция правдоподобия в модели (3) имеет вид Ь =11 Ну (у1), а логарифмическая функ-

п ,=]

ция правдоподобия есть I = ^ 1п Иу (yi).

Для получения оценки максимума правдоподобия (ОМП) надо максимизировать функцию I. Необходимым условием для этого в общем случае является решение системы из (р + 5) уравнений о равенстве нулю производных функции I по оцениваемым параметрам

(^о, АЬ р, т ои, sv, г).

Выразим Н7(у,.) через функцию стандартного нормального распределения:

0

Ну (у,) = /£ (у, — г,) = / / ^ + у, — ^, t) dt =

V1-r2 2pouOv 1-Ф(-т/ои)

/

exp

(t + У, - г,-)2 2о2

X

Xexp

(t-m)2

2о2

1

2(1-r2 )

(t+yr -г,)2r2 - 2r(t+y,,-г,) ф-1

о

о

Ф((Г-тУ s и )-Ф(-тК ) 1-ф(-тК )

+

2

/

0

1

+r2

хф-1 f ф( -mV s и )-Q(-m/g и)xX

2 \

1-ф(-т/сти)

//

(9)

I

<u

I

Производные, необходимые для исследования функции hY (y. ), приведены в Приложе- ^

нии 2. о

Получить оценки параметров b0, bi,..., b p, m,°и, °V,r в явном виде затруднительно, по- <S

этому систему уравнений следует решать численными методами. Для этого в случае p = 3 _Ц

использовался макрос, написанный в программе Excel, который позволяет вычислить с за- ^

данной точностью (при необходимости до е-300) значения функции hY(yi ) в окрестности ® некоторой начальной точки (Ь°, b°, b2, b°, , °U, °v,r 0). Он также позволяет вычислить

n <5

SO

ln Ну (yt) и убедиться, что в найденной точ- >s

ке максимума они равны нулю. При этом начальную точку вида (b°, b°, b2, b0, , °U , °V ,0) ^ будем искать в пакете Stata 10.0 с помощью методологии, описанной в работе (Айвазян и др., 2012). Созданный макрос позволяет перебором с заданным шагом найти наибольшее зна-

n

чение функции l = ^ ln HY (y i).

i=1

Напомним, что целью данного исследования является проверка справедливости предпосылки о независимости случайных величин V и U. Для ее достижения будем проверять статистическую гипотезу о равенстве нулю коэффициента r — Hr: r = 0 при альтернативе HAr : r * 0.

Для проверки воспользуемся статистикой Lr = 2[lnL(HA)-ln L(Hr)], где lnL(HA) — логарифм максимального значения функции правдоподобия при альтернативной гипотезе, ln L(H r) — логарифм максимального значения функции правдоподобия при нулевой гипо-

n

тезе (т. е. максимальное значение функции l = ^ ln HY (y i) при фиксированном r = 0).

i=i

Распределение статистики Lr, см. (Self, Liang, 1987; Айвазян и др., 2012), в условиях справедливости гипотезы Hr является (асимптотически при n ^ оо) смесью распределений случайных величин % (0) и % (1) с весами 1/2 и 1/2, где под распределением % (0) понимается вырожденное распределение, при котором случайная величина равна нулю с вероятностью единица. Таким образом, гипотезу Hr следует отвергнуть, если при заданном уровне значимости a значение тестовой статистики Lr окажется больше a-квантиля Lr (a) = %2-2а (1) упомянутого распределения, т. е. квантиля уровня 1 -2a распределения % (1)3.

3. Проверка независимости случайных величин V, и Ц на смоделированных данных

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

3 Можно также провести классический тест на независимость компонент копулы, описанный в (Genest, Remillard, 2004).

4 Данные могут быть предоставлены авторами ([email protected]).

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

1. По набору из 80 векторов данных ((, x), xf, xf), построенных по реальным наблюдениям для американских компаний из отрасли «Healthcare» за 2009-2012 годы, некоторые из которых уже использовались в работе (Айвазян и др., 2012)5, в соответствии со спецификацией, описанной в (Айвазян и др., 2012), были получены оценки параметров

((0, А , b2, b3 , SV , ^ SU ).

2. Были смоделированы зависимые величины Ui ~ N+ (m;sf ) и Vi ~ N(0;oV ) с коэффициентом корреляции p = 0.945.

3. Были построены yt = ln b0 + b, x) + b2 xf + b3 xf + V — Ut.

Найдем и сравним оценки, полученные при построении модели:

1) в предположении, что Ui и V являются независимыми случайными величинами, т. е. r = 0;

2) в предположении, что Ui и Vi являются зависимыми случайными величинами с ненулевой корреляцией и, соответственно, r Ф 0.

При моделировании были взяты следующие параметры:

b0 = e089, b, = 0.69, b2 = 0.2), b3 = 0.), of = o.)6,m = 0, sf = 0.04, p = 0.945, J, = г,. + e = г, +V — Ut, z,. = ln be + b)x1 + bvxv + bf x3.

Диаграмма разброса для смоделированных Ut и V представлена на рис. 1.

0.6 - U ♦

0.5 - ♦

0.4 - 4 *

0.3 - 0.2 -> * / /

а,- 00— V

-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5

Рис. 1. Диаграмма разброса U i и Vi

В случае ошибочного предположения о независимости Uf и V будем действовать согласно схеме спецификации трехфакторных моделей, описанной в (Айвазян и др., 2012). Вычисления проводим в пакете Stata 10.0. Получим следующую цепочку рассуждений:

5 При этом x = (x¡, xf, xi) = (ln Kt, ln Lt ,ln I¡).

Согласно (Айвазян и др., 2012), в качестве оценки интеллектуального капитала была использована: Ii = р — р, где р = (l + ROAil 100) K¡ — оценка рыночной стоимости компании, без учета интеллектуального капитала, ROAi — показатель отдачи активов (в %), P¡ — рыночная стоимость компании, равная числу акций, умноженному на стоимость одной акции.

1м2;М + ; Яо; Я—; Н2; Я—;М1;М+; Я0; Я—; Я1; Я-; Я1,2; Н+1

I М1 1.

Подробный протокол приведен в Приложении 3.

В Приложении 2 показано, что в случае г = 0 аналитический вид функции правдоподо-

п

бия Ь = ^ (yi) совпадает с видом функции правдоподобия с независимыми компонен-

1=1

тами Ц и V, максимизируемой в пакете Stata 10.0.

В предположении г = 0 полученная модель обозначается М1, наибольшее значение 1пЬ(Нг) = -20.7811. В предположении г Ф 0 полученная модель обозначаетсяМг, наибольшее значение 1п Ь(Н) = —19.5869. Значения оцененных параметров этих моделей приведены в табл. 1.

Таблица 1. Оценки параметров моделей М1 (при г = 0) и Мг (при г Ф 0)

М1 Mr

Оценки факторов производства

lnK 0.648*** 0.66

ln L 0.227*** 0.217

in/ 0.143*** 0.144

const 0 94*** 0.99

Оценки параметров компонентов ошибки

m 0 -0.871

°2V 0.018 0.436

°2U 0.257 0.781

r 0 0.966

I

ф

I

оа

(о Ф

О

1

3

п и

О >S

о

Р-значение при проверке гипотезы об отсутствии неэффективности

Число наблюдений

Логарифм функции правдоподобия

0.000

80 -20.7811

80 -19.5869

Примечание. В скобках указаны стандартные ошибки; *** — значимость оценок коэффициентов на 1%-ном уровне. Значимость коэффициентов в модели Мг и параметров компонентов ошибки не оценивалась.

Статистика Ьг = 2[1п Ь(Н)—1п Ь(Нг)] = 2.3884. При уровне значимости а = 0.1 имеем X2- 2а (1) = 1 6424. Гипотезу о равенстве нулю коэффициента корреляции для данной выборки компаний следует отвергнуть в пользу альтернативной Нг.

Отметим, что созданный макрос не осуществляет проверку значимости оценок параметров. Но целью данной работы является проверка статистической гипотезы о независимости случайных величин Ц и V, кроме того, полученные оценки коэффициентов при основных факторах производства оказались достаточно близки в обеих моделях. Поэтому сделан выбор в пользу модели Rl = е099К066Ь021710144/--1и, Ц ~ N+(—0.871; 0.781), V- ~ N(0; 0.436),

и Vi являются зависимыми случайными величинами с коэффициентом г = 0.966. Соответствующ V ~ N(0; 0.436).

Соответствующая стохастическая граничная функция равна Ri = е099 К10ббЬ°'2171°Л44е,/'

В условиях предположений г = 0 и г Ф 0 рассчитаны (см. Приложение 4) оценки технической эффективности и их ранги.

exp(-U)

♦ lu ♦

1.2 п exp(-U)

>Г* ♦

а) б)

Рис. 2. Диаграммы разброса: а) е~и и ТЕН для М1; б) е~и и TEri для Мг

На рисунке 2а приведена диаграмма разброса истинных е~и' и оцененных технических эффективностей по моделиМ1. На рисунке 2б приведена диаграмма разброса истинных в~ i и оцененных технических эффективностей по модели Мг. Как видно по графикам, именно оценки ТЕП, рассчитанные в предположении г Ф 0, согласуются со значениями i , рассчитанными непосредственно по смоделированным данным.

Таблица 2. Матрицы коэффициентов корреляции и ранговой корреляции

TE1 TEr e-U rk(TE 1) rk(TE) rk(e-U)

TE1 1 rk(TE 1) 1

TEr -0.9930 1 rk(TE) -0.9985 1

e-U -0.9179 0.9283 1 rk(e-U) -0.9874 0.9885 1

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

4. Выводы

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

1.0

0

0.6

0.4

0.2 -

TE 1

0.0

0

0

0.2

0.4

0.6

0

1.2

6 Этот результат верен также в случае, если вместо оценок ТЕи использовать оценки технической эффективности, полученные по модели м2 (см. Приложение 3).

2. Введение показателя тесноты связи г между случайными величинами Ui и V при по- § строении нормальной копула-функции является адекватным инструментом выявления их ¡5 зависимости, позволяющим обнаружить наличие сильной корреляции. ^

3. Аналитический вид найденной функции плотности совместного распределения ве- ^ личин Ц и V может быть использован в дальнейших исследованиях со схожей тематикой. о Вычисленные производные логарифмической функции правдоподобия I в модели (3) могут <5 быть применены для нахождения точки максимума I при создании программного обеспе- Л чения, использующего другой подход к поиску точек экстремума (например градиентный ^ метод). ®

4. Предположение о независимости случайных величин Ц и Vi в условиях, когда эти слу- ^ чайные величины фактически являются статистически зависимыми, может приводить к су- <3

■ж—г ®

щественным искажениям оценок технической эффективности. При отсутствии «разумного»

обоснования независимости величин Vi и Ц целесообразность использования предпосылки ч об их независимости должна определяться целями исследования.

о

Список литературы

Айвазян С. А., Афанасьев М. Ю., Руденко В. А. (2012). Некоторые вопросы спецификации трех-факторных моделей производственного потенциала компании, учитывающих интеллектуальный капитал. Прикладная эконометрика, 27 (3), 36-69.

Айвазян С. А., Фантаццини Д. (2014). Методы эконометрики. Продвинутый курс с приложениями в финансах. М.: Магистр.

Благовещенский Ю. Н. (2012). Основные элементы теории копул. Прикладная эконометрика, 26 (2), 113-130.

Кучмент Л. С., Демидов В. Н. (2013). Об использовании теории копул для определения вероятностных характеристик весеннего половодья. Метеорология и гидрология, 4, 68-78.

Пеникас Г. И., Симакова В. Г. (2009). Управление процентным риском на основе копулы^ЛБСИ моделей. Прикладная эконометрика, 13 (1), 3-36.

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

Aigner D. J., Amemiya T., Poirier D. J. (1976). On the estimation of production frontiers: Maximum likelihood estimation of the parameters of discontinuous density function. International Economic Review, 17 (2), 377-396.

Amsler Ch., Prokhorov A., Schmidt P. (2014). Using copulas to model time dependence in stochastic frontier models. Econometric Reviews, 33 (5-6), 497-522.

Carta A., Steel M. F. J. (2012). Modelling multi-output stochastic frontiers using copulas. Computational Statistics & Data Analysis, 56 (11), 3757-3773.

Frees E., Valdez E. (1998). Understanding relationships using copulas. North American Actuarial Journal, 2 (1), 1-25.

Genest Ch., Remillard B. (2004). Test of independence and randomness based on the empirical copula process. TEST. An Official Journal of the Spanish Society of Statistics and Operations Research, 13 (2), 335-369.

Joe H. (1997). Multivariate models and dependence concepts. London: Chapman Hall.

Kumbhakar S. C., Lovell C. A. K. (2000). Stochastic frontier analysis. Cambridge University Press. Lai H. P., Huang C. (2013). Maximum likelihood estimation of seemingly unrelated stochastic frontier regressions. Journal of Productivity Analysis, 40 (1), 1-14.

Meeusen W., van den Broeck J. (1977). Efficiency estimation from Cobb-Douglas production functions with composed error. International Economic Review, 18, 435-444.

Rosenberg J., Schuermann T. A. (2006). General approach to integrated risk management with skewed, fat-tailed risks. Journal of Financial Economics, 79, 569-614.

Self S. G., Liang K.-Y. (1987). Asymptotic properties of maximum likelihood estimators and likelihood ratio tests under nonstandard conditions. Journal of the American Statistical Association, 82, 605-610.

Sklar A. (1996). Random variables, distribution functions, and copulas: Personal look backward and forward. Lecture Notes. Monograph series, 28, 1-14.

Smith M. D. (2008). Stochastic frontier models with dependent error components. The Econometrics Journal, 11 (1), 172-192.

Приложение 1

Вывод формулы (8).

c(ui, u2)=jsp exp(-2 (s-1 -

где £ = (Ф-1 (u1),Ф-1 (u2)) , S — ковариационная матрица, I — единичная матрица.

S =

/1 r\

r 1, V V

s12 =л/1-7, s-1 = —

1-r2

1 1 - rx

-r 1

V ' 1 /

тг r ^ 1 Г z,, ч ф((y-mVQy)-Ф(-тК)

u1 = Fv(x) = .— J e vdt, U2 = Fy(y) = ^- ' v-,

V2pqf r„ 1-Ф(-т^)

где Ф(у) — функция стандартного нормального распределения. c( Fv (x), Fy (y)) =

= Vrbrexp(-2 (Ф-1 (Fv (x)), Ф-1 (Fy (y)))T-(S-1 -I)( (Fv (x)), Ф-1 (Fy (y))) =

(r2 (Ф-1 (Fv (x)))2 - 2rФ-1 (Fv (х))Ф-1 (Fy (y)) + r2 (Ф-1 (Fy (y)))2)

r

1

V1-

2 exp r2

Заметим, что Ф 1 (Fv (х)) = х/ оу, f (х, у) = с( Fv (х), Fu (у))- ^ (х)- /и (у).

Подставив в плотность копула-функции вышеуказанные выражения, получим плотность совместного распределения V и и:

f (x, y) =

X exp

1

1

1

V1-r2 2psusv 1-Ф(-т/^и) 1

•exp

1 x2 1 (y-m)2

2 s2 2 sU

X

2(1- r2)

I 2 2 0 / x r 2rx -1 —:---Ф 1

V

\

Ф(( y -m)/Su )-Ф(-m/gu)

1—ф(—m/s u )

+

+r2

/ / Ф-

\ v

Ф(( y -m)/ s u )-Ф(-т^ u) 1^(-m/su)

w

2 \

I

Ф

I

«i

ai

(0 Ф

О

1

Следует отметить, что при г = 0 данная формула совпадает с совместной плотностью | распределения независимых случайных величин:

fHe3ae (Х,У)

1

1

2pSu Sv ^Ф^^)

•exp

1 x2 1 (y-m)2

2 sV 2 sU

QQ 'S «i

«i

d

Приложение 2

Введем дополнительные обозначения:

a(t) = Ф-

Ф((t -mV su )-Ф(-m|su) 1^(-m/su)

g,(t )=exp

1

2(1-r2)

(t+yt-z,)2r2 2ra(t)(t+yt -z,) 2 2/

s

s

-+r2 a2 (t)

(t+y -z,)2 (t-m)2

2s2

2s 2

ТогДа hj (y t) = fe(y t - z,.) =

VI- г2 2РСТи о> 1-Ф(-и/о-и )

Составим систему из (р + 5) уравнений:

f g, (t) dt.

1 n 00 / 00

^ = bs2 (1-r 2 ) 2 f g. (t)'(t + У - Z - ra(t)sv )dt f g,(t) dt Pos V(1 ' ) 1=1 L 0 /0

= 0,

1 no / 00

hk =s2 (1- 2) 2 x\ f gt(t)-(t + y, - Z, - ra(t)sv )t f gt(t)dt

sV(1 r ) ,=1 L 0 / 0

= 0 , j = 1,2,...,p,

n

с =2

"Vsv +f g,(t)-((t + y - z,)2 - ra(t)sv(t + y, - z,))dt /К(1-r2)f g(t)dt

= 0.

0

0

=1

Ku. =2

- +

jU-e

-u

,2 al

Or

JlPal (1-0(-u/ou ))

+

j g, (t )■

ra'au (t)(t + J, - г,.) r2a(t)aOu (t) it - u)

+

aF (1-r2)

1-r2

+

a

j g , (t) dt

ГДе aL(t) = -

U a2

(t-u)2

= 0,

•e2°u -(l-0((t-uVau)) + ((t-uVaU) 2°U -(l-0(-u/Ou))

_a(t)' 2

e 2 -(l-®(-u/a„))

k.=2

—e

-u

2 a2

/ g (t

raU (t)(t + J, - г,) r2a(t)aU (t) t

+

av (l-r2)

l-r2

UV/+ _U

+

a

u /

ГДе aU (t) = -

n

*;=2

,/2Pa„(l—ф(—uM,)) j g, (t) dt

0

-u2 (_-u)2

e2°u -(ф((_-uVau)-l) + e 24 '(l —Ф(—u/au))

a (t )2

Ou -e" 2 '(1 —Ф(—u/Ou ))

l-r2

\

— f (t) f r(t + J, - г)2 — a(t)(1 + r2)(t + j, - г,) + ra2 (t) j gi(t)'l aV(l-r2) - +

av (l-r2)

l-r2

(l-r2) J gt (t) dt

= 0,

= 0.

Нетрудно убедиться, что в случае г = 0 вид логарифма функции правдоподобия

п

I = ^1п ^ (у;) совпадает с видом, используемым в пакете Stаtа 10.0 для независимых ком-

1=1

понентов и и V:

п

Кеэас =Е1П ^ (У ) =

1=1 /

= — n ln

/ \\ u

+aV -Ф — | +2lnФ

\

\°u // 1=1

! 2 2 \ uav- (J - г)au

\ auavVa2 + aV /

i n —2 2

\2

j, - г, +u

vV°u+aV /

Подробный вывод логарифма функции правдоподобия для случая г = 0 можно найти, например, в (Kumbhаkаr, Lovell, 2000).

l

i=i

0

u

I

\

0

1=1

r

,=1

0

n

Приложение 3 §

ф

5

Приведем необходимые обозначения моделей из (Айвазян и др., 2012): ^

М,: R = Ь0Кь№Iье*-, где V ~ N(0;а^), и ~ (0;а2), 2

М2: R = ЬоКь 1Ь еу-и, где V ~ N(0; а2у), и ~ (т; а2), |

н-

а также гипотезы:

Н0: bi — 0 для некоторого i £ {1,2,3} (среди факторов производства существует незначимый фактор или фактор с отрицательным коэффициентом); «з

Н : > 0 для всех i £ {1,2,3} (все факторы производства значимы и имеют положи- |

00

тельные коэффициенты);

Н, : а2и = 0 (неэффективности нет); ^

Н : а^ > 0 (неэффективность присутствует);

\т = 0 ,

Н2: \ (оба параметрат и а2иравны нулю, т. е. в моделиМ2 нет неэффективности);

К =0

т* 0

2 (хотя бы один из параметров т и а2и не равен нулю, т. е. в модели М2 при-аи * 0

сутствует неэффективность);

Н, 2: т = 0 в модели М2 (неэффективность в моделяхМ, и М2 неразличима);

Н^2: т*0 в модели М2 (неэффективность в моделях М, и М2 различима). Следуя схеме спецификации, описанной в (Айвазян и др., 2012), получим оценки модели М,, которые приведены в соответствующем столбце табл. П3.1. Гипотеза Н0 отвергается. При проверке гипотезы Н2 об отсутствии неэффективности в моделиМ2 было получено значение статистики г = -3.074. Следовательно, гипотеза Н2 отвергается.

Переходим к построению модели М,, результаты оценки которой также можно увидеть в табл. П3.1. Гипотеза Н0 о незначимости отвергается на 10%-ном уровне. Гипотеза Н, об отсутствии неэффективности в модели М, отвергается.

Для проверки гипотезы Н12 необходимо вычислить значения исправленных критериев Акаике:

Л1Сс(М,) = 2 • 6 - 21пL(M1) + 2'6'(6 +1) = 54.713 ,

80 - 6 -,

Л1Сс(М7) = 2• 7 -21пL(М2) + 2 •7•(7 +1) = 55.98 .

80 - 7 -,

Таким образом, гипотеза Н12 не отвергается. Итоговой моделью является М,:

^ = е094К^^227110Л43е^-и', V ~ N(0; 0.0,8), и ~ N+ (0; 0.257).

Таблица П3.1. Результаты оценок моделей на изучаемых смоделированных данных

м2 m1

Оценки факторов производства

lnK 0.657*** 0.648***

lnL 0.219*** 0.227***

Inl 0 144*** 0.143***

const 0 87*** 0 94***

Оценки параметров компонентов ошибки

m -1.569 0

oVv 0.024 0.018

oU 0.73 0.257

Р-значение при проверке 0.001 0.000

отсутствия неэффективности

Число наблюдений 80 80

Логарифм функции правдоподобия - 20.3072 - 20.7811

Примечание. В скобках указаны стандартные ошибки; *** — значимость оценок коэффициентов на 1%-ном уровне. Значимости параметров компонентов ошибки в модели M1 и дисперсий в модели M2 не оценивались.

Приложение 4

Согласно обозначениям, введенным в Приложении 2,

1-Ф

TEU = E (e-u|e, ) =

оU ovjo2v +о2и j

1-Ф

\ovsToV + oU J

•exp

'oV oU+Ve, oUx v(o V +o U)

(более подробный вывод TE1i можно найти в (Kumbhakar, Love11, 2000)). Тогда

/(*+ У, - г,, *К

TE, = E (e-u|e, ) = / e-t • f (t|e,) dt = f e

1__1_ 1

л/l- r2 Vnovot

fe (J, - ^ ) 1 00 ■Г-ФЫ^ f (t) dt_ f e"g,(t) dt

л/1- r2 V pou ov 1-Ф (-m/o U )

f g, (t) dt f g, (t) dt

I

\

1

0