2) Амплитуда колебаний индивида вокруг своего центра масс (mob_2).
Для расчета данного показателя мобильности для начала необходимо найти координаты центра масс индивида по следующей формуле:
(10)
где: и - широта и долгота центра;
и - широта и долгота базовой станции i;
- вес базовой станции i (доля издержек, приходящихся на базовую станцию i в общих расходах индивида).
Далее рассчитывается сам показатель мобильности индивида j:
(11)
где: - показатель мобильности индивида j;
- расстояние от базовой станции i до центра масс в метрах;
- вес базовой станции i.
Для нахождения расстояния от базовой станции i до центра масс будет использоваться формула гаверсинуса (Миловзоров и др., 2014):
(12)
где: - широта и долгота центра масс в радианах;
- широта и долгота базовой станции i в радианах.
Вышеописанные два способа оценивают, насколько мобилен индивид в плане территориального передвижения. Помимо них также будет рассчитан индекс Херфиндаля-Хиршмана в качестве показателя мобильности.
3) Индекс Херфиндаля-Хиршмана (HHI).
Индекс Херфиндаля-Хиршмана - это показатель уровня концентрации рынка, который определяется как сумма квадратов долей рынка, занимаемых каждым его участником (Rhoades, 1993). В данной работе индекс Херфиндаля-Хиршмана будет рассчитан как сумма квадратов процентных долей затрат, приходящихся на каждый контакт абонента:
(13)
где: - процентная доля затрат абонента, приходящаяся на контакт i.
Чем выше данный показатель, тем более сконцентрированы разговоры абонента в пространстве, то есть он менее мобилен. Все три показателя мобильности интерпретируются в зависимости от базовой альтернативы: при увеличении мобильности индивид более или менее вероятно выберет альтернативу j по сравнению с базовой альтернативой. Следовательно, при одинаковом распределении мобильности на разных тарифных планах включение данной переменной нецелесообразно. Для начала взглянем на гистограммы распределений трех показателей мобильности.
Из гистограммы распределения переменной mob_1 видно (Рисунок 3), что анализируемые тарифные планы действительно отличаются друг от друга и не возникает сомнений в целесообразности включения данной переменной в модель.
По гистограмме распределения переменной mob_2 (Рисунок 4) сложно сделать вывод об однородности выбранных тарифов: распределения переменной на разных тарифах практически дублируют друг друга с небольшими различиями.
Рис. 3 Распределение переменной mob_1 по тарифам
Рис. 4 Распределение переменной mob_2 по тарифам
Рис. 5 Распределение переменной HHI по тарифам
Из гистограммы распределения переменной HHI явно видно, что выбранные для исследования тарифы не являются однородными по данной переменной. Тем не менее, для точного ответа на вопрос об однородности тарифов далее необходимо использовать критерий . Ниже представлены результаты проведенного теста.
Таблица 2
Тест на однородность
|
Показатель |
Вычисленное значение |
p-value |
Табличное значение |
|
|
mob_1 |
1999,122 |
0,1 |
161,827 |
|
|
mob_2 |
307,5969 |
0,05 |
168,613 |
|
|
HHI |
1360,9 |
0,01 |
181,84 |
Нулевая гипотеза об однородности распределений показателей мобильности по тарифам отвергается на всех уровнях значимости, что говорит о целесообразности включения данных переменных в модель.
Таким образом, после оценки моделей и сравнения их между собой представляется возможным проверить выдвинутые в работе гипотезы и сделать выводы о структуре выбора тарифного плана и влиянии рассматриваемых переменных на вероятность переключения на другой тариф.
Данные
В выборке изначально было 2059 абонента восьми различных тарифных планов за период с января по ноябрь 2012 года включительно: «Знатоки общения» (tariff1), «Капитал Platinum» (tariff2), «Кнопкины точка ру» (tariff3)", «Пионер» (tariff4), «Пять звезд» (tariff5), «Радуга Фристайл» (tariff6), «Своя страна +» (tariff7), «Союз» (tariff8), - но после удаления выбросов осталось 1154 абонента (7134 наблюдения выбора конкретного тарифа), среди которых 101 абонент сменил тарифных план в рассматриваемый период. Ниже представлено распределение абонентов по тарифам.
Рис. 6 Распределение абонентов по тарифам
Как было уже сказано, в качестве alternative-specific факторов, влияющих на вероятность смены тарифного плана, были выбраны наиболее популярные услуги:
- стоимость звонков внутри сети;
- стоимость звонков на номера других сотовых операторов;
- стоимость звонков на номера фиксированной телефонной связи;
- стоимость SMS сообщений на номера других сотовых операторов;
- стоимость SMS сообщений внутри сети.
В силу того, что в предоставленных данных отсутствуют демографические характеристики абонентов, то помимо характеристик тарифного плана в модель также будет включена такая характеристика индивида, как его мобильность в качестве объясняющей переменной. Таким образом, функция полезности индивида i на тарифе j будет выражаться следующим образом:
(14)
где:- неизвестные параметры объясняющих переменных;
- стоимость звонков внутри сети индивида i на тарифе j;
- стоимость звонков на номера других операторов индивида i на тарифе j;
- стоимость звонков на номера фиксированной связи индивида i на тарифе j;
- стоимость исходящих SMS сообщений на номера других операторов индивида i на тарифе j;
- стоимость исходящих SMS сообщений внутри сети индивида i на тарифе j;
- мобильность индивида i.
Так как выбранные для исследования модели предполагают специальную структуру данных, то для каждого абонента необходимо рассчитать его действительные месячные издержки по каждой выбранной услуге, а также альтернативные издержки, то есть сколько бы он потратил на каждую услугу, если бы выбрал любой другой доступный тарифный план.
Данный этап работы является одним из наиболее сложных и трудоемких, так как на основе информации о прейскурантах тарифов (Приложение 1) для каждой альтернативы нужно сконструировать формулы для расчета издержек. Если тарификация на тарифе является линейной, то есть стоимость любой единицы услуги одинакова, то сложностей с расчетом ежемесячных затрат не возникает. Однако для тарифов, предусматривающих дифференциацию стоимости по накопленной сумме минут, на первом шаге необходимо рассчитать сумму накопленной в течение суток или месяца длительности как по всем направлениям звонков, так и для каждого направления отдельно. Аналогично необходимо поступить с SMS сообщениями, и далее, на втором шаге, при расчете затрат использовать в формуле рассчитанные накопленные суммы. Что качается тарифов с абонентской платой, то фиксированную плату надо распределить по каждой услуге в зависимости от доли совокупных затрат, приходящуюся на услугу. Это можно объяснить тем, что, если индивид в основном совершает звонки и редко пользуется SMS сообщениями, то он рассматривает абонентскую плату как оплату за голосовые сообщения. Описательные статистики выбранных переменных представлены в таблице:
Таблица 3
Описательные статистики alternative-specific переменных
|
Переменная |
Среднее значение |
Стандартное отклонение |
Минимальное значение |
Максимальное значение |
|
|
стоимость звонков внутри сети (руб./мес.) |
196,4385 |
158,9844 |
0 |
1188,75 |
|
|
стоимость звонков на номера других сотовых операторов (руб./мес.) |
144,179 |
163,9143 |
0 |
1277,604 |
|
|
стоимость звонков на номера фиксированной телефонной связи (руб./мес.) |
24,79228 |
36,67959 |
0 |
428 |
|
|
стоимость SMS сообщений на номера других сотовых операторов (руб./мес.) |
33,36061 |
64,5484 |
0 |
663 |
|
|
стоимость SMS сообщений внутри сети (руб./мес.) |
23,80022 |
52,15155 |
0 |
622,8 |
На основе формул, описанных в разделе 3, был произведен расчет всех трех показателей мобильности, описательные статистики которых представлены ниже.
Таблица 4
Описательные статистики case-specific переменных
|
Переменная |
Среднее значение |
Стандартное отклонение |
Минимальное значение |
Максимальное значение |
|
|
mob_1 |
24,95522 |
23,30009 |
1 |
209 |
|
|
mob_2 |
13455,78 |
40792,69 |
3,72e-10 |
937888,1 |
|
|
HHI |
1462,952 |
1153,535 |
84,24309 |
10000 |
Также следует отметить, что не все выбранные тарифные планы были доступны для подключения на протяжении всего рассматриваемого периода, что тоже должно быть учтено в структуре данных. Так, например, тарифный план «Знатоки общения» доступен лишь с января по май включительно. Единственным тарифом, который доступен на протяжении всего периода, является «Капитал Platinum», именно поэтому он будет выбран в качестве базовой альтернативы при построении модели, то есть интерпретация коэффициентов модели будет осуществлена в сравнении с данным тарифным планом. Даты начала и закрытия тарифов представлены на диаграмме ниже.
Рис.7 Период действия тарифов
Результаты
В работе были оценены три типа модели: Conditional, Mixed и Nested logit. Conditional logit модель позволяет понять, как изменение стоимости услуг на тарифном плане и изменение личностных характеристик (мобильности) влияет на вероятность смены абонентом текущего тарифного плана.
Далее была оценена Mixed logit модель, которая учитывает гетерогенность людей по чувствительности к изменению стоимости услуг. Это позволяет коэффициентам при alternative-specific регрессорах в данной модели быть не фиксированными, что является реалистичным, так как чувствительность к изменению цены звонков внутри сети для двух абонентов, один из которых в основном совершает внутрисетевые звонки, а другой, как правило, отправляет SMS сообщения и редко совершает звонки, будет разной.
Наконец, последняя модель - Nested logit, позволяет проверить гипотезу о структуре выбора тарифного плана. Логично предположить, что в процессе подбора тарифа индивид вначале неким образом классифицирует для себя множество тарифов и лишь затем делает выбор из той группы тарифов, которая является для него наиболее оптимальной. В данной работе была смоделирована следующая классификация тарифов: индивид группирует тарифы по их особенностям (наличие или отсутствие дифференциации стоимости, наличие абонентской платы), предполагая, что минимальный выгодный уровень потребления в этих группах отличается.
В конечном итоге, оценив все три типа моделей и сравнив их между собой, представляется возможным опровергнуть или подтвердить гипотезу о двухшаговой структуре выбора, а также проверить другие выдвинутые гипотезы. Для оценивания моделей использовался статистический пакет Stata. Оценивание Conditional и Nested logit моделей возможно лишь на основе кросс-секционной структуры данных, поэтому вначале было осуществлено оценивание моделей и последующее их сравнение на кросс-секционных данных, затем - оценивание Mixed logit модели, которая позволяет учесть панельную структуру данных. Оценки моделей, где за мобильность было взято количество уникальных базовых станций, представлены в таблице:
Таблица 5
Оценки моделей с включением переменной mob_1
|
Переменная |
Кросс-секционные данные |
Панельные данные |
|||
|
Conditional |
Mixed |
Nested |
Mixed |
||
|
Price1 |
-0,0053470*** |
-0,0055223*** |
-0,0073411*** |
-0,0141344*** |
|
|
Price2 |
-0,0039514*** |
-0,0050717*** |
-0,0073411*** |
0,0007602 |
|
|
Price3 |
-0,0107695*** |
-0,0147341*** |
-0,0223226*** |
-0,0117466*** |
|
|
Price4 |
-0,0058523*** |
-0,0057923*** |
-0,0150540*** |
-0,0050489 |
|
|
Price5 |
-0,0116916*** |
-0,0163181*** |
-0,0199513*** |
-0,0219243*** |
|
|
tariff1 |
2,9530000*** |
6,8641540*** |
|||
|
tariff3 |
-1,4137140*** |
-0,3012990 |
|||
|
tariff4 |
3,2412100*** |
6,5654170*** |
|||
|
tariff5 |
2,1419920*** |
6,1297310*** |
|||
|
tariff6 |
3,4161040*** |
7,3364490*** |
|||
|
tariff7 |
1,5828770*** |
6,2751710*** |
|||
|
tariff8 |
0,8715456*** |
1,5408780*** |
|||
|
mob_1*tariff1 |
-0,0666786*** |
-0,0680429*** |
-0,0923735*** |
||
|
mob_1*tariff3 |
-0,0296767 |
-0,0441697 |
-1,0848410*** |
||
|
mob_1*tariff4 |
-0,0241532*** |
-0,0256336*** |
-0,1041617*** |
||
|
mob_1*tariff5 |
-0,0316000*** |
-0,0328275*** |
-0,0466105*** |
||
|
mob_1*tariff6 |
-0,0560851*** |
-0,0575262*** |
-0,0487691*** |
||
|
mob_1*tariff7 |
-0,1389196*** |
-0,1391655*** |
-0,5409510*** |
||
|
mob_1*tariff8 |
-0,0249343*** |
-0,0281077*** |
-0,0382091*** |
||
|
SD Price1 |
0,0011947 |
0,0492535*** |
|||
|
SD Price2 |
0,0059649*** |
0,0716980*** |
|||
|
SD Price3 |
0,0141959*** |
0,2408872*** |
|||
|
SD Price4 |
0,0005480 |
0,1323511*** |
|||
|
SD Price5 |
0,0152434*** |
0,2074835*** |
|||
|
one-part_with_diff_mob_1 |
-0,0353631*** |
||||
|
two-part_mob_1 |
-0,0221268*** |
||||
|
one-part _tau |
-0,0221268*** |
||||
|
with_diff_tau |
-16,598660*** |
||||
|
with_flat_rate_tau |
1,2141080*** |
||||
|
p-value: * p < 0,10; ** p < 0,05; *** p < 0,01 |