Заглавная страница Избранные статьи Случайная статья Познавательные статьи Новые добавления Обратная связь FAQ Написать работу КАТЕГОРИИ: АрхеологияБиология Генетика География Информатика История Логика Маркетинг Математика Менеджмент Механика Педагогика Религия Социология Технологии Физика Философия Финансы Химия Экология ТОП 10 на сайте Приготовление дезинфицирующих растворов различной концентрацииТехника нижней прямой подачи мяча. Франко-прусская война (причины и последствия) Организация работы процедурного кабинета Смысловое и механическое запоминание, их место и роль в усвоении знаний Коммуникативные барьеры и пути их преодоления Обработка изделий медицинского назначения многократного применения Образцы текста публицистического стиля Четыре типа изменения баланса Задачи с ответами для Всероссийской олимпиады по праву Мы поможем в написании ваших работ! ЗНАЕТЕ ЛИ ВЫ?
Влияние общества на человека
Приготовление дезинфицирующих растворов различной концентрации Практические работы по географии для 6 класса Организация работы процедурного кабинета Изменения в неживой природе осенью Уборка процедурного кабинета Сольфеджио. Все правила по сольфеджио Балочные системы. Определение реакций опор и моментов защемления |
Расчет графика плотности вероятности нормального законаСодержание книги
Поиск на нашем сайте
Сумма плотностей вероятности, равная 2,5037, близка к среднеквадратичному отклонению s у = 2,53. Если продолжить кривую плотности вероятности за пределы графика (рис.2.13), то сумма будет точно равна s у. Чтобы перейти от плотности вероятности к теоретической частоте n т, применяется формула (2.36) Здесь N = 147 (суммарное число измерений); h – размер (шаг) класса. Поскольку в качестве аргумента взят номер класса у, шаг h = 1, s у = 2,53. Следовательно,
Такое же соотношение получится, если в качестве аргумента взять случайную величину х. Из табл.2.5 имеем h = 2, s = 5,06, и отношение h /s не изменится. Сумма теоретических частот 145,4 близка к общему числу значений N = 147. Продолжив кривую за пределы графика, можно убедиться, что суммы теоретических и фактических частот совпадут. Сравнение фактических и теоретических частот в табл.2.13 показывает их сходство. Можно совместить частоты на одном графике, построив гистограмму по фактическим частотам n, а кривую плотности вероятности – по теоретическим частотам n т (рис.2.13). Последняя графа табл.2.13 позволяет рассчитать критерий c2: (2.40) где Nk – число классов сравнения частот. Чтобы избежать грубых случайных расхождений, рекомендуется объединять соседние классы с малым числом фактических частот. В табл.2.13 объединены два первых и два последних класса. В результате число классов Nk = 11. Вычисленное значение c2 = 2,338. Приняв вероятность a = 0,05 и зная число степеней свободы k = 11 – 3 = 8, из табл.2.11 найдем предельное значение = 15,51. Так как c2 < , то следует признать, что с вероятностью p > (1 – a) = 0,95 распределение фактических частот не противоречит нормальному закону. Такая же проверка соответствия фактических частот теоретическим может быть выполнена для других законов распределения. Возможен и другой способ проверки соответствия, но только нормальному закону. Он основан на том, что асимметрия и эксцесс нормального закона равны нулю. Оценивая степень отклонения фактических значений асимметрии и эксцесса от нуля с помощью какого-либо критерия, можно сделать заключение о соответствии или несоответствии распределения случайной величины нормальному закону. Обычно используется критерий распределения Стьюдента (см. табл.2.10), а при большом числе исходных данных – критерий нормального закона (см. табл.2.7). Проверяют две гипотезы: при числе степеней свободы k = n – 1 асимметрия А = 0 и эксцесс Е = 0. Согласно формуле (2.24) имеем (2.41) или проще , , где s А = , s Е = . Если tA и tЕ будут меньше предельного значения t пр, то распределение случайной величины не противоречит нормальному закону. Если tA или tЕ больше предельного значения t пр, то распределение противоречит нормальному закону. В качестве предельного значения можно брать t пр = 3, что соответствует вероятности q = 0,997 (см. табл.2.6). Иногда вероятность принятия решения задается (например, q = 0,95), тогда t пр определяют по табл.2.7 (t пр = 1,96). В рассматриваемом примере асимметрия и эксцесс рассчитаны по табл.2.5: А = 0,166; Е = –0,269 (см. пример 2.2). Вычислим s А = = 0,202; s Е = = 0,404, тогда tA = 0,166/0,202 = 0,82; tЕ = 0,269/0,404 = 0,67. При вероятности q = 0,95 имеем t пр = 1,96. Так как tA и tЕ меньше t пр, то еще раз получаем подтверждение того, что распределение содержаний не противоречит нормальному закону. Аналогичным способом можно проверить соответствие распределения случайной величины логнормальному распределению, оперируя с логарифмами значений случайной величины. Критерии tA и tЕ не требуют построения гистограммы, их удобно использовать при любом числе значений случайной величины, но область их применения ограничена нормальным и логнормальным законами. Критерий c2 более универсален, его можно использовать для сравнения гистограмм с любыми законами распределения. Возникает вопрос, с какой целью производится проверка гипотез о законах распределения. Ответ заключается в том, что при статистической обработке значений случайной величины нужно знать вероятность (т.е. надежность) принятия решений, а вероятность можно определить лишь тогда, когда известен закон распределения случайной величины.
2.2.8. Преобразование случайной величины
Большинство решений, принимаемых на базе статистических закономерностей, основано на нормальном законе распределения, играющем универсальную роль. Как было отмечено, при определенных условиях к нему приближаются логнормальный закон, распределение Стьюдента, распределение χ2 и многие другие. Однако реальное распределение свойств геологических объектов часто отличается от нормального, что вызывает затруднения в принятии решений и в оценке достоверности получаемых выводов. Поэтому принятию решений обычно предшествует проверка соответствия распределения случайной величины нормальному закону, и, если соответствия нет, то можно попытаться преобразовать случайную величину, приведя ее распределение к нормальному. Подобное преобразование применялось выше, когда вместо случайной величины х вводилась новая случайная величина z = ln x. В результате асимметричное логнормальное распределение преобразовывалось в симметричное нормальное. Представляют интерес такие преобразования, которые превращают произвольно распределенную случайную величину х в случайную величину z, распределение которой близко к нормальному. Задача заключается в подборе наилучшей функции преобразования. Преобразование обычно меняет область существования случайной величины. Например, если случайная величина х меняется в пределах от нуля до +¥, то преобразованная случайная величина z = ln x имеет область существования от –¥ до +¥. Поэтому учет области существования случайной величины может помочь в выборе наилучшего преобразования. Если случайная величина x имеет область существования от a до b, то преобразование (2.42) меняет пределы ее существования от –¥ до +¥, что во многих случаях эффективно. Частным случаем является ситуация, когда а = 0, b = 1 = 100 % (например, содержание химических элементов не может быть меньше нуля и больше 100 %), и формула преобразования имеет вид . (2.43) Если значения случайной величины х очень малы, то ею в знаменателе можно пренебречь, и получается формула z = ln x, лежащая в основе логнормального распределения. Наоборот, если значения х близки к единице, то получается формула преобразования в правоасимметричное логнормальное распределение z = –ln(1 – x). Если случайная величина колеблется в пределах от –1 до +1 (например, коэффициент корреляции или многие тригонометрические функции), то эффективным является преобразование (2.44) или преобразование, предложенное Фишером (2.45) Для преобразования могут быть использованы также степенные функции вида z = xa или z = x – a, где а может принимать значения от 1/2 до 3 [2].
8 Пример 2.5. Рассмотрим подбор функции преобразования для случайной величины х с асимметричным распределением (табл.2.14).
Таблица 2.14 Подбор функции преобразования случайной величины х (содержание TiO2 в магнетите)
Применим формулы преобразования: и z = ln x. Они дали положительные результаты – распределение преобразованной величины по критериям асимметрии и эксцесса не противоречит нормальному закону. Лучший результат дало преобразование z = ln x, особенно по критерию асимметрии.7
Подбор функции преобразования удобно контролировать с помощью графика пробит-функции [12]. По оси абсцисс откладывают значения случайной величины x или z, а по оси ординат – квантили нормального распределения t (рис.2.14), которые соответствуют вероятностям p = (2 n – 1)/2 N, где n – порядковый номер случайной величины в упорядоченном ряду наблюдений (табл.2.14), N = 20 – число наблюдений. На пересечении абсциссы и ординаты ставят точки. Если точки на графике расположены вдоль прямой линии, то распределение случайной величины близко к нормальному. В противном случае распределение отличается от нормального. В рассматриваемом примере график пробит-функции для исходной случайной величины явно нелинейный, а для преобразованной случайной величины z = ln x – близкий к линейному, что подтверждает эффективность выбранного преобразования. Для построения графика пробит-функции рассчитаны значения вероятностей p, а соответствующие им квантили t получены интерполяцией данных табл.2.7.
Геологические приложения
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Последнее изменение этой страницы: 2016-04-26; просмотров: 387; Нарушение авторского права страницы; Мы поможем в написании вашей работы! infopedia.su Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Обратная связь - 3.141.21.106 (0.009 с.) |