Saya ingin menggunakan BIC untuk pemilihan model HMM: BIC = -2*logLike + num_of_params * log(num_of_data) Jadi bagaimana cara menghitung jumlah parameter dalam model HMM. Pertimbangkan HMM 2-negara sederhana, di mana kami memiliki data berikut: data = [1 2 1 1 2 2 2 1 2 3 3 2 3 2 1 2 2 3 4 5 5...
Bootstrap, dalam bentuk standarnya, dapat digunakan untuk menghitung interval kepercayaan dari perkiraan statistik asalkan pengamatan itu benar. I. Visser et al. di " Interval Keyakinan untuk Parameter Model Hidden Markov ," menggunakan bootstrap parametrik untuk menghitung CI untuk parameter HMM....
Saya pikir fungsi tautan kanonik berasal dari parameter alami keluarga eksponensial. Katakan, pertimbangkan keluarga lalu adalah fungsi tautan kanonik. Ambil distribusi Bernoulli sebagai contoh, kami memiliki Jadi, fungsi tautan
Saya mencoba menganalisis data tindakan berulang dan berusaha untuk membuatnya bekerja R. Data saya pada dasarnya adalah sebagai berikut, saya memiliki dua kelompok perlakuan. Setiap subjek dalam setiap kelompok diuji setiap hari dan diberi skor (persentase benar pada suatu tes). Data dalam format...
Catatan: pertanyaan ini adalah repost, karena pertanyaan saya sebelumnya harus dihapus karena alasan hukum. Sambil membandingkan PROC CAMPURAN dari SAS dengan fungsi lmedari nlmepaket di R, saya menemukan beberapa perbedaan yang agak membingungkan. Lebih khusus lagi, derajat kebebasan dalam...
Saya baru-baru ini kebetulan membaca tentang Bayes empiris (Casella, 1985, Pengantar analisis data Bayes empiris) dan itu sangat mirip dengan model efek acak; karena keduanya memiliki estimasi yang menyusut menjadi rata-rata global. Tapi saya belum membacanya secara menyeluruh ... Adakah yang...
Saya mengalami masalah kinerja menggunakan MCMCglmmpaket di R untuk menjalankan model efek campuran. Kode ini terlihat seperti ini: MC1<-MCMCglmm(bull~1,random=~school,data=dt,family="categorical" , prior=list(R=list(V=1,fix=1), G=list(G1=list(V=1, nu=0))) , slice=T, nitt=iter, ,burnin=burn,...
Saya tidak dapat menemukan cara melakukan regresi linier dalam R dalam untuk desain ukuran berulang. Dalam pertanyaan sebelumnya (masih belum dijawab) disarankan kepada saya untuk tidak menggunakan lmmelainkan menggunakan model campuran. Saya menggunakan lmcara berikut:
Dengan rangkaian kovariat dan keluarga distribusi yang sama, bagaimana saya bisa membandingkan model yang memiliki fungsi tautan yang berbeda? Saya pikir jawaban yang benar di sini adalah "AIC / BIC", tetapi saya tidak 100% yakin. Apakah mungkin untuk memiliki model bersarang jika mereka memiliki...
Saya memiliki kumpulan data yang mewakili peluruhan eksponensial. Saya ingin mencocokkan fungsi eksponensial untuk data ini. Saya sudah mencoba log mentransformasikan variabel respons dan kemudian menggunakan kuadrat terkecil agar sesuai dengan garis; menggunakan model linier umum dengan fungsi log...
Saya telah diberikan satu set data yang berisi jumlah penghargaan yang diterima oleh siswa di satu sekolah menengah di mana prediktor jumlah penghargaan yang diperoleh termasuk jenis program di mana siswa terdaftar dan skor pada ujian akhir mereka dalam matematika. Saya bertanya-tanya apakah ada...
Saya kenal dengan notasi seperti: manaβ0j=β0+uj, danysaya j= β0+ βsayaxsaya j+ uj+ esaya j= β0 j+ βsayaxsaya j+ esaya jyij=β0+βixij+uj+eij=β0j+βixij+eij\begin{align} y_{ij} &= \beta_0 + \beta_i x_{ij} + u_j + e_{ij}\\ &= \beta_{0j} + \beta_i x_{ij} + e_{ij} \end{align}β0 j= β0+...
Saya memiliki dataset dengan 8000 cluster dan 4 juta pengamatan. Sayangnya perangkat lunak statistik saya, Stata, berjalan agak lambat ketika menggunakan fungsi data panel untuk regresi logistik:, xtlogitbahkan dengan subsampel 10%. Namun, ketika menggunakan logitfungsi fungsi nonpanel muncul...
Apa yang menjelaskan perbedaan dalam nilai-p dalam hal berikut aovdan lmpanggilan? Apakah perbedaannya hanya karena berbagai jenis perhitungan jumlah
Saya melakukan model linear umum, di mana saya harus menentukan keluarga yang berbeda dari yang normal. Berapa distribusi residu yang diharapkan? Misalnya, apakah residu harus didistribusikan secara
Saya ingin memahami apa yang dilakukan kode berikut. Orang yang menulis kode tidak lagi berfungsi di sini dan hampir tidak berdokumen. Saya diminta untuk menyelidikinya oleh seseorang yang berpikir " ini model regresi logistik bayesian " bglm <- function(Y,X) { # Y is a vector of binary...
Parameter apa start, etastart, mustartdi GLM () function ? Saya telah mencari di dokumen dan internet tetapi saya belum menemukan penjelasan yang jelas apa artinya ini. Ini menyerupai "nilai awal" bayesian untuk rantai, tapi saya ragu ini terkait, karena fungsi glm () dalam R adalah statistik...
Pertimbangkan model efek tak teramati linier dari tipe: mana adalah karakteristik yang tidak diobservasi tetapi waktu-invarian dan adalah kesalahan, dan indeks masing-masing pengamatan dan waktu. Pendekatan tipikal dalam regresi efek tetap (FE) adalah dengan menghapus melalui individu dummies...
Saya memiliki kumpulan data deret waktu yang saya coba paskan dengan Hidden Markov Model (HMM) untuk memperkirakan jumlah status laten dalam data. Kode pseudo saya untuk melakukan ini adalah sebagai berikut: for( i in 2 : max_number_of_states ){ ... calculate HMM with i states ......
Saya menyesuaikan model efek acak dengan glmerbeberapa data bisnis. Tujuannya adalah untuk menganalisis kinerja penjualan oleh distributor, dengan mempertimbangkan variasi regional. Saya memiliki variabel berikut: distcode: ID distributor, dengan sekitar 800 level region: ID geografis tingkat...