Лекции статистика. Учебнометодическое обеспечение курса
![]()
|
2.3. Статистическое изучение вариацииПризнаки, изучаемые статистикой, варьируются (отличаются друг от друга) у различных единиц совокупности в один и тот же период или момент времени. Например, величина внешнеторгового оборота (ВО) варьируется по подразделениям Федеральной таможенной службы; величина экспорта (импорта) варьируется по направлениям экспорта (по разным странам-партнерам по внешней торговле), по видам товаров и т.п. Причиной вариации являются разные условия существования разных единиц совокупности. Например, огромное число причин влияет на масштабы внешней торговли различных стран мира. Для управления и изучения вариации статистикой разработаны специальные методы исследования вариации, система показателей, с помощью которой вариация измеряется, характеризуются ее свойства. Первым этапом статистического изучения вариации является построение рядараспределения(иливариационногоряда) – упорядоченного распределения единиц совокупности по возрастающим (чаще) или по убывающим (реже) значениям признака и подсчет числа единиц с тем или иным значением признака. Существует 3 вида ряда распределения:
Построим ряд распределения внешнеторгового оборота (ВО) по таможенным постам России, для чего необходимо провести статистическое наблюдение, то есть собрать первичный статистический материал, который представляет собой величину ВО по таможенным постам. Результаты наблюдения ВО по 35 таможенным постам региона за отчетный период представим в виде ранжированного по возрастанию величины ВО ряда распределения (таблица 0). Таблица 0. Внешнеторговый оборот (ВО) по 35 таможенным постам, млн.долл. № постаВО№постаВО№постаВО124,161354,122565,31227,061454,912669,24329,121555,742771,39431,171655,912877,12537,081756,072979,12639,111856,803084,34741,581956,933186,89844,842057,073291,74946,802158,393396,011048,372259,6134106,841151,442359,9535111,161252,562462,05Итого2100,00 Построим интервальный ряд распределения ВО по таможенным постам, для чего необходимо выбрать оптимальное число групп (интервалов признака) и установить длину (размах) интервала. Поскольку при анализе ряда распределения сравнивают частоты в разных интервалах, необходимо, чтобы длина интервалов была постоянной11. Оптимальное число групп выбирается так, чтобы достаточной мере отразилось разнообразие значений признака в совокупности и в то же время закономерность распределении, его форма не искажалась случайными колебаниями частот. Если групп будет слишком мало, не проявится закономерность вариации; если групп будет чрезмерно много, случайные скачки частот исказят форму распределения. Чаще всего число групп в ряду распределения определяют по формуле Стерждесса (5) или , (6) или (5) или , (6): ![]() ![]() где k – число групп (округляемое до ближайшего целого числа); N – численность совокупности. Из формулы Стерджесса видно, что число групп – функция объема данных (N). Зная число групп, рассчитывают длину (размах) интервала12 по формуле (7): ![]() где Xмax и Xmin — максимальное и минимальное значения в совокупности. В нашем примере про ВО по формуле Стерждесса (5) или , (6) определим число групп: k = 1 + 3,322lg35 = 1+ 3,322*1,544 = 6,129 ≈ 6. Рассчитаем длину (размах) интервала по формуле (7): h = (111,16 – 24,16)/6 = 87/6 = 14,5 (млн.долл.). Теперь построим интервальный ряд с 6 группами с интервалом 14,5 млн.долл. (см. первые 3 столбца табл. 0). Таблица 0. Интервальный ряд распределения ВО по таможенным постам, млн.долл. iГруппы постов по величине ВО XiЧисло постов fiСередина интервала Хi’Хi’fiНакопл. частота fi’|Хi’ - ![]() ![]() ![]() ![]() ![]() ![]() Рис. 0. Гистограмма распределения Рис. 0. Полигон распределения Данные табл. 0 и рис. 0 показывают характерную для многих признаков форму распределения: чаще встречаются значения средних интервалов признака, реже – крайние (малые и большие) значения признака. Форма этого распределения близка к нормальному закону распределения, которое образуется, если на варьирующую переменную влияет большое число факторов, ни один из которых не имеет преобладающего значения. Если имеется дискретный ряд распределения или используются середины интервалов (как в нашем примере про ВО – в таблице 0 в 4-м столбце рассчитаны середины интервалов как полусумма значений начала и конца интервала), то графическое изображение такого ряда называется полигоном (см. рис. 0)14, которое получается соединением прямыми точек с координатами Xi и fi. Вторым этапом статистического изучения вариации является расчет характеристик ряда распределения, которые описывают количественно его структуру, строение. Такова, например, медиана – величина варьирующего признака, делящая совокупность на две равные части – со значением признака меньше медианы и со значением признака больше медианы15. В нашем примере про ВО (табл. 0) медиана – это 18-й таможенный пост из 35 с величиной ВО 56,8 млн.долл. Из этого примера видно принципиальное различие между медианой и средней величиной: медиана не зависит от значений на краях ранжированного ряда. Даже если бы ВО 35-го таможенного поста был в 10 раз больше, величина медианы не изменилась бы. Поэтому медиану часто используют как более надежный показатель типичного значения признака, нежели средняя арифметическая, если ряд значений неоднороден, включает резкие отклонения от средней. В интервальном ряду распределения для нахождения медианы применяется формула: ![]() где Ме – медиана; X0 – нижняя граница интервала, в котором находится медиана; h – величина (размах) интервала; ![]() fMe – частота в медианном интервале. В табл. 0 медианным является среднее из 35 значений, т.е. 18-е от начала значение ВО. Как видно из столбца накопленных частот (6-й столбец), оно находится в третьем интервале. Тогда по формуле ( ): ![]() Аналогично медиане вычисляются значения признака, делящие совокупность на 4 равные по численности части – квартили, которые обозначаются заглавной латинской буквой Q с подписным значком номера квартиля. Ясно, что Q2 совпадает с Ме. Для первого и третьего квартилей приводим формулы и расчет по данным табл. 0: ![]() ![]() Так как Q2 = Ме = 59,30 млн.долл., видно, что различие между первым квартилем и медианой (–15,87) больше, чем между медианой и третьим квартилем (12,89). Этот факт свидетельствует о наличии некоторой несимметричности в средней области распределения, что заметно и на рис. 0. Значения признака, делящие ряд на 5 равных частей, называются квинтилями, на 10 частей – децилями, на 100 частей – перцентилями. Эти характеристики применяются при необходимости подробного изучения структуры ряда распределения16. Безусловно, важное значение имеет такая величина признака, которая встречается в изучаемом ряду распределения чаще всего. Такую величину принято называть модой. В дискретном ряду мода определяется без вычисления как значение признака с наибольшей частотой. Обычно встречаются ряды с одним модальным значением признака. Если в ряду распределения встречаются 2 или несколько равных (и даже несколько различных, но больших чем соседние) значений признака, то он считается соответственно бимодальным или мультимодальным. Это свидетельствует о неоднородности совокупности, возможно, представляющей собой агрегат нескольких совокупностей с разными модами. В интервальном ряду распределения интервал с наибольшей частотой является модальным. Внутри этого интервала находят условное значение признака, вблизи которого плотностьраспределения (число единиц совокупности, приходящихся на единицу измерения варьирующего признака) достигает максимума. Это условное значение и считается точечноймодой. Логично предположить, что такая точечная мода располагается ближе к той из границ интервала, за которой частота в соседнем интервале больше частоты в интервале за другой границей модального интервала. Отсюда получаем обычно применяемую формулу ( ): ![]() где Мо – мода; Х0– нижнее значение модального интервала; fMo – частота в модальном интервале; fMo-1 – частота в предыдущем интервале; fMo+1 – частота в следующем интервале за модальным; h – величина интервала. По данным табл. 0 рассчитаем точечную моду по формуле ( ): ![]() К изучению структуры ряда распределения средняя арифметическая величина также имеет отношение, хотя основное значение этого обобщающего показателя другое. В интервальном ряду распределения ВО по таможенным постам средняя арифметическая рассчитывается как взвешенная по частоте середина интервалов X (расчет числителя – в 5-м столбце табл. 0) по формуле ( ): ![]() ![]() Различие между средней арифметической величиной (60,82), медианой (59,30) и модой (58,96) в нашем примере невелико. Чем ближе распределение по форме к нормальному закону, тем ближе значения медианы, моды и средней величины между собой. Третьим этапом статистического изучения вариации является расчет показателей размера и интенсивности вариации. Простейшим показателем является размах вариации – абсолютная разность между максимальным и минимальным значениями признака из имеющихся в изучаемой совокупности значений (8): ![]() Поскольку величина размаха характеризует лишь максимальное различие значений признака, она не может измерять закономерную силу его вариации во всей совокупности. Предназначенный для данной цели показатель должен учитывать и обобщать все различия значений признака в совокупности без исключения. Число таких различий равно числу сочетаний по два из всех единиц совокупности (в нашем примере про ВО число сочетаний составит ![]() ![]() В нашем примере про ВО по данным табл. 0 среднее линейное отклонение вычисляется как взвешенное по частоте отклонение по модулю середин интервалов от средней арифметической величины (расчет числителя произведен в 7-м столбце табл. 0), т.е. по формуле (9): ![]() Это означает, что в среднем величина ВО в изучаемой совокупности таможенных постов отклонялась от средней величины ВО в РФ на 14,678 млн.долл. Простота расчета и интерпретации составляют положительные стороны показателя Л, однако математические свойства модулей «плохие»: их нельзя поставить в соответствие с каким-либо вероятностным законом, в том числе и с нормальным распределением, параметром которого является не средний модуль отклонений, а среднееквадратическоеотклонение, обозначаемое малой греческой буквой сигма ( ![]() ![]() ![]() В нашем примере про ВО по данным табл. 0 среднее квадратическое отклонение величины ВО по формуле (12) составило (расчет числителя произведен в 8-м столбце табл. 0): ![]() Среднее квадратическое отклонение по величине в реальных совокупностях всегда больше среднего модуля отклонений. Разница между ними тем больше, чем больше в изучаемой совокупности резких, выделяющихся отклонений, что служит индикатором «засоренности» совокупности неоднородными с основной массой элементами. Для нормального закона распределения отношение ![]() ![]() Квадрат среднего квадратического отклонения представляет собой дисперсию отклонений, на использовании которой основаны практически все методы математической статистики, ее формула имеет вид (13) – для несгруппированных данных (простая дисперсия) и (14) – для сгруппированных (взвешенная дисперсия): ![]() ![]() Еще одним показателем силы вариации, характеризующим ее не по всей совокупности, а лишь в ее центральной части, служит среднееквартильноерасстояние (отклонение), т.е. средняя величина разности между квартилями, определяемая по формуле (15): ![]() В нашем примере про ВО по формуле (15): ![]() Сила вариации в центральной части совокупности, как правило, меньше, чем в целом по всей совокупности. Соотношение между средним линейным отклонением и средним квартильным расстоянием служит для изучения структуры вариации: большое значение такого соотношения свидетельствует о наличии слабоварьирующего «ядра» и сильно рассеянного вокруг него окружения в изучаемой совокупности. Для нашего примера про ВО соотношение Л/q = 1,021, что говорит о совсем незначительном различии силы вариации в центральной части совокупности и на ее периферии. Для оценки интенсивности вариации и для сравнения ее в разных совокупностях и тем более для разных признаков необходимы относительныепоказателивариации, которые вычисляются как отношение абсолютных показателей силы вариации, рассмотренных ранее, к средней арифметической величине признака, то есть показатели (16) – (19):
В нашем примере про ВО эти показатели составляют: ![]() ![]() ![]() Оценка степени интенсивности вариации возможна только для каждого отдельного признака и совокупности определенного состава, она состоит в сравнении наблюдаемой вариации с некоторой обычной ее интенсивностью, принимаемой за норматив17. Так, для совокупности таможенных постов вариация величины ВО может быть определена как слабая, если ![]() ![]() ![]() Различная сила, интенсивность вариации обусловлены объективными причинами, поэтому нельзя говорить о каком-либо универсальном критерии вариации (например, 33%), так как для разных явлений и признаков этот критерий различен. Например, цена продажи американского доллара в коммерческих банках Н.Новгорода 26 июля 2007 года варьировала от 25,45 до 26,00 при средней цене 25,595 руб., тогда по формуле (16) ![]() Четвертым этапом статистического изучения вариации является расчет моментов распределения и показателей его формы. Для дальнейшего изучения характера вариации используются средние значения разных степеней отклонений отдельных величин признака от его средней арифметической величины. Эти показатели называются центральныемоментыраспределения порядка, соответствующего степени, в которую возводятся отклонения (табл. 0) или просто моментов (нецентральные моменты в таможенной статистике практически не используются). Таблица 0. Центральные моменты ПорядокмоментаФормулапо несгруппированным даннымпо сгруппированным даннымПервый μ1 ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() В нашем примере про ВО показатель асимметрии по формуле (20) составил (расчет числителя произведен в 9-м столбце табл. 0): ![]() Английский статистик К.Пирсон на основе разности между средней арифметической величиной и модой предложил другой показатель асимметрии (21): ![]() В нашем примере по данным табл. 0 показатель асимметрии по формуле (21) составил: ![]() Показатель асимметрии Пирсона (21) зависит от степени асимметричности в средней части ряда распределения, а показатель асимметрии (20) – от крайних значений признака. Таким образом, в нашем примере про ВО в средней части распределения наблюдается меньшая асимметрия, чем по краям, что видно и по графику (рис. 0). Распределения с сильной правосторонней и левосторонней асимметрией показаны на рис. 0. Левосторонняя As < 0 Правосторонняя As > 0 ![]() Мо ![]() ![]() ![]() Рис. 0. Асимметрия распределения С помощью момента четвертого порядка характеризуется еще более сложное свойство рядов распределения – эксцесс(от англ. «излишество»). Показатель эксцесса рассчитывается по формуле (22): ![]() Чаще всего эксцесс интерпретируется как «крутизна» распределения, что не совсем верно. График распределения может выглядеть сколь угодно крутым в зависимости от силы вариации признака: чем слабее вариация, тем круче кривая распределения при данном масштабе. Не говоря уже о том, что, изменяя масштабы по осям абсцисс и ординат, любое распределение можно искусственно сделать «крутым» и «пологим». Чтобы показать, в чем состоит эксцесс распределения, и правильно его интерпретировать, нужно сравнить ряды с одинаковой силой вариации (одной и той же величиной σ) и разными показателями эксцесса. Чтобы не смешать эксцесс с асимметрией, все сравниваемые ряды должны быть симметричными. Такое сравнение изображено на рис. 0. Ex > 0 Нормальное распределение Ex = 0 Ex < 0 ![]() Рис. 0. Эксцесс распределения Наличие положительного эксцесса означает наличие слабоварьирующего «ядра» и сильно рассеянного вокруг него окружения в изучаемой совокупности. Отрицательный эксцесс означает отсутствие такого «ядра». В нашем примере по формуле (22) эксцесс составил (расчет числителя произведен в 10-м столбце табл. 0): ![]() Пятым этапом статистического изучения вариации является проверка соответствия ряда распределения теоретическому (нормальному, логнормальному, биномиальному, распределению Руассона и др.) с помощью критериев согласия, среди которых чаще всего применяют критерии Пирсона χ2, Колмогорова и Романовского. Данный этап не входит в программу изучения для студентов заочного отделения (в случае необходимости – см. конспект лекций для дневного отделения). |