Pertanyaan yang diberi tag «r»

Saya menggunakan paket randomForest di R (R versi 2.13.1, randomForest versi 4.6-2) untuk regresi dan melihat bias yang signifikan dalam hasil saya: kesalahan prediksi tergantung pada nilai variabel respon. Nilai tinggi di bawah prediksi dan nilai rendah di atas prediksi. Pada awalnya saya menduga ini adalah konsekuensi dari data saya, …

9 r  regression  random-forest 

Saya membaca A. Agresti (2007), Pengantar Analisis Data Kategorikal , 2. edisi, dan saya tidak yakin apakah saya memahami paragraf ini (hal.106, 4.2.1) dengan benar (walaupun seharusnya mudah): Dalam Tabel 3.1 tentang mendengkur dan penyakit jantung pada bab sebelumnya, 254 subjek melaporkan mendengkur setiap malam, di antaranya 30 menderita penyakit …

9 r  generalized-linear-model  logistic 

Model Bradley – Terry – Luce (BTL) menyatakan bahwa , di mana adalah probabilitas objek yang dinilai "lebih baik", lebih berat, dll, daripada objek , dan , dan adalah parameter.p i j j i δ i δ jhalj i= L o gi t- 1( δj- δsaya)pji=logit−1(δj−δi)p_{ji} = logit^{-1}(\delta_j - \delta_i)halsaya …

9 r  logistic  categorical-data  paired-data  bradley-terry-model 

Data saya adalah serangkaian waktu populasi yang dipekerjakan, L, dan rentang waktu, tahun. n.auto=auto.arima(log(L),xreg=year) summary(n.auto) Series: log(L) ARIMA(2,0,2) with non-zero mean Coefficients: ar1 ar2 ma1 ma2 intercept year 1.9122 -0.9567 -0.3082 0.0254 -3.5904 0.0074 s.e. NaN NaN NaN NaN 1.6058 0.0008 sigma^2 estimated as 1.503e-06: log likelihood=107.55 AIC=-201.1 AICc=-192.49 BIC=-193.79 …

9 r  regression  arima 

@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? Misalnya, kami berulang kali mengukur , , dan pada …

9 r  repeated-measures  simulation 

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 dipertimbangkan: model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003), offset=(log(population)), data=data, subset=28:36, family=poisson()) predict (model_glm, type="response") Saya mendapatkan case bukan rate ... Saya sudah mencoba juga model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003)+ offset(log(population)), data=data, …

9 r  generalized-linear-model  prediction  offset 

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 #Offset Namun, saya tidak yakin bagaimana e1071::svm()memilih kelas positif dan negatif, jadi …

9 r  svm  e1071 

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 korelasi untuk semua sel. Saya hanya ingin satu nomor sebagai …

9 r  correlation 

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 variabel yang berbeda. Saya ingin menerapkan algoritma genetika …

9 r  machine-learning  svm  genetic-algorithms 

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 menggunakan survdiff(Surv(time, status) ~ group, data = obs, rho = 0) …

9 r  survival  wilcoxon-signed-rank 

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 pada Senin / Selasa), tetapi secara drastis lebih sedikit …

9 r  time-series  seasonality  forecasting 

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 laporan berulang dari kode R saya dan visualisasi ggplot2. Saya mengerti bahwa LaTeX sepertinya …

9 r 

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 dalam objek pomp. Contoh di sini: http://cran.r-project.org/web/packages/pomp/vignettes/intro_to_pomp.pdf Katakanlah saya …

9 r 

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 dalam versi 1.7.3 dari paket R glmnet. Saya mengalami beberapa masalah menjalankan glmnet dengan …

9 r  multinomial  glmnet 

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 Proportional daripada dengan pengkodean yang sebenarnya. …

9 r  survival  cox-model