Berikut adalah sedikit latar belakang tentang situasi saya: data saya merujuk pada jumlah mangsa yang berhasil dimakan oleh pemangsa. Karena jumlah mangsa terbatas (25 tersedia) di setiap percobaan, saya memiliki kolom "Sampel" yang mewakili jumlah mangsa yang tersedia (jadi, 25 di setiap percobaan), dan yang lain disebut "Hitungan" yang merupakan jumlah keberhasilan ( berapa banyak mangsa yang dimakan). Saya mendasarkan analisis saya pada contoh dari buku R pada data proporsi (halaman 578). Variabel penjelas adalah Temperatur (4 level, yang saya perlakukan sebagai faktor), dan Jenis Kelamin pemangsa (jelas, pria atau wanita). Jadi saya berakhir dengan model ini:
model <- glm(y ~ Temperature+Sex+Temperature*Sex data=predator, family=quasibinomial)
Setelah mendapatkan tabel Analysis of Deviance, ternyata Suhu dan Jenis Kelamin (tetapi bukan interaksi) memiliki pengaruh yang signifikan terhadap konsumsi mangsa. Sekarang, masalah saya: Saya perlu tahu suhu mana yang berbeda, yaitu, saya harus membandingkan 4 suhu satu sama lain. Jika saya memiliki model linier, saya akan menggunakan fungsi TukeyHSD, tetapi karena saya menggunakan GLM saya tidak bisa. Saya telah mencari melalui paket MASS dan mencoba untuk mengatur matriks kontras tetapi untuk beberapa alasan itu tidak berhasil. Ada saran atau referensi?
Inilah ringkasan yang saya dapatkan dari model saya, jika itu membantu membuatnya lebih jelas ...
y <- cbind(data$Count, data$Sample-data$Count)
model <- glm(y ~ Temperature+Sex+Temperature*Sex data=predator, family=quasibinomial)
> summary(model)
# Call:
# glm(formula = y ~ Temperature + Sex + Temperature * Sex, family=quasibinomial, data=data)
# Deviance Residuals:
# Min 1Q Median 3Q Max
# -3.7926 -1.4308 -0.3098 0.9438 3.6831
# Coefficients:
# Estimate Std. Error t value Pr(>|t|)
# (Intercept) -1.6094 0.2672 -6.024 3.86e-08 ***
# Temperature8 0.3438 0.3594 0.957 0.3414
# Temperature11 -1.0296 0.4803 -2.144 0.0348 *
# Temperature15 -1.2669 0.5174 -2.449 0.0163 *
# SexMale 0.3822 0.3577 1.069 0.2882
# Temperature8:SexMale -0.2152 0.4884 -0.441 0.6606
# Temperature11:SexMale 0.4136 0.6093 0.679 0.4990
# Temperature15:SexMale 0.4370 0.6503 0.672 0.5033
# ---
# Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
# (Dispersion parameter for quasibinomial family taken to be 2.97372)
# Null deviance: 384.54 on 95 degrees of freedom
# Residual deviance: 289.45 on 88 degrees of freedom
# AIC: NA
# Number of Fisher Scoring iterations: 5
glht
fungsi dalammultcomp
paket . Untuk melakukan tes TukeyHSD untuk suhu, gunakan seperti ituglht(my.glm, mcp(Temperature="Tukey"))
. Dan btw: Rumus Model Anda dapat disingkat:model<-glm(y ~ Temperature*Sex data=predator, family=quasibinomial)
. Dengan tanda bintang ( ) interaksi dan efek utama dipasang.Temperature
sebagai faktor? Apakah Anda tidak memiliki nilai numerik yang sebenarnya? Saya akan menggunakannya sebagai variabel kontinu & kemudian seluruh masalah ini diperdebatkan.Jawaban:
Anne, saya akan menjelaskan bagaimana melakukan beberapa perbandingan secara umum. Mengapa ini tidak berhasil dalam kasus spesifik Anda, saya tidak tahu; Maafkan saya.
Tetapi biasanya, Anda bisa melakukannya dengan
multcomp
paket dan fungsinyaglht
. Berikut ini sebuah contoh:Jika Anda ingin menghitung perbandingan berpasangan antara
rank
menggunakan Tukey's HSD, Anda bisa melakukannya dengan cara ini:Catatan: Seperti yang dicatat @gung dalam komentar, Anda harus - bila memungkinkan - memasukkan suhu sebagai variabel kontinu daripada variabel kategorikal. Mengenai interaksi: Anda dapat melakukan tes rasio kemungkinan untuk memeriksa apakah istilah interaksi secara signifikan meningkatkan kesesuaian model. Dalam kasus Anda, kode akan terlihat seperti itu:
Jika tes ini tidak signifikan, Anda dapat menghapus interaksi dari model Anda. Mungkin
glht
akan berhasil?sumber