Saya telah meneliti paket tikus , dan saya belum menemukan cara untuk menggunakan beberapa imputasi untuk membuat model Cox, kemudian memvalidasi model itu dengan fungsi paket rmsvalidate() . Berikut adalah beberapa contoh kode dari apa yang saya miliki sejauh ini, menggunakan kumpulan data...
Saya mencoba melakukan analisis pertumbuhan kelas laten (LCGA) dan / atau model campuran pertumbuhan (GMM) di R. Data yang saya gunakan adalah peningkatan jumlah garpu repositori git (variabel diskrit, bukan kategorikal), yang Anda bisa lihat di dataset ini . Saya mencoba lavaan, yang membantu...
Bagaimana saya bisa menghitung sensitivitas dan spesifisitas (atau tindakan analog) dari tes diagnostik berkelanjutan dalam memprediksi hasil yang berkelanjutan (misalnya, tekanan darah) tanpa mendikotomisasi hasilnya? Ada ide? Tampaknya para peneliti telah melakukan ini menggunakan pemodelan efek...
Saya memiliki masalah dalam menemukan solusi tentang cara menjalankan tes post-hoc (Tukey HSD) setelah 2-faktor (keduanya dalam-mata pelajaran) mengukur ANOVA berulang dalam R. Untuk ANOVA, saya telah menggunakan fungsi aov: summary(aov(dv ~ x1 * x2 + Error(subject/(x1*x2)), data=df1)) Setelah...
Saya menggunakan lme4paket dalam R untuk melakukan beberapa pemodelan efek campuran logistik. Pemahaman saya adalah bahwa jumlah setiap efek acak harus nol. Ketika saya membuat mainan model linier campuran menggunakan lmer, efek acak biasanya <10−1010−1010^{-10}mengkonfirmasikan keyakinan saya...
Saya mencoba memahami apa perbedaan antara HMM standar dan HMM Bayesian. Wikipedia hanya menyebutkan secara singkat bagaimana modelnya tetapi saya perlu tutorial yang lebih rinci. Apakah ada yang tahu tentang makalah atau implementasi yang bisa saya lihat? Saya juga punya masalah dengan...
Saya mencoba menjalankan regresi Cox pada dataset 2.000.000 sampel baris sebagai berikut hanya menggunakan R. Ini adalah terjemahan langsung dari PHREG di SAS. Sampel tersebut mewakili struktur dataset asli. ## library(survival) ### Replace 100000 by 2,000,000 test <-
Dari Wikipedia General linear model (GLM) adalah model linear statistik. Ini dapat ditulis sebagai 1 mana adalah matriks dengan serangkaian pengukuran multivariat, adalah matriks yang mungkin merupakan matriks desain , adalah matriks yang berisi parameter yang biasanya diperkirakan dan adalah...
Jadi, bagi saya kelihatannya fungsi bobot dalam lm memberikan bobot pengamatan lebih besar, semakin besar nilai 'bobot' pengamatan, sedangkan fungsi lme di lme justru melakukan sebaliknya. Ini dapat diverifikasi dengan simulasi sederhana. #make 3 vectors- c is used as an uninformative random...
Ini adalah spin-off dari pertanyaan ini: Bagaimana membandingkan dua kelompok dengan beberapa pengukuran untuk setiap individu dengan R? Dalam jawaban di sana (jika saya mengerti dengan benar) saya belajar bahwa varians dalam subjek tidak mempengaruhi kesimpulan yang dibuat tentang mean kelompok...
Pertanyaan keseluruhan saya adalah: mengapa menggunakan bayesglmmetode klasifikasi lainnya? catatan: Saya hanya tertarik pada prediksi. Saya memiliki jumlah data yang layak (~ 100.000 obs). Saya merasa ukuran sampel cukup besar sehingga parameter regresi logistik reguler akan terdistribusi...
Subitisasi adalah penghitungan cepat, akurat untuk tampilan dengan angka rendah, dibedakan dari penghitungan dengan non-linearitas yang tajam dalam plot waktu respons. Di bawah ini adalah plot yang representatif, dari Watson, DG, Maylor, EA, & Bruce, LAM (2007). Perhatikan bahwa waktu...
Saya mencoba untuk membuat model dengan menggunakan regresi binomial negatif (GLM binomial negatif). Saya memiliki ukuran sampel yang relatif kecil (lebih dari 300), dan data tidak diskalakan. Saya perhatikan bahwa ada dua cara untuk mengukur good of fit - satu penyimpangan dan yang lainnya adalah...
Apakah ada yang bisa saya lakukan LASSO dengan Regresi Binomial Negatif pada R? Saya melakukan regresi binomial negatif pada dataset saya karena data terlalu tersebar untuk memaksakan regresi poisson. Sementara itu, saya juga menghadapi beberapa masalah multikolinieritas. Saya sudah mencoba...
Saya memiliki tabel berikut ini di R df <- structure(list(x = structure(c(12458, 12633, 12692, 12830, 13369, 13455, 13458, 13515), class = "Date"), y = c(6080, 6949, 7076, 7818, 0, 0, 10765, 11153)), .Names = c("x", "y"), row.names = c("1", "2", "3", "4", "5", "6", "8", "9"), class =...
Saya memiliki model linier umum berikut. Objek glmDVdimodelkan sebagai proporsi keberhasilan atas uji coba total. Objek x_iadalah variabel kontinu. Seperti apakah ini dalam notasi matematika? winp.glm = glm(glmDV ~ x1 + x2 + x3 + x4 + x5 + x6 + x7, data=myData,
Saya memiliki satu set data yang saya coba lakukan regresi, dan gagal. Situasi: Ribuan operator robot pertempuran saling bertarung satu sama lain menggunakan robot pertempuran. Beberapa robot pertempuran kuat dan kuat, dan yang lainnya lemah; yang kuat menang lebih sering dan memberikan lebih...
Proyek saya saat ini mungkin mengharuskan saya untuk membuat model untuk memprediksi perilaku sekelompok orang tertentu. set data pelatihan hanya berisi 6 variabel (id hanya untuk tujuan identifikasi): id, age, income, gender, job category, monthly spend di mana monthly spendadalah variabel...
Saya sedang mengerjakan model data hitung nol-naik menggunakan psclpaket. Saya hanya ingin tahu mengapa tidak ada pengembangan model untuk model data hitungan satu kali lipat! Juga mengapa tidak ada pengembangan model data bimodal, katakanlah zero-and-2-inflated! Setelah saya menghasilkan data...
Membiarkan E( h )E(h)\mathcal{E(h)} sebuah fungsi yang diberi beberapa hipotesis hhh mengembalikan kesalahan generalisasi untuk perbaikan itu hhh. Saya membaca beberapa catatan tentang pemilihan model dan kesalahan generalisasi dan dikatakan: "Jika kita memiliki akses ke E( h