Pertanyaan yang diberi tag «r»

R in Action (Kabacoff, 2011) menyarankan rutinitas berikut untuk menguji penyebaran berlebihan dalam regresi logistik: Fit regresi logistik menggunakan distribusi binomial: model_binom <- glm(Species=="versicolor" ~ Sepal.Width, family=binomial(), data=iris) Fit regresi logistik menggunakan distribusi quasibinomial: model_overdispersed <- glm(Species=="versicolor" ~ Sepal.Width, family=quasibinomial(), data=iris) Gunakan chi-squared untuk menguji penyebaran berlebihan: pchisq(summary(model_overdispersed)$dispersion * model_binom$df.residual, …

9 r  regression  distributions  logistic  overdispersion 

Saya melakukan regresi linier menggunakan fungsi Rm: x = log(errors) plot(x,y) lm.result = lm(formula = y ~ x) abline(lm.result, col="blue") # showing the "fit" in blue tapi itu tidak pas. Sayangnya saya tidak bisa memahami manualnya. Bisakah seseorang mengarahkan saya ke arah yang benar agar lebih cocok? Maksud saya, saya …

9 r  regression 

Apakah mungkin untuk memplot plot dependensi parsial untuk menampilkan probabilitas kelas dan memperkirakan efek prediktor untuk model GBM ? Sesuatu yang mirip dengan partialPlotdari randomForestpaket. Menurut artikel ini , sebagian plot dapat dilakukan dengan gbm. Terima kasih sebelumnya atas bantuan Anda.

9 r  gbm 

Saya memiliki data kecelakaan kendaraan bermotor menurut jam sehari. Seperti yang Anda harapkan, mereka tinggi di tengah hari dan puncaknya pada jam sibuk. Geom_density default ggplot2 memuluskannya dengan baik Subset dari data, untuk crash terkait drive-minum, tinggi di kedua ujung hari (malam hari dan pagi hari) dan tertinggi di ekstrem. …

9 r  pdf  kernel-smoothing  splines 

Katakanlah kita memiliki data berikut: set.seed(123) data <- data.frame(x = c(rnorm(50, 1, 1), rnorm(50, 5, 2)), y = c(rep('A', 50), rep('B', 50))) Yang menghasilkan boxplot berikut ( boxplot(data$x ~ data$y)): Sekarang katakanlah saya ingin menguji apakah dua sampel memiliki parameter lokasi yang sama (median dan / atau rata-rata). Dalam kasus …

9 r  nonparametric  wilcoxon-mann-whitney 

Saya bukan ahli, jadi maafkan saya jika beberapa terminologi sedikit canggung. Senang memberikan informasi lebih lanjut di mana diperlukan. Saya memiliki dua vektor dari 50 nilai numerik berpasangan di R. Saya ingin melakukan uji pengacakan atau permutasi dua sisi untuk menentukan apakah perbedaannya disebabkan oleh kebetulan atau tidak. Tes permutasi …

9 r  statistical-significance  permutation-test 

Saya bertanya-tanya apakah ada cara yang baik untuk menghitung kriteria pengelompokan berdasarkan formula BIC, untuk keluaran k-means dalam R? Saya agak bingung bagaimana cara menghitung BIC sehingga saya bisa membandingkannya dengan model clustering lainnya. Saat ini saya menggunakan implementasi paket statistik dari k-means.

9 r  clustering  k-means  bic 

Saya memproduksi skrip untuk membuat sampel bootstrap dari catsdataset (dari -MASS-paket). Mengikuti buku teks Davidson dan Hinkley [1] Saya menjalankan regresi linier sederhana dan mengadopsi prosedur dasar non-parametrik untuk bootstrap dari pengamatan awal, yaitu pemasangan kembali pasangan . Sampel asli dalam bentuk: Bwt Hwt 2.0 7.0 2.1 7.2 ... 1.9 …

9 r  bootstrap  fixed-effects-model 

Saya memiliki sejumlah besar prediktor (lebih dari 43.000) untuk memprediksi variabel dependen yang dapat mengambil 2 nilai (0 atau 1). Jumlah pengamatan lebih dari 45.000. Sebagian besar prediktor adalah unigrams, bigrams, dan trigram kata-kata, sehingga ada tingkat kolinearitas yang tinggi di antara mereka. Ada banyak sparsity di dataset saya juga. …

9 r  multiple-regression  lasso  glmnet 

Menggunakan plot.rqdalam quantregpaket di R, kita dapat memplot distribusi estimasi koefisien, dan mendapatkan sesuatu seperti ini: Apa garis merah putus-putus? Googling ekstensif telah mengungkapkan bahwa yang di tengah adalah rata-rata dari semua 99 perkiraan, tetapi masih tidak tahu tentang garis merah putus-putus.

9 r  data-visualization  quantile-regression 

Saya menggunakan rank(a, ties.method="max")peringkat a. Tetapi saya tidak begitu yakin apa yang ties.method="max"dilakukannya. Bisakah Anda membantu?

9 r  ranks  ties 

Saya ingin memperoleh interval kepercayaan 95% untuk centroid berdasarkan pada kesamaan Gower antara beberapa sampel mulivariat (data komunitas dari inti sedimen). Saya sejauh ini menggunakan vegan{}paket dalam R untuk mendapatkan kesamaan Gower yang dimodifikasi antara core (berdasarkan Anderson 2006; sekarang termasuk dalam R sebagai bagian dari vegdist()). Adakah yang tahu …

9 r  confidence-interval 

Saya memiliki dataset medis dengan sekitar 200 variabel. Salah satu variabel adalah penanda-bio (konsentrasi enzim tertentu). Distribusinya condong ke kanan, dan masalahnya adalah bahwa nilai di atas level tertentu disensor / terputus pada level itu. Jadi, sementara rata-rata variabel sekitar 10, nilai apa pun yang lebih besar dari 50 dicatat …

9 r  epidemiology  data-imputation  censoring 

Saya bertanya-tanya apakah ada yang punya pengalaman menggunakan fungsi mouse, seperti yang dijelaskan dalam mouse: Imputasi Multivarian oleh Chained Equations di R (JSS 2011 45 (3))? Saya memiliki dataset dengan sejumlah variabel, masing-masing dengan berbagai tingkat data yang hilang. Pertanyaan utama saya adalah: katakanlah saya menggunakan regresi linier Bayesian untuk …

9 r  regression  bayesian  data-imputation  mice 

Kami memiliki kumpulan data dengan dua kovariat dan variabel pengelompokan kategoris dan ingin tahu apakah ada perbedaan yang signifikan antara kemiringan atau intersep di antara kovariat yang terkait dengan variabel pengelompokan yang berbeda. Kami telah menggunakan anova () dan lm () untuk membandingkan kecocokan dari tiga model yang berbeda: 1) …

9 r  regression  anova  confidence-interval