Katakanlah saya menghasilkan probabilitas hasil berdasarkan faktor tertentu dan plot kurva hasil itu. Apakah ada cara untuk mengekstrak persamaan untuk kurva itu dari R? > mod = glm(winner~our_bid, data=mydat, family=binomial(link="logit")) > summary(mod) Call: glm(formula = winner ~...
Apa perbedaan antara model GLM (regresi logistik) dengan variabel respon biner yang mencakup subjek dan waktu sebagai kovariat dan model GEE analog yang memperhitungkan korelasi antara pengukuran pada berbagai titik waktu? GLM saya terlihat seperti: Y(binary) ~ A + B1X1(subject id) + B2X2(time)...
Bagaimana saya menggunakan output dari analisis komponen utama (PCA) dalam model linier umum (GLM), dengan asumsi PCA digunakan untuk pemilihan variabel untuk GLM? Klarifikasi: Saya ingin menggunakan PCA untuk menghindari penggunaan variabel yang berhubungan dalam GLM. Namun, PCA memberi saya...
Saya tahu ini mungkin pertanyaan dasar ... Tapi sepertinya saya tidak menemukan jawabannya. Saya pas GLM dengan keluarga Poisson, dan kemudian mencoba untuk melihat prediksi, namun offset tampaknya
Saya membaca A. Agresti (2007), Pengantar Analisis Data Kategorikal , 2. edisi, dan saya tidak yakin apakah saya memahami paragraf ini (hal.106, 4.2.1) dengan benar (walaupun seharusnya mudah): Dalam Tabel 3.1 tentang mendengkur dan penyakit jantung pada bab sebelumnya, 254 subjek melaporkan...
Saat ini saya sedang mengerjakan serangkaian model deret waktu Poisson yang mencoba memperkirakan efek perubahan dalam bagaimana penghitungan diperoleh (beralih dari satu uji diagnostik ke uji lain) sambil mengendalikan tren lain dari waktu ke waktu (katakanlah peningkatan umum pada kejadian...
Tutup. Pertanyaan ini di luar topik . Saat ini tidak menerima jawaban. Ingin meningkatkan pertanyaan ini? Perbarui pertanyaan sehingga sesuai topik untuk Cross Validated. Ditutup 2 tahun yang lalu . Kode berikut PredictNew <- predict (glm.fit,...
Dengan data kontinu, regresi linier mengasumsikan bahwa istilah kesalahan didistribusikan N (0, )σ 2Y=β1+β2X2+uY=β1+β2X2+uY=\beta_1+\beta_2X_2+uσ2σ2\sigma^2 1) Apakah kita berasumsi bahwa Var (Y | x) juga ~ N (0, )?σ2σ2\sigma^2 2) Apa distribusi kesalahan ini dalam regresi logistik? Ketika data...
Saya mendapatkan koefisien yang sangat besar selama regresi logistik, lihat koefisien dengan krajULKV: > summary(m5) Call: glm(formula = cbind(ml, ad) ~ rok + obdobi + kraj + resid_usili2 + rok:obdobi + rok:kraj + obdobi:kraj + kraj:resid_usili2 + rok:obdobi:kraj, family =...
Pada halaman 232 dari "Teman R untuk menerapkan regresi" catatan Fox dan Weisberg Hanya keluarga Gaussian yang memiliki varian konstan, dan dalam semua GLM lainnya, varian bersyarat dari y pada tergantung pada μ ( x )xx\bf{x}μ ( x )μ(x)\mu(x) Sebelumnya, mereka mencatat bahwa varian bersyarat...
Secara default ketika kami menggunakan glmfungsi dalam R, ia menggunakan metode iteratively reweighted least square (IWLS) untuk menemukan estimasi kemungkinan maksimum parameter. Sekarang saya punya dua pertanyaan. Apakah estimasi IWLS menjamin maksimum global dari fungsi kemungkinan?...
Saya punya data ini: set.seed(1) predictor <- rnorm(20) set.seed(1) counts <- c(sample(1:1000, 20)) df <- data.frame(counts, predictor) Saya menjalankan regresi poisson poisson_counts <- glm(counts ~ predictor, data = df, family = "poisson") Dan regresi binomial
Saya memiliki kumpulan data sekitar 5000 fitur. Untuk data itu saya pertama kali menggunakan uji Chi Square untuk pemilihan fitur; setelah itu, saya mendapat sekitar 1500 variabel yang menunjukkan hubungan signifikan dengan variabel respon. Sekarang saya harus menyesuaikan regresi logistik. Saya...
Saya memiliki data dari desain eksperimental berikut: pengamatan saya adalah jumlah dari jumlah keberhasilan ( K) dari jumlah percobaan yang sesuai ( N), diukur untuk dua kelompok yang masing-masing terdiri dari Iindividu, dari Tperawatan, di mana dalam setiap kombinasi faktor tersebut terdapat...
Saya pas beberapa model formulir .. glm(DV ~ I(1/IV), family = Gamma(link = "log") .. dan saya mencari cara untuk membandingkan model yang diperoleh untuk variabel yang berbeda. Saya bertanya-tanya apakah nilai alfa adalah penggunaan praktis? Untuk tiga plot di bawah ini nilai alpha adalah...
Saya mencoba membantu seorang mahasiswa dari seorang kolega. Siswa mengamati dan menghitung perilaku burung (jumlah panggilan) dalam pengaturan eksperimental. Jumlah panggilan yang dikaitkan dengan burung yang diamati secara spesifik selama setiap percobaan tidak dapat ditentukan tetapi menghitung...
Saya mensimulasikan uji coba Bernoulli dengan acak logitθ∼N(logitθ0,12)logitθ∼N(logitθ0,12)\text{logit}\, \theta \sim {\cal N}(\text{logit}\, \theta_0, 1^2)antara kelompok dan kemudian saya cocokkan dengan model yang sesuai dengan lme4paket: library(lme4) library(data.table) I <- 30 # number of...
Model logistik bertingkat berikut dengan satu variabel penjelas pada level 1 (level individu) dan satu variabel penjelas pada level 2 (level grup): π 0 j = γ 00 + γ 01 z j + u 0 j … ( 2 ) π 1 j = γ 10 + γ 11 z j + u 1 j … ( 3 )logit ( hlmsaya j) = π0 j+ π1 jxsaya j… ( 1...
Di bidang saya, cara yang biasa untuk memplot data berpasangan adalah sebagai serangkaian segmen garis miring tipis, menimpanya dengan median dan CI median untuk dua kelompok: Namun, plot semacam ini menjadi lebih sulit dibaca karena jumlah titik data menjadi sangat besar (dalam kasus saya yang...
Saya sedang mempelajari perbedaan antara regularisasi dalam regresi RKHS dan regresi linier, tetapi saya kesulitan memahami perbedaan penting antara keduanya. (xi,yi)(xi,yi)(x_i,y_i)f(⋅)f(⋅)f(\cdot)f(x)≈u(x)=∑i=1mαiK(x,xi),f(x)≈u(x)=∑i=1mαiK(x,xi),\begin{equation}f(x)\approx u(x)=\sum_{i=1}^m...