Dalam menguji asumsi regresi paralel dalam regresi logistik ordinal saya menemukan ada beberapa pendekatan. Saya telah menggunakan pendekatan grafis (seperti yang dijelaskan dalam buku Harrell) dan pendekatan yang dirinci menggunakan paket ordinal dalam R.
Namun saya juga ingin menjalankan tes Brant (dari Stata) untuk variabel individu dan juga untuk model total. Saya telah melihat sekeliling tetapi tidak dapat menemukannya diimplementasikan dalam R.
Apakah ada implementasi tes Brant di R?
Jawaban:
Saya menerapkan tes brant di R. Paket dan fungsinya disebut brant dan sekarang tersedia di CRAN.
Tes brant didefinisikan oleh Rollin Brant untuk menguji asumsi regresi paralel (Brant, R. (1990) Menilai proporsionalitas dalam model odds proporsional untuk regresi logistik ordinal. Biometrik , 46 , 1171-1178).
Berikut ini contoh kode:
Dalam contoh, asumsi regresi paralel berlaku, karena semua nilai p di atas 0,05. Omnibus adalah untuk seluruh model, sisanya untuk koefisien indvidual.
sumber
Ya - sebenarnya paket ordinal yang Anda tautkan dapat melakukannya (meskipun mereka tidak menyebutnya tes Brant). Lihatlah halaman 6 dan 7 dari tautan Anda, yang menunjukkan "uji rasio kemungkinan dari kemiringan yang sama atau asumsi peluang proporsional," yang persis seperti yang Anda cari.
sumber
Beberapa catatan tentang topik tersebut
Paket R
VGAM
dalamCumulative
perintah (Regresi Ordinal dengan Peluang Kumulatif) memungkinkan untuk mengubah asumsi peluang proporsional, dengan opsiparallel=FALSE
.Ini dikenal sebagai masalah umum (dari buku: Model Regresi untuk Variabel Ketergantungan Kategori Menggunakan Stata, Edisi Kedua, Oleh J. Scott Long, Jeremy Freese)
Makalah ini masuk jauh dalam topik ini, yang jelas dan ditulis dengan baik, tetapi tidak mempertimbangkan paket VGAM atau perintah "kumulatif": Regresi logistik ordinal dalam studi epidemiologi
sumber
Ini tutorial tentang regresi logistik ordinal di R meliputi pengujian asumsi kemungkinan proporsional.
sumber