Pertanyaan yang diberi tag «r»

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 ~ our_bid, family = binomial(link = "logit"), data = mydat) Deviance Residuals: Min 1Q …

9 r  logistic  generalized-linear-model 

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 meletakkan a, b, c, dll di atas kotak yang …

9 r  statistical-significance  boxplot 

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 cocok …

9 r  mixed-model  lme4-nlme 

Saya memiliki tipe data berikut. Saya telah mengevaluasi 10 orang yang masing-masing diulang 10 kali. Saya memiliki 10x10 matriks relasi (hubungan antara semua kombinasi individu). set.seed(1234) mydata <- data.frame (gen = factor(rep(1:10, each = 10)), repl = factor(rep(1:10, 10)), yld = rnorm(10, 5, 0.5)) Gen ini adalah varietas tanaman yang …

9 r  mixed-model 

Membangun posting Bagaimana mengelola proyek analisis statistik dan ProjectTemplatepaket dalam R ... T: Bagaimana Anda membangun struktur direktori proyek statistik Anda ketika banyak bahasa banyak fitur (misalnya, R DAN Splus)? Sebagian besar diskusi tentang topik ini terbatas pada proyek yang terutama menggunakan satu bahasa. Saya prihatin dengan cara meminimalkan kecerobohan, …

9 r  project-management  splus 

Saya pikir ini adalah pertanyaan yang agak mendasar, tetapi saya baru sadar bahwa saya tidak begitu mengerti istilah koreksi kontinuitas. Saya menggunakan R dan menemukan sintaks yang sama correct=TRUEuntuk keduanya chisq.testdan mcnemar.test. Apakah mereka mengacu pada metode koreksi kontinuitas yang berbeda? Saya telah mendengar tentang kesinambungan yate untuk uji pearson …

9 r  chi-squared  yates-correction 

Latar belakang : Israel (dan timur tengah pada umumnya) dipenuhi dengan protes. Saya ingin tahu, ketika diberi gambar, untuk memperkirakan berapa banyak orang dalam gambar itu (sering kali gambar kerumunan besar). Pemodelan apa yang dapat menawarkan beberapa solusi untuk masalah ini? (dan, tentu saja, dapatkah dilakukan dengan paket open source. …

8 r  image-processing  pattern-recognition 

Saya membantu seorang rekan untuk mem - bootstrap model efek-campuran meta-analisis menggunakan kerangka paket metafor R yang ditulis oleh @Wolfgang. Menarik dan mengkhawatirkan, untuk salah satu koefisien model saya mendapatkan distribusi bimodal ketika bootstrap (lihat panel kanan bawah pada gambar di bawah). Saya kira salah satu penyebab utama bisa menjadi …

8 r  mixed-model  bootstrap  meta-analysis  bimodal 

Pertanyaan berikut didasarkan pada diskusi yang ditemukan di halaman ini . Diberikan variabel respons y, variabel penjelas kontinu xdan faktor fac, dimungkinkan untuk mendefinisikan Model Aditif Umum (GAM) dengan interaksi antara xdan facmenggunakan argumen by=. Menurut file bantuan ?gam.models dalam paket R mgcv, ini dapat dilakukan sebagai berikut: gam1 <- …

8 r  interaction  gam  mgcv 

Saya mengalami kesulitan untuk masuk ke halaman yang sama dengan penyelia saya ketika datang untuk memvalidasi model saya. Saya telah menganalisis residu (diamati terhadap nilai yang dipasang) dan saya menggunakan ini sebagai argumen untuk membahas hasil yang diperoleh oleh model saya, namun penyelia saya bersikeras bahwa satu-satunya cara untuk memvalidasi …

8 r  references  glmm  validation  glmmtmb 

Saya memperkirakan beberapa tindakan berulang model Fixed Effects, dengan komponen kesalahan bersarang, berdasarkan pada variabel pengelompokan, yaitu model non-bersarang, menggunakan plm. Saya sekarang tertarik uji apakah model lengkap berbeda secara signifikan, yaitu mana adalah model lengkap untuk dan adalah model lengkap untuk danHHai:βFe m a l e=βM.a l eHHai:βFemSebuahle=βM.SebuahleH_o: \beta_{Female} …

8 r  hypothesis-testing  repeated-measures  nested-data  plm 

Saya memiliki dataset dengan 5.818.446 baris dan 51 kolom, di mana 50 di antaranya adalah prediktor. Respons saya kuantitatif, jadi saya tertarik dengan model regresi. Saya mencoba menyesuaikan hutan acak dengan data saya menggunakan paket tanda sisipan. Namun, saya tidak memiliki cukup RAM untuk melakukannya. Saya sudah mencari solusi untuk …

8 r  random-forest  large-data 

Saya melihat bagaimana jarak Euclidean minimum yang diharapkan antara titik-titik seragam acak dan titik asal berubah saat kami meningkatkan kepadatan titik-titik acak ( titik per unit persegi ) di sekitar titik asal. Saya telah berhasil menemukan hubungan antara keduanya yang digambarkan sebagai berikut: Expected Min Distance=12Density−−−−−−√Expected Min Distance=12Massa jenis\text{Expected Min …

8 r  expected-value  monte-carlo  uniform  minimum 

Saya punya dataset dengan 338 prediktor dan 570 instance (sayangnya tidak dapat diunggah) di mana saya menggunakan Lasso untuk melakukan pemilihan fitur. Secara khusus, saya menggunakan cv.glmnetfungsi dari glmnetsebagai berikut, di mana mydata_matrixmatriks biner 570 x 339 dan hasilnya adalah biner: library(glmnet) x_dat <- mydata_matrix[, -ncol(mydata_matrix)] y <- mydata_matrix[, ncol(mydata_matrix)] …

8 r  classification  lasso  glmnet 

Kita dapat mensimulasikan regresi linier tanpa keacakan, yang berarti kita membuat daripada . Maka jika kita cocok dengan model linier koefisien akan identik dengan "kebenaran dasar". Berikut ini sebuah contoh.y=Xβy=Xβy=X\betay=Xβ+ϵy=Xβ+ϵy=X\beta+\epsilon set.seed(0) n <- 1e5 p <- 3 X <- matrix(rnorm(n*p), ncol=p) beta <- runif(p) # y <- X %*% beta …

8 r  logistic  regularization  linear  separation