Apakah stan (khususnya, rstan) memiliki fasilitas bawaan untuk menghasilkan distribusi posterior prediktif? Tidak sulit untuk menghasilkan distribusi dari stan fit, tetapi saya lebih suka tidak menemukan kembali
Saya bertanya-tanya apakah ada cara yang baik untuk menghitung kriteria pengelompokan berdasarkan formula BIC, untuk keluaran k-means dalam R? Saya agak bingung bagaimana cara menghitung BIC sehingga saya bisa membandingkannya dengan model clustering lainnya. Saat ini saya menggunakan implementasi...
Saya mengeksplorasi sifat psikometrik dari ukuran laporan diri 10 item. Saya memiliki sekitar 400 kasus dalam dua sampel independen. Item selesai pada skala Likert 4 poin. EFA jelas mendukung solusi satu faktor (misalnya, nilai eigen pertama di atas 6, semua lainnya di bawah 1) dan alpha Cronbach...
karena saya melangkah ke peramalan dengan model ARIMA, saya mencoba memahami bagaimana saya dapat meningkatkan ramalan berdasarkan ARIMA yang sesuai dengan musiman dan pergeseran. Data saya adalah seri waktu berikut (lebih dari 3 tahun, dengan tren yang jelas ke atas dan musim yang terlihat, yang...
Jadi saya bekerja dengan model regresi logistik di R. Meskipun saya masih baru dalam statistik, saya merasa seperti saya mendapatkan sedikit pemahaman untuk model regresi sekarang, tetapi masih ada sesuatu yang mengganggu saya: Melihat gambar yang ditautkan, Anda melihat ringkasan R yang dicetak...
Di R, jika saya memanggil lm()fungsi dengan cara berikut: lm.1 = lm(response ~ var1 + var2 + var1 * var2) summary(lm.1) Ini memberi saya model linear dari variabel respons dengan var1, var2dan interaksi di antara mereka. Namun, bagaimana tepatnya kita menafsirkan istilah interaksi secara...
Saya sedang mencari persamaan yang tepat untuk menghitung kovarians sampel yang tidak bias tertimbang. Sumber-sumber internet sangat jarang pada tema ini dan mereka semua menggunakan persamaan yang berbeda. Persamaan yang paling mungkin saya temukan adalah persamaan
EDIT: Sejak membuat posting ini, saya telah menindaklanjuti dengan posting tambahan di sini . Ringkasan teks di bawah ini: Saya sedang mengerjakan model dan telah mencoba regresi linier, transformasi Box Cox dan GAM tetapi belum membuat banyak kemajuan Dengan menggunakan R, saya saat ini sedang...
Situasi Saya memiliki dataset dengan satu variabel dependen dan satu variabel bebas x . Saya ingin mencocokkan regresi linear piecewise berkelanjutan dengan k breakpoint diketahui / diperbaiki terjadi di ( a 1 , a 2 , ... , a k ) . Breakpoins dikenal tanpa ketidakpastian, jadi saya tidak ingin...
Saya memiliki randomForestmodel klasifikasi yang baik yang ingin saya gunakan dalam aplikasi yang memprediksi kelas kasus baru. Kasing baru memiliki nilai yang pasti hilang. Prediksi tidak akan berfungsi seperti itu untuk NAS. Bagaimana saya harus melakukan ini? data(iris) # create first the new...
Saya ingin memahami penggunaan simulasi Monte Carlo dalam chisq.test()fungsi di R. Saya memiliki variabel kualitatif yang memiliki 128 level / kelas. Ukuran sampel saya adalah 26 (saya tidak dapat mengambil sampel lebih banyak "individu"). Jadi jelas, saya akan memiliki beberapa level dengan 0...
Latar Belakang Dalam sebuah makalah dari Epstein (1991): Pada memperoleh nilai-nilai klimatologis harian dari rata-rata bulanan , formulasi dan algoritma untuk menghitung interpolasi Fourier untuk nilai periodik dan genap spasi diberikan. Di koran, tujuannya adalah untuk mendapatkan nilai harian...
Saya bukan ahli, jadi maafkan saya jika beberapa terminologi sedikit canggung. Senang memberikan informasi lebih lanjut di mana diperlukan. Saya memiliki dua vektor dari 50 nilai numerik berpasangan di R. Saya ingin melakukan uji pengacakan atau permutasi dua sisi untuk menentukan apakah...
Saya melakukan meta-analisis ukuran efek d dalam R menggunakan paket metafor. d mewakili perbedaan dalam skor memori antara pasien dan sehat. Namun beberapa penelitian melaporkan hanya subskala ukuran minat d (misalnya beberapa skor memori yang berbeda atau skor dari tiga blok pengujian memori yang...
Ini mungkin pertanyaan yang konyol, tetapi ketika membuat model dengan tanda sisipan dan menggunakan sesuatu seperti LOOCVatau (bahkan lebih tepatnya) LGOCV, apa manfaat membagi data ke dalam kereta dan set tes jika ini pada dasarnya merupakan langkah validasi silang ngomong-ngomong Saya membaca...
Saya punya model seperti ini: require(nlme) set.seed(123) n <- 100 k <- 5 cat <- as.factor(rep(1:k, n)) cat_i <- 1:k # intercept per kategorie x <- rep(1:n, each = k) sigma <- 0.2 alpha <- 0.001 y <- cat_i[cat] + alpha * x + rnorm(n*k, 0, sigma) plot(x, y) m1 <- lm(y ~...
Saya berpikir tentang masalah yang memprediksi log (menghabiskan) pelanggan menggunakan regresi linier. Saya sedang mempertimbangkan fitur apa yang akan digunakan sebagai input dan bertanya-tanya apakah akan baik-baik saja untuk menggunakan persentil variabel sebagai input. Misalnya saya bisa...
Katakanlah kita memiliki data berikut: set.seed(123) data <- data.frame(x = c(rnorm(50, 1, 1), rnorm(50, 5, 2)), y = c(rep('A', 50), rep('B', 50))) Yang menghasilkan boxplot berikut ( boxplot(data$x ~ data$y)): Sekarang katakanlah saya ingin menguji apakah dua sampel memiliki parameter...
Saya memiliki data kecelakaan kendaraan bermotor menurut jam sehari. Seperti yang Anda harapkan, mereka tinggi di tengah hari dan puncaknya pada jam sibuk. Geom_density default ggplot2 memuluskannya dengan baik Subset dari data, untuk crash terkait drive-minum, tinggi di kedua ujung hari (malam...
Saya punya contoh data yang diatur sebagai berikut: Volume <- seq(1,20,0.1) var1 <- 100 x2 <- 1000000 x3 <- 30 x4 = sqrt(x2/pi) H = x3 - Volume r = (x4*H)/(H + Volume) Power = (var1*x2)/(100*(pi*Volume/3)*(x4*x4 + x4*r + r*r)) Power <- jitter(Power, factor = 1, amount =...