Fövqəladə vəziyyət analizi: Fəsil 9 nümunələri üçün R kodu

Göndərildi: 08.09.2021
Məqalənin müəllifi Adəm Quliyev

Mənasını qısa təsvir etmək üçün funksiya arqumentinin üzərinə sürün. Dəyişən adlar aşağıda göstərilən nümunələrdən götürülür.

Epitools paketindən istifadə edərək əmsal nisbəti:

\ (\ chi ^ 2 \) fövqəladə vəziyyət testi:

Fisherin dəqiq sınağı:

Digər yeni metodlar:

  • Nisbi risk.
  • G-test.
  • Bir vəziyyətə nəzarət işi üçün mozaika sahəsi.

Şəkil 9.1-1. Titanikdə seks və yaşamaq

Gəmi batan Titanikdə cinsi əlaqə ilə sağ qalma arasındakı əlaqəni göstərən mozaika planı.

Məlumatları oxuyun və yoxlayın.

Cinsiyyət və sağ qalma arasındakı əlaqənin fövqəladə cədvəli. Addmargins funksiyası fövqəladə cədvəl ilə satır və sütun cəmlərini əhatə edir.

Dərnəyin mozaika sahəsi.

Nümunə 9.2. Aspirin və xərçəng

2x2 fövqəladə cədvəldə aspirin müalicəsi ilə xərçəng insidansı arasındakı əlaqəni qiymətləndirmək üçün nisbət nisbətinisbi risk.

Məlumatları oxuyun və yoxlayın.

Aspirin müalicəsi ilə xərçəng insidansı arasında əlaqəni göstərən 2 x 2 fövqəladə cədvəl(Cədvəl 9.2-1).

Birliyi göstərən mozaika sahəsi(Şəkil 9.2-1).

Epitools paketini istifadə edərək əmsal nisbətini hesablayın. Əgər hələ etməmisinizsə, paketi quraşdırın (bunun R-yə qoşulduqdan sonra hər kompüter üçün bir dəfə edilməsi lazımdır). İstifadəsi üçün epitools paketini yükləyin (bunun hər R seansında yalnız bir dəfə edilməsi lazımdır).

Kitabdakı ehtimal nisbəti və nisbi risk düsturları ehtimal olunan cədvəlin xərçəng cədvəlinə bənzər bir tərtibə sahib olduğunu, yəni

Bu tərtib oddsratio funksiyası üçün işləyir.

Nəticədə oddsratio çıxışı bir çox təsadüfi hesablamaları əhatə edir. Yalnız 95% güvən intervalı daxil olmaqla nisbət nisbəti üçün təmiz nəticə əldə etmək üçün bunun əvəzinə aşağıdakılardan istifadə edin.

Oran nisbəti və güvən intervalı üzrə addım-addım hesablamalar üçün aşağıdakıları edirik:

Epitools paketini istifadə edərək nisbi riski hesablayın. Riskratio funksiyası ilə gözlənilən plan kitabda istifadə edilmiş düzənin tam əksidir. Əmri kitab tərzində (məsələn, xəritələr cədvəli) bir fövqəladə cədvəllə istifadə etmək üçün cədvəli çevirməli (yer dəyişdirməli) və sıraların sırasını tərs çevirməliyik. Riskratio funksiyasına dair aşağıdakı dəlillərlə bunu bir anda edə bilərik.

Diqqət yetirin ki, nəticənin eyni məlumatlar üçün kitabda verilən nisbi risk dəyərindən bir qədər fərqlənir (1.007), çünki burada yuvarlaqlaşdırma xətası azalır.

Yalnız 95% güvən intervalı daxil olmaqla nisbi risk üçün təmiz nəticə üçün aşağıdakı əmri istifadə edin:

Nisbi riskin addım-addım hesablanması üçün aşağıdakıları edirik.

Nümunə 9.3. Toksoplazma və sürücülük qəzaları

2 x 2 vəziyyətə nəzarət işindəToxoplasma infeksiyası ilə sürücülük tezliyi arasındakı əlaqəni təxmin etmək nisbət nisbəti.

Məlumatları oxuyun və yoxlayın.

2 x 2 fövqəladə vəziyyət cədvəli(Cədvəl 9.3-1).

Bir vəziyyətə nəzarət işi üçün mozaika sahəsi(Şəkil 9.3-1).

Epitools istifadə edərək əmsal nisbəti.

Nümunə 9.4. Qurd quş alır

\ (\ chi ^ 2 \)yırtıcı quşların iştirakı ilə killifişin trematod infeksiyası vəziyyəti ilə onların aqibəti (yeyib-yemədiyi) arasındakı əlaqəni testetmək üçün fövqəladə vəziyyət testi.

Məlumatları oxuyun və yoxlayın.

Cədvəllərdə və qrafiklərdə infeksiya kateqoriyalarının üstünlük verilən qaydasını təyin edin.

Fövqəladə vəziyyət cədvəli(Cədvəl 9.4-1). Addmargins funksiyası cədvəlin ekranına sətir və sütun cəmlərini əlavə edir.

Mozaika sahəsi(Şəkil 9.4-1).

\ (\ chi ^ 2 \) fövqəladə vəziyyət testi. Yatesin düzəldilməsinin qarşısını almaq üçün arqumenti düzgün = YALAN daxil edirik. Bunun 2 x 2 cədvəl xaricində heç bir təsiri yoxdur, ancaq nümayiş məqsədi ilə burada saxlayırıq.

Null hipotez altında gözlənilir tezliklərinnəticələri ilə daxil edilir, lakin adətən nümayiş deyil. Bunları çıxartması gözlənilən saveChiTest $ yazın. Sətir və sütun cəmlərini də görmək istəyirsinizsə addmargins funksiyasını daxil edin.

\ (G \) -testqurdu-quş məlumatlarına tətbiq edilir (Bölmə 9.6). R-in uyğunluq məlumatları ilə \ (G \) testini həyata keçirmək üçün sadə, daxili bir funksiyası yoxdur. Brent Larget tərəfindən hazırlanan bir komanda üçün g.test kodu burada mövcuddur. Aşağıda, bu kodu g.test funksiyasından istifadə etməyə imkan verən mənbəyi təqdim edirik.

Nümunə 9.5. Vampir yarasa yemi

Fişerininəklərin estrus statusu ilə inəklərin vampir yarasaları tərəfindən dişlənib-vurulmadığı arasındakı dəqiq əlaqə testi.

Məlumatları oxuyun və yoxlayın.

Fövqəladə vəziyyət cədvəli(Cədvəl 9.5-1).

Sıfır müstəqillik fərziyyəsi altında gözlənilən tezliklər. R, fərziyyələrin pozulduğundan şikayətlənir. Yalnız fərq etməmisiniz.

Fişerin dəqiq sınağı. Çıxışda bahis nisbətinin qiymətləndirilməsi daxildir, lakin kitabımızda istifadə olunandan fərqli bir metoddan istifadə olunur.