Pertanyaan yang diberi tag «r»

Apakah "residu pelajar" dan "residu standar" sama dalam model regresi? Saya membangun model regresi linier dalam R dan ingin memplot grafik nilai residu Studentized v / s, tetapi tidak menemukan cara otomatis untuk melakukan ini di R. Misalkan saya punya model library(MASS) lm.fit <- lm(Boston$medv~(Boston$lstat)) kemudian menggunakan plot(lm.fit)tidak memberikan plot …

10 r  regression  residuals  terminology 

Situasi: Dua burung (jantan dan betina) melindungi telur mereka di sarang dari pengganggu. Setiap burung dapat menggunakan serangan atau ancaman untuk perlindungan, dan hadir atau tidak ada. Ada pola yang muncul dari data bahwa perilaku mungkin saling melengkapi - serangan laki-laki sementara perempuan menggunakan tampilan ancaman dan sebaliknya. Pertanyaan saya …

10 r  sequence-analysis 

Muat paket yang dibutuhkan. library(ggplot2) library(MASS) Hasilkan 10.000 angka yang pas untuk distribusi gamma. x <- round(rgamma(100000,shape = 2,rate = 0.2),1) x <- x[which(x>0)] Gambarkan fungsi kerapatan probabilitas, seandainya kita tidak tahu distribusi x mana yang cocok. t1 <- as.data.frame(table(x)) names(t1) <- c("x","y") t1 <- transform(t1,x=as.numeric(as.character(x))) t1$y <- t1$y/sum(t1[,2]) ggplot() …

10 r  mathematical-statistics  goodness-of-fit  gamma-distribution  ggplot2 

Konteks Saya mencoba memahami bagaimana R's coxph () menerima dan menangani entri berulang untuk subjek (atau pasien / pelanggan jika Anda suka). Beberapa menyebutnya format Panjang ini, yang lain menyebutnya 'tindakan berulang'. Lihat misalnya kumpulan data yang menyertakan kolom ID di bagian Jawaban di: Paket terbaik untuk model Cox dengan …

10 r  repeated-measures  survival  cox-model  frailty 

Menurut dokumentasi removeSparseTermsfungsi dari tmpaket, inilah yang diperlukan sparsity: A term-document matrix where those terms from x are removed which have at least a sparse percentage of empty (i.e., terms occurring 0 times in a document) elements. I.e., the resulting matrix contains only terms with a sparse factor of less …

10 r  text-mining  natural-language 

Fungsi R cv.glm (library: boot) menghitung perkiraan kesalahan prediksi validasi silang K-fold untuk model linier umum dan mengembalikan delta. Apakah masuk akal untuk menggunakan fungsi ini untuk regresi laso (library: glmnet) dan jika demikian, bagaimana hal itu dapat dilakukan? Pustaka glmnet menggunakan cross-validation untuk mendapatkan parameter belokan terbaik, tapi saya …

10 r  regression  cross-validation  lasso  glmnet 

Saya menganalisis data pelacakan mata dari percobaan yang dirancang. Versi sederhana dari data saya terlihat seperti ini (Anda bisa mendapatkan data dput () di sini ), head(lookDATA) participant fixationImage fixationCount 1 9 Automobile 81 2 9 Bird 63 3 9 Chair 82 4 9 Dog 64 5 9 Face 90 …

10 r  mixed-model  poisson-regression  lsmeans  lme4-nlme 

Apakah mungkin untuk melakukan analisis kekuatan untuk uji Kolmogorov Smirnov 2-sisi pada R? Saya menguji apakah dua distribusi empiris berbeda menggunakan ks.test (), dan saya ingin menambahkan analisis daya. Saya belum dapat menemukan analisis daya bawaan untuk tes KS di R. Ada saran? Sunting : Ini adalah distribusi yang dihasilkan …

10 r  power-analysis  kolmogorov-smirnov 

Saya menggunakan rumus dari pengujian ab Bayesian untuk menghitung hasil tes AB menggunakan metodologi Bayesian. Pr(pB>pA)=∑i=0αB−1B(αA+i,βB+βA)(βB+i)B(1+i,βB)B(αA,βA)Pr(pB>pA)=∑i=0αB−1B(αA+i,βB+βA)(βB+i)B(1+i,βB)B(αA,βA) \Pr(p_B > p_A) = \sum^{\alpha_B-1}_{i=0} \frac{B(\alpha_A+i,\beta_B+\beta_A)}{(\beta_B+i)B(1+i,\beta_B)B(\alpha_A, \beta_A)} dimana αAαA\alpha_A dalam satu ditambah jumlah keberhasilan untuk A βAβA\beta_A dalam satu ditambah jumlah kegagalan untuk A αBαB\alpha_B dalam satu ditambah jumlah keberhasilan untuk B βBβB\beta_B dalam …

10 r  bayesian  ab-test 

Saya menganalisa beberapa data perilaku menggunakan lme4di R, sebagian besar mengikuti tutorial Bodo Winter baik , tapi saya tidak mengerti jika saya menangani interaksi dengan benar. Lebih buruk lagi, tidak ada orang lain yang terlibat dalam penelitian ini menggunakan model campuran, jadi saya sedikit terombang-ambing ketika datang untuk memastikan semuanya …

10 r  mixed-model  p-value  lme4-nlme 

Saya melakukan analisis deret waktu menggunakan R. Saya harus menguraikan data saya menjadi komponen tren, musiman dan acak. Saya memiliki data mingguan selama 3 tahun. Saya telah menemukan dua fungsi di R - stl()dan decompose(). Saya telah membaca bahwa stl()ini tidak baik untuk dekomposisi multiplikasi. Adakah yang bisa memberi tahu …

10 r  time-series 

Dari Quick-R Robert Kabacoff yang saya miliki # Bootstrap 95% CI for regression coefficients library(boot) # function to obtain regression weights bs <- function(formula, data, indices) { d <- data[indices,] # allows boot to select sample fit <- lm(formula, data=d) return(coef(fit)) } # bootstrapping with 1000 replications results <- boot(data=mtcars, …

10 r  regression  p-value  bootstrap 

Perbandingan dua sampel proporsi, estimasi ukuran sampel: R vs Stata Saya mendapat hasil berbeda untuk ukuran sampel, sebagai berikut: Dalam R power.prop.test(p1 = 0.70, p2 = 0.85, power = 0.90, sig.level = 0.05) Hasil: (jadi 161) untuk setiap kelompok.n=160.7777n=160.7777n = 160.7777 Di Stata sampsi 0.70 0.85, power(0.90) alpha(0.05) Hasil: untuk …

10 r  stata  sample-size  power-analysis  jmp 

Saya belum menemukan literatur tentang penerapan Hutan Acak untuk MNIST, CIFAR, STL-10, dll. Jadi saya pikir saya akan mencobanya sendiri dengan permutasi MNIST invarian . Di R , saya mencoba: randomForest(train$x, factor(train$y), test$x, factor(test$y), ntree=500) Ini berjalan selama 2 jam dan mendapat 2,8% kesalahan pengujian. Saya juga mencoba scikit-belajar , …

10 r  machine-learning  classification  random-forest  scikit-learn 

Dalam kode berikut ini saya melakukan regresi logistik pada data yang dikelompokkan menggunakan glm dan "dengan tangan" menggunakan mle2. Mengapa fungsi logLik di R memberi saya kemungkinan logLik (fit.glm) = - 2.336 yang berbeda dari yang logLik (fit.ml) = - 5.514 saya dapatkan dengan tangan? library(bbmle) #successes in first column, …

10 r  self-study  generalized-linear-model