Статистичний аналіз тенденцій захворюваності в Україні - дипломная работа

Міністерство освіти і науки України

Дніпропетровський національний університет

ДИПЛОМНА РОБОТА

СТАТИСТИЧНИЙ АНАЛІЗ ТЕНДЕНЦІЙ ЗАХВОРЮВАНОСТІ В УКРАЇНІ

Виконавець: студентка групи МС-02-1

Власова В.М.

Дніпропетровськ

2007

РЕФЕРАТ

Дипломна робота: 177 с., 10 джерел, 239 рисунків, 1 додаток.

Перелік ключових слів: регресія, МНК-оцінка, паралельність лінійних регресій, збіг лінійних регресій, гіпотеза, рівень значущості, критерій, захворюваність, рік спостереження, довірчий інтервал.

Об’єкт дослідження: захворюваності на туберкульоз, рак, СНІД, гепатити А та Б в України з 1990 по 2005 роки.

Мета роботи: дослідження тенденцій захворюваності на туберкульоз, рак, СНІД, гепатити А та Б в двадцяти чотирьох областях України, А.Р.Крим, містах Київ та Севастополь та Україні в цілому у період з 1990 по 2005 роки.

Annotation

The graduation research of the fifth year student Vlasova Violetta (Dnipropetrovsk national university, mathematics and mechanics faculty, department of probability theory and statistics) deals with statistical analysis of tendencies of morbidity on tuberculosis, cancer, AIDS, hepatitis A and B in Ukraine with 1990 for 2005 years.

ЗМІСТ

ВСТУП

РОЗДІЛ 1. ЛІНІЙНА РЕГРЕСІЯ

1.1 Метод найменших квадратів. Властивості оцінок найменших квадратів. Оцінювання . Теорія розподілів. Оцінювання при наявності лінійних обмежень

1.2 F-критерій

1.3 Лінійна одновимірна регресія

1.4 Порівняння прямих регресії. Критерій паралельності прямих. Критерій збігу прямих

РОЗДІЛ 2. ДОСЛІДЖЕННЯ ЗАХВОРЮВАНОСТІ НА ТУБЕРКУЛЬОЗ, РАК, СНІД, ГЕПАТИТ А, ГЕПАТИТ Б

2.1 Дослідження захворюваності на туберкульоз (всі форми), туберкульоз легенів

2.2 Дослідження захворюваності на рак

2.3 Дослідження захворюваності на СНІД

2.4 Дослідження захворюваності на гепатит А

2.5 Дослідження захворюваності на гепатит Б

РОЗДІЛ 3. ПОРІВНЯННЯ ЗАХВОРЮВАННОСТІ СЕРЕД ОБЛАСТЕЙ УКРАЇНИ, А.Р. КРИМ, МІСТАХ КИЇВ ТА СЕВАСТОПОЛЬ

3.1 Порівняння тенденцій росту захворюваності на туберкульоз серед областей України, А.Р. Крим, м. Київ та Севастополь

3.2 Порівняння тенденцій росту захворюваності на рак серед областей України, А.Р. Крим, м. Київ та Севастополь

3.3 Порівняння тенденцій росту захворюваності на СНІД серед областей України, А.Р. Крим, м. Київ та Севастополь

3.4 Порівняння тенденцій зменшення захворюваності на гепатит А серед областей України, А.Р. Крим, м. Київ та Севастополь

3.5 Порівняння тенденцій зменшення захворюваності на гепатит Б серед областей України, А.Р. Крим, м. Київ та Севастополь

ВИСНОВКИ

СПИСОК ВИКОРИСТАНИХ ДЖЕРЕЛ

ДОДАТОК

ВСТУП

Актуальність теми. У дипломній роботі досліджується захворюваність на туберкульоз, рак, СНІД, гепатити А та Б. На сучасний момент по цим захворюванням в Україні склалась наступна ситуація.

Епідемію туберкульозу в Україні було оголошено ВООЗ в 1995 році. Щоденно в країні виявляють приблизно 82 нових випадки захворювання на туберкульоз і 30 хворих помирає від нього. За офіційними даними, кількість хворих на туберкульоз в Україні складає 1,4% населення, тобто 700 тисяч чоловік, а за експертними оцінками – близько 1,5 мільйона. З них близько 120 тисяч – з активною, найбільш заразною формою туберкульозу. Основними причинами, що призвели до масштабного поширення туберкульозу в Україні, є, в першу чергу, поглиблення соціально-економічної кризи в країні, різке зниження життєвого рівня та добробуту переважної більшості населення, безробіття, незбалансоване, недостатнє і неякісне харчування значних прошарків населення, збільшення числа соціально-дезадаптованих осіб в країні (бідних, жебраків, бродяг, мігрантів, переселенців, алкоголіків, наркоманів та інших), тих, що мають незадовільні житлові умови, занепад системи охорони здоров'я та низька санітарна культура населення. Чимало недоліків мають місце в організації та координації протитуберкульозних заходів, а також у роботі органів і закладів охорони здоров'я.

Погіршенню епідеміологічної ситуації також сприяє епідемія ВІЛ/СНІДу в Україні. Понад 30% ВІЛ-інфікованих хворіє на туберкульоз і близько 40% з них помирає від нього. Наявність ВІЛ-інфекції є потужним чинником, який сприяє розвитку активної форми туберкульозу в носіїв латентної інфекції та підвищує ймовірність розвитку рецидивуючого туберкульозу.

Україна займає одне з перших місць в Європі за темпами поширення ВІЛ-інфекції/СНІДу. За період з 1987 по 2006 роки офіційно було зареєстровано 104 645 випадків ВІЛ-інфекції, в тому числі 17 851 випадок захворювання на СНІД та 9 983 випадків смертей від захворювань, обумовлених СНІДом. Високі показники захворюваності відзначаються у Миколаївській області (86,3), Донецькій - (78,4), Дніпропетровській - (71,8), Одеській – (61,3), у м. Севастополь – (70,9) на 100 тисяч населення.

Як вважають експерти, якщо епідемію СНІДу в Україні не зупинити, зовсім незабаром країна зможе зайняти місце поряд з африканськими країнами. Також як і на африканському континенті, для більшості українців вартість противірусного лікування занадто висока. У лікарнях не вистачає навіть найнеобхідніших ліків, тому на лікування можуть розраховувати лише діти.

На початок 2006 року в онкологічних закладах системи МОЗ України на обліку перебувало 864 273 хворих, що на 3,7% більше, ніж в попередньому році.

Щорічно в Україні реєструються понад 150 тис. нововиявлень злоякісних захворювань. Найвищі показники захворюваності за останні 5 років спостерігаються у південних областях, а також Кіровоградській, Полтавській областях, АР Крим.

Майже 90 тисяч жителів України щорічно помирають від раку, при цьому 35 відсотків померлих - особи працездатного віку. Онкопатологія є причиною понад 15% усіх випадків смерті і поступається лише серцево-судинним захворюванням. У 2005 році смертність від злоякісних новоутворень в Україні склала 90 997 випадків. 38-40% онкологічних хворих помирають протягом першого року після встановлення діагнозу, що свідчить про запізнілу діагностику. У розвинених країнах світу цей показник не перевищує 30 %.

Вірусний гепатит А залишається однією з найпоширеніших інфекцій в Україні. З 1995 року має місце тенденція до зниження захворюваності на вірусний гепатит А, однак за даними ВООЗ (Всесвітня організація охорони здоров’я) щорічно за медичною допомогою звертаються від 70 до 150 тисяч осіб. В Україні за останні 5 років гепатитом перехворіло понад 375 тисяч осіб, кожний третій хворий - дитина до 15 років.

За мірою шкідливості для здоров’я населення та масштабами захворюваності вірусні гепатити в Україні займають домінуюче становище в структурі інфекційної патології. На цю групу інфекційних хвороб щорічно припадає до 20-30% захворювань, за винятком грипу та гострих інфекційних захворювань верхніх дихальних шляхів.

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

Профілактику гепатиту A та гепатиту В шляхом вакцинації в Україні розпочали нещодавно. Разом з тим усім відомо, що такий метод профілактики є найбільш ефективним.

З метою забезпечення епідемічного благополуччя населення України та попередження інфекцій, щеплення з метою запобігання захворюванню на гепатит В включене в календар обов’язкових профілактичних щеплень за віком, а вакцинація з метою профілактики гепатиту А рекомендована на ендемічних територіях та за епідемічними показаннями.

Туберкульоз - це повітряно-крапельна, а точніше повітряно-пилова інфекція. Збудником туберкульозу є мікобактерія туберкульозу. Вона відкрита у 1882 році німецьким вченим Робертом Кохом. Однак, туберкульоз може бути спричинений і мікобактеріями африканською, бичачою, мишачою, пташиною та іншими.

До туберкульозу сприйнятливі приблизно 55 видів домашніх і диких тварин, риб та біля 25 видів птахів, поміж ними та людиною можливе перехресне зараження. Зараження туберкульозом відбувається від хворих на туберкульоз людей і можливе від великої рогатої худоби, свиней, овець, кіз, коней, кролів, птиці та інших. Немає такого органу людини чи тварини, який би не уражався туберкульозом. Залежно від того, який орган уражений мікобактерії можуть виділятися від хворого на туберкульоз з харкотинням, калом, сечею, спермою, із нориць тощо. У 95 % випадків мікобактерія туберкульозу потрапляє в організм людини при вдиханні запорошеного повітря. Заразитися при вживанні заражених продуктів харчування, зокрема сирого молока, погано провареного м'яса від хворих на туберкульоз тварин. Рознощиками туберкульозу можуть бути мухи, таргани, деякі кліщі; зараження можливе також при поцілунках, допалюванні чужих цигарок, через книжки, пошкоджену шкіру і т. ін.

Мікобактерії дуже стійкі в навколишньому середовищі. Так, мікобактерії туберкульозу в молочних продуктах живуть до 240-300 днів, в пилу, в тому числі вуличному, до 3 місяців. Вони добре витримують нагрівання до 85^о С і охолодження до мінус 200^о С. При температурі мінус 23^о С мікобактерії туберкульозу зберігають життєздатність протягом 7 років, а прямі сонячні промені їх вбивають через 2 - 6 годин. Збудники туберкульозу зберігаються у річній воді до 5 місяців, у ґрунті - 1 - 2 роки, у приміщеннях при розсіяному світлі - до 1,5 місяців, у фекаліях і на пасовищах - до 1 року. Вегетативні, або бактеріальні, форми мікобактерій можуть трансформуватися у вірусні, ультрадрібні фільтруючі чи L-форми, і навпаки.

Мікобактерії туберкульозу (органи дихання, харчовий канал, шкіра тощо) спричиняють туберкульозне запальне вогнище (первинний афект). Туберкульозне запалення розвивається і водночас у запальний процес залучаються довколишні лімфатичні судини, найближчі до вогнища запалення лімфатичні вузли. Потім мікобактерії туберкульозу потрапляють в кров (бактеріємія). Мікобактерії туберкульозу розносяться лімфою і кров'ю в різні здорові органи і тканини і осідають там і спричиняють там туберкульозне запалення. Це явище називають лімфогенною й гематогенною дисемінацією туберкульозного процесу. Ця дисемінація характеризується утворенням первинних туберкульозних вогнищ у різних органах (легені, кістки, нирки тощо). Ці вогнища можуть загоюватися при сприятливій імунній опірності організму. Тоді організм набуває імунітету до туберкульозу

Рак. В останній час чітко простежується зростання захворюваності населення в зв'язку забрудненням навколишнього середовища. Особливо явно ця тенденція простежується в зростанні онкологічних захворювань на територіях, що постраждали від впливу аварії на Чорнобильській АЕС.

Ріст захворюваності злоякісними пухлинами обумовлює високу долю останніх у структурі інвалідності й смертності, а збільшення кількості осіб молодого віку і дітей серед хворих притягає до цієї проблеми все більшу увагу спеціалістів. Захворюваність, взагалі, й ураження злоякісними пухлинами особливо залежить, в першу чергу, від толерантності організму. Тому головною метою протиракової боротьби варто вважати організацію моніторингу стану здоров'я населення на основі методів раннього виявлення онкологічних захворювань.

СНІД. ВІЛ-інфекція - захворювання, яке послаблює імунітет людини - вражає переважно молодих людей віком від 18 до 27 років. Головними чинниками розповсюдження ВІЛ серед молоді є вживання ін'єкційних наркотиків і незахищені статеві стосунки.

Гепатит А. Вірус гепатиту А має кислотостійку оболонку. Це допомагає вірусам, які потрапили до організму із забрудненими продуктами і водою, пройти кислий захисний бар'єр шлунку. Вірус гепатиту А стійкий у водному середовищі, тому епідемії гепатиту А часто мають водний шлях передачі. Вірус гепатиту А відрізняється високою імуногенністю, після перенесеного захворювання формується стійкий довічний імунітет.

Гепатит А є однією з найпоширеніших інфекцій людини. У країнах з теплим кліматом і незадовільними санітарними умовами гепатитом А хворіє дуже багато людей.

Джерелом інфекції є хвора на гепатит А людина, яка з фекаліями виділяє в навколишнє середовище мільярди вірусів. При вживанні зараженої вірусом гепатиту А води або харчових продуктів (особливо погано термічно оброблених морепродуктів) віруси проникають в кишечник, потім, всмоктуючись, з потоком крові потрапляють в печінку і упроваджуються в її клітки - гепатоцити. Вірусні частинки-віріони розмножуються в цитоплазмі кліток печінки. Після виходу з кліток печінки вони потрапляють в жовчні протоки і виділяються з жовчю в кишечник.

Запальний процес в печінці, що приводить до пошкодження гепатоцитов, має імунологічну основу. Клітки імунної системи людини, Т-лімфоцити розпізнають уражені вірусом гепатоцити і атакують їх. Це приводить до загибелі інфікованих гепатоцитів, розвитку запалення (гепатит) і порушенню функції печінки.

Гепатит Б. Вірусний гепатит Б - це інфекційне захворювання печінки, обумовлене вірусом. Інфекція, обумовлена вірусом гепатиту Б, поширена повсюдно і кожен може піддатися зараженню.

Вірус гепатиту Б передається, в основному, через кров. Хронічний вірусний гепатит Б - небезпечне захворювання, яке може призвести до важких наслідків і загибелі хворого. Тому важливо своєчасно розпізнати і почати лікування гепатиту Б. Існують вакцини проти вірусу гепатиту Б, що дають надійний захист від розвитку захворювання.

Гепатит Б зустрічається найчастіше у дорослих людей (пік захворюваності доводиться на вікові групи 20-49 років). Зниження захворюваності у дітей і підлітків в розвинених країнах вдалося добитися за допомогою регулярної вакцинації.

Перенесення вірусу здійснюється з кров'ю або іншими біологічними рідинами хворої людини, які потрапляють безпосередньо в кров людини, яка інфікується. Це відбувається при спільному використовуванні колючо - ріжучих предметів (манікюрний набір, верстати для гоління), одного шприца для введення наркотиків, при пірсинге, татуаже з використанням погано обробленого інструментарію, при проведенні медичних маніпуляцій, статевим шляхом і від інфікованої матері дитині під час проходження його через родові шляхи.

Вихідні дані. Вихідними даними дослідження є офіційні дані центра медичної статистики Міністерства охорони здоров’я України.

Постановка задачі. На основі вихідних даних необхідно провести дослідження тенденцій захворюваності на туберкульоз, рак, СНІД, гепатит А та гепатит Б в Україні на протязі останніх 15 років.

В роботі необхідно розв’язати наступні задачі.

І. Методами лінійного регресійного аналізу дослідити

1) залежність захворюваності на туберкульоз (всі форми) та туберкульоз легенів від року спостереження,

2) залежність захворюваності на рак від року спостереження,

3) залежність захворюваності на СНІД від року спостереження,

4) залежність захворюваності на гепатит А від року спостереження,

5) залежність захворюваності на гепатит Б від року спостереження

(спостереження захворюваності відбувалося з 1990 по 2005 роки в кожній з 24 областей України, А.Р. Крим, м. Київ, м. Севастополь та Україні в цілому). Зробити висновки.

ІІ. Методами лінійного регресійного аналізу провести порівняння захворюваності на туберкульоз (всі форми), туберкульоз легенів, рак, СНІД, гепатит А, гепатит Б серед областей України, А.Р. Крим, м. Київ та Севастополь з метою виявлення регіонів України, в яких темпи росту або спадання захворюваності однакові або захворюваність кількісно однокова. Зробити висновки.

Методи дослідження. В роботі використовуються методи лінійного регресійного аналізу.

РОЗДІЛ 1. ЛІНІЙНА РЕГРЕСІЯ

1.1 Метод найменших квадратів. Властивості оцінок найменших квадратів. Оцінювання σ² . Теорія розподілів. Оцінювання при наявності лінійних обмежень

Нехай Y - випадкова величина, яка флуктуює навколо деякого невідомого параметра η, тобто Y = η + ε, де ε - флюктуація або „помилка". Наприклад, ε може бути „природною" флуктуацією, яка властива самому експерименту, або може бути помилкою у вимірюванні значення η.

Припустимо, що η можна подати у вигляді

η = β₀ + β₁ x₁ + … + β_{p -1} x_{p -1} ,

де х₁ , х₂ , ..., x_{p -1} - відомі постійні величини, а β_j (j = 0, 1, .., p - 1) - невідомі параметри, які підлягають оцінюванню. Якщо значення х_j , j = 0, 1, .., p – 1 змінюються і при цьому спостерігається n значень Y₁ , Y₂ , ...,Y_n змінною Y, то

Y_i = β₀ + β₁ x_{i 1} + … + β_{p -1} x _{i , p -1} + ε_i , i = 1, 2, ..., n, (1.1.1)

де x_ij i-те значенням для х_j . В матричному вигляді (1.1.1) запишеться

або

Y = Xβ + ε, (1.1.2)

де x₁₀ = x₂₀ = ... = x_{n 0} = 1.

Означення. Матриця X = Х^{( n p )} називається регресійною матрицею. При цьому значення x_ij зазвичай вибираються так, щоб стовпці цієї матриці були лінійно незалежними, тобто ранг матриці X дорівнював р. Проте в деяких випадках при плануванні експерименту елементи матриці X обираються рівними тільки нулю і одиниці, і її стовпці можуть виявитися лінійно залежними. В цьому випадку матрицю X називають матрицею плану.

Далі х_j називатимемо регресором, а Y – відкликом.

Модель (1.1.1) або (1.1.2) лінійна по відношенню до невідомих параметрів β_j , тому її називають лінійною моделлю.

Одним з методів знаходження оцінки вектора β є метод найменших квадратів. Цей метод полягає в мінімізації суми по відношенню до вектора β. Точніше, вважаючи θ = Xβ, мінімізуємо величину ε'ε = ||Y- θ||² по відношенню до θ [Х] = Ω, де Ω - образ оператора X, тобто Ω = {у: у = Хх} для деякого х. Якщо змінювати значення вектора θ в межах Ω, то ||Y- θ||² (квадрат довжини вектора Y- θ) досягає мінімуму при тому значенні θ = , для якого (Y - ) Ω (рис.1.1.1). Тому

X'(Y - ) = 0,

Або

Х' = Х'Y. (1.1.3)

Вектор визначається однозначно, оскільки він є ортогональною проекцією вектора Y на Ω. Якщо тепер стовпці матриці X лінійно незалежні, то існує

Рис. 1.1.1 Метод найменших квадратів полягає у знаходженні такої точки А, для якої відстань АВ мінімальна

єдиний вектор , для якого = X . Підставлячи в (1.1.3), одержуємо нормальне рівняння

Х'Х = Х'Y. (1.1.4)

Оскільки ми припускаємо, що матриця X має ранг р, то матриця Х'Х додатньо визначена і, отже, не вироджена. Тому рівняння (1.1.4) має єдиний розв’язок, а саме

= ( Х'Х)^-1 Х'Y

Цей розв’язок називається оцінкою найменших квадратів вектора β.

Оцінку для β можна одержати й в інший спосіб.

ε'ε = (Y-Хβ)'(Y-Хβ) = Y'Y - 2β'Х'Y+ β'Х'Хβ

(використовуємо той факт, що β'Х'Y = (β'Х'Y)' = Y'Хβ). Продиференцюємо ε'ε по β. Прирівнюючи одержану похідну ε'ε/ β нулю, приходимо до рівняння

- 2Х'Y +2Х'Хβ = 0, (1.1.5)

Або

Х'Хβ = Х'Y.

Звідки

= ( Х'Х)^-1 Х'Y

Покажемо, що знайдена стаціонарна точка є мінімумом функції ε’ε. Перепишемо (Y-Хβ)’(Y-Хβ) у вигляді

(Y-Хβ)'(Y-Хβ) = (Y-Х )'(Y-Х ) + ( - β)'Х'Х( - β). (1.1.6)

Розпишемо

(Y-Х )'(Y-Х ) + ( - β)'Х'Х( - β) = (Y'-Х' ')(Y-Х) +

+ ( ' - β')(Х'Х - Х'Хβ) = Y'Y - Y'X - 'X'Y + 'X'X +

+ 'X'X - 'X'X - 'X'X + 'X'X =

= {X'X = X'Y, оскільки - розв’язок нормального рівняння} =

= Y'Y - Y'X - 'X'Y + 'X'Y + 'X'Y - 'X'X β – β'X'Y + β'X'Xβ =

= Y'Y - Y'Xβ – β'X'Y + β'X'X β = (Y - Xβ)'(Y - Xβ)

Ліва частина в (1.1.6) досягає мінімуму при β = .

Далі позначимо = Х . Елементи вектора

e = Y – = Y – Х = (I_n - Х(Х'Х)^-1 Х')Y = (I_n - Р)Y (1.1.7)

називаються залишками (ми позначили тут скорочено Х(Х'Х)^-1 Х' через Р). Мінімальне значення ε'ε називається залишковою сумою квадратів (RSS)).

RSS = (Y - Х )'(Y - Х )= Y'Y - 2 Х' Y + 'Х'Х =

= Y’Y - 'Х' Y + '[Х'Х - Х'Y] =

= Y'Y - 'Х'Y (1.1.8)

Або

RSS = Y'Y - 'Х'Х (1.1.9)

Відмітимо, що і е єдині.

Оскільки = Х = Х(Х'Х)^-1 Х'Y = РY, то Р є матрицею лінійного перетворення, яке є ортогональним проектуванням n-мірного евклідова простору Е_n на Ω. Аналогічно I_n - Р є матрицею ортогонального проектування Е_n на - ортогональне доповнення до Ω в Е_n . Тому вираз Y = РY + (I_n - Р)Y є єдиним ортогональним розкладом вектора Y на дві складові, одна з яких лежить в Ω, а інша - в . Деякі основні властивості матриць Р і (I_n - Р) наведено в теоремі 1.1.1. Спочатку сформулюємо деякі означення.

Означення. Слідом trX матриці Х називають суму її діагональних елементів

trX = 1 + x₂₁ + x₃₂ + … + x_{np -1}

Означення. Матриця Р називається ідемпотентною, якщо Р² = Р. Симетрична ідемпотентна матриця називається проекційною. Якщо Р – проекційна матриця, то trР = rankР.

Теорема 1.1.1.

(I) Матриці Р і I_n - Р симетричні та ідемпотентнi.

(II) rank[I_n - Р] = tr[I_n - Р] = n - р.

(III) (I_n - Р)Х = 0.

Доведення.

(I) Р' = (X(X'X)^-1 X')' = X((X'X)^-1 )'X' = X(X'X)^-1 X' = P

Отже, матриця Р є симетричною і (I_n - Р)' = I_n - Р' = I_n - Р. Крім того,

Р² = X(Х'Х)^-1 Х'Х(Х'Х) ^-1 X' = XI_p (Х'Х)^-1 X' = Р,

і (I_n – Р)² = I_n - 2Р + P² = I_n – Р.

(II) Оскільки матриця I_n - Р симетрична та ідемпотентна, то вона проекційна і tr(I_n – Р) = rank(I_n – Р). Тоді

rank[I_n - Р] = tr[I_n - Р] = n - trР,

де

trР = tr[X (Х'Х)^-1 X'] = tr[Х'Х (Х'Х)^-1 ] = trI_p = р.

(III) (I_n - Р)Х = Х - Х(Х'Х)^-1 Х'Х = Х - Х = 0.

Теорема доведена.

Теорема 1.1.2.

Нехай Р = X(Х'Х)^-1 X', тоді R(P) = R(X), тобто простір, породжений стовпцями матриці P є простором, породженим стовпцями матриці Х.

Доведення.

R(P) = {z: z = Pα} для деякого α, R(X) = {Y: Y = Xγ} для деякого γ.

Вибираємо z R(P), тоді z = Pα. Отже,

z = Pα = X(X'X)^-1 X'α = Xβ,

отже z R(X).

Вибираємо Y R(X), тоді Y = Xγ

Y = Xγ = X(X'X)^-1 X'Xγ = X(X'X)^-1 X'Xγ = PY,

отже Y R(P).

Теорема доведена.

Теорема 1.1.3.

(Y - ) = 0 або

Доведення.

(Y - ) = { = X = X(X'X)^-1 X'Y = PY} = (PY)'(Y – PY) = Y'P'(1 – P)Y = = Y'P(1 – P)Y = Y'(P – P² )Y = Y'(P – P)Y = 0.

Теорема доведена.

Якщо припустити, що помилки ε такі, що , то

M[ ] = (X’X)^-1 X’M[Y] = (X’X)^-1 X’X β = β (1.1.9)

тобто є незміщеною оцінкою вектора β. Якщо, окрім того, припустити, що всi ε_i , і = 1, …, n - некорельовані і мають однакову дисперсію, тобто

соv[ε_i , ε_j ] = ,

то D[ε] = σ² I_n ,

D[Y] = D[Y - Xβ] = D[ε], отже D[Y] = σ² I_n .

Звідси одержуємо

D[ ] = D[(Х'Х)^-1 Х'Y] = сov((Х'Х)^-1 X'Y, (Х'Х)^-1 X'Y) =

= (X'X)^-1 X'cov(Y,Y)((X'X)^-1 X')' = (X'X)^-1 X'DYX(X'X)^-1 =

= (X'X)^-1 X'σ² IX(X'X)^-1 = σ² (X'X)^-1 (X'X) (X'X)^-1 = σ² (X'X)^-1 (1.1.10)

Виникає таке питання: чому за оцінку вектора β ми вибираємо саме (оцінку найменших квадратів), а не будь – яку іншу оцінку? Далі покажемо, що в деякому розумному класі оцінок _j , є оцінкою параметра β_j з найменшою дисперсією. Цю оцінку _j можна „виділити" з вектора = ( ₀ , ₁ , ..., _{p -1} )' множенням зліва на вектор-рядок c', у якого (j +1)-й елемент рівний одиниці, а всі інші елементи дорівнюють нулю. Таку специфічну властивість оцінки _j , можна узагальнити на випадок довільної лінійної комбінації а' . Для цього використовуємо наступну теорему.

Теорема 1.1.4.

Нехай - оцінка найменших квадратів вектора = Хβ. Тоді в класі всіх лінійних незміщених оцінок лінійної комбінації c'θ оцінка c' є єдиною оцінкою, яка має мінімальну дисперсію. (Будемо говорити, що c' є найкращою лінійною незміщеною оцінкою (НЛНО) для c'θ)

Доведення.

Оцінку найменших квадратів вектора = Хβ представимо у вигляді

= X = X(Х'Х)^-1 X'Y = X(Х'Х)^-1 X'Y = PY,

при цьому

PX = X(Х'Х)^-1 X'X = X(Х'Х)^-1 X'X = XI = X .

Перевіримо, що c' - лінійна незміщена оцінка для c'θ. Дійсно,

M[c' ] = Mc'РY = c'P MY = c'Pθ = c'PXβ = c'Xβ = c'θ

для всіх θ Ω = [Х] і c' = c'PY = (P'c)'Y = (Рс)'Y. Розглянемо іншу лінійну незміщену оцінку для c'θ. Тоді M[d'Y] = c'θ з одного боку, а з іншого

M[d'Y] = d'MY = d'θ,

Тоді

c'θ = d'θ (с' - d')θ = 0 (с- d)'θ = 0, тобто (c - d) Ω = R(X).

Оскільки R(X) = R(P) в силу теореми 1.1.2, то

(c – d) R(P), (c – d)'P = 0 ((c – d)'P)' = 0' P(c – d) = 0

Pc = Pd

Порахуємо дисперсію оцінки c' :

Dc' = D[(Рd)'Y] = D[(Рd)'Y] = Dd'P'Y = cov(d'P'Y, d'P'Y) =

= d'P'cov(Y, Y)(d'P')' = d'PDYPd = d'Pσ² IPd = σ² d'Р² d = σ² d'Рd,

Тоді

D[d'Y] - D[c' ] = D[d'Y] - D[(Рd)' Y] =

= d'DYd - σ² d'Pd = σ² d'd - σ² d'Pd =

= σ² (d'd - d'Рd) = σ² d'(I_n - Р)d = {I_n – P = (I_n – P)² } =

= σ² d'(I_n - Р)(I_n - Р)d = {I_n – P = (I_n – P)'} =

= σ² d'(I_n - Р)'(I_n - Р)d = σ² [(I_n - Р)d]'[(I_n - Р)d] ≥ 0

Рівність нулю досягається тоді й тільки тоді, коли

(I_n - Р)d = 0

d – Pd = 0

d = Рd = Рс

Тоді D(d'Y) ≥ D(c' ), при цьому c'θ = d'θ. Це і означає, що c' має мінімальну дисперсію і є єдиною оцінкою з такою властивістю в класі всіх лінійних незміщених оцінок лінійних комбінацій c'θ.

Теорема доведена.

Теорема доведена в припущенні, що матриця X має ранг p, так що Р = X (Х'Х)^-1 X', і θ =Хβ випливає, що β = (Х'Х)^-1 Х'θ.

Нехай с' = а'(Х'Х)^-1 X', тоді звідси оцінка а'β = a'(X’X)^-1 X' = с' є НЛНО з мінімальною дисперсією для а'β при кожному а.

Зауваження. Якщо похибки ε_і незалежні й однаково розподілені ε ~ або, в еквівалентній формі, Y ~ , то a' має мінімальну дисперсію серед усіх незміщених оцінок, а не тільки в класі лінійних незміщених оцінок.

Зокрема, МНК – оцінка _і , і = 0, …, p – 1 є також оцінкою максимальної правдоподібності, і вона ефективна оцінка для β_і .

Якщо ж розподіл ε_i не є нормальним, то МНК – оцінка _і відрізняється від оцінки максимальної правдоподібності. В цьому випадку МНК – оцінка _і асимптотично ефективна для β_і .

Оцінимо параметр σ² = Dε_i , але спочатку сформулюємо низку лем.

Лема 1.1.1. Нехай Y = Y^{( n ×1)} – випадковий вектор, А^{( n × n )} = A – симетрична матриця. Якщо MY = θ, DY = ∑, тоді математичне сподівання квадратичної форми Y'AY дорівнює

M(Y'AY) = tr(A∑) + θ'Aθ

.Наслідок

Якщо ∑ = σ² I, то tr(A∑) = σ² trA.

Лема 1.1.2.

Нехай маємо n незалежних випадкових величин Y₁ , Y₂ , …, Y_n з середніми θ₁ , θ₂ , …, θ_n , однаковими дисперсіями μ₂ та однаковими третіми та четвертими центральними моментами μ₃ та μ₄ відповідно (μ_r = M(Y_i – θ_i )^r ). Якщо A = = А^{( n × n )} – симетрична матриця, а a – вектор – стовпець, утворений її діагональними елементами, тоді дисперсія квадратичної форми Y'AY дорівнює

D(Y'AY) = (μ₄ – 3(μ₂ )² )a'a + 2(μ₂ )² trA² + 4(μ₂ )² θ'A² θ + 4μ₃ θ'Aa

Теорема 1.1.4.

Якщо

М[Y] = Xβ, де Х = X^{( n × p )} , rangX = p, D[Y] = σ² I_n ,

тоді оцінка

є незміщеною оцінкою для σ² .

Доведення.

Похибку ε запишемо у вигляді:

ε = Y - = Y - Х = { = (X'X)^-1 X'Y } = Y – X(X'X)^-1 X'Y =

= (I_n – X(X'X)^-1 X')Y = (I_n - Р)Y.

Тоді

(n - p)S² = (Y - X )'(Y - X ) = ((I_n – P)Y)'((I_n – P)Y) = Y'(I_n – P)'(I_n – P)Y = {(I_n – P)' = I_n – P – симетрична} =Y'(I_n – P)² Y = Y'(I_n – P)Y.

Виразимо Y'(I_n – P)Y з рівності:

(Y – Xβ)'(I_n – P)(Y – Xβ) = Y'(I_n – P)Y – Y'(I_n – P)Xβ – (Xβ)'(I_n – P)Y + (Xβ)'(I_n – P)Xβ;

Y'(I_n – P)Y = (Y – Xβ)'(I_n – P)(Y – Xβ) + Y'(I_n – P)Xβ + (Xβ)'(I_n – P)Y - (Xβ)'(I_n – P)Xβ.

Порахуємо M(n – p)S²

M(n – p)S² = MY'(I_n – P)Y = {лема 1.1.1} = M(Y – Xβ)'(I_n – P)(Y – Xβ) +

+ MY'(I_n – P)Xβ + M(Xβ)'(I_n – P)Y – M(Xβ)'(I_n – P)Xβ =

= M(Y – Xβ)'(I_n – P)(Y – Xβ) + (Xβ)'(I_n – P)Xβ + (Xβ)'(I_n – P)Xβ –

- (Xβ)'(I_n – P)Xβ = M(Y – MY)'(I_n – P)(Y – MY) =

= + (Xβ)'(I_n – P)Xβ =

= + (Xβ)'(I_n – P)Xβ =

= + (Xβ)'(I_n – P)Xβ =

= σ² (p₁₁ + p₂₂ + … + p_nn ) + β'X'(I_n – P)Xβ =

= σ² tr(I_n – P) + β'X'(I_n – P)Xβ = =

= σ² (n – p) + 0 = σ² (n – p)

Отже,

M(n – p)S² = σ² (n – p) MS² = σ² .

Теорема доведена.

Виявляється, що S² , подібно до , має певні властивості оптимальності, які наведено в наступній теоремі.

Теорема 1.1.5.

Нехай Y₁ , Y_2, …, Y_n – незалежні випадкові величини, які мають однакові дисперсії μ₂ = 3σ² і однакові треті та четверті моменти μ₃ і μ₄ . Якщо M[Y] = Xβ, де матриця Х = Х^{( n × p )} , rangX = p, то DY = σ² I і (n – p)S² є єдиною невід’ємною квадратичною незміщеною оцінкою для (n – p)σ² , яка має мінімальну дисперсію при μ₄ = 3σ⁴ або при рівності всіх діагональних елементів матриці P.

Доведення.

Оскільки σ² > 0, то будемо розглядати тільки невід’ємні оцінки.

Нехай Y'АY незміщена квадратична оцінка для (n - р)σ² . Порахуємо математичне сподівання та дисперсію оцінки Y'АY

(n - р)σ² = M[Y'АY] = σ² trА + β'Х'АХβ

для всіх β, тоді trА = n - р і β'Х'АХβ = 0 для всіх β. Отже, Х'АХ = 0 А- додатньо напіввизначена симетрична матриця з Х'АХ = 0 випливає, що АХ = 0.

Позначимо а – вектор, утворений діагональними елементами матриці А і γ₂ = (μ₄ - 3σ⁴ )/σ⁴ , тоді згідно з лемою 1.1.2,

D[Y'АY] = (μ₄ – 3(μ₂ )² )a'a + 2(μ₂ )² trA² + 4(μ₂ )² (Xβ)'A² (Xβ) + 4μ₃ (Xβ)'Aa =

= = (μ₄ – 3(μ₂ )² )a'a + 2(σ² )² trA² + 4(σ² )² β'X'AXβ +

+ 4μ₃ β'(AX)'a = σ⁴ γ₂ а'а + 2σ⁴ trА² . (1.1.11)

Далі розглянемо оцінку (n - р)S² , яка належить класу незміщених квадратичних оцінок для (n - р)σ² згідно з теоремою 1.1.4

(n - р)S² = (Y - X )’(Y - X ) = Y(I_n - Р)Y = Y'RY

(де для стислості, введене позначення I_n - Р = R), trR² = trR = n - р.

Розглянемо D[Y'RY]:

D[Y'RY] = σ⁴ γ₂ r'r + 2σ⁴ trR² = σ⁴ γ₂ r'r + 2σ⁴ (n - р). (1.1.12)

де r – вектор, утворений діагональними елементами матриці R.

Для того, щоб знайти достатні умови для мінімальності дисперсії оцінки Y'АY, покладемо А = R + D. Оскільки A та R симетричні, то матриця D також симетрична і trА = trR + trD.

Підставляємо: (n – p) = (n – p) + 0 таким чином, trD = 0. Оскільки АХ = 0, то АР = АХ(Х'Х)^-1 X' = 0, тоді

A = R + D

AP = RP + DP

AP = P – P² + DP

0 = P – P + DP

DP = 0

Тоді

DR = D – DP = D – 0 = D

(останнє рівне також D = D' = RD, так як D симетрична).

Позначимо a = r + d, r – вектор діагональних елементів матриці R, d– вектор діагональних елементів матриці D.

A² = (R + D)² = R² + DR + RD + D² = R + 2D + D²

tr A² = trR + 2trD + trD² = (n - р) + trD² .

Підставляючи а = r + d і tr A² в (1.1.11), одержуємо

D[Y'АY] = σ⁴ γ₂ a'а + 2σ⁴ trA² = σ⁴ γ₂ (r + d)'(r + d) + 2σ⁴ (n – p + trD² ) =

= σ⁴ γ₂ (r' + d')(r + d) + 2σ⁴ (n – p + trD² ) =

= σ⁴ γ₂ (d'r + d'd + r'r + r'd) + 2σ⁴ (n – p + trD² ) =

= σ⁴ γ₂ r'r + 2σ⁴ (n – p) + 2σ⁴ =

= D[Y'RY] + 2σ⁴ .

Щоб знайти оцінку з мінімальною дисперсією, потрібно мінімізувати D[Y'АY] за умов tr D = 0 і DR = D. У загальному випадку виконати таку мінімізацію досить важкою. Проте в двох важливих окремих випадках ця мінімізація виконується не важко. Перший випадок - це ситуація, коли γ₂ = 0 При цьому

D[Y'AY] = D[Y'RY] + 2σ²

Остання ж величина досягає мінімуму, коли d_ij = 0 для всіх i, j, тобто коли D = 0 і А = R. Другий випадок - це випадок рівності всіх діагональних елементів матриці Р. При цьому всі вони рівні р₁₁ = p₂₂ = … = p_nn

trR = trI – trP = n – p tr Р = р.

Тому

р₁₁ + p₂₂ + … + p_nn r_ii = p

np_ii = p p_ii = p/n

Тоді діагональні елементи матриці R = (I – P) дорівнюють r_ii = 1 – p_ii = 1 – p/n = (n - р)/n для кожного і

D[Y'AY] = D[Y'RY] + 2σ⁴ ( =

= =

= D(Y'RY) + 2σ⁴ =

= D[Y'RY] + 2σ⁴ , (1.1.13)

Далі для будь–якої випадкової величини ξ виконується нерівність γ₂ ≥-2. Дійсно,

0 ≤ D(ξ – Mξ)² = M(ξ – Mξ)⁴ – (M(ξ - Mξ)² )² = μ₄ – (μ₂ )² =

= μ₄ – 3(μ₂ )² + 2(μ₂ )² = (μ₂ )² (μ₄ / (μ₂ )² – 3 + 2) =

= = (μ₂ )² (γ₂ + 2), отже γ₂ ≥ -2

отже D[Y'АY] досягає мінімуму, коли d_ij = 0 для всіх i, j. Таким чином, в обох випадках дисперсія виявляється мінімальною тоді і тільки тоді, коли А = R. Теорема доведена. Доведена теорема говорить про те, що незміщена квадратична оцінка для σ² , з мінімальною дисперсією існує тільки при певних обмеженнях, наведених в теоремі. У припущенні нормальності, тобто при γ₂ = 0, оцінка S² є незміщеною оцінкою для σ² , яка має мінімальну дисперсією в класі всіх незміщених оцінок, а не тільки в класі квадратичних незміщених оцінок. Раніше ми припускали відносно похибок ε_i , що M[ε] = 0 і D[ε] = σ² I_n . Якщо додатково припустити, що похибки ε_i розподілені нормально, тобто ε ~ N_n (0, σ² I_n ) (отже Y ~ N_n (Xβ, σ² I_n )), то можна одержати низку наступних результатів, пов'язаних з розподілами.

Теорема 1.1.6. Якщо Y ~ N_n (Xβ, σ² I_n ), де Х = Х^{( n × p )} , rangX = p, тоді

(I) ~ N_p (β, σ² (X'X)^-1 );

(II) ( - β)'X'X( - β)/σ² ~ ;

(III) не залежить від S² ;

(IV) RSS/σ² = (n – p)S² /σ² ~ .

Доведення. (I) МНК – оцінка вектора β має вигляд = (Х'Х)^-1 Х'Y, тоді = СY, де C = (Х'Х)^-1 Х' - матриця розміру р×n, для якої rangС = rang(Х'Х)^-1 Х' = rangХ^-1 (Х')^-1 X' = rangХ^-1 = p. Вектор Y ~ N_n (Xβ, σ² I_n ). Генератриса моментів для вектора дорівнює

M = M .

M(t) = M = M = = M = = =

- генератриса моментів , де cXβ = (X’X)^-1 β = β,

cσ² Ic’ = (X'X)^-1 X'σ² I((X'X)^-1 X')' = σ² (X'X)^-1 X'X(X'X)^-1 = σ² (X'X)^-1 .

Генератриса функції моментів нормального розподілу ξ ~ N(a; σ² ):

M(t) = Me^tξ = ,

Генератриса моментів для вектора однозначно визначає щільність розподілу вектора і дорівнює M(t) = Me^{t '} , , t = (t₁ , t₂ , …, t_p )'

(II) ( - β)'Х'Х( - β)/σ² = =

= ( - β)'(D )^-1 ( - β) = ( ₁ – β₁ , …, _p – β_p )(D )^-1 =

= (D )^-1

~ N(β; σ² (X'X)^-1 ),

- β ~ N(0; σ² (X'X)^-1 ),

, тоді . Отже, .

(III) Необхідно довести, не залежить від S² . Порахуємо cov( ,Y-X )

cov( , Y - X ) = cov((X'X)^-1 X'Y, Y – X(X'X)^-1 X'Y) =

= cov((X'X)^-1 X'Y, Y - PY) = cov((X'X)^-1 X'Y, (I – P)Y) =

= (X'X)^-1 X'cov(Y, Y)(I – P)' = {(I – P)' = I – P} =

= (X'X)^-1 X'DY(I – P) = {DY = σ² } = (X'X)^-1 X'σ² I(I – P) =

= σ² (X'X)^-1 X'(I – P) = = 0

Залишилось скористатись наступною теоремою:

Нехай Y ~ N(Xβ; σ² I), U = AY, V = BY, матриця А₁ складена з лінійно незалежних рядків матриці А, U₁ = A₁ Y. Якщо cov(U, V) = 0, то

1) випадковий вектор U₁ не залежить від V'V;

2) випадкові величини U'U та V'V незалежні.

Позначимо

U₁ = , V = Y - X , U = U₁ =

U₁ = (X'X)^-1 X'Y, V = Y - X = (I – P)Y.

Оскільки cov(U₁ , V) = 0, тоді U₁ = не залежить від V'V=(Y - X )'(Y - X ) = = (n – p)S² .

(IV) Розглянемо

Q₁ = (Y – Xβ)'(Y – Xβ) = (Y - X + Х( - β))'(Y - X + X( - β)) =

= (Y – X )'(Y – X ) + (Y – X )'X ( - β) + ( - β)'X'(Y - X) +

+ ( - β)'X'X ( - β) =

= =

= (Y – X )'(Y – X ) + ( - β)'X'X ( - β) = Q + Q₂ . (1.1.15)

Тут ми позначили

(Y - X )'(Y - Х ) = Q, ( - β)'Х'Х( - β) = Q₂ .

При цьому відношення

Q₁ /σ² = = (ε_i ~ N(0; σ² ), ε_i /σ ~ N(0; 1)), Q₂ /σ² ~ .

Отже, Q = Q₁ + Q₂ , Q₁ ~ , Q₂ ~ (n > p). Тому Q/σ² = Q₁ /σ² – Q₂ /σ ~ ~ .

Теорема доведена.

Нехай лінійна модель регресії має вигляд Y = Xβ + ε, X = X^{( n × p )} , rangX = p, ε ~ N(0; σ² I).

Необхідно оцінити параметр β, при лінійних обмеженнях H: Aβ = c,

де А = А^{( q × p )} – відома матриця, c = c^{( q ×1)} – відомий вектор. (1.1.16)

Обмеження (1.1.16) можна переписати у вигляді:

H: Aβ = c

H: β = ,

де a'_i – i-тий рядок матриці А

H: a'_i β = c_i , i = 1, 2, …, q.

Використаємо метод множників Лагранжа для розв’язання цієї задачі.

В подальшому будемо використовувати такий вираз:

λ₁ (a'₁ β – с₁ ) + λ₂ (a'₂ β – с₂ ) + … + λ_q (a'_q β – с_q ) =

= (λ₁ , λ₂ , …, λ_q ) = λ'(Aβ – c) = (λ'(Aβ – c))' =

= (Aβ – c)'λ = (β'A' - c')λ (1.1.17)

Мінімізуємо суму квадратів залишків ε'ε при лінійних обмеженнях H:

Aβ = c.

r = ε'ε + λ₁ (a'₁ β – с₁ ) + … + λ_q (a'_q β – с_q ) = ε'ε + (β'A' - c')λ = (Y – Xβ)'(Y – Xβ) + (β'A' - c')λ = (Y' – X'β') (Y – Xβ) + (β'A' - c')λ = Y'Y - Y'Xβ - β'X'Y + β'X'Xβ + (β'A' - c')λ = Y'Y - 2β'X'Y + β'X'Xβ + β'A'λ - c'λ

З (1.1.18) випливає, що

X'Xβ = X'Y - A'λ

= (X'X)^-1 X'Y - (X'X)^-1 A' (1.1.20)

= - (X'X)^-1 A' (1.1.21)

Формулу (1.1.21) підставляємо в (1.1.19)

c = A = A - (X'X)^-1 A'

c - A = - (X'X)^-1 A'

(A(X'X)^-1 A')^-1 (c - A ) = -

Останнє підставляємо в (1.1.21)

= + (X'X)^-1 A'(A(X'X)^-1 A')^-1 (c - A )

мінімізує ε'ε при обмеженнях Aβ = c.

1.2 F-критерій

Розглянемо лінійну модель Y =Хβ + ε, в якій матриця X має розмір n р і ранг р, ε ~ N_n (0, σ² I_n ). Нехай ми хочемо перевірити гіпотезу H: Аβ = c, де А - відома (q p) - матриця рангу q, а с - відомий (q 1) - вектор. Позначимо

RSS = (Y –X )'(Y-X ) = (n – p)S²

RSS_H = (Y –X _H )'(Y-X _H )

Де _H = + (Х'Х)^-1 А'[А(Х'Х)^-1 А']^-1 (с-А ), (1.2.1)

і RSS_H - мінімальне значення ε'ε при обмеженнях Аβ = с.

Теорема 1.2.1.

(I) RSS_H - RSS = (А - c)' [А (Х'Х)^-1 А']^-1 (А - c),

(II) М [RSS_H - RSS] = σ² q + (Аβ -с)' [А(Х'Х)^-1 А']^{- 1} (Аβ - с).

(III) Якщо гіпотеза Н: Аβ = с справедлива, то статистика

F =

має розподіл Фішера F_{q , n - p} (F-розподіл з q і n - p ступенями вільності відповідно).

(IV) Якщо с = 0, то статистика F приймає вигляд

F = ,

де Р_H - симетрична і ідемпотентна матриця і Р_Н P = PР_Н = Р_Н

Доведення.

(I) Спочатку доведемо тотожність:

||Y - X _H ||² = ||Y - X ||² + ||X( - _H )||²

Розглянемо

||X( - )||² = (X( - β))'X( - β) = ( - β)'X'X ( - β) = ( - _H + _H - β)'X'X ( - _H + _H - β) =

= ( - _H )'X'X ( - _H ) + ( _H - β)'X'X ( _H - β) =

= 2 ((X'X)^-1 A' )'X'X( _H - β) = A(X'X)^-1 X'X( _H - β) = A( _H - β) = (A _H - Aβ) = (c – c) = 0 = (X( - _H ))'X( - _H ) + (X( _H - β))'X( _H - β) = ||X( - _H )||² + ||X( _H - β)||² .

Далі,

ε'ε = (Y – Xβ)'(Y – Xβ) = ||Y – Xβ||² = (Y - X )'(Y - X ) +

+ ( - β)'X'X( - β) = ||Y - X ||² + ||X( - β)||²

Підставляємо

||X( - β)||² :

ε'ε = ||Y - X ||² + ||X( - _H )||² + ||X( - β)||²

ε'ε досягає мінімального значення при ||X( - β)||² = 0, тобто

X( - β) = 0

β = , Х ≠ 0 (оскільки стовпці Х лінійно незалежні)

Покладаючи в ε'ε β = , знаходимо

||Y - X _H ||² = ||Y - X ||² + ||X( - _H )||²

Тоді

RSS_H – RSS = (Y - X _H )'(Y - X _H ) – (Y - X )'(Y - X ) =

= ||Y - X _H ||² - ||Y - X ||² = ||X( - _H )||² = (X( - _H ))'(X( - _H )) =

= ( - _H )'X'X( - _H ) =

= =

= ((X'X)^-1 A'(A(X'X)^-1 A')^-1 (A - c))'X'X((X'X)^-1 A'(A(X'X)^-1 A')^-1 (A - c)) =

= (A - c)'(A(X'X)^-1 A')^-1 A(X'X)^-1 (X'X)(X'X)^-1 A'(A(X'X)^-1 A')^-1 (A - c) =

= (A - c)'(A(X'X)^-1 A')^-1 (A - c).

(II) Скористаємось лемою.

Нехай Y = Y^{( n ×1)} - випадковий вектор, A^{( n × n )} = A - симетрична матриця. Якщо MY = θ, DY = ∑, тоді

M(Y'AY) = tr(A∑) + θ'Aθ.

Раніше, доведено, що

~ N_p (β, σ² (Х'Х)^-1 ), A ~ N_q (Aβ, σ² A(Х'Х)^-1 A').

Позначимо Z = А - c і В = А(Х'Х)^-1 А'. Тоді

M[Z] = M(А – c) = A - c = = Аβ – c і

D[Z] = D(А – c) = D[A ] = σ² B

Тоді

M[RSS_H - RSS] = M[Z'В^{- 1} Z] = tr[σ² В^{- 1} В] + (Аβ - с)' В^{- 1} (Аβ - с) =

= tr[σ² I_q ] + (Aβ – c)'B^-1 (Aβ – c) =

= σ² q + (Aβ – c)'B^-1 (Aβ – c). (1.2.2)

(III) Відомо, що ~ N_q (β,σ² А(Х'Х)^-1 ), тоді

A ~ N_q (Aβ, σ² A(Х'Х)^-1 A') і

А - с ~ N_q (Aβ - c, σ² A(Х'Х)^-1 A'),

, тоді .

Розглянемо (RSS_H – RSS)/σ²

= (А - с)' (D[А ])^-1 (А - с),

Раніше доведено, що RSS/σ² ~ (теорема 1.1.6 (IV)), тоді статистика

при справедливій гіпотезі Н має вигляд [ /q]/[ /(n - р)]. Отже, якщо гіпотеза Н справедлива, то F ~ F_{q , n - p} .

(IV) Нехай у виразі (1.2.1) c = 0, тоді маємо

= X( _H - (Х'Х)^-1 А'[А(Х'Х)^-1 А']^-1 А ) = X -

- X(Х'Х)^-1 А'[А(Х'Х)^-1 А']^-1 А =

= X(Х'Х)^-1 Х'Y - X(Х'Х)^-1 А'[А(Х'Х)^-1 А']^-1 А(Х'Х)^-1 Х'Y =

={X(Х'Х)^-1 X' - X(Х'Х)^-1 А'[А(Х'Х)^-1 А']^-1 А(Х'Х)^-1 Х'}Y= (Р-Р₁ )Y, (1.2.3)

Тобто

(1.2.4)

де Р_H - симетрична матриця. Спростивши вираз для матриці Р₁ , знаходимо, що Р₁ симетрична і ідемпотентна і Р₁ Р = РР₁ = Р₁ . Звідси одержуємо

= Р² - Р₁ P – РP₁ + = P - 2P₁ + P₁ = P - P₁ = P_H (1.2.5)

P_H P = (P - P₁ )P = P - P₁ = P_H (1.2.6)

і РР_H = Р_H (останнє одержуємо транспонуванням).

Y - X = Y - X(Х'Х)^-1 Х'Y = Y(I - X(Х'Х)^-1 Х') = (I – P)Y.

Тоді

RSS = (Y - X )'(Y - X ) = ((I – P)Y)'(I – P)Y =

= Y'(I – P)'(I – P)Y = Y'(I_n - Р)Y

Aналогічно

RSS_H = (Y - X _H )'(Y - X _H ) = Y'(I_n – Р_H )Y. (1.2.7)

Таким чином,

RSS_H – RSS = Y'(I_n – Р_H )Y - Y'(I_n - Р)Y = Y'(I – Р_H – I + P)Y = Y'(P – Р_H )Y.

Отже,

Теорема доведена.

F – критерій для перевірки гіпотези H: Aβ = c.

Гіпотезу H: Aβ = c відхиляють при

і не відхиляють в супротивному разі. Рівень значущості критерію α.

1. 3 Лінійна одновимірна регресія

Нехай Y_i = β₀ + β₁ x_i + ε_i (i = 1,2, …, n) і ми хочемо перевірити гіпотезу H: β₁ = 0. Тоді X = (I_n , x),

, ,

Підставляючи ці вирази у формулу = (X’X)^-1 X’Y, після деяких спрощень одержуємо

₀ = (1.3.1)

= ₀ + ₁ x_i = (1.3.2)

Нарешті, знаходимо вираз для F- статистики

(1.3.3)

Де

Помітимо, що з (1.3.4) випливає, що

Де

є квадратом вибіркового коефіцієнта кореляції між Y і х. Відношення r є мірою ступені лінійності зв'язку меж Y і х, оскільки, згідно з (1.3.5),

RSS = (1.3.6)

Отже, чим більше значення г² , тим менше RSS і, тим краще підібрана пряма відповідає спостереженням.

1.4 Порівняння прямих регресій. Критерій паралельності прямих. Критерій збігу прямих

Нехай необхідно порівняти K ліній регресій

Y = α_k + β_k x_k + ε (k =1, 2, ..., K),

де M[ε] = 0 і дисперсії D[ε] = σ² однакові для всіх K ліній. Якщо для k-й лінії є n_k пар спостережень (x_ki , Y_ki ) (i = 1, 2, ..., n_k ), то модель приймає вигляд

Y_ki = α_k + β_k x_ki + ε_ki (k =1, 2, ..., n_k ), (1.4.1)

де ε_ki - незалежні випадкові величини з розподілом N(0, σ² ).

Введемо позначення Y' = (Y₁₁ , Y₁₂ , …, Y_{1 n 1} , …, Y_{Kn 1} , …, Y_Knk ) запишемо модель у вигляді Y = Xγ + ε, де

Тут X -матріца розміру N×2K рангу 2К, а N = .

Використовуючи загальну теорію підрозділу 1.2, можна перевірити будь-яку гіпотезу вигляду Н:Аγ = с. Дві гіпотези такого роду розглядаються нижче.

Критерій паралельності прямих

Розглянемо задачу перевірки паралельності всіх K ліній. Тоді гіпотеза

Н:Аγ = с має вигляд H₁ : β₁ = β₂ = . . . = β_K = β, або β₁ - β_K = β₂ – β_K = ... = β_{K -1} - - β_K = 0. У матричній формі H₁ приймає вигляд

або Аγ = 0, де А-матрица розміру (К- 1)×2K рангу K-1. Використовуючи загальну теорію регресії з q = K-1, n = N і р = 2К, одержуємо, що статистика критерія для перевірки гіпотези H₁ , має вигляд

(1.4.2)

Для знаходження RSS необхідно мінімізувати

S = ε'ε = .

Продиференціюємо S по α_k та β_k

З перших K рівнянь системи знаходимо, що

Підставляємо α_k в (1.4.4)

(1.4.5)

Тоді Підставляємо оцінки в S. Знаходимо

RSS =

(1.4.6)

Для знаходження мінімізуємо суму S = відносно α і β. Продиференціюємо S по α_k та β

(k = 1, 2, …, K)

З перших K рівнянь системи знаходимо

Підставляємо знайдені α_k в друге рівняння системи

Тоді

Підставляючи оцінки в S, знаходимо

Отже, гіпотезу H₁ : β₁ = β₂ = . . . = β_K = β відхиляємо при

і не відхиляємо в супротивному разі. Рівень значущості критерію α.

Критерій збігу прямих

Для вирішення питання про збіг всіх K ліній регресії розглянемо гіпотезу H₂ : α₁ = α₂ = ... = α_K і β₁ = β₁ = ... = β_K , або, що те ж саме, H₂ : α₁ - α_K = = α₂ - α_K = ... = α_{K -1} - α_K = 0 і β₁ - β_K = ... = β_{K -1} - β_K = 0.

В матричній формі гіпотеза H₂ приймає вигляд:Aγ = 0, зокрема

Матриця А має розмір (2K-2) рядків, 2K стовпців (ранг 2М – 2).

F – статистика для перевірки гіпотези H₂ має вигляд

(1.4.7)

Знайдемо RSS_{H 2} . Для знаходження RSS_{H 2} необхідно мінімізувати суму по α і β, одержуємо

Із першого рівняння знаходимо

Підставляючи α в перше рівняння системи одержуємо β

Тоді

Отже,

Отже, гіпотезу H₂ :α₁ = α₂ = ... = α_K і β₁ = β₁ = ... = β_K відхиляємо при

> F_{α ,(2 K -2),( N -2 K )}

і не відхиляємо в супротивному разі. Рівень значущості критерію α.

На практиці спочатку застосовують критерій паралельності прямих, а потім вже, якщо гіпотеза Н₁ не відхиляється, перевіряють гіпотезу H₂ (вважаючи, що гіпотеза H₁ справелива), використовуючи статистику

.

РОЗДІЛ 2. ДОСЛІДЖЕННЯ ЗАХВОРЮВАНОСТІ НА ТУБЕРКУЛЬОЗ, РАК, СНІД, ГЕПАТИТ А, ГЕПАТИТ Б

2.1 Дослідження захворюваності на туберкульоз (всі форми), туберкульоз легенів

Необхідно перевірити гіпотезу про наявність лінійної залежності між y- захворюваністю на туберкульоз на 100 000 населення та х – роком проведення спостереження, оцінити невідомі параметри лінійної одновимірної регресії , побудувати довірчі інтервали для параметрів .

Захворюваність на туберкульоз в кожній з 24 областей, А.Р.Крим, м.Київ та м. Севастополі, а також в Україні в цілому розглядається окремо.

Для кожної області дослідження залежності захворюваності на туберкульоз від року спостереження проводиться так:

1) знаходяться МНК- оцінки коефіцієнтів регресії ;

2) перевіряється гіпотеза H₀ : про значущість вільного члена регресії за допомогою t – критерію;

3) перевіряється гіпотеза H₀ : про значущість коефіцієнта при незалежній змінній x за допомогою F- критерію та t- критерію;

4) перевіряється гіпотеза H₀ : про значущість коефіцієнта кореляції за допомогою t- критерію;

5) будується довірчий інтервал для з коефіцієнтом довіри 0,95;

6) будується довірчий інтервал для з коефіцієнтом довіри 0,95;

7) будується довірчий інтервал для з коефіцієнтом довіри 0,95.

А. Р. Крим

Рис.2.1.1 Захворюваність на туберкульоз на 100000 населення в А.Р. Крим (всі форми)

Лінійна регресія 4,2862х – 8512,2 значуща.

F = 124,76 > = 4,6; t₁ = 11,17 > = 2,145, гіпотеза H₀ :b₁ = 0 відхиляється на рівні значущості α = 0,05.

t₀ = |-11,105| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 11,17 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Вінницька область

Рис. 2.1.2 Захворюваність на туберкульоз на 100000 населення в Вінницькій області, (всі форми)

Лінійна регресія 2,5399х – 5024 значуща.

F = 82,203 > = 4,6, t₁ = 9,07 > = 2,145, гіпотеза H₀ :b₁ = 0 відхиляється на рівні значущості α = 0,05.

t₀ = |-8,978| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 9,07 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Волинська область

Рис. 2.1.3 Захворюваність на туберкульоз на 100000 населення в Волинській області, (всі форми)

Лінійна регресія 2,396х – 4735,2 значуща.

F = 129,62 > = 4,6, t₁ = 11,385 > = 2,145, гіпотеза H₀ :b₁ = 0 відхиляється на рівні значущості α = 0,05.

t₀ = |-11,2639| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 11,385 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Дніпропетровська область

Рис. 2.1.4 Захворюваність на туберкульоз на 100000 населення в Дніпропетровській області, (всі форми)

Лінійна регресія 4,3928х – 8720,9 значуща.

F = 183,16 > = 4,6, t₁ = 13,53 > = 2,145 – гіпотеза H₀ :b₁ = 0 відхиляється на рівні значущості α = 0,05.

t₀ = |-13,45| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 13,53 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Донецька область

Рис. 2.1.5 Захворюваність на туберкульоз на 100000 населення в Донецькій області, (всі форми)

Лінійна регресія 5,3024х – 10530 значуща.

F = 665,16 > = 4,6, t₁ = 25,791 > = 2,145 – гіпотеза H₀ :b₁ = 0 відхиляється на рівні значущості α = 0,05.

t₀ = |-25,64| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 25,791 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Житомирська область

Рис. 2.1.6 Захворюваність на туберкульоз на 100000 населення в Житомирській області, (всі форми)

Лінійна регресія 3,6297х – 7186,3 значуща.

F = 197,78 > = 4,6, t₁ = 14,063 > = 2,145 – гіпотеза H₀ :b₁ = 0 відхиляється на рівні значущості α = 0,05.

t₀ = |-13,94| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 14,063 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Закарпатська область

Рис. 2.1.7 Захворюваність на туберкульоз на 100000 населення в Закарпатській області, (всі форми)

Лінійна регресія 1,8896х – 3730,7 значуща.

F = 63,471 > = 4,6, t₁ = 7,967 > = 2,145 – гіпотеза H₀ :b₁ = 0 відхиляється на рівні значущості α = 0,05.

t₀ = |-7,875| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 7,967 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Запорізька область

Рис. 2.1.8 Захворюваність на туберкульоз на 100000 населення в Запорізькій області, (всі форми)

Лінійна регресія 4,6251х – 9178,1 значуща.

F = 173,619 > = 4,6, t₁ = 13,176 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-13,09| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 13,176 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Ів.-Франківська область

Рис. 2.1.9 Захворюваність на туберкульоз на 100000 населення в Ів.-Франківській області, (всі форми)

Лінійна регресія 2,9951х – 5930,4 значуща.

F = 222,779 > = 4,6, t₁ = 14,926 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-14,795| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 14,926 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Київська область

Рис. 2.1.10 Захворюваність на туберкульоз на 100000 населення в Київській області, (всі форми)

Лінійна регресія 2,4616х – 4867,2 значуща.

F = 255,154 > = 4,6, t₁ = 15,974 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-15,811| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 15,974 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Кіровоградська область

Рис. 2.1.11 Захворюваність на туберкульоз на 100000 населення в Кіровоградській області, (всі форми)

Лінійна регресія 4,6137х – 9156,1 значуща.

F = 120,162 > = 4,6, t₁ = 10,962 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-10,891| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 10,962 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Луганська область

Рис. 2.1.12 Захворюваність на туберкульоз на 100000 населення в Луганській області, (всі форми)

Лінійна регресія 5,9235х – 11767 значуща.

F = 261,454 > = 4,6, t₁ = 16,170 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-16,080| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 16,170 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Львівська область

Рис. 2.1.13 Захворюваність на туберкульоз на 100000 населення в Львівській області, (всі форми)

Лінійна регресія 3,2575х – 6452,1 значуща.

F = 582,815 > = 4,6, t₁ = 24,142 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-23,938| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 24,142 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Миколаївська область

Рис. 2.1.14 Захворюваність на туберкульоз на 100000 населення в Миколаївській області, (всі форми)

Лінійна регресія 6,0119х – 11943 значуща.

F = 114,02 > = 4,6, t₁ = 10,678 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-10,619| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 10,678 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Одеська область

Рис. 2.1.15 Захворюваність на туберкульоз на 100000 населення в Одеській області, (всі форми)

Лінійна регресія 3,8169х – 7563,1 значуща.

F = 180,966 > = 4,6, t₁ = 13,452 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-13,344| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 13,452 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Полтавська область

Рис. 2.1.16 Захворюваність на туберкульоз на 100000 населення в Полтавській області, (всі форми)

Лінійна регресія 2,0362х – 4026,1 значуща.

F = 144,948 > = 4,6, t₁ = 12,039 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-11,918| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 12,039 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Рівненська область

Рис. 2.1.17 Захворюваність на туберкульоз на 100000 населення в Рівненській області, (всі форми)

Лінійна регресія 2,7975х – 5531,4 значуща.

F = 246,950 > = 4,6, t₁ = 15,715 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-15,555| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 15,715 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Сумська область

Рис. 2.1.18 Захворюваність на туберкульоз на 100000 населення в Сумській області, (всі форми)

Лінійна регресія 1,7218х – 3392,3 значуща.

F = 86,422 > = 4,6, t₁ = 9,296 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-9,169| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 9,296 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Тернопільська область

Рис. 2.1.19 Захворюваність на туберкульоз на 100000 населення в Тернопільській області, (всі форми)

Лінійна регресія 2,1978х – 4342,8 значуща.

F = 92,134 > = 4,6, t₁ = 9,599 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-9,495| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 9,599 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Харківська область

Рис. 2.1.20 Захворюваність на туберкульоз на 100000 населення в Харківській області, (всі форми)

Лінійна регресія 4,7826х – 9496,6 значуща.

F = 90,884 > = 4,6, t₁ = 9,533 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-9,477| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 9,533 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Херсонська область

Рис. 2.1.21 Захворюваність на туберкульоз на 100000 населення в Херсонській області, (всі форми)

Лінійна регресія 9,5147х – 18920 значуща.

F = 82,923 > = 4,6, t₁ = 9,106 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-9,065| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 9,106 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Хмельницька область

Рис. 2.1.22 Захворюваність на туберкульоз на 100000 населення в Хмельницькій області, (всі форми)

Лінійна регресія 2,1478х – 4244,2 значуща.

F = 185,139 > = 4,6, t₁ = 13,607 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-13,461| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 13,607 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Черкаська область

Рис. 2.1.23 Захворюваність на туберкульоз на 100000 населення в Черкаській області, (всі форми)

Лінійна регресія 2,4782х – 4905,3 значуща.

F = 302,667 > = 4,6, t₁ = 17,397 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-17,239| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 17,397 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Чернівецька область

Рис. 2.1.24 Захворюваність на туберкульоз на 100000 населення в Чернівецькій області, (всі форми)

Лінійна регресія 2,0654х – 4080,4 значуща.

F = 96,164 > = 4,6, t₁ = 9,806 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-9,699| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 9,806 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Чернігівська область

Рис. 2.1.25 Захворюваність на туберкульоз на 100000 населення в Чернігівській області, (всі форми)

Лінійна регресія 3,1078х – 6151,6 значуща.

F = 102,457 > = 4,6, t₁ = 10,122 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-10,030| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 10,122 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

м. Київ

Рис. 2.1.26 Захворюваність на туберкульоз на 100000 населення в м. Київ, (всі форми)

Лінійна регресія 1,6359х – 3232,8 значуща.

F = 153,807 > = 4,6, t₁ = 12,402 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-12,269| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 12,402 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

м. Севастополь

Рис. 2.1.27 Захворюваність на туберкульоз на 100000 населення в м. Севастополь, (всі форми)

Лінійна регресія 3,0962х – 6141,8 значуща.

F = 452,541 > = 4,6, t₁ = 21,273 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-21,126| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 21,273 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Україна

Рис. 2.1.28 Захворюваність на туберкульоз на 100000 населення в Україні, (всі форми)

Лінійна регресія 3,7551х – 7446,5 значуща.

F = 457,416 > = 4,6, t₁ = 21,387 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості α = 0,05.

t₀ = |-21,232| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості α = 0,05.

t_R = 21,387 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості α = 0,05.

Таблиця 2.1.1 Захворюваність на туберкульоз, всі форми

Захворюваність на туберкульоз легенів

А.Р. Крим

Рис.2.1.29 Захворюваність на туберкульоз легень на 100000 населення в А.Р. Крим

Лінійна регресія 3,6518х – 7254,1 значуща.

F = 95,027 > = 4,6, t₁ = 9,748 > = 2,145, гіпотеза H₀ :b₁ = 0 відхиляється на рівні значущості 0,05.

t₀ = |-9,697| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05. t_R = 9,748 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Вінницька область

Рис. 2.1.30 Захворюваність на туберкульоз легень на 100000 населення в Вінницькій області

Лінійна регресія 2,5818х – 5115 значуща.

F = 70,976 > = 4,6, t₁ = 8,425 > = 2,145, гіпотеза H₀ :b₁ = 0 відхиляється на рівні значущості 0,05.

t₀ = |-8,358| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 8,425 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Волинська область

Рис. 2.1.31 Захворюваність на туберкульоз легень на 100000 населення в Волинській області

Лінійна регресія 2,0864х – 4127,4 значуща.

F = 162,232 > = 4,6, t₁ = 12,737 > = 2,145, гіпотеза H₀ :b₁ = 0 відхиляється на рівні значущості 0,05.

t₀ = |-12,617| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 12,737 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Дніпропетровська область

Рис. 2.1.32 Захворюваність на туберкульоз легень на 100000 населення в Дніпропетровській області

Лінійна регресія 4,1368х – 8216,5 значуща.

F = 127,344 > = 4,6, t₁ = 11,285 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-11,224| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 11,285 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Донецька область

Рис. 2.1.33 Захворюваність на туберкульоз легень на 100000 населення в Донецькій області

Лінійна регресія 4,9575х – 9849,5 значуща.

F = 416,123 > = 4,6, t₁ = 20,399 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-20,295| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 20,399 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Житомирська область

Рис. 2.1.34 Захворюваність на туберкульоз легень на 100000 населення в Житомирській області

Лінійна регресія 3,1461х – 6232,4 значуща.

F = 91,352 > = 4,6, t₁ = 9,558 > = 2,145 – гіпотеза H₀ :b₁ = 0 відхиляється на рівні значущості 0,05.

t₀ = |-9,481| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 9,558 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Закарпатська область

Рис. 2.1.35 Захворюваність на туберкульоз легень на 100000 населення в Закарпатській області

Лінійна регресія 1,7039х – 3363,5 значуща.

F = 39,867 > = 4,6, t₁ = 6,314 > = 2,145 – гіпотеза H₀ :b₁ = 0 відхиляється на рівні значущості 0,05.

t₀ = |-6,241| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 6,314 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Запорізька область

Рис. 2.1.36 Захворюваність натуберкульоз легень на 100000 населення в Запорізькій області

Лінійна регресія 4,4161х – 8768,4 значуща.

F = 119,702 > = 4,6, t₁ = 10,941 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-10,878| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 10,941 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Ів.-Франківська область

Рис. 2.1.37 Захворюваність на туберкульоз легень на 100000 населення в Ів.-Франківській області

Лінійна регресія 2,4707х – 4892 значуща.

F = 114,795 > = 4,6, t₁ = 10,714 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-10,623| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 10,714 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Київська область

Рис. 2.1.38 Захворюваність на туберкульоз легень на 100000 населення в Київській області

Лінійна регресія 2,4418х – 4833,2 значуща.

F = 275,108 > = 4,6, t₁ = 16,586 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-16,44| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 16,586 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Кіровоградській області

Рис. 2.1.39 Захворюваність на туберкульоз легень на 100000 населення в Кіровоградській області

Лінійна регресія 4,255х – 8449,4 значуща.

F = 128,668 > = 4,6, t₁ = 11,343 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-11,279| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 11,343 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Луганська область

Рис. 2.1.40 Захворюваність на туберкульоз легень на 100000 населення в Луганській області

Лінійна регресія 5,3836х – 10697 значуща.

F = 145,512 > = 4,6, t₁ = 12,063 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-12,003| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 12,063 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Львівська область

Рис. 2.1.41 Захворюваність на туберкульоз легень на 100000 населення в Львівській області

Лінійна регресія 2,9554х – 5858,2 значуща.

F = 393,732 > = 4,6, t₁ = 19,843 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-19,696| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 19,843 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Миколаївська область

Рис. 2.1.42 Захворюваність на туберкульоз легень на 100000 населення в Миколаївській області

Лінійна регресія 5,52х – 10969 значуща.

F = 86,279 > = 4,6, t₁ = 9,289 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-9,243| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 9,289 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Одеська область

Рис. 2.1.43 Захворюваність на туберкульоз легень на 100000 населення в Одеській області

Лінійна регресія 3,4589х – 6856,7 значуща.

F = 113,536 > = 4,6, t₁ = 10,655 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-10,577| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 10,655 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Полтавська область

Рис. 2.1.44 Захворюваність на туберкульоз легень на 100000 населення в Полтавській області

Лінійна регресія 2,0893х – 4136,9 значуща.

F = 141,154 > = 4,6, t₁ = 11,881 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-11,780| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 11,881 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Рівненська область

Рис. 2.1.45 Захворюваність на туберкульоз легень на 100000 населення в Рівненській області

Лінійна регресія 2,4664х – 4882,3 значуща.

F = 81,36 > = 4,6, t₁ = 9,02 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-8,941| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 9,02 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Сумська область

Рис. 2.1.46 Захворюваність на туберкульоз легень на 100000 населення в Сумській області

Лінійна регресія 1,63х – 3220,3 значуща.

F = 98,003 > = 4,6, t₁ = 9,9 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-9,794| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 9,9 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Тернопільська область

Рис. 2.1.47 Захворюваність на туберкульоз легень на 100000 населення в Тернопільській області

Лінійна регресія 1,8989х – 3756,4 значуща.

F = 92,908 > = 4,6, t₁ = 9,639 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-9,548| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 9,639 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Харківська область

Рис. 2.1.48 Захворюваність на туберкульоз легень на 100000 населення в Харківській області

Лінійна регресія 4,6232х – 9185,8 значуща.

F = 70,08 > = 4,6, t₁ = 8,371 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-8,329| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 8,371 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Херсонська область

Рис. 2.1.49 Захворюваність на туберкульоз легень на 100000 населення в Херсонській області

Лінійна регресія 8,4086х – 16724 значуща.

F = 46,062 > = 4,6, t₁ = 6,787 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-6,76| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 6,787 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Хмельницька область

Рис. 2.1.50 Захворюваність на туберкульоз легень на 100000 населення в Хмельницькій області

Лінійна регресія 1,9579х – 3871,8 значуща.

F = 160,873 > = 4,6, t₁ = 12,684 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-12,560| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 12,684 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Черкаська область

Рис. 2.1.51 Захворюваність на туберкульоз легень на 100000 населення в Черкаській області

Лінійна регресія 1,9132х – 3787,2 значуща.

F = 273,318 > = 4,6, t₁ = 16,532 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-16,387| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 16,532 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Чернівецька область

Рис. 2.1.52 Захворюваність на туберкульоз легень на 100000 населення в Чернівецькій області

Лінійна регресія 2,1461х – 4251 значуща.

F = 102,363 > = 4,6, t₁ = 10,117 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-10,035| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 10,117 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Чернігівській області

Рис. 2.1.53 Захворюваність на туберкульоз легень на 100000 населення в Чернігівській області

Лінійна регресія 2,9232х – 5792,9 значуща.

F = 143,748 > = 4,6, t₁ = 11,989 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-11,898| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 11,989 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

м. Київ

Рис. 2.1.54 Захворюваність на туберкульоз легень на 100000 населення в м. Київ

Лінійна регресія 1,8314х – 3628,4 значуща.

F = 248,39 > = 4,6, t₁ = 15,76 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-15,636| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 15,76 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

м. Севастополь

Рис. 2.1.55 Захворюваність на туберкульоз легень на 100000 населення в м. Севастополь

Лінійна регресія 2,9225х – 5802,9 значуща.

F = 508,098 > = 4,6, t₁ = 22,541 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-22,412| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 22,541 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Україна

Рис. 2.1.56 Захворюваність на туберкульоз легень на 100000 населення в Україні

Лінійна регресія 3,4857х – 6916,8 значуща.

F = 251,374 > = 4,6, t₁ = 15,855 > = 2,145 – гіпотеза H₀ :b₁ =0 відхиляється на рівні значущості 0,05.

t₀ = |-15,754| > = 2,145 – гіпотеза H₀ :b₀ = 0 відхиляється на рівні значущості 0,05.

t_R = 15,855 > = 2,145 – гіпотеза H₀ :R = 0 відхиляється на рівні значущості 0,05.

Таблиця 2.1.2 Захворюваність на туберкульоз легенів