Bagaimana cara melakukan beberapa tes post-hoc chi-square pada tabel 2 X 3?

9

Kumpulan data saya terdiri dari angka kematian total atau kelangsungan hidup suatu organisme pada tiga tipe lokasi, inshore, midchannel, dan offshore. Angka-angka dalam tabel di bawah ini mewakili jumlah situs.

              100% Mortality            100% Survival
Inshore             30                       31 
Midchannel          10                       20 
Offshore             1                       10

Saya ingin tahu apakah # situs tempat kematian 100% terjadi signifikan berdasarkan jenis situs. Jika saya menjalankan 2 x 3 chi-square, saya mendapatkan hasil yang signifikan. Apakah ada perbandingan berpasangan post-hoc yang dapat saya jalankan atau haruskah saya benar-benar menggunakan ANOVA logistik atau regresi dengan distribusi binomial? Terima kasih!

chl
sumber

Jawaban:

7

Tabel kontingensi harus berisi semua kategori yang saling eksklusif pada kedua sumbu. Dalam negeri / Midchannel / Offshore terlihat baik-baik saja, namun kecuali "mortalitas kurang dari 100%" berarti "100% bertahan hidup" dalam pengaturan biologis ini, Anda mungkin perlu membuat tabel yang menjelaskan semua kasus yang diamati atau menjelaskan mengapa Anda membatasi analisis Anda secara ekstrem. ujung sampel.

Karena 100% angka bertahan berarti 0% kematian, Anda dapat memiliki tabel dengan kolom 100% = kematian / 100%> kematian> 0% / kematian = 0%. Dalam hal ini Anda tidak akan lagi membandingkan persentase, tetapi membandingkan ukuran kematian ordinal di tiga kategori tipe situs. (Bagaimana dengan menggunakan nilai persentase asli alih-alih kategori?) Versi uji Kruskal-Wallis mungkin sesuai di sini yang mempertimbangkan kaitan secara tepat (mungkin tes permutasi).

Ada tes post hoc yang ditetapkan untuk tes Kruskal-Wallis: 1 , 2, 3 . (Pendekatan resampling dapat membantu mengatasi ikatan.)

Regresi logistik dan regresi binomial mungkin lebih baik karena mereka tidak hanya memberi Anda nilai p, tetapi juga estimasi yang berguna dan interval kepercayaan dari ukuran efek. Namun untuk mengatur model-model itu dibutuhkan lebih banyak detail mengenai 100%> mortalitas> 0% lokasi.

GaBorgulya
sumber
4

Saya akan menganggap bahwa "kelangsungan hidup 100%" berarti bahwa situs Anda hanya berisi satu organisme. jadi 30 berarti 30 organisme mati, dan 31 berarti 31 organisme tidak. Berdasarkan hal ini chi-square harus baik-baik saja, tetapi hanya akan memberi tahu hipotesis mana yang tidak didukung oleh data - itu tidak akan memberi tahu Anda apakah dua hipotesis masuk akal lebih baik atau tidak. Saya menyajikan analisis probabilitas yang mengekstraksi informasi ini - itu setuju dengan uji chi-square, tetapi memberi Anda lebih banyak informasi daripada uji chi-square, dan cara yang lebih baik untuk menyajikan hasilnya.

Model adalah model bernouli untuk indikator "kematian", ( menunjukkan sel dari tabel , dan menunjukkan unit individu dalam sel).i 2 × 3 jYsayajBsayan(1,θsayaj)saya2×3j

Ada dua asumsi global yang mendasari uji chi-square:

  1. dalam sel tabel yang diberikan, semuanya sama, yaitu θ i j = θ i k = θ iθsayajθsayaj=θsayak=θsaya
  2. yang secara statistik independen, mengingat . Ini berarti bahwa parameter probabilitas memberi tahu Anda segala sesuatu tentang - semua informasi lain tidak relevan jika Anda tahu θ i Y i j θ iYsayajθsayaYsayajθsaya

Nyatakan sebagai jumlah , (jadi ) dan biarkan menjadi ukuran grup (jadi ). Sekarang kami memiliki hipotesis untuk diuji: Y i j X 1 = 30 , X 2 = 10 , X 3 = 1 N i N 1 = 61 , N 2 = 30 , N 3 = 11XsayaYsayajX1=30,X2=10,X3=1NsayaN1=61,N2=30,N3=11

HSEBUAH:θ1=θ2,θ1=θ3,θ2=θ3

Tapi apa alternatifnya? Saya akan mengatakan kombinasi lain yang mungkin sama atau tidak sama.

H B 2 : θ 1θ 2 , θ 1 = θ 3 , θ 2θ 3 H B 3 : θ 1 = θ 2 , θ 1θ 3 , θ 2

HB1:θ1θ2,θ1θ3,θ2=θ3
HB2:θ1θ2,θ1=θ3,θ2θ3
HB3:θ1=θ2,θ1θ3,θ2θ3
HC:θ1θ2,θ1θ3,θ2θ3

Salah satu hipotesis ini harus benar, mengingat asumsi "global" di atas. Tetapi perhatikan bahwa tidak satupun dari ini yang menentukan nilai spesifik untuk tarif - sehingga mereka harus diintegrasikan. Sekarang mengingat benar, kami hanya memiliki satu parameter (karena semuanya sama), dan seragam sebelumnya adalah pilihan konservatif, menyatakan ini dan asumsi global oleh . jadi kita punya:HSEBUAHsaya0

P(X1,X2,X3|N1,N2,N3,HSEBUAH,saya0)=01P(X1,X2,X3,θ|N1,N2,N3,HSEBUAH,saya0)dθ
=(N1X1)(N2X2)(N3X3)01θX1+X2+X3(1-θ)N1+N2+N3-X1-X2-X3dθ
=(N1X1)(N2X2)(N3X3)(N1+N2+N3+1)(N1+N2+N3X1+X2+X3)

Yang merupakan distribusi hypergeometrik dibagi dengan konstanta. Demikian pula untuk kita akan memiliki: HB1

P(X1,X2,X3|N1,N2,N3,HB1,saya0)=01P(X1,X2,X3,θ1θ2|N1,N2,N3,HB1,saya0)dθ1dθ2
=(N2X2)(N3X3)(N1+1)(N2+N3+1)(N2+N3X2+X3)

Anda dapat melihat pola untuk yang lain. Kita dapat menghitung peluang untuk mengatakan dengan hanya membagi dua ekspresi di atas. Jawabannya adalah sekitar , yang berarti data mendukung lebih dari sekitar faktor - bukti yang cukup lemah dalam mendukung tingkat yang sama. Probabilitas lain diberikan di bawah ini.HSEBUAHvsHB14HSEBUAHHB14

HyhalHaithessayashalrHaibSebuahbsayalsayaty(HSEBUAH|D)0,018982265(HB1|D)0,004790669(HB2|D)0,051620022(HB3|D)0.484155874(HC|D)0.440451171

Ini menunjukkan bukti kuat terhadap tingkat yang sama, tetapi tidak dalam bukti kuat mendukung alternatif defintie. Sepertinya ada bukti kuat bahwa laju "lepas pantai" berbeda dengan dua tarif lainnya, tetapi bukti yang tidak meyakinkan tentang apakah tarif "dalam negeri" dan "saluran tengah" berbeda. Inilah yang tes chi-square tidak akan memberi tahu Anda - itu hanya memberi tahu Anda bahwa hipotesis adalah "omong kosong" tetapi bukan alternatif apa yang harus ditempatkan pada tempatnya.SEBUAH

probabilityislogic
sumber
1

Berikut ini adalah kode untuk melakukan tes chi square serta menghasilkan berbagai statistik uji. Namun, uji statistik dari asosiasi margin tabel tidak berguna di sini; jawabannya jelas. Tidak ada yang melakukan tes statistik untuk melihat apakah musim panas lebih panas daripada musim dingin.

Chompy<-matrix(c(30,10,1,31,20,10), 3, 2)
Chompy
chisq.test(Chompy)
chisq.test(Chompy, simulate.p.value = TRUE, B = 10000)
chompy2<-data.frame(matrix(c(30,10,1,31,20,10,1,2,1,2,1,2,1,2,3,1,2,3), 6,3))
chompy2
chompy2$X2<-factor(chompy2$X2) 
chompy2$X3<-factor(chompy2$X3)
summary(fit1<-glm(X1~X2+X3, data=chompy2, family=poisson))
summary(fit2<-glm(X1~X2*X3, data=chompy2, family=poisson)) #oversaturated
summary(fit3<-glm(X1~1, data=chompy2, family=poisson)) #null
anova(fit3,fit1)
library(lmtest)
waldtest(fit1)
waldtest(fit2) #oversaturated
kruskal.test(X1~X2+X3, data=chompy2)
kruskal.test(X1~X2*X3, data=chompy2)
Patrick McCann
sumber
3
Akan menarik bagi pembaca (dan OP) jika Anda dapat memberikan rincian tentang sintaks R yang berbeda (dan tes yang mendasarinya) yang Anda berikan, dan terutama bagaimana tes Kruskal-Wallis dibandingkan dengan model log-linear.
chl
Anda dapat melihat ini dengan menyalin dan menempelkan kode ke konsol R.
Patrick McCann
1
Tentu. Tanggapan datang dari diri mereka sendiri dengan menjalankan kode, tentu saja.
chl
0

Saya yakin Anda bisa menggunakan "interval kepercayaan simultan" untuk melakukan beberapa perbandingan. Rujukannya adalah Agresti et al. 2008 Interval kepercayaan simultan untuk membandingkan parameter binomial. Biometrik 64 1270-1275.

Anda dapat menemukan kode R yang sesuai di http://www.stat.ufl.edu/~aa/cda/software.html

Tu.2
sumber