Dalam R, bagaimana cara menentukan model lmer tanpa efek tetap global? Misalnya, jika saya mengatakan sesuatu seperti lmer(y ~ (1 | group) + (0 + x | group), data = my_df) model yang pas akan ysaya j= a +αsaya+βsayaxsaya jyij=a+αi+βixijy_{ij} = a + \alpha_{i} + \beta_{i} x_{ij} Bagaimana saya...
Saya mencoba memvisualisasikan beberapa data konsumen, yang memiliki 4 kategori. Pengguna bebas untuk beralih di antara berbagai kategori. Saya ingin memvisualisasikan tiga atau empat sakelar terakhir untuk setiap individu. Jadi kita akan mulai dengan plot dengan kolom dengan 4 proporsi ditumpuk....
Saya memiliki plot kotak dari 13 grup yang saya tampilkan dalam satu plot. Kelompok-kelompok tersebut memiliki populasi yang tidak seimbang dan tidak terdistribusi secara normal. Saya ingin menunjukkan pasangan mana yang secara statistik serupa (yaitu, memiliki nilai p kruskal.test <0,05) dengan...
Jika saya memiliki variabel dengan 4 level, secara teori saya perlu menggunakan 3 variabel dummy. Dalam praktiknya, bagaimana ini sebenarnya dilakukan? Apakah saya menggunakan 0-3, apakah saya menggunakan 1-3 dan membiarkan 4 kosong? Ada saran? CATATAN: Saya akan bekerja di R. UPDATE: Apa yang...
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 ~...
Saya punya beberapa data deret waktu di mana variabel yang diukur adalah bilangan bulat positif (jumlah). Saya ingin menguji apakah ada tren naik dari waktu ke waktu (atau tidak). Variabel independen (x) berada dalam kisaran 0-500 dan variabel dependen (y) berada dalam kisaran 0-8. Saya pikir saya...
Tutup. Pertanyaan ini di luar topik . Saat ini tidak menerima jawaban. Ingin meningkatkan pertanyaan ini? Perbarui pertanyaan sehingga sesuai topik untuk Cross Validated. Ditutup 9 bulan lalu . Masalah yang disebutkan dalam pertanyaan ini diperbaiki...
Saya telah membuat versi termplot saya sendiri yang sedikit ditingkatkan yang saya gunakan dalam contoh ini, Anda dapat menemukannya di sini . Saya sebelumnya telah diposting di SO tetapi semakin saya memikirkannya, saya percaya bahwa ini mungkin lebih terkait dengan interpretasi model bahaya Cox...
Saya mencari contoh kode sederhana tentang cara menjalankan Partikel Filter di R. Paket pomp tampaknya mendukung bit matematika ruang keadaan, tetapi contoh-contoh agak sulit untuk mengikuti secara pemrograman untuk pengembang OO sederhana seperti saya, khususnya cara memuat data yang diamati ke...
Apakah ada prinsip umum tentang apakah seseorang harus menghitung korelasi pearson untuk dua variabel acak X dan Y sebelum mengambil transformasi log atau setelahnya? Apakah ada prosedur untuk menguji mana yang lebih tepat? Mereka menghasilkan nilai yang sama tetapi berbeda, karena transformasi log...
Tutup. Pertanyaan ini di luar topik . Saat ini tidak menerima jawaban. Ingin meningkatkan pertanyaan ini? Perbarui pertanyaan sehingga sesuai topik untuk Cross Validated. Ditutup 3 tahun yang lalu . Saya sangat tertarik mempelajari cara membuat
Sebagai hobi sampingan, saya telah menjelajahi rangkaian waktu perkiraan (khususnya, menggunakan R). Untuk data saya, saya memiliki jumlah kunjungan per hari, untuk setiap hari kembali hampir 4 tahun. Dalam data ini ada beberapa pola yang berbeda: Senin-Jumat memiliki banyak kunjungan (tertinggi...
Katakanlah Anda memiliki data kelangsungan hidup seperti ini: obs <- data.frame( time = c(floor(runif(100) * 30), floor((runif(100)^2) * 30)), status = c(rbinom(100, 1, 0.2), rbinom(100, 1, 0.7)), group = gl(2,100) ) Untuk melakukan tes peringkat log standar, seseorang dapat...
Saya menggunakan paket kernlab di R untuk membangun SVM untuk mengklasifikasikan beberapa data. SVM bekerja dengan baik karena memberikan 'prediksi' keakuratan yang layak, namun daftar variabel input saya lebih besar daripada yang saya inginkan dan saya tidak yakin akan pentingnya relatif dari...
Saya memiliki masalah dalam menggunakan cor()dan cor.test()fungsi. Saya hanya memiliki dua matriks (hanya nilai numerik, dan jumlah baris dan kolom yang sama) dan saya ingin memiliki nomor korelasi dan nilai p yang sesuai. Ketika saya menggunakan cor(matrix1, matrix2)saya mendapatkan koefisien...
Saya membutuhkan paket yang dapat memberi saya persamaan untuk model SVM linier. Saat ini saya menggunakan e1071 seperti: library(e1071) m = svm(data, labels, type='C', kernel='linear', cost=cost, probability=FALSE, scale=scale) w = t(m$coefs) %*% data[m$index,] #Weight vector b = -model$rho...
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
Kami memiliki kumpulan data dengan dua kovariat dan variabel pengelompokan kategoris dan ingin tahu apakah ada perbedaan yang signifikan antara kemiringan atau intersep di antara kovariat yang terkait dengan variabel pengelompokan yang berbeda. Kami telah menggunakan anova () dan lm () untuk...
@whuber telah menunjukkan cara mensimulasikan hasil multivarian ( , , dan ) untuk satu titik waktu.y 2 y 3y1y1y_1y2y2y_2y3y3y_3 Seperti yang kita ketahui, data longitudinal sering terjadi dalam studi medis. Pertanyaan saya adalah bagaimana mensimulasikan tindakan berulang hasil multivarian di R?...
Saya telah menggunakan Root Mean Squared Error(RMSE) untuk mengukur akurasi nilai yang diprediksi menggunakan model. Saya mengerti bahwa nilai yang dikembalikan menggunakan satuan ukuran saya (bukan persentase). Namun, saya ingin mengutip nilai saya sebagai persentase. Pendekatan yang saya ambil...