Pertanyaan yang diberi tag «r»

Masalah klasifikasi dengan batas nonlinier tidak dapat diselesaikan dengan perceptron sederhana . Kode R berikut adalah untuk tujuan ilustrasi dan didasarkan pada contoh ini dengan Python): nonlin <- function(x, deriv = F) { if (deriv) x*(1-x) else 1/(1+exp(-x)) } X <- matrix(c(-3,1, -2,1, -1,1, 0,1, 1,1, 2,1, 3,1), ncol=2, byrow=T) …

10 r  neural-networks  svm  kernel-trick  perceptron 

Apakah ada paket R (atau bahkan fungsi R dasar) yang mengimplementasikan metode Fisher atau Stouffer untuk menggabungkan nilai-p? Pengkodean ini seharusnya hampir sepele, tetapi saya lebih suka menggunakan (dan mengutip) sebuah paket. Contoh kode dalam pertanyaan ini: Metode Fisher untuk menyisir nilai-p - bagaimana dengan ekor bawah?

10 r  combining-p-values 

Saya ingin membandingkan 2 sampel rata-rata untuk pengembalian stok 1 menit. Saya menganggap mereka terdistribusi Laplace (sudah diperiksa) dan saya membagi pengembalian menjadi 2 kelompok. Bagaimana saya dapat mengecek apakah keduanya berbeda secara signifikan? Saya pikir saya tidak bisa memperlakukan mereka seperti distribusi Normal, karena meskipun mereka lebih dari 300 …

10 r  confidence-interval  laplace-distribution 

Dalam lmerfungsi di lme4dalam Rada panggilan untuk membangun model matriks efek acak, , seperti yang dijelaskan di sini , halaman 7 - 9.ZZZ Menghitung mensyaratkan KhatriRao dan / atau produk Kronecker dari dua matriks, dan X_i . J i X iZZZJiJiJ_iXiXiX_i Matriks JiJiJ_i adalah suap: "Matriks indikator indeks faktor pengelompokan", …

10 r  mixed-model  lme4-nlme  matrix 

Saya telah mengikuti vektor X dan Y sederhana: > X [1] 1.000 0.063 0.031 0.012 0.005 0.000 > Y [1] 1.000 1.000 1.000 0.961 0.884 0.000 > > plot(X,Y) Saya ingin melakukan regresi menggunakan log X. Untuk menghindari log (0), saya mencoba untuk memberi +1 atau +0.1 atau +0.00001 atau …

10 r  regression  lognormal 

Bayangkan sebuah situasi: Kami memiliki catatan sejarah (20 tahun) dari tiga tambang. Apakah kehadiran perak meningkatkan kemungkinan menemukan emas di tahun depan? Bagaimana cara menguji pertanyaan seperti itu? Berikut ini contoh data: mine_A <- c("silver","rock","gold","gold","gold","gold","gold", "rock","rock","rock","rock","silver","rock","rock", "rock","rock","rock","silver","rock","rock") mine_B <- c("rock","rock","rock","rock","silver","rock","rock", "silver","gold","gold","gold","gold","gold","rock", "silver","rock","rock","rock","rock","rock") mine_C <- c("rock","rock","silver","rock","rock","rock","rock", "rock","silver","rock","rock","rock","rock","silver", "gold","gold","gold","gold","gold","gold") time <- seq(from …

10 r  time-series  hypothesis-testing  stochastic-processes 

Diberikan tiga variabel, ydan x, yang positif kontinu, dan z, yang kategorikal, saya memiliki dua model kandidat yang diberikan oleh: fit.me <- lmer( y ~ 1 + x + ( 1 + x | factor(z) ) ) dan fit.fe <- lm( y ~ 1 + x ) Saya berharap untuk …

10 r  hypothesis-testing  mixed-model  lme4-nlme 

Saya mencari untuk mengelompokkan / menggabungkan node dalam grafik menggunakan pengelompokan grafik di 'r'. Ini adalah variasi mainan yang menakjubkan dari masalah saya. Ada dua "cluster" Ada "jembatan" yang menghubungkan cluster Berikut ini adalah jaringan kandidat: Ketika saya melihat jarak koneksi, "hopcount", jika Anda mau, maka saya bisa mendapatkan matriks …

10 r  clustering  data-visualization  numerics 

Saya memiliki beberapa variabel dan saya tertarik untuk menemukan hubungan non-linear di antara mereka. Jadi saya memutuskan untuk memasang spline atau loess, dan mencetak plot yang bagus (lihat kode di bawah). Tapi, saya juga ingin memiliki beberapa statistik yang memberi saya ide seberapa besar kemungkinan hubungan itu adalah masalah keacakan …

10 r  regression  splines  loess 

Saya melakukan bootstrap dengan model campuran (beberapa variabel dengan interaksi dan satu variabel acak). Saya mendapat hasil ini (hanya sebagian): > boot_out ORDINARY NONPARAMETRIC BOOTSTRAP Call: boot(data = a001a1, statistic = bootReg, R = 1000) Bootstrap Statistics : original bias std. error t1* 4.887383e+01 -1.677061e+00 4.362948e-01 t2* 3.066825e+01 1.264024e+00 5.328387e-01 …

10 r  confidence-interval  bootstrap  bias 

Saya menjalankan regresi OLS yang dikumpulkan menggunakan paket plm di R. Meskipun, pertanyaan saya lebih lanjut tentang statistik dasar, jadi saya coba posting di sini dulu;) Karena hasil regresi saya menghasilkan residu heteroskedastik, saya ingin mencoba menggunakan kesalahan standar kuat heteroskedastisitas. Sebagai hasil dari coeftest(mod, vcov.=vcovHC(mod, type="HC0"))saya mendapatkan tabel yang …

10 r  anova  multiple-regression  heteroscedasticity  robust-standard-error 

Saya menjalankan glms di R (model linear umum). Saya pikir saya tahu nilai-nilai - sampai saya melihat bahwa memanggil ringkasan untuk GLM tidak memberi Anda perwakilan pvalue utama dari model secara keseluruhan - setidaknya tidak di tempat di mana model linear melakukannya. Saya bertanya-tanya apakah ini diberikan sebagai nilai untuk …

10 r  statistical-significance  generalized-linear-model  p-value  descriptive-statistics 

Saya bermain-main dengan dataset kanker payudara dan membuat sebar semua atribut untuk mendapatkan ide yang paling berpengaruh dalam memprediksi kelas malignant(biru) dari benign(merah). Saya mengerti bahwa baris mewakili sumbu x dan kolom mewakili sumbu y tetapi saya tidak dapat melihat pengamatan apa yang dapat saya lakukan tentang data atau atribut …

10 r  data-visualization  interpretation  scatterplot 

Semua literatur Pemodelan Distribusi Spesies menunjukkan bahwa ketika memprediksi ada / tidaknya suatu spesies menggunakan model yang menghasilkan probabilitas (misalnya, RandomForests), pilihan ambang batas probabilitas yang digunakan untuk benar-benar mengklasifikasikan suatu spesies sebagai ada atau tidaknya adalah penting dan kita harus tidak selalu mengandalkan default 0,5. Saya butuh bantuan dengan …

10 r  classification  random-forest  threshold 

Saya mengalami kesulitan memahami bagaimana varImpfungsi ini bekerja untuk model randomForest dengan caretpaket. Pada contoh di bawah ini, fitur var3 sama pentingnya dengan varImpfungsi caret , tetapi model final randomForest yang mendasarinya memiliki kepentingan nol untuk fitur var3. Mengapa demikian? require(randomForest) require(caret) rf <- train(x, y, method = "rf", trControl …

10 r  caret  random-forest