Pertanyaan yang diberi tag «r»

Saya mencoba mensimulasikan dataset yang cocok dengan data empiris yang saya miliki, tetapi saya tidak yakin bagaimana memperkirakan kesalahan dalam data asli. Data empiris mencakup heteroskedastisitas, tetapi saya tidak tertarik untuk mengubahnya, tetapi menggunakan model linier dengan istilah kesalahan untuk mereproduksi simulasi data empiris. Sebagai contoh, katakanlah saya memiliki beberapa …

9 r  simulation  heteroscedasticity  gamlss  dglm 

Saya mencoba membuat model yang diperkecil untuk memprediksi banyak variabel dependen (DV) (~ 450) yang sangat berkorelasi. Variabel independen saya (IV) juga banyak (~ 2000) dan sangat berkorelasi. Jika saya menggunakan laso untuk memilih model yang direduksi untuk setiap output secara terpisah, saya tidak dijamin mendapatkan subset variabel independen yang …

9 r  lasso  multivariate-regression 

Saya menyesuaikan model efek campuran dengan istilah spline dalam aplikasi di mana tren dari waktu ke waktu dikenal sebagai kurva-linier. Namun, apa yang ingin saya nilai adalah apakah tren kurva-linier terjadi karena penyimpangan individu dari linieritas, atau apakah itu efek pada tingkat kelompok yang membuat tingkat kelompok yang sesuai tampak …

9 r  splines  lme4-nlme 

Saya baru mendukung mesin vektor. Penjelasan singkat The svmfungsi dari e1071paket di R menawarkan berbagai pilihan: Klasifikasi-C nu-klasifikasi satu-klasifikasi (untuk deteksi kebaruan) regresi eps nu-regresi Apa perbedaan intuitif antara lima jenis? Yang mana yang harus diterapkan dalam situasi apa?

9 r  classification  svm  e1071 

Katakanlah saya memiliki data yang memiliki ketidakpastian. Sebagai contoh: X Y 1 10±4 2 50±3 3 80±7 4 105±1 5 120±9 Sifat ketidakpastian dapat berupa pengukuran ulang atau eksperimen, atau misalnya ketidakpastian instrumen pengukuran. Saya ingin menyesuaikan kurva menggunakan R, sesuatu yang biasanya saya lakukan lm. Namun, ini tidak memperhitungkan …

9 r  least-squares  error-propagation 

Karena standard error dari regresi linier biasanya diberikan untuk variabel respon, saya bertanya-tanya bagaimana cara mendapatkan interval kepercayaan di arah lain - misalnya untuk x-intersep. Saya dapat memvisualisasikan apa yang mungkin terjadi, tetapi saya yakin pasti ada cara mudah untuk melakukan ini. Di bawah ini adalah contoh dalam R bagaimana …

9 r  regression  confidence-interval  bootstrap 

Adakah yang bisa memberi tahu saya mengapa saya mendapatkan hasil yang berbeda dari Rkuadrat terkecil tertimbang dan solusi manual dengan operasi matriks ? Secara khusus, saya mencoba untuk secara manual menyelesaikan , di mana W adalah matriks diagonal pada bobot, A adalah matriks data, b adalah vektor respons. WAx=WbWAx=Wb\mathbf W …

9 r  regression  least-squares  weighted-regression  weighted-data 

Singkatnya pertanyaan saya: apakah ada metode untuk meningkatkan waktu berjalan R MICE (data imputasi)? Saya bekerja dengan satu set data (30 variabel, 1,3 juta baris) yang berisi (secara acak) data yang hilang. Sekitar 8% dari pengamatan di sekitar 15 dari 30 variabel mengandung NAs. Untuk memasukkan data yang hilang, saya …

9 r  multiple-imputation  mice 

Bagaimana cara memberi bobot lebih pada pengamatan terbaru di R? Saya menganggap ini sebagai pertanyaan atau keinginan yang umum diajukan, tetapi saya memiliki kesulitan untuk mencari tahu bagaimana menerapkannya. Saya telah mencoba mencari banyak untuk ini tetapi saya tidak dapat menemukan contoh praktis yang baik. Dalam contoh saya, saya akan …

9 r  time-series  regression  random-forest 

Menggunakan R pada beberapa data dan mencoba melihat apakah data saya heteroscedastic atau tidak, saya telah menemukan dua implementasi tes Breusch-Pagan, bptest (paket lmtest) dan ncvTest (paket mobil). Namun, ini menghasilkan hasil yang berbeda. Apa perbedaan keduanya? Kapan sebaiknya Anda memilih untuk menggunakan satu atau yang lain? > model <- …

9 r  regression  heteroscedasticity  assumptions  breusch-pagan 

Saya mencoba untuk membuat kelompok kumpulan data yang berbeda menggunakan algoritma yang tidak diawasi (clustering). Masalahnya adalah saya memiliki banyak fitur (~ 500) dan sejumlah kecil case (200-300). Sejauh ini saya hanya melakukan masalah klasifikasi yang selalu saya beri label data sebagai set pelatihan. Di sana saya menggunakan beberapa kriteria …

9 r  clustering  feature-selection  unsupervised-learning 

Saya mencoba mengonversi kolom faktor menjadi variabel dummy: str(cards$pointsBin) # Factor w/ 5 levels ".lte100",".lte150",..: 3 2 3 1 4 4 2 2 4 4 ... labels <- model.matrix(~ pointsBin, data=cards) head(labels) # (Intercept) pointsBin.lte150 pointsBin.lte200 pointsBin.lte250 pointsBin.lte300 # 741 1 0 0 0 0 # 407 1 1 0 …

9 r  categorical-data  categorical-encoding 

Masalah Saya ingin melakukan beberapa kesimpulan tentang sistem analog dengan mati dengan jumlah sisi yang tidak diketahui. Die digulung beberapa kali, setelah itu saya ingin menyimpulkan distribusi probabilitas atas parameter yang sesuai dengan jumlah sisi yang dimiliki die, θ. Intuisi Jika setelah 40 gulungan Anda telah mengamati 10 merah, 10 …

9 r  probability  bayesian  jags 

Saya menggunakan lmer di R untuk memeriksa efek kondisi ( cond) pada beberapa hasil. Berikut adalah beberapa data yang dibuat, di mana s adalah pengidentifikasi subjek dan a, bdan ckondisi. library("tidyr") library("dplyr") set.seed(123) temp <- data.frame(s = paste0("S", 1:30), a = rnorm(30, -2, 1), b = rnorm(30, -3, 1), c …

9 r  lme4-nlme  contrasts 

Saya sudah mendapatkan data bulanan dari tahun 1993 hingga 2015 dan ingin melakukan perkiraan data ini. Saya menggunakan paket tsoutliers untuk mendeteksi outliers, tetapi saya tidak tahu bagaimana cara melanjutkan perkiraan dengan set data saya. Ini kode saya: product.outlier<-tso(product,types=c("AO","LS","TC")) plot(product.outlier) Ini adalah output saya dari paket tsoutliers ARIMA(0,1,0)(0,0,1)[12] Coefficients: sma1 …

9 r  time-series  forecasting  arima  outliers