Saya meninjau makalah yang memiliki eksperimen biologis berikut. Alat digunakan untuk mengekspos sel pada berbagai jumlah tegangan geser cairan. Saat tegangan geser yang lebih besar diterapkan ke sel, lebih banyak sel mulai terlepas dari substrat. Pada setiap tingkat tegangan geser, mereka menghitung sel-sel yang tetap melekat, dan karena mereka tahu jumlah total sel yang terpasang di awal, mereka dapat menghitung perlekatan fraksional (atau detasemen).
Jika Anda memplot fraksi patuh vs tegangan geser, hasilnya adalah kurva logistik. Secara teori, setiap sel individu adalah pengamatan tunggal, tetapi jelas ada ribuan atau puluhan ribu sel, sehingga kumpulan data akan menjadi raksasa, jika itu diatur dengan cara biasa (dengan setiap baris menjadi pengamatan).
Jadi, tentu saja, pertanyaan saya (sebagaimana dinyatakan dalam judul) harus masuk akal sekarang. Bagaimana kita melakukan regresi logistik menggunakan hasil fraksional sebagai DV? Apakah ada transformasi otomatis yang dapat dilakukan di glm?
Sejalan dengan hal yang sama, jika ada potensi pengukuran 3 (atau lebih), bagaimana seseorang melakukan ini untuk regresi logistik multinomial?
sumber
http://www.ats.ucla.edu/stat/r/dae/mlogit.htm
Jawaban:
The
glm
fungsi dalamR
memungkinkan 3 cara untuk menentukan formula untuk model regresi logistik.Yang paling umum adalah bahwa setiap baris bingkai data mewakili pengamatan tunggal dan variabel responsnya adalah 0 atau 1 (atau faktor dengan 2 level, atau varibale lainnya dengan hanya 2 nilai unik).
Pilihan lain adalah dengan menggunakan matriks 2 kolom sebagai variabel respon dengan kolom pertama menjadi jumlah 'keberhasilan' dan kolom kedua adalah jumlah 'kegagalan'.
Anda juga dapat menentukan respons sebagai proporsi antara 0 dan 1, lalu tentukan kolom lain sebagai 'bobot' yang memberikan jumlah total dari mana proporsinya (jadi respons 0,3 dan bobot 10 sama dengan 3 ' sukses 'dan 7' kegagalan ').
Salah satu dari 2 cara terakhir akan sesuai dengan apa yang Anda coba lakukan, yang terakhir tampaknya paling langsung untuk bagaimana Anda menggambarkan data Anda.
sumber
Sebagai permulaan, jika Anda memiliki variabel dependen yang proporsional, Anda dapat menggunakan Regresi Beta. Ini tidak meluas (dengan pengetahuan saya yang terbatas) hingga beberapa proporsi.
Untuk ikhtisar Regresi Beta dan implementasi R, periksa betareg .
sumber
Saya sudah menggunakan
nnet::multinom
(paket nnet adalah bagian dari MASS) untuk tujuan yang sama, ia menerima input kontinu dalam [0, 1].Jika Anda memerlukan referensi: C. Beleites et.al .: Raman penilaian spektroskopi jaringan astrositoma: menggunakan informasi referensi lunak. Anal Bioanal Chem, 2011, Vol. 400 (9), hlm. 2801-2816
sumber