86309 (Логарифмічно-лінійний аналіз)

КОНТРОЛЬНА РОБОТА

ЛОГАРИФМІЧНО ЛІНІЙНИЙ АНАЛІЗ

Зміст

1. Основні поняття логлінійного аналізу

2. Аналіз зв’язку категоризованих змінних

3. Канонічна кореляція при аналізі таблиць спряженості ознак

4. Побудова логарифмічної моделі

Висновки

Список використаних джерел

2. Основні поняття логлінійного аналіза

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

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

Одним з основних питань, яке виникає при обробці даних, є питання про наявність взаємозв'язку між ознаками. Методом багатовимірного статистичного аналізу, який використовується для обробки багатовимірних таблиць спряженості, є логлінійний аналіз.

Основна ідея логлінійного аналізу полягає в тому, що в кожній комірці таблиці спряженості записується розклад натурального логарифма частоти на суму ефектів всіх взаємозв'язків досліджуваних ознак. Модель к-го порядку називається насиченою, якщо вона крім головних ефектів містить будь-які взаємозв'язки величин від другого до к-го порядків включно. Так, наприклад, для 3-х факторів A, B, C насичена модель має вигляд:

(1.1)

де μ - загальне середнє, - головні ефекти досліджуваних ознак. Так, наприклад, можна інтерпретувати як додаток (спад, якщо її знак від’ємний) i-ой категорії фактора A в порівнянні із загальним середнім, що є мірою того, наскільки ймовірна i-а категорія ознаки А в порівнянні з іншими; - ефекти відповідних взаємодій. Насичена модель містить стільки ж параметрів, скільки і комірок в таблиці спряженості. Для того, щоб модель була не надлишковою, тобто, щоб число параметрів не перевищило число елементів ( ) досліджуваної таблиці спряженості, на значення в моделі накладаються наступні обмеження:

(1.2)

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

1) - такі моделі називаються моделями незалежності. Вони містять лише головні ефекти, відсутні ефекти взаємодії чинників.

2) - категорії факторів В і С – рівноймовірні.

3) - категорії факторів А і С – рівноймовірні.

4) - категорії факторів В і А – рівноймовірні.

5) - всі категорії всіх факторів – рівноймовірні.

6) - має місце взаємозв'язок двох ознак А і В. Така модель не відноситься до ієрархічних.

7) - має місце взаємозв'язок двох ознак В і С.

2. Аналіз зв’язку категоризованих змінних

Дуже часто вивчення взаємозв'язків між якісними змінними зводитися до аналізу двох змінних х і у, які набувають низку значень – х₁, …, х_m і y₁,…, y_p. Позначимо номер рядка двовимірної таблиці через і (і = 1,…, m), номер стовпця – через j (j = 1,…,р). Кожен об'єкт характеризується значеннями двох змінних х_i і y_j. У комірках таблиці записується число об'єктів у вибірці, що володіють даним поєднанням значень ознак (коміркова частота) – n_ij.

Як правило, аналіз двовимірних таблиць спряженості обмежується висновками про наявність зв'язку і її тісноти. Класичним тестом, який використовується для встановлення факту наявності зв'язку, є критерій χ². Величина χ² звичайно визначається як сума квадратів між фактичними (n_ij) і теоретичними ( ) частотами двовимірного розподілу, поділена на теоретичні частоти.

(2.1)

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

(2.2)

де - сума за і-тим рядком таблиці, - сума за j-м стовпцем таблиці (так звані маргінальні частоти). Обираючи рівень значущості λ, визначають відповідне критичне значення критерію χ²_λ;df при числі ступенів вільності df = (m - 1)(p - 1). Якщо розрахункова величина χ² перевищує критичне значення χ²_λ;df, то на даному рівні значущості нульова гіпотеза Н₀ може бути відхилена.

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

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

При цьому не виявляється дуже важлива інформація про зміну інтенсивності зв'язку змінних із зміною їх значень: при одних значеннях зв'язок ознак може бути щільним, при інших – слабким і навіть взагалі не спостерігатися. Узагальнюючи різний ступінь інтенсивності зв'язку в єдиному показнику, ми можемо дійти парадоксальних висновків, загубити практично значущі випадки. Щоб уникнути цього, корисно аналізувати окремі фрагменти початкової таблиці. Перехід від таблиці в цілому до окремих її частин доцільний як у разі багатовимірного аналізу, так і за наявності великого числа категорій двовимірного розподілу. На виділення частинних зв'язків заснована ціла низка статистичних методів, наприклад дисперсійний аналіз, логлінійний аналіз.

Використання розкладання статистики χ² припускає перетворення вихідної таблиці спряженості в безліч таблиць розмірністю 2×2, кожна з яких характеризує особливий аспект зв'язку, що вивчається, – зв'язок між певними значеннями змінних. Число таких таблиць повинне дорівнювати числу ступенів вільності критерію χ², при його обчисленні за вихідною таблицею. Кожний з локальних зв'язків у виділеному фрагменті загальної таблиці оцінюється за допомогою критерію χ², що має одну степінь вільності. Зв'язок загальної величини статистики χ² із значеннями, знайденими по виділених частинах таблиці, виражається рівністю:

(2.3)

де χ_k² – k –та компонента загальної величини χ², знайдена за загальною таблицею 2х2.

Рівність (2.3) справедлива при означенні оцінки χ² методом максимальної правдоподібності (χ_ML²), яка має вигляд:

(2.4)

Зазвичай вираз (2.4), знайдений відповідно до загального методу найменших квадратів, не має властивості адитивності: в цьому випадку сума окремих χ_k² не обов'язково буде точно дорівнювати загальному χ² для всієї таблиці. Відмітимо, що властивість адитивності χ_ML² використовується в одному з сучасних методів багатовимірного аналізу якісних змінних – логлінійному аналізі.

3. Канонічна кореляція в аналізі таблиць спряженості

Один з напрямів аналізу таблиці спряженості пов'язаний з «оцифруванням» якісних ознак – з приписуванням градаціям якісних змінних числових міток. Такий підхід дозволяє розповсюдити на якісні дані методи багатовимірного статистичного аналізу, розроблені відносно кількісних змінних.

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

Дослідити це питання дозволяє критерій χ² через адитивні компоненти, які відповідають лінійним і нелінійним ефектам в структурі зв'язку між змінними. При цьому лінійні ефекти пов’язують з першим перетворенням змінних, нелінійні – з другим і т.д. перетвореннями. Адитивність ефектів випливає з ортогональності канонічних змінних.

Покажемо, що канонічний аналіз таблиці зв'язаності відповідає розкладанню статистики χ² на ряд доданків, число яких залежить від розмірності таблиці. Перетворимо вираз таким чином:

(3.1)

(3.2)

Остання формула може бути подана за допомогою суми діагональних елементів, тобто слід симетричної матриці , де N - матриця розмірності (m × p) з елементами (m – число рядків таблиці, p – число стовпців):

(3.3)

Якщо число рядків таблиці не дорівнює числу стовбців, то, як правило, матрицю С формують так, щоб її розмірність була мінімальною (min(m,p)). Оскільки слід матриці дорівнює сумі її власних чисел, то вираз (3.3) приймає наступний вигляд:

(3.4)

де λ_k - k-е власне число матриці С.

Враховуючи, що власні числа є показниками кореляції (R²) між кожною парою канонічних змінних, виділених з вихідних наборів даних, запишемо рівність (3.4) у вигляді:

(3.5)

З m (або p) власних чисел матриці С максимальне завжди дорівнює одиниці, йому відповідає вироджений набір міток 1 = (1,...,1). Тому вираз (3.5) доцільно переписати так :

(3.6)

Найбільше з чисел, що залишилися (m - 1) або (p - 1) власних чисел відповідає гіпотезі лінійності зв'язку між категоризованими змінними; наступне за величиною значення λ_k відповідає гіпотезі про складніший характер взаємозв’язку змінних. Така інтерпретація компонент χ² представляється можливою з причини того, що кожна подальша пара канонічних змінних є функцією першої перетвореної пари, а все розкладання χ² є спадаючою послідовністю.