Saya mengerti model AR (p): inputnya adalah deret waktu yang dimodelkan. Saya benar-benar terjebak ketika membaca tentang model MA (q): inputnya adalah inovasi atau kejutan acak seperti yang sering dirumuskan. Masalahnya adalah saya tidak bisa membayangkan bagaimana mendapatkan komponen inovasi...
Ini adalah pertanyaan yang cukup umum: anggap saya ingin membangun model untuk memprediksi pengamatan berikutnya berdasarkan sebelumnya NNN pengamatan (NNNdapat menjadi parameter untuk mengoptimalkan secara eksperimental). Jadi pada dasarnya kami memiliki jendela geser fitur input untuk memprediksi...
Saya mencoba mencari tahu spesifikasi terbaik untuk dataset saya. Saya mencoba menyelidiki efektivitas zona ekonomi khusus di Polandia dalam arti pertumbuhan ekonomi dalam tiga model panel data yang sama untuk variabel yang dijelaskan: a) tingkat pengangguran terdaftar b) PDB per kapita c)...
Sejumlah sumber menunjukkan bahwa ada banyak konsekuensi negatif dari diskritisasi (kategorisasi) variabel kontinu sebelum analisis statistik (sampel referensi [1] - [4] di bawah). Sebaliknya [5] menunjukkan bahwa beberapa teknik pembelajaran mesin diketahui menghasilkan hasil yang lebih baik...
Saya melakukan analisis survival dalam R dengan survivalpaket. Saya pikir saya bekerja dengan data terpotong kiri, tapi saya tidak sepenuhnya yakin bagaimana menanganinya. Saya memiliki kohort pasien yang didiagnosis antara 1990 dan 2012. Semua pasien memiliki waktu diagnosis yang jelas (waktu...
Proses menghasilkan data adalah:y=sin(x+I(d=0))+sin(x+4∗I(d=1))+I(d=0)z2+3I(d=1)z2+N(0,1)y=sin(x+I(d=0))+sin(x+4∗I(d=1))+I(d=0)z2+3I(d=1)z2+N(0,1)y = \text{sin}\Big(x+I(d=0)\Big) + \text{sin}\Big(x+4*I(d=1)\Big) + I(d=0)z^2 + 3I(d=1)z^2 + \mathbb{N}\left(0,1\right) Misalkan menjadi urutan dari ke...
Pertama mari kita mensimulasikan beberapa data untuk regresi logistik dengan bagian-bagian tetap dan acak: set.seed(1) n <- 100 x <- runif(n) z <- sample(c(0,1), n, replace=TRUE) b <- rnorm(2) beta <- c(0.4, 0.8) X <- model.matrix(~x) Z <- cbind(z, 1-z) eta <- X%*%beta +...
Saya mengalami masalah dengan memeriksa asumsi untuk variabel kontinu dalam model bahaya proporsional. Jika variabel adalah faktor dengan banyak level, maka saya bisa menggunakan tes logrank atau memeriksa apakah log (-log) transformasi kurva survival adalah paralel. Tetapi bagaimana jika suatu...
Saya berjuang untuk sepenuhnya memahami beberapa notasi dalam buku di mana mereka menggunakan simbol "crosshair" - pertama seperti mana adalah matriks dan yang kedua seperti mana dan keduanya adalah matriks.⨁i=1nZj⨁i=1nZj\bigoplus\limits_{i=1}^n{} Z_j ZjZjZ_jIn⊗ΦIn⊗ΦI_n \otimes...
Beberapa dari Anda mungkin telah membaca makalah yang bagus ini: O'Hara RB, Kotze DJ (2010) Jangan log-transform data hitungan. Metode dalam Ekologi dan Evolusi 1: 118-122. klick . Saat ini saya sedang membandingkan model binomial negatif dengan model gaussian pada data yang diubah. Tidak seperti...
Saya mencari ukuran yang sesuai dari "varians yang dijelaskan" dari Poisson GLM (menggunakan fungsi log-link). Saya telah menemukan sejumlah sumber daya yang berbeda (baik di situs ini dan di tempat lain) yang membahas sejumlah pseudo- berbedaR2R2R^2 tindakan, tetapi hampir setiap situs...
Dalam buku mereka "Analisis Multilevel: Pengantar Pemodelan Multilevel Dasar dan Lanjutan" (1999), Snijders & Bosker (bab 8, bagian 8.2, halaman 119) mengatakan bahwa korelasi intersep-slope, dihitung sebagai kovarians intersep-slope yang dibagi oleh akar kuadrat dari produk intercept variance...
Saya ingin menentukan model regresi logistik di mana saya memiliki hubungan berikut: E[Yi|Xi]=f(βxi1+β2xi2)E[Yi|Xi]=f(βxi1+β2xi2)E[Y_i|X_i] = f(\beta x_{i1} + \beta^2x_{i2}) dimana fff adalah fungsi invit logit. Apakah ada cara "cepat" untuk melakukan ini dengan fungsi R yang sudah ada atau...
Saya mencoba memperkirakan kesenjangan upah gender untuk pekerja kantoran pria dan wanita di sebuah perusahaan besar Swedia untuk menguji apakah ada diskriminasi gender. Tes Hausman menolak nol bahwa efek tetap individu adalah acak dan karena itu saya tidak dapat mengandalkan OLS gabungan atau efek...
Saya telah menggunakan fungsi glm.fit di R untuk menyesuaikan parameter dengan model regresi logistik. Secara default, glm.fit menggunakan kuadrat terkecil yang diulang-ulang secara berulang agar sesuai dengan parameter. Apa alasan mengapa algoritma ini gagal untuk bertemu, ketika digunakan untuk...
Saya memiliki data yang menggambarkan seberapa sering suatu peristiwa berlangsung selama satu jam ("angka per jam", nph) dan berapa lama acara berlangsung ("durasi dalam detik per jam", dph). Ini adalah data asli: nph <- c(2.50000000003638, 3.78947368414551, 1.51456310682008, 5.84686774940732,...
Saya ingin mencocokkan model linier dengan R family=binomial(link="identity"), namun, keluarga binomial tidak memiliki tautan identitas. Apa yang harus saya
Saya ingin tahu tentang bagaimana paket lmerTest di R, khususnya fungsi "rand", menangani tes efek acak. Pertimbangkan contoh dari pdf lmerTest tentang CRAN yang menggunakan kumpulan data "wortel" bawaan : #import lme4 package and lmerTest package library(lmerTest) #lmer model with correlation...
Jadi saya bermain-main dengan ide menulis algoritma yang menumbuhkan dan memangkas pohon regresi dari data dan kemudian, di terminal node pohon, cocok dengan GLM. Saya sudah mencoba untuk membaca ide tetapi saya tidak bisa menemukan nama yang konsisten untuk teknik ini. Saya sudah membacanya...
Katakanlah kami tertarik pada bagaimana nilai ujian siswa dipengaruhi oleh jumlah jam belajar siswa tersebut. Kami mengambil sampel siswa dari beberapa sekolah yang berbeda. Kami menjalankan model efek campuran berikut: nilai ujiansaya= a +β1×jam. belajarsaya+sekolahj+esayanilai...