CC BY
79
СТОХАСТИЧЕСКАЯ МОДЕЛЬ
ЦЕНООБРАЗОВАНИЯ НА РЫНКЕ НЕДВИЖИМОСТИ: ФОРМИРОВАНИЕ ЛОГНОРМАЛЬНОЙ ГЕНЕРАЛЬНОЙ СОВОКУПНОСТИ
УДК 339.13.017
Олег Витальевич Русаков,
к.ф-мат.н., доцент кафедры теории вероятностей и математической статистики Санкт-Петербургского государственного университета Тел.: (812) 428-42-12 Эл. почта: [email protected]
Михаил Борисович Ласкин,
к.ф-мат.н., доцент, директор ООО «Инвест-Строй» Тел.: (812) 315-50-92 Эл. почта: [email protected]
Ольга Ильинична Джаксумбаева,
к.э.н., Санкт-Петербургский государственный университет, Экономический факультет
Тел.: (812) 273-02-76
Эл. почта: [email protected]
Мы строим стохастическую модель ценообразования на рынке вторичной недвижимости. Метод построения модели основывается на последовательных сравнениях цен. Мы доказываем, что при стандартных условиях на коэффициенты сравнения единственный невырожденный предел по распределению для цены имеет закон логнормального распределения вероятностей. Соответствие эмпирических распределений цен теоретически доказанному логнор-мальному закону мы подтверждаем на многочисленных данных статистики по недвижимости Санкт-Петербурга. При этом мы в основном используем мощный и чувствительный статистический критерий согласия Колмогорова-Смирнова. Опираясь на «Федеральный стандарт оценки №2», мы делаем вывод, что рекомендуемая наиболее вероятная цена - мода распределения корректно и однозначно определяется при аппроксимации логнормальным законом. Так как у логнормального распределения среднее превосходит наиболее вероятное значение - моду, то цены, вычисленные по математическому ожиданию, систематически завышены.
Ключевые слова: рыночная стоимость недвижимости, стохастическая модель ценообразования, геометрическое броуновское движение, мода логарифмически нормального закона распределения, применение критерия согласия Колмогорова-Смирнова.
Введение
При определении рыночной стоимости объектов недвижимости методом сравнительного подхода возникают следующие проблемы как теоретического, так и практического характера: выбора объектов сравнения, формирования достаточно представительной выборки, построения оценки рыночной стоимости на основе статистических данных.
Однако, в практике оценивания, как правило, выборки не являются представительными, содержат малое количество объектов сравнения, а проверка статистических гипотез, как справедливо отмечалось в [1], часто не производится.
В Федеральном Стандарте Оценки № 2 [2] в п.6 указано следующее положение: «При определении рыночной стоимости определяется наиболее вероятная цена, по которой объект оценки может быть отчужден на дату оценки на открытом рынке в условиях конкуренции, когда стороны сделки действуют разумно, располагая всей необходимой информацией, а на величине сделки не отражаются какие-либо чрезвычайные обстоятельства».
Аналогичные формулировки, в основе которых фигурирует «наиболее вероятная цена», содержат и многие признанные зарубежные стандарты оценки, такие как: IVS (p.30 a) [3], TEGOVA (EVS p. 5.3.1) [4], USPAP (Standard rule 6-2, p. c.) [5], RICS (p. 3.2.1) [6].
С математической точки зрения это означает, что определение рыночной стоимости, заданное стандартами, опирается на факт вероятностной природы цены: цена - величина случайная. При этом сама рыночная стоимость является числовой характеристикой распределения случайной величины (цены), её модой и показывает наиболее вероятное значение случайной величины.
Абсолютное большинство используемых в современной практике методов оценивания рыночной стоимости при статистической обработке имеющегося материала ориентированы на нормальное распределение, что плохо согласуется с эмпирическими распределениями цен. Так, в частности, одной из особенностей нормального распределения являются его симметричность, а значит, в частности, и равенство характеристик: моды, математического ожидания и медианы. Между тем, цены на однородные активы никогда не распределены нормально. Это следует хотя бы из того, что цены положительны, а нормально распределенная случайная величина со строго положительной вероятностью принимает отрицательные значения. Более того, распределение цен в большинстве случаев асимметрично влево.
В настоящей статье, при достаточно простых и естественных предположениях, устанавливается теоретический факт слабой сходимости (сходимости по распределению: в одномерном случае - слабая сходимость равносильна сходимости по распределению; однако, в данной статье в некоторых местах мы выходим за рамки одномерного случая, например, - в ключевой Лемме 1) процесса последовательных сравнений цен к логарифмически нормальному распределению. Получены статистические результаты, подтверждающие данный факт. Ряд ключевых примеров - приводится. Устанавливается ряд особенностей рыночной стоимости, которые вытекают из логарифмически нормального распределения цен. Рассмотрены примеры, основанные на статистическом материале по ценам предложений вторичного рынка жилой недвижимости в г. Санкт-Петербурге (СПб) по состоянию на 31.03.2014 г.
Oleg V. Rusakov,
PhD in physic and mathematic, Associate professor, Theory Probability and Mathematical Statistics department, Saint-Petersburg State University Tel.: (812) 428-42-12 E-mail: [email protected]
Michael B. Laskin,
PhD in physic and mathematic, Associate-professor, Director OOO "Invest-Stroy", Saint-Petersburg Tel.: (812) 315-50-92 E-mail [email protected]
Olga I. Jaksumbaeva,
PhD in economics, Economics department, Saint-Petersburg State University Tel.: (812) 273-02-76 E-mail: [email protected]
STOCHASTIC PRICING MODEL FOR THE REAL ESTATE MARKET: FORMATION OF LOG-NORMAL GENERAL POPULATION
We construct a stochastic model of real estate pricing. The method of the pricing construction is based on a sequential comparison of the supply prices. We proof that under standard assumptions imposed upon the comparison coefficients there exists an unique non-degenerated limit in distribution and this limit has the lognormal law of distribution. The accordance of empirical distributions of prices to the theoretically obtained log-normal distribution we verify by numerous statistical data of real estate prices from Saint-Petersburg (Russia). For establishing this accordance we essentially apply the efficient and sensitive test of fit of Kolmogorov-Smirnov. Basing on "The Russian Federal Estimation Standard N2", we conclude that the most probable price, i.e. mode of distribution, is correctly and uniquely defined under the log-normal approximation. Since the mean value of log-normal distribution exceeds the mode - most probable value, it follows that the prices valued by the mathematical expectation are systematically overstated.
Keywords: real estate market value, stochastic model of pricing, geometric Brownian motion, mode of the logarithmically normal law distribution, applications of the Kolmogorov-Smirnov test of fit.
1. Стохастическая модель ценообразования
Мы строим стохастическую модель ценообразования на рынке жилой недвижимости, основанную на последовательном сравнении цен объектов недвижимости. Коэффициенты сравнения мультипликативны, малы, случайны, имеют простое распределение (элементарны и однородны в определенном смысле). Мы показываем, что при достаточно большом числе сравнений N (начиная с 500, примерно, для ошибки аппроксимации порядка 5%, т.к. таковая ошибка обратно пропорциональна корню из N) при стандартных предположениях о малости и однородности коэффициентов сравнения с высокой точностью формируется логнормальное распределение цены (имеется в виду удельная цена - за 1 кв. м.). Мы полагаем, что выполнены следующие простые и стандартные предположения:
а) об однородности коэффициентов сравнения то есть величин, показывающих во сколько раз цена следующего (ближайшего) объекта сравнения больше или меньше цены предыдущего;
б) о независимости от предыдущих шагов сравнения увеличения/ уменьшения цены. При выполнении данных предположений мы получаем приближение для распределения удельной цены объекта недвижимости к распределению вероятностей логарифмически нормального закона. Полученное приближение следует понимать в смысле предела по распределению (слабой сходимости) цен при последовательном случайном сравнении, когда число объектов сравнения N ^ да.
Определение 1. Мы говорим, что случайная величина X имеет логнормальное распределение (lognormal distribution), если ее логарифм нормально распределен, то есть она допускает следующее представление (в виде равенства по распределению) d
X = exp^ + aN(0,1)}, (1)
здесь и далее:
вещественный параметр ц - математическое ожидание случай-
ной величины lnX; положительный параметр а1 - волатильность (volatility) - дисперсия распределения lnX, а - соответствующее стандартное отклонение; и в то же время: exp{/} - медиана распределения величиныX; exp{/ + а1 / 1} - математическое ожидание величины X; exp{/ - а1} - мода распределения величины X!
Обозначение N(/, а1) мы используем как для закона нормального распределения со средним /, дисперсией а1, так и для соответствующей
d
случайной величины. Значок = обозначает равенство по распределению.
Введем неслучайную величину S0 > 0, которой мы придаем смысл цены (цена за один квадратный метр) начального объекта сравнения. Запишем случайную величину X в виде
X = Soexp{(M - lnSo) + (1) + oN(0, 1)}.
Введем параметры «отношение «сигнал/шум», которые также будем называть относительным риском. По сути - это принятые модификации известного коэффициента Шарпа (Sharpe parameters) [7].
R = R(ß, S0,a) =
A ß- ln S0
a
?ч A Ц- ln S0 R2 = R2(p,S0,a2) =И 2 0 . (3) а
Величину ц - S0 мы интерпретируем, как полезный сигнал.
2. Процесс последовательных сравнений и его пределы по распределению
Процесс последовательных сравнений реализуется по следующей схеме. Выбирается некоторый начальный объект сравнения с ценой ЗЬ. Выбор предполагается возможно случайным, а если неслучайным, то не из крайних значений спектра цен совокупности объектов сравнения. Переход к следующему объекту сравнения предполагается таким, что цена может увеличиться или уменьшится незначительно. Мы рассматриваем бинарную схему, лежащую в основе последовательности испытаний Бернулли. Цена следующего объекта сравнения определяется по формуле
= ЗД,
(4)
где У1 - случайная величина, принимающая два значения U > 1 (Up) и0 < D < 1 (Down) с соответствующими вероятностями 0 < p < 1 и q = 1 - p. При этом предполагается, что величины U и D незначительно отличаются от единицы (в большую и меньшую стороны, соответственно). Переход к следующим объектам происходит по рекуррентному правилу, причем каждый раз умножение происходит на независимый коэффициент, и все коэффициенты одинаково распределены. Для того, чтобы доказать логнормальность предельного распределения данной рекуррентной процедуры, мы будем доказывать нормальность предела сумм логарифмов независимых одинаково распределенных коэффициентов сравнения (У). Вследствие свойства непрерывности экспоненты и ограниченности траекторий броуновского движения на замкнутом интервале [8] будет достаточно этих оснований для установления необходимого факта логнормальности предельного закона распределения цен. Важно отметить, что вероятности p и q зависят от выбора начальной цены сравнения S0. Мы изначально предполагаем, что существует некоторый спектр (интервал) цен объектов сравнения, который зависит от количества объектов N (вовлеченных в процесс последовательных сравнений). Число объектов, вовлеченных в процесс сравнения, должно быть соизмеримо с числом объектов всей существующей, или предполагаемой совокупности и они должны быть распределены однородно по совокупности.
В случае, когда цена So находится в начале спектра цен объектов сравнения, то вероятность взять следующий объект более дорогим (то есть вероятность p события {Si > S0}) превосходит вероятность q взять следующий объект более дешевым (события {S1 < S0}). При этом - чем ниже (левее в спектре) цена S0, тем больше вероятность p превосходит вероятность q. Таким образом, мы предполагаем зависимость p = p(S0), q = q(S0). Также заметим, что в случае, когда выбор начальной
цены сравнения 50 случаен, возникает ее зависимость от количества объектов Ы, и, следовательно, появляется зависимость вероятностей р и q от N. Следующее предположение для упрощения построения конструкции процесса сравнения - это «логарифмически симметричные» коэффициенты сравнения, то есть
D = 1.
U
(5)
j=i
+ (n -[nt])£м+1 + tC,
(6)
Разумеется, что существует взаимно-однозначное соответствие между ломаными (£„) и последовательностями <2, ..., <п. Это соответствие задается формулой
£t+1 = g„\ — I-g„ I - I-C„
п ) V „) „
k = 0, 1,
, n - 1.
(7)
Для дальнейших построений будем опираться на следующую стохастическую лемму «общего характера» [9]. Здесь эта лемма приводится в более слабом виде, чем в первоисточнике, но в достаточном для поставленных в настоящей работе целей. Приводится также и необходимое доказательство.
Используем следующие стандартные обозначения математического анализа: при п ^ да функция о -«о малое», функция О - «О большое».
Пусть случайные величины {<}, ] е К, независимы, одинаково распределены и имеют распределение, зависящее от некоторого натурального п, произвольного, но фиксированного. Случайная величина < вводится с целью задания закона распределения.
Дальнейшая интерпретация -следующая: случайные величины вида (<) представляют собой соответствующие центрированные (то есть, после вычитания соответствующих математических ожиданий) логарифмы коэффициентов сравнения вида 1^, и этих коэффициентов сравнения рассматривается п штук. Константа СП, введенная в следующей ниже формуле (6) - есть сумма математических ожиданий логарифмов (накопленное математическое ожидание логарифмов) коэффициентов сравнения вида 1пУ
По последовательности <1, <2, ..., <п построим непрерывную случайную ломаную следующего вида:
д Н () = ^ +
где t е [0, 1], а Cn - некоторая константа, зависящая только от n.
Лемма 1. Пусть для всех натуральных п и для некоторого 8 > 0 конечно математическое ожидание Е(п)| <да. Пусть при п ^ да случайные величины < = ¿¡(п),} = 1, ..., п; и константы Сп, ап > 0удовлетворяют следующим условиям: (О Е<1(п) = 0;
а1 гт2 =
(О) Б%(п) = = —(1 + а(1)); п п
(Ш) Е\^(п)\г = 0(пг/2), г > 0;
(/V) Сп = С (1 + О (1)), где С и а > 0 - некоторые постоянные, не зависящие от п. Тогда конечномерные распределения случайных ломаных gn(t) при п ^ да слабо сходятся к конечномерным распределениям случайного процесса
аWt + а, t е [0, 1], (8)
где Wt - стандартное броуновское движение. В пункте (т) не предполагается существования момента произвольного порядка. Следует понимать так, что для всякого г > 1, такого что Е\<5(п)|г < да имеет место условие (Ш). Естественно, здесь подойдут все г < 2 + 8.
О слабой сходимости более подробно смотри напр.: [10] для случайных величин, [8] для векторов и функций.
Доказательство Леммы 1. Обозначим /(х) характеристическую функцию (преобразование Фурье меры) для распределения случайной величины Прежде всего заметим, что в силу условий (/) - (Ш) при всяком фиксированном х справедливо разложение
1п I (х)=-1 - х2+;(11, 2 п V п )
п ^ да.
Отметим, что подобного рода разложение используется, как правило, при доказательстве методом
характеристических функций центральной предельной теоремы в форме Леви (см. [10]). Рассмотрим некоторое произвольное (не обязательно равномерное) разбиение временного промежутка [Ь, 1], заданное посредством точек Ь = /Ь < ••• < ^ < 1. Вычислим распределение случайного вектора, составленного из значений в этих точках случайной ломаной gn: д
^ = (т) =
= ( £п (Ч ) > £п (*2 ) >•••> £п (т )) > (1Ь)
где gn(t) определены посредством равенства (6). В дальнейшем, при доказательстве настоящей леммы, опустим индекс п при обозначении случайной функции (случайной ломаной) gn(t) = g(t). Обозначим через ближайший слева к точке ^ узел равномерного разбиения на п частей интервала времени [Ь, 1]:
* Л Г 1 1 I 11
о]-=-К-}-' (11)
L (tj) - g (t* )|——0
n ^ да.
(13)
* * а / * * \а Ф = p'„(m) = ( g (AX-, g (С) ) =
= [gI ki-I'•••'gI km -
где к = [п/],] = 1, ..., т. Рассмотрим случайный вектор у, составленный из приращений вектора ф*:
у(п) = у = (У1, ..., Ут) (15) с компонентами
= * [j J"* [j 1
Очевидно, что
(16)
VJ
= i $ + Сп (kj - j) -. (17)
=к- +1
здесь {а} обозначает дробную часть числа а. Из (6) нетрудно вывести равенство
g (tj) - g (tj) =
Так как первый сомножитель в правой части этого равенства равномерно по п ограничен единицей, а второй сомножитель имеет стремящиеся к нулю (при п ^ да) дисперсию и среднее, то использование классического неравенства Чебышева легко устанавливает следующую сходимость по вероятности
Кроме того, (у1, ..., ум) - независимые случайные величины, и даже более - каждая компонента вектора у состоит из суммы независимых случайных величин. Поэтому для всякого j = 1, ..., м характеристическая функция fj случайной величины yj допускает представление
fj (х; п) = fj (х) = = expi^ixCn(kj -kj-1)П} f'x),
x e R, (18)
где f(x) - характеристическая функция случайной величины . Так как
lim( kj - kj -)X- = tj - tj,
j = 1, ., m (19)
то при каждом фиксированном вещественном x для характеристической функции fj(x) при каждом j = 1, ., m справедливо предельное соотношение
lim fj (x; n) =
= exp {ixC(tj - tj-1)}x
а
x exp —(t -t -i)x
(20)
Таким образом, если у вектора ф существует предельное распределение, то оно совпадает с предельным распределением вектора
и ковариационной матрицей
8ij6L(tj - t;_i), i, j, ..., m, (H)
где dij = 0 при i ф j и Sy = 1 при i = j -символ Кронекера.
Что же касается векторов ф
*
и ф , то их распределения совпадают с точностью до бесконечно малого
члена (при n ^ да). В свою очередь,
*
вектор ф получается из вектора у суммированием независимых компонент последнего, поэтому нетрудно
вычислить общее предельное рас*
пределение векторов ф и ф :
гауссовское распределение с вектором средних (ti, ..., tm)C и ковариационной матрицей CT1min{ti, tj}, i, j, ., m.
Эти характеристики полностью соответствуют известным характеристикам конечномерных распределений винеровского процесса с коэффициентом масштаба а и со сдвигом Ct. Таким образом, мы установили необходимый нам факт сходимости конечномерных распределений:
Конечномерные распределения процесса, заданного случайными ломаными gn(t), сходятся при n ^ да к конечномерным распределениям процесса
aWt + Ct, t £ [0, 1],
(23)
(14)
Здесь мы использовали предельное соотношение для последовательности {Сп} из условия Леммы 1, а также представление для характеристической функцииДх) случайной величины . Таким образом, по теореме о непрерывном соответствии между характеристическими функциями и распределениями следует, что распределения случайного вектора у сходятся к т-мерному гауссовскому распределению с вектором средних
(/1 - /ь, /2 - /1, ..., /] - /}-1)С (21)
где - стандартный винеровский процесс. Лемма доказана.
Теперь заметим, что распределение винеровского процесса в единице - стандартно нормально N(0, 1). Таким образом, нам достаточно в формулировке леммы рассмотреть момент времени / = 1, когда количество наших объектов сравнения N совпадает с параметром п из предыдущей леммы. Тогда мы получим сходимость (слабую сходимость) распределений логарифмов цен к распределению нормального закона Щ(С, а2), так как распределение случайной величины аЖ\ + С нормально со средним С и дисперсией ст2. Лемма 1 устанавливает факт сходимости по распределению к нормальному закону для логарифмов цен. Однако, нам нужна сходимость распределений собственно цен к логнормальному закону. Следующее Следствие 1 доказывает, что экспоненты от сумм (■ сходятся по распределению к логнормальному
2
п
п
закону. Рассмотрим последовательность при натуральном п случайных величин (Яи(1)), которые являются концами случайных ломаных, что определены равенством (6).
Следствие 1 к Лемме 1. Последовательность (ехр^п(1)}) сходится по распределению при п ^ да к логнормальной случайной величине ехр{С + аИ(00, 1)}.
Доказательство Следствия 1. Введем обозначения Оп = gn(1), при всех п е К; Ф(и) =
^ехр<|-^и2 |<и - функция распределения стандартного нормального закона, и е Я. Заметим, что сходимость по распределению равносильна сходимости функций распределений во всех точках непрерывности функции распределения предельного закона (см. [8]). При этом отметим, что функции распределения как нормального закона, так и логнормального, непрерывны на всей оси и на правой полуоси, соответственно. Из Леммы 1 имеем для любого вещественного а при п ^ да, Р(Оп < а) ^ Ф(а).
Так как экспонента строго монотонно растущая функция и применение ее к обеим частям сохранит неравенство, а вероятность эквивалентных неравенств равна, то Р(Оп < а) = Р(ехр(Оп) < ехр (а)) ^ ^ Ф(а). Обозначим хо = ехр(а).
Тогда при п ^ да, Р(ехр(Оп) < < х0) ^ Ф(1пх0). При повторном использовании свойства экспоненты строго монотонно расти, легко видеть, что Ф(1пхо) = Р(МР,1) < 1пхо) = = Р(ехр(МР,1)) < хо). В силу произвольности а е Я и взаимно-одно-значности логарифма, как функции из правой полуоси в Я, окончательно имеем требуемое в Следствии 1 утверждение. Следствие 1 доказано.
Замечание к Лемме 1. Параметр t е [0, 1] в формулировке Леммы естественно интерпретируется, как (непрерывное) время действия процесса сравнения. Этот параметр легко обобщается на случай t е [0, Т] для произвольного Т > 0. При вовлечении всё большего числа объектов сравнения в процесс сравнения при выполнении условия однородности мы получаем процесс
Броуновского движения для логарифма цены с параметром: процессом линейно накапливаемой вола-тильности а2t = D{аWt}, t е [0, Т]. В случае же неоднородности, зависимости волатильности от времени а2 = аг(Г) (но ограниченности в смысле среднего-квадратичного), и зависимости среднего от времени С = С(Г) мы получаем сходимость логарифмов цены к процессу
0 < t < Т, (24)
где dWu - стандартный гауссовский «белый шум». Процессом накопленной волатильности в этом случае является интеграл jt0^2(u)du.
3. Процесс последовательных сравнений и стохастическая модель ценообразования
Запишем схему сравнения по шагам
SJ+1 = SjYj+1, у < N, (25)
- в рекуррентном виде;
^ = 50 , I < М, (26)
1=1
- в явном виде. Здесь {у}, ] < N -последовательность независимых одинаково распределенных случайных величин с распределением (задаваемым случайной величиной У)
у А , Р 1А д
и, р(5о) 1 (27)
и д(в0)
Для применения Леммы 1 перейдём к логарифмам
^ = 1п У},] < N, (28)
и, затем, проведём центрирование
- = 1п У} - Е 1п У},
1 ± N.
(29)
= (и- 1п + о (1) - с. =
n
= С. (n);
(30)
(дисперсия от сдвига не зависит, а случайные величины ] < N, получены из «штрихованных» величин £,) путем сдвига). Видим, что условия Леммы 1 выполнены; в формулировке Леммы 1 мы полагает t = 1 и утверждаем необходимый нам факт сходимости распределения логарифма цены к нормальному закону распределения при N ^ да
Ь {5о П Ъ } =
= 1п5о + £^ + ^ 1=1 1=1 ^^N(0,1) + р
(31)
Положим (в обозначениях Леммы 1) п = N. Положим первый и второй центральный момент равными следующим выражениям ^ ^ да)
= = ^ + о(1)
(Логарифмы 1п5Ь сократились!). Обозначение <»» мы полагаем общепринятым для сходимости по распределению, слабой сходимости распределений.
Замечание 1. Отметим, что в предельном выражении (31) распределение коэффициента сравнения вида У вовсе не обязано иметь бинарное распределение. Единственное требование - условия Леммы 1, связанные с независимостью и моментами. Поэтому мы можем говорить об инвариантности метода доказательства относительно выбора закона распределения коэффициентов сравнения.
Замечание 2. Вследствие факта сходимости в Лемме 1 накопленных сумм к процессу броуновского движения мы сможем получить сформированное логнормальное распределение и в случае, когда в процесс сравнения будет вовлечена подвыборка или увеличенная выборка объема [М], 0 < t < Т < да, [•] -целая часть. Важно лишь, чтоб коэффициенты сравнения были независимыми и однородно малыми в определенном смысле. В таких случаях у нас лишь будет меняться дисперсия нормального распределения от логарифма в соответствии со свойством броуновского движения DW(t) = ^ t > 0. В случае неоднородности у нас возникнет нетривиальный процесс накопленной во-латильности (Смотри Замечание к Лемме 1.)
Замечание 3. Известно, что сходимость по распределению еще
2
не влечет сходимость моментов (математического ожидания, дисперсии) и других характеристик распределений. Однако, характеристики будут сходиться при слабой сходимости распределений в случае, когда они равномерно ограничены. Данное условие в нашей модели автоматически выполняется, т.к. логарифмы всех цен ограничены одной константой, возможно, и весьма большой.
4. Инвариантность распределения относительно выбора начального объекта сравнения и анализ коэффициентов
Рассмотрим вопрос зависимости вероятностей у коэффициентов сравнения от начального значения 50. Применим известные элементарные формулы для дисперсии и математического ожидания бинарных случайных величин к случайным величинам вида £' с распределением, заданным (29), ] = 1, ..., Ы,
Е$ = рЫи + сНпО =
= (2р - 1)1пи, (32)
= D£J = (1пи - 1п£)2рд =
= (2 1пЦ)2рд = 4р(1 -р)(1пи)2 (33)
Приравниваем вычисленные математическое ожидание (32) и дисперсию (33) к выражениям в (30) и получаем, пренебрегая членами малого порядка, систему из двух уравнений (относительно неизвестных р и 50)
(2 р - 1)1п и = (р- 1п 5-0) N
у1р(1 - Р) 21п и = -£= V!
Поделим в (34) первое уравнение на второе уравнение и получим следующее Утверждение 1, устанавливающее соотношение между вероятностью р и коэффициентом относительного риска Я, определенного (3), (а, следовательно, - между р и 50 при известных значениях параметров ¡л и а2).
Утверждение 1. Вероятность подъема цены р в процессе последовательных сравнений зависит только от коэффициента относительного
(34)
риска и количеством объектов сравнения следующим образом:
1 Р-(1 - Р)_
Ч Р(1 - р)
1 р- 1п 5"0
а
л/ы
="Ттт■ (35)
Комментарий к Утверждению
1. Уравнение (35) имеет в правой части член, убывающий к нулю со скорость порядка ЧЫ. Это является следствием того факта, что в следствие механизма работы Центральной предельной теоремы для последовательности испытаний Бернулли отклонение количества успехов от ожидаемого количества успехов попадает в интервал, пропорциональный ЧК. Другими словами, при случайном выборе логарифма начального значения 1п50 из N данных возможных значений мы с вероятностью, стремящейся к единице при N ^ ю, выберем случайное число 1п5о, отличающееся от среднего Л на величину пропорциональную ЧК, когда логарифм коэффициента (одного) сравнения имеет порядок малости 1 / ЧК.
Введем обозначение 1(р) для функции в левой части (35), заданной
на интервале (0, 1) и затем рассмотрим поведение этой функции,
А 1
i (р=2
1р 1 - р
N 1 - р \ р )
(36)
Нетрудно убедиться, что функция 1(р) непрерывная (даже гладкая - бесконечно-дифференцируема), строго монотонно возрастающая, стремится к минус бесконечности при р ^ 0, равна нулю при р = 1/2 (то есть при симметричном случайном блуждании) и стремится к плюс бесконечности при р ^ 1. Из этих свойств следует, что: в уравнении (36) существует и единственное решение р е (0, 1) при всякой вещественной левой части. Более того, при известных значениях параметров ц и а в объеме выборки N мы по заданному значению 50 можем вычислить р = р(50), - и наоборот. Так как мы работаем с выборкой (цен на объекты недвижимости), имеем объем этой выборки Ы, то мы всегда можем вычислить выборочные среднее ¡л и дисперсию а2, подставить их в (35) вместо ц и а, соответственно, и получить прямое (вычисляемое) соотношение между 50 и р = р(50). Простой анализ построенной сто-
Рис. 1. График функции 1(р).
хастической модели и полученных формул позволяет сделать следующие важные практические выводы.
Вывод 1. Так как выполняется равенство по распределению (1), то среднее логнормального закона распределения равно ехр{л + а2 / 2}, а мода равна ехр{л - а2}, и легко видеть, что отношение среднего к моде не зависит от л, всегда больше единицы и равно ехр{(3 / 2)а2}. Разлагая экспоненту, получаем при малых а, что это отношение приблизительно равно 1 + (3 / 2)а2. Это означает, что применение средних арифметических для оценки рыночной стоимости (в качестве рекомендуемого наиболее вероятного значения - моды) приводит к систематической ошибке в сторону завышения, выраженную в относительной погрешности ехр{(3 / 2)а2} - 1(х 100%). Например, при а = 0,258 систематическая ошибка составляет 10,5%.
Вывод 2. Иногда в практике оценки вместо средних арифметических (оценка математического ожидания) применяют средние геометрические (оценка медианы). Так как медиана логнормального закона равна ехр{м}, то отношение медианы к моде не зависит от ¡, всегда больше единицы и равно ехр{а2}. Разлагая экспоненту, получаем при малых а, что это значение приблизительно равно 1 + а2. Это означает, что применение средних геометрических для оценки рыночной стоимости (наиболее вероятного значения - моды) приводит к систематической ошибке в сторону превышения, выраженную в относительной погрешности ехр{а2} - 1(х 100%). В частности, при а = 0,258 систематическая ошибка составляет 6,9%.
Вывод 3. Из свойств логнор-мального распределения следует, что при постоянном ц, с ростом а2 увеличивается разница между математическим ожиданием, медианой и модой. Таким образом, при постоянном ц и растущей а2 среднее арифметическое растет, а рыночная стоимость (мода - наиболее вероятное значение) уменьшается.
Итак, некоторый начальный спектр (интервал) цен достаточно
общего вида (такого, чтобы величины р и 1 - р, зависящие от Ы, не вырождались с ростом Ы) определяет параметры л и а формирующегося с помощью последовательных сравнений логнормального закона цены. Мы полагаем, что участники рынка действуют разумно и выбирают объекты сравнения, руководствуясь определенными принципами однородности - объекты имеют примерно одинаковый набор ценооб-разующих факторов. Тогда для разных групп объектов, обладающих одинаковой однородностью, будут формироваться логнормальные законы распределения с различными параметрами л и а. Такая модель автоматически приведет к смеси логнормальных распределений, что мы кратко и обсудим в следующем разделе.
5. Множество выборок объектов сравнения - формирование смеси распределений
Из обычной практики следует, что объекты сравнения далеко не всегда однородны. В этом случае элементы для сравнения мы берем из выборок, сформированных по тому или иному принципу однородности. Для каждой такой выборки формируется в соответствии с изложенным выше процессом сравнения свое логнормальное распределение, со своими параметрами ц и а2. При этом совокупное распределение получается путем смешивания (о смеси распределений см. [11]) полученных логнормальных распределений, когда смешивающие коэффициенты (веса) получаются пропорциональными количеству единиц сравнения в выборке, сформированной по данному принципу однородности. Отметим, что может быть принят альтернативный принцип формирования смешивающих коэффициентов - пропорционально совокупной денежной массе по каждой зоне ценовых предпочтений (по каждому району). Случай, когда веса пропорциональны денежной массе согласуется с тем, что логнор-мальное распределение получается для цены - денежной единицы.
Определение смеси. Дадим определение для конечной смеси, то
есть, когда смешивается конечное число распределений. Пусть имеется т функций распределения ^1, ..., Ет - «элементов смеси», и т неотрицательных чисел, называемых «весами», а1, ..., ат, сумма которых равна единице. Тогда функция Н, определенная посредством равенства
Н(х) = а^(х) + ... + ат^т(х),
х е Я (37)
является функцией распределения, которое называется смесью распределений Fl, ..., Ет с соответствующими весами а1, ..., ат.
Необходимо отметить, что если у распределений существуют плотности f, ..., /т, то и у распределения Н существует плотность h, равная
Н(х) = аф(х) + ... + ат/т(х),
х е Я (38)
Обращаем внимание читателя на то, что определение смеси исходно никак не связано со случайными величинами (можно было бы построить аналог (37) для любых функций), но в прикладном смысле распределения без случайных величин практически не рассматриваются. Введем случайные величины Х1, ..., Хт, имеющие соответствующие функции распределения ., ^т, и, так называемую, «смешивающую» случайную величину V, принимающую значения 1, ., т с соответствующими вероятностями а1, ..., ат (а числа а по построению естественно трактуются, как вероятности). Важно, чтоб случайная величина V не зависела от Х1, ..., Хт, при этом случайные величины (Х1, ..., Хт) между собой могут как угодно зависеть. Тогда смесь реализуется по следующему алгоритму: сначала реализуем V, а затем реализуем Ху с тем номером, на который указала V. Полученную в результате такого алгоритма случайную величину обозначим У. С помощью формулы полной вероятности легко видеть, что функция распределения для У есть Н. Отметим, что моменты распределений «смешиваются» при формировании смеси, то есть
ЕУ = а1ЕХ1к + ... + а„ЕХтк (39)
для любого натурального k. Также заметим, что ничего принципиаль-
ного не изменится, если рассматривать счетный набор весов. Говоря языком чистой математики, смесь распределений - это выпуклая аффинная комбинация распределений.
Для дальнейших применений нам потребуется следующий факт, который, образно говоря, заключается в том, что «смесь нормальных распределений после потенцирования превращается в смесь логнор-мальных распределений» с теми же весами.
Утверждение 2:
A. Пусть функция распределения Н случайной величины У задается в виде, соответствующем (37), причем функции распределения ^1, ..., Гт являются функциями распределения нормального закона с соответствующими параметрами Щ^-, а2-), ] = 1, .••, т. Тогда случайная величина ехр(У) имеет функцию распределения вида а\0\ + ... + атОт, где Gj - функция распределения логнормального закона с параметрами: - медиана, а2-- волатильность; ] = 1, ..., т.
B. Обратно, пусть распределение случайной величины 2 предста-вимо в виде смеси логнормальных распределений с соответствующими параметрами: - медиана, О- -волатильность; ] = 1, ..., т. Тогда распределение случайной величины У = 1п(2) является смесью с теми же весами нормальных распределений
а;), 1 = 1, т.
Доказательство Утверждения
2. Докажем сначала первую часть А Утверждения 2. Имеем для любого положительного у, используя определение функции распределения для случайных величин,
Р(ехр(У) < у) = Р(У < 1п(у)) =
= ^арХ] < 1п(у)) = .
]
= ^а]Р(вХ' < у)
]
Мы использовали тот факт, что при применении к обеим частям неравенства монотонно растущей функции - неравенство не меняется. Логарифм и экспонента очевидно таковые. Часть В доказывается аналогично при эксплуатации тех же
свойств логарифма и экспоненты. Утверждение 2 доказано.
Задача выделения из смеси распределений ее компонент относится к классу некорректных задач. Для решения таких задач часто используются факты и приемы, относящиеся к предметной области. В нашем случае - это разбиение общей выборки на естественные зоны, кластеры, районы и т.п.
Доказанное Утверждение 2. позволяет нам переходить от единичного нормального распределения к смеси нормальных распределений, и, в конечном итоге, к смеси логнормальных распределений, используя универсальную Лемму 1 и метод последовательных сравнений для каждой зоны, кластера и т.п.
Примеры в Приложении 2 иллюстрируют тот факт, что в зонах ценовых предпочтений (приведены только некоторые из 12 зон по СПб), в кластерах «тип дома» и даже в районах наблюдается хорошее соответствие логнормаль-ному закону, как элементу общей смеси по СПб логнормальных распределений.
Важно заметить, что даже несмотря на тот факт, что смесь логнор-мальных распределений не является логнормальным распределением, в случае, когда множество совокупных значений параметров л и а2 (взятых по всем зонам СПб) имеет небольшой разброс, компактно, полученное в результате смешивания [теоретическое] распределение может остаться унимодальным (то есть иметь одну моду - одно значение с максимальной вероятностью). Этот факт имеет свое подтверждение на реальных данных, см. Приложение 1, Рис. 2.
В любом случае, когда значения Л и о2 в смешивающих распределениях имеют небольшой разброс, и смешивающие распределения имеют веса примерно одинакового порядка, то итоговое распределение имеет, по меньшей мере, хорошо выделяемый, узкий интервал «максимальной вероятности». Данный тезис находит свою иллюстрацию даже в статистике, взятой по всему городу СПб (см. Приложение 1, Рис. 7).
6. Критерий согласия, статистическая оценка и калибровка параметров логнормального распределения цены
Все приведенные в настоящей статье диаграммы и аппроксимации получены с помощью статистического пакета R. Модельным параметрам Л и о2 присваиваются значения оценочных. Оцениваются (стандартным методом вычисления выборочного среднего и среднеквадратичного отклонения (ско)) параметры среднего и масштаба для логарифмов цены. Модельные параметры л и о2 мы полагаем равными тем значениям, что получены в результате оценивания.
Оценки параметров л и о2 (и и о2, соответственно) могут проводиться отдельно, но могут и быть взяты из работы тестов, представленных ниже. Статистические выводы начинаются с оценки согласия [fit] и качества согласия с логнормальным распределением цены по отдельно взятой выборке. Проверяется согласие с логнормальным распределением универсальным методом, предлагаемым пакетом R, - функцияfitdistr пакета R [12]. Алгоритм работы функции fitdistr основан на применении классического метода максимального правдоподобия для оценки параметра предлагаемого распределения (предлагается логнормальное распределение) с использованием меры «качества», основанной на среднеквадратичной метрике. Очевидно, что среднеквадратичная метрика чересчур чувствительна к «выбросам» и мало чувствительна к тяжести хвостов распределения, поэтому представляется оправданной возможность дальнейшей «корректировки» параметров, в том числе и визуальным способом. Далее проводится визуальный анализ полученной теоретической плотности [соответствующего выборке логнормального закона] с эмпирической плотностью, представленной в виде гистограммы. На этом шаге параметры л и о2 калибруются путем «подгонки», с тем, чтобы добиться максимального визуального соответствия с эмпирической плотностью. Относительное изменение про-
исходит с порядком не более 10%. Пусть полученные после подгонки параметры обозначаются ц* и о*2, соответственно. Для установления факта значимого (или не значимого) статистического соответствия анализируемой выборки закону логнор-мального распределения используется мощный критерий согласия Колмогорова-Смирнова (КС-тест) (см., [13], [14], [15], [16], [17]). При удовлетворении наперед заданному значению р^аШена основе КС-теста принимается гипотеза о логнор-мальности с параметрами ц* и а*2 для распределения цены случайно взятого объекта недвижимости. После того, как распределение прошло КС-тест наиболее вероятной ценой берется мода (как и предписывается, собственно, самим определением понятия «моды») логнормального распределения с параметрами ц* и о*2, то есть
ехр(ы* и a*2).
(40)
Для согласия выборки с законом распределения смеси логнормаль-ных распределений применялись функции пакета mixtools в рамках пакета Я. С помощью этих функций формально выражалась смесь, построенная на основе разделения (разбиения) на районы, зоны ценовых предпочтений и т.п. Затем формировалась гипотеза о логнор-мальности элементов построенной смеси и пакетом mixtools оценивались параметры логнормальных распределений - элементов смеси. После этого исследуемая выборка проходила КС-тест на согласие с построенной смесью. Обработанные этой функцией данные показали (см. Приложение 1, Рис. 2), что выборка «масс-маркет» по СПб дает хорошо обусловленную моду для смеси, когда элементы смеси - районы города. Более того, все районы СПб, включая премиум сегмент, дают эмпирическое распределение, вполне удовлетворительно согласующееся посредством применения КС-теста со смесью логнормальных распределений. Представляется перспективным использовать в дальнейшем пакет fitdistrplus в рамках пакета Я. В пакете fitdistrplus есть опции, позволяющие улучшать оценки
параметров, в частности, снижать чувствительность к «выбросам», лучше обрабатывать тяжелые хвосты и прочее.
Приложение 1
Рассмотрим ряд примеров. Обработан статистический материал - объявления о продажах объектов вторичной жилой недвижимости, опубликованных в [18]. Из выборки исключены некорректные объявления: нет цены или площади, другие ошибки. Рассматривается следующий масштаб распределения цен: тыс. руб. за 1 кв. м.
Разделение смеси по административным районам города. Масс-маркет
Пример статистического анализа выборки: «масс-маркет» Санкт-Петербурга. На приведенном ниже рис. 2 показаны гистограмма, кривая плотности логарифмически нормального распределения, - при оценивании «одиночным» логнормаль-ным распределением, и кривая смеси при оценивании смесью логнор-мальных распределений. Каждый элемент смеси представляет из себя «масс-маркет» одного района города. Представлены все 18 районов СПб. Объем всей выборки: 12772 объявления о продаже. Веса смеси пропорциональны количествам объявлений. Легенда для диаграммы: пунктир - логнормальное распределение, сплошная линия - смесь логнормальных распределений. Вертикальной линией обозначена мода смеси; вертикальным пунктиром обозначена мода одного логнормального распределения. Полученные результаты: наиболее вероятные значения оказались равными
88,0669 тыс. руб. за кв. м. и 88,2176 тыс. руб. за кв. м., соответственно для единичного логнормального распределения и для смеси. Как видно, разница весьма незначительная, и при этом, распределение, полученное в результате смеси, сохраняет свойство унимодальности. Полученное в результате КС-теста, примененного к данной смеси, зна-чениеp - value = 1,258e-6.
Необходимо отметить, что несмотря на сравнительно небольшую величину a значение p-value мало (специфичной особенностью КС-теста при больших выборках является чувствительность к небольшим компонентам смести, что чрезмерно снижает значение p-value), что свидетельствует о крайне малой значимости модели одного логнормального распределения для решения задачи аппроксимации распределения данной выборки: p - value < 1,517e-12, - на несколько порядков меньше, чем при использовании значимости для всего города СПб одного логнормального распределения (Для модели смеси), - и это является важным статистическим аргументом для разбиения города на подвыборки, формирующие отдельно взятые логнормальные распределения.
Следующий пример показывает, что выделение некоторых отдельно взятых административных районов города СПб дает хорошо обусловленное с точки зрения p-value КС-теста логнормальное распределение.
Разделение смеси
по административным районам
города
Разделение смеси по административным районам города показывает,
Рис. 2. Масс-маркет. Смесь
Рис. 3. Адмиралтейский район СПб
Рис. 4. Зона ценовых предпочтений «Золотой треугольник»
Рис. 5. Зона ценовых предпочтений № 4
что большинство районов города (17 из 18) являются районами застройки, достаточно однородной, чтобы распределение цены 1 кв. м. в них могло быть представлено логарифмически нормальным распределением. Центральный район недостаточно хорошо аппроксимируется логнор-мальным распределением, но легко может быть разделен на зоны, аппроксимируемые логнормальным распределением. На приведенном ниже рис. 3 показаны гистограмма и кривая плотности логарифмически нормального распределения для Адмиралтейского района г Санкт-Петербурга (за исключением премиум-сегмента - объектов на Конногвардейском бульваре).
Выборка цен предложений вторичной жилой недвижимости в Адмиралтейском районе г. Санкт-Петербурга на 31.03.2014 г. и ее аппроксимация логнормальным распределением. Объем выборки: N = 403. Параметры логнормального закона: ц = 4,57, а = 0,226. Среднее арифметическое по выборке равно 99,799 тыс. руб. за кв. м. Оценка рыночной стоимости 1 кв. м., мода по логнормальной аппроксимации равна 91,983 тыс. руб. за кв. м. Результат КС-теста: р^а1ые = 0,3245
Разделение исходной выборки по зонам ценовых предпочтений
Поскольку районы могут быть недостаточно однородными, в КУГИ г. Санкт-Петербурга в 2012 году был разработан принцип разбиения города по зонам ценовых предпочтений. Всего было выделено 12 зон (1 зона - наиболее дорогие адреса, увеличение номера зоны означает сдвиг к периферии) см. [19]. Разбиение исходной выборки по зонам ценовых предпочтений дает 12 подвыборок хорошо приближающихся логарифмически нормальным законом распределения. Например, на рис. 4 и рис. 5, представлены гистограммы и кривые плотности логнормального закона распределения для зоны 1 (Золотой треугольник) и зоны 4.
Выборка цен предложений вторичной жилой недвижимости в г. Санкт-Петербурге для зоны ценовых предпочтений №1 на 31.03.2014 г. и ее аппроксимация логнормальным
распределением. Объем выборки: n = 47. Параметры логнормального закона: ц = 5,3, а = 0,46. Среднее арифметическое по выборке равно 213,358 тыс. руб. за кв. м. Оценка рыночной стоимости 1 кв. м., то есть мода по логнормальной аппроксимации равна 162,130 тыс. руб. за кв. м. Результат КС-теста: р^аЫе = 0,713.
Выборка цен предложений вторичной жилой недвижимости в г. Санкт-Петербурге для зоны ценовых предпочтений №4 на 31.03.2014 г. и ее аппроксимация логнормальным распределением. Объем выборки: n = 734. Параметры логнормального закона: ^ = 4,631, а = 0,191.Среднее
арифметическое по выборке равно 106,339 тыс. руб. за кв. м. Оценка рыночной стоимости 1 кв. м., то есть мода по логнормальной аппроксимации равна 98,940 тыс. руб. за кв. м. Результат теста Колмогорова-Смирнова: р^аЫе = 0,05332.
Разделение исходной выборки по типам домов
Разделение выборки по типам домов дает 12 подвыборок (с учетом объединения зданий разных номерных серий в группы панельной и блочной застройки). Все подвыборки хорошо аппроксимируются логнормальным распределением. На Рис. 6 представ-
Рис. 6. Застройка сталинского периода
Рис. 7. Весь город СПб
лены гистограмма и кривая плотности логнормального распределения для «сталинских домов».
Выборка цен предложений вторичной жилой недвижимости в г. Санкт-Петербурге для застройки сталинского периода на 31.03.2014 г. и ее аппроксимация логнормальным распределением. Объем выборки: N = 546. Параметры логнормального закона р = 4,612, а = 0,185.Среднее арифметическое по выборке равно 104,035 тыс. руб. за кв. м. Оценка рыночной стоимости 1 кв. м., то есть мода по логнормальной аппроксимации, равна 97,298 тыс. руб. за кв. м. Результат теста Колмогорова-Смирнова: р^аЫе = 0,1575 .
Весь город. Статистика, идентифицированная по зонам ценовых предпочтений
Рассмотрим выборку по всему городу СПб - статистику объявлений по всему городу, которая попала в идентификаторы зон ценовых предпочтений (см. [19]). Такая выборка, безусловно, является смесью не вполне однородных выборок, различающихся по целому спектру ценообразующих факторов, таких
как место (административные, территориальные и т.п. деления), типам домов, этажности, типам квартир и т.д. На рис. 7 видно, что эмпирическое распределение имеет «тяжелый» правый «хвост», указывающий на необходимость дальнейшего разделения представленной смеси, видимо выделение такого/таких элементов смеси, как премиум сегмент/сегменты. Легенда для диаграммы: пунктир -ядерная оценка плотности, полученная с помощью пакета Я. Объем выборки: 10627. Параметры логнор-мального закона: р = 4,539, а = 0,20. Среднее арифметическое по выборке равно 97,462 тыс. руб. за кв. м. Оценка рыночной стоимости 1 кв. м., то есть мода по логнормальной аппроксимации равна 89,579 тыс. руб. за кв. м. Результат теста Колмогорова-Смирнова: р^аЫе < 2,861е-8.
Малое значение р-т1ие, как и в случае Рис. 2 свидетельствует о малой значимости для всего города одного логнормального распределения, -и, как отмечалось выше, это является важным статистическим аргументом для разбиения города на подвыбор-ки, формирующие отдельно взятые логнормальные распределения.
Литература
1. Грибовский С.В. «Математические методы оценки стоимости недвижимого имущества.» Учебное пособие. М.: Финансы и Статистика, 2008
2. Федеральный стандарт оценки №2 «Цель оценки и виды стоимости.» Приказ Минэкономразвития РФ № 255 от 20.07.2007. Зарегистрирован в Минюсте РФ 23.08.2007
3. International Valuation Standard Council, London, 2013. URL:http://www.valuersinstitute. com.au/docs/professional_practice/ International%20Valuation%20 Stand
4. European Valuation Standards, 7-th edition, Brussels, 2012. URL: http:// www.tegova.org/en/p4fe1fcee0b1db
5. Uniform Standards of Professional Appraisal Practice, 2014-2015 edition, Annapolis, Maryland, USA, 2014. URL: http://www.uspap.org/
6. Royal Institution of Chartered Surveyors Valuation Professional Standard, London, 2014. URL:http:// www.rics.org/uk/shop/RICS-Valuation-Professional-Standards-Red-Book-2013-19750.aspx
7. Sharpe, William R., "The Sharpe Ratio." Journal of Portfolio Management21(1): 49-58, 1994.
8. Биллингсли, П. «Сходимость вероятностных мер.» М.: Наука, 1977
9. Русаков О.В., Солев ВН. «Предельный переход для классических моделей в стохастической финансовой математике (учебно-методическое пособие)» Издательство Санкт-Петербургского государственного университета, 1999.
10. Ширяев А.Н. «Вероятность» М.: Наука, 1989.
11. Феллер, В. «Введение в теорию вероятностей и ее приложения.» т. II, М.: Мир, 1984.
12. Venables, W.N. and Ripley, B.D. "Modern Applied Statistics with S." 4th edition. Springer, 2002.
13. Birnbaum, Z.W. and Fred, H.T. "One-sided confidence contours for probability distribution functions." The Annals of Mathematical Statistics, 22/4: 592-596, 1951.
14. William, J. Conover "Practical Nonparametric Statistics." New York: John Wiley & Sons. Pages 295-301 (one-sample Kolmogorov
test), 309-314 (two-sample Smirnov test), 1971.
15. Durbin, J. "Distribution theory for tests based on the sample distribution function." SIAM, 1973.
16. Большев, Л.Н., Смирнов, Н.В., «Таблицы математической статистики.» М.: Наука, 1983.
17. Marsaglia, G., Tsang, W.W., and Wang, J. "Evaluating Kolmogorov's Distribution." Journal of Statistical Software, 8/18, 2003. URL: http:// www.jstatsoft.org/v08/i18/.
18. Бюллетень Недвижимости №13/1605 Санкт-Петербург, 24.03.2014.
19. Отчет об определении кадастровой стоимости объектов недвижимости (за исключением земельных участков), расположенных на территории Санкт-Петербурга. КУГИ Правительства Санкт-Петербурга, том 2, раздел 2.3, СПб, 2012 http:// rosreestr.ru/wps/portal/p/cc_ib_portal_ services/cc_ib_ais_fdg/co/ Субъект РФ: С-Петербург. Виды объектов: Недвижимость, помещения. Отчет №32-1-0733/2012(2), Санкт-Петербург, 30.11.2012
References
1. Gribovsky, S. V. "Mathematical methods of survey (estimation) of real estate value" Textbook. M.: Finance and Statistics, (in Russian), 2008.
2. Federal Standard of Survey (Valuation) № 2 "The purpose of valuation and types of prices". Decree of the Ministry of Economic Development of
the Russian Federation №255 dated at 20.07.2007. Registered in the Ministry of Justice of RF 23.08.2007
3. International Valuation Standard Council, London, 2013. URL: http:// www.valuersinstitute.com.au/docs/pro-fessional_practice/International%20 Valuation%20Standards%202013.pdf
4. European Valuation Standards, 7thedition, Brussels, 2012. URL: http:// www.tegova.org/en/p4fe1fcee0b1db
5. Uniform Standards of Professional Appraisal Practice, 2014-2015 edition, Annapolis, Maryland, USA, 2014. URL: http://www.uspap.org/
6. Royal Institution of Chartered Surveyors Valuation Professional Standard, London, 2014. URL: http:// www. ric s. org/uk/shop/RICS-Valu-ation-Professional-Standards-Red-Book-2013-19750.aspx
7. Sharpe, William R., "The Sharpe Ratio."Journal of Portfolio Management 21(1):49-58, 1994.
8. Billingsley, P. "Convergence of Probability Measures" John Willey & Sons, New York, 1968.
9. Rusakov, O. V., Solev, V. N. "Passage to the Limit for Classical Models in Stochastic Financial Mathematics (study guide)" Saint-Petersburg State University Press, (in Russian) 1999.
10. Shiryaev, A. N. "Probability" M.: Nauka (in Russian), 1989.
11. Feller, W. "An Introduction to Probability Theory and Its Applications" Volume II, 3rd ed., John Willey & Sons, New York, 1971.
12. Venables, W.N. and Ripley, B.D. "Modern Applied Statistics with S."4th ed., Springer, 2002.
13. Birnbaum, Z.W. and Fred, H.T. "One-sided confidence contours for probability distribution functions." The Annals of Mathematical Statistics, 22/4: 592-596, 1951.
14. William, J. Conover "Practical Nonparametric Statistics." New York: John Wiley & Sons. Pages 295-301 (one-sample Kolmogorov test), 309314 (two-sample Smirnov test),1971.
15. Durbin, J. "Distribution theory for tests based on the sample distribution function." SIAM, 1973.
16. Bolshev, L. N., Smirnov N. V., "Tables of Mathematical Statistics" M.: Nauka (in Russian), 1983 Real Estate Bulletin № 13/1605 Saint-Petersburg, 24.03.2014.
17. Marsaglia, G., Tsang, W.W., and Wang, J. "Evaluating Kolmogorov's Distribution." Journal of Statistical Software, 8/18, 2003. URL: http:// www.jstatsoft.org/v08/i18/.
18. Real Estate Bulletin № 13/1605 Saint-Petersburg, 24.03.2014.
19. Report on determination of the real estate cadastre value (excluding land costs), situated in Saint-Petersburg. CASP Government of Saint Petersburg, edition 2, section 2.3, Saint-Petersburg, 2012. URL: http://rosreestr.ru/wps/portal/p/ cc_ib_portal_services/cc_ib_ais_fdg/co/ Federal Subject of RF: Saint-Petersburg. Types of property: Real Estate, Buildings. Report №32-1-0733/2012(2), Saint-Petersburg, 30.11.2012.
