Pertanyaan yang diberi tag «r»

Saya mencari untuk menghasilkan 450 titik data dalam R. Ada tiga set yang berbeda 150 masing-masing didistribusikan dalam sebuah band lingkaran dengan jari-jari yang berbeda (pada 1, 2,8 dan 5). Secara khusus, saya ingin mereproduksi grafik pertama pada hal. 466 dari The Elements of Statistics Learning. Saya akan sangat berterima …

11 r  self-study  random-generation 

Saya mengalami masalah menggunakan 2l.normmetode imputasi bertingkat di mice. Sayangnya saya tidak dapat memposting contoh yang dapat direproduksi karena ukuran data saya - ketika saya mengurangi ukurannya, masalahnya hilang. Untuk variabel tertentu, micebuat kesalahan dan peringatan berikut: Error in chol.default(inv.sigma2[class] * X.SS[[class]] + inv.psi) : the leading minor of order …

11 r  missing-data  multiple-imputation  mice 

Dalam R, stepperintah ini dimaksudkan untuk membantu Anda memilih variabel input ke model Anda, bukan? Berikut ini berasal dari example(step)#-> swiss& step(lm1) > step(lm1) Start: AIC=190.69 Fertility ~ Agriculture + Examination + Education + Catholic + Infant.Mortality Df Sum of Sq RSS AIC - Examination 1 53.03 2158.1 189.86 <none> …

11 r  self-study  stepwise-regression 

Saya ingin menguraikan data deret waktu berikut ke dalam komponen musiman, tren, dan residual. Data tersebut adalah Profil Energi Pendinginan setiap jam dari sebuah bangunan komersial: TotalCoolingForDecompose.ts <- ts(TotalCoolingForDecompose, start=c(2012,3,18), freq=8765.81) plot(TotalCoolingForDecompose.ts) Oleh karena itu ada efek musiman harian dan mingguan yang didasarkan pada saran dari: Bagaimana cara menguraikan rangkaian …

10 r  time-series  forecasting  multiple-seasonalities  tbats 

tl; dr: Dimulai dengan dataset yang dihasilkan di bawah nol, saya melakukan resampled case dengan penggantian dan melakukan tes hipotesis pada setiap dataset yang di-resampled. Tes hipotesis ini menolak nol lebih dari 5% dari waktu. Di bawah ini, simulasi yang sangat sederhana, saya menghasilkan dataset dengan , dan saya memasangkan …

10 r  bootstrap  simulation  resampling 

Mengapa saya mendapatkan prediksi berbeda untuk ekspansi polinomial manual dan menggunakan polyfungsi R ? set.seed(0) x <- rnorm(10) y <- runif(10) plot(x,y,ylim=c(-0.5,1.5)) grid() # xp is a grid variable for ploting xp <- seq(-3,3,by=0.01) x_exp <- data.frame(f1=x,f2=x^2) fit <- lm(y~.-1,data=x_exp) xp_exp <- data.frame(f1=xp,f2=xp^2) yp <- predict(fit,xp_exp) lines(xp,yp) # using poly …

10 r  regression  polynomial 

Saya bertanya-tanya bagaimana cara pendekatan dengan benar pelatihan dan pengujian model LASSO menggunakan glmnet di R? Secara khusus, saya bertanya-tanya bagaimana cara melakukannya jika kurangnya set data uji eksternal mengharuskan saya gunakan validasi silang (atau pendekatan serupa lainnya) untuk menguji model LASSO saya. Biarkan saya memecah skenario saya: Saya hanya …

10 r  cross-validation  glmnet  lasso 

Seperti yang dinyatakan judul, saya mencoba mereplikasi hasil dari glmnet linear menggunakan pengoptimal LBFGS dari perpustakaan lbfgs. Pengoptimal ini memungkinkan kita untuk menambahkan istilah regularizer L1 tanpa harus khawatir tentang diferensiabilitas, selama fungsi objektif kami (tanpa istilah regularizer L1) adalah cembung. minβ∈Rp12n∥β0+Xβ−y∥22+αλ∥β∥1+12(1−α)λ∥β∥22minβ∈Rp12n‖β0+Xβ−y‖22+αλ‖β‖1+12(1−α)λ‖β‖22\min_{\beta \in \mathbb{R}^p} \frac{1}{2n}\Vert \beta_0 + X\beta - y …

10 r  regression  regularization  glmnet  elastic-net 

Jadi pertama-tama saya melakukan riset di forum ini, dan saya tahu pertanyaan yang sangat mirip telah ditanyakan tetapi biasanya mereka belum dijawab dengan benar atau kadang-kadang jawabannya tidak cukup detail untuk saya pahami. Jadi kali ini pertanyaan saya adalah: Saya punya dua set data, masing-masing, saya melakukan regresi polinomial seperti: …

10 r  regression  statistical-significance  regression-coefficients  group-differences 

Kemungkinan ada lebih dari satu kesalahpahaman serius dalam pertanyaan ini, tetapi ini tidak dimaksudkan untuk membuat perhitungannya benar, tetapi lebih untuk memotivasi pembelajaran deret waktu dengan beberapa fokus dalam pikiran. Dalam mencoba memahami penerapan deret waktu, tampaknya seolah-olah tren data membuat prediksi nilai masa depan menjadi tidak masuk akal. Misalnya, …

10 r  time-series  data-visualization 

Saya ingin memahami bagaimana saya bisa mendapatkan persentase varians dari kumpulan data, bukan di ruang koordinat yang disediakan oleh PCA, tetapi terhadap serangkaian vektor (rotasi) yang sedikit berbeda. set.seed(1234) xx <- rnorm(1000) yy <- xx * 0.5 + rnorm(1000, sd = 0.6) vecs <- cbind(xx, yy) plot(vecs, xlim = c(-4, …

10 r  variance  pca  linear-algebra 

Saya memiliki satu set data yang terdiri dari 482 observasi. data=Populationfull Saya akan melakukan analisis asosiasi genotipe untuk 3 SNP. Saya mencoba membangun model untuk analisis saya dan saya menggunakan aov (y ~ x, data = ...). Untuk satu sifat saya memiliki beberapa efek tetap dan kovariat yang telah saya …

10 r  anova 

Apa yang Anda sebut jenis plot ini, dan apakah mungkin membuatnya di R? EDIT: terima kasih banyak semua - sangat membantu. Judul terbaik sejauh ini: plot biola quantised!

10 r  data-visualization 

Saya melakukan penelitian tentang game edukasi, dan beberapa proyek saya saat ini melibatkan menggunakan data dari BoardGameGeek ( BGG ) dan VideoGameGeek (VGG) untuk menguji hubungan antara elemen desain game (yaitu, "diatur dalam Perang Dunia II", "melibatkan rolling dadu" ) dan peringkat pemain dari game-game tersebut (yaitu skor dari 10). …

10 r  categorical-data  factor-analysis  binary-data 

Bagaimana saya bisa mendapatkan p,d and qnilai dalam ARIMA(p,d,q)model yang diperkirakan oleh auto.arima(mytimeseries)? arima_model <- auto.arima (mytimeseries, ic = 'bic') Jika kita melihat output dari arima_model $ arma kita mendapatkan, [1] 1 0 0 0 1 2 0 Apa arti dari angka-angka yang muncul dalam urutan di atas?

10 r  arima