Pertanyaan yang diberi tag «r»

Latar Belakang Saya memiliki data dari studi lapangan di mana ada empat tingkat pengobatan dan enam ulangan di masing-masing dua blok. (4x6x2 = 48 pengamatan) Blok berjarak sekitar 1 mil, dan di dalam blok, ada grid 42, 2m x 4m plot dan jalan setapak lebar 1m; studi saya hanya menggunakan …

10 r  spatial  linear-model  covariance 

Saya ingin mengetahui nilai-nilai yang (x, y)digunakan dalam merencanakan plot(b, seWithMean=TRUE)di mgcv paket. Adakah yang tahu bagaimana saya bisa mengekstraksi atau menghitung nilai-nilai ini? Berikut ini sebuah contoh: library(mgcv) set.seed(0) dat <- gamSim(1, n=400, dist="normal", scale=2) b <- gam(y~s(x0), data=dat) plot(b, seWithMean=TRUE)

10 r  time-series  smoothing  mgcv 

Pertanyaan: Dengan rantai MCMC 10 dimensi, katakanlah saya siap memberikan Anda sebuah matriks undian: 100.000 iterasi (baris) dengan 10 parameter (kolom), bagaimana cara terbaik saya mengidentifikasi mode posterior? Saya terutama memperhatikan beberapa mode. Latar Belakang:Saya menganggap diri saya seorang ahli statistik yang mengerti secara komputasi, tetapi ketika seorang kolega bertanya …

10 r  bayesian  mcmc  k-nearest-neighbour  mode 

Buku statistik yang saya baca merekomendasikan omega kuadrat untuk mengukur efek percobaan saya. Saya telah membuktikan menggunakan desain plot terpisah (campuran antara subyek dan desain antara subyek) bahwa faktor dalam subyek saya signifikan secara statistik dengan p <0,001 dan F = 17. Sekarang saya ingin melihat seberapa besar perbedaannya ... …

10 r  anova  effect-size  split-plot 

Saya mencoba menafsirkan tipe model logistik berikut: mdl <- glm(c(suc,fail) ~ fac1 + fac2, data=df, family=binomial) Apakah output dari predict(mdl)peluang keberhasilan yang diharapkan untuk setiap titik data? Apakah ada cara sederhana untuk menabulasi peluang untuk setiap tingkat faktor model, daripada semua titik data?

10 r  logistic  generalized-linear-model 

Latar Belakang Jadi pertama-tama latar belakang untuk mengukur tingkat pemahaman yang mungkin saya miliki. Saat ini menyelesaikan tesis MSc, statistik telah menjadi bagian yang dapat diabaikan dari ini meskipun saya memiliki pemahaman dasar. Pertanyaan saya saat ini membuat saya ragu apa yang bisa / harus saya lakukan dalam praktik, membaca …

10 r  time-series  variance  descriptive-statistics  structural-change 

Dari apa yang saya baca, antara lain di situs statistik UCLA, kelompok konsultasi tes rasio kemungkinan dan tes wald cukup mirip dalam menguji apakah dua model glm menunjukkan perbedaan yang signifikan dalam kesesuaian untuk dataset (permisi jika kata-kata saya mungkin agak off). Intinya saya dapat membandingkan dua model dan menguji …

10 r  regression  likelihood-ratio 

Saya menguji kointegrasi menggunakan tes Johansen. Saya telah melihat pertanyaan seperti bagaimana menafsirkan hasil tes, tetapi ketika saya menafsirkan milik saya, saya memiliki beberapa keraguan. Dalam hasil saya r = 3sejak itu 4.10 < 10.49, jadi saya tidak bisa membentuk seri stasioner. Ini adalah sama untuk r = 2 dan …

10 r  hypothesis-testing  cointegration 

Tutup. Pertanyaan ini di luar topik . Saat ini tidak menerima jawaban. Ingin meningkatkan pertanyaan ini? Perbarui pertanyaan sehingga sesuai topik untuk Cross Validated. Ditutup 4 tahun yang lalu . Saya melatih model menggunakan caretpaket dalam R selama hampir 3 hari. Perhitungan berjalan secara paralel (beberapa proses). Sayangnya tidak ada …

10 r  machine-learning  svm  caret 

Bagaimana saya bisa mendapatkan efek acak dikumpulkan untuk lmer setelah beberapa kali imputasi? Saya menggunakan mouse untuk menghubungkan beberapa frame data. Dan lme4 untuk model campuran dengan intersep acak dan kemiringan acak. Pooling lmer baik-baik saja, kecuali bahwa itu tidak mengumpulkan efek acak. Saya telah mencari banyak solusi tanpa keberuntungan. …

10 r  lme4-nlme  multiple-imputation 

Saya menggunakan ccfuntuk menemukan korelasi antara 2 seri waktu. Saya mendapatkan plot yang terlihat seperti itu: Perhatikan bahwa saya terutama tertarik pada korelasi untuk lag = 0. Pertanyaan: Apakah mengartikannya dengan benar bahwa ada korelasi silang untuk lag = 0, karena untuk lag ini korelasi silang berada di atas garis …

10 r  interpretation  cross-correlation 

Gambar di bawah ini (Gambar 1 dari hal. 646 tulisan ini ) membandingkan nilai yang diamati dengan nilai yang diharapkan di bawah distribusi Poisson. Kemudian menjalankan uji chi-squared untuk melihat apakah nilai yang diamati berbeda dari nilai yang diharapkan di bawah distribusi Poisson. Dengan menggunakan R, bagaimana mungkin untuk menghasilkan …

10 r  distributions  chi-squared  poisson-distribution 

Saya punya sampel yang merupakan vektor dengan 220 angka. Berikut ini tautan ke histogram data saya. . Dan saya ingin memeriksa apakah data saya cocok dengan distribusi Pareto, tetapi saya tidak ingin melihat plot QQ dengan distribusi itu, tetapi saya membutuhkan jawaban yang tepat dengan nilai-p dalam R, seperti uji …

10 r  distributions  pareto-distribution 

Saya ingin menghasilkan data acak dari distribusi normal terbatas menggunakan R. Sebagai contoh saya mungkin ingin mensimulasikan variabel dari distribusi normal dengan mean=3, sd= 2dan nilai apa pun yang lebih besar dari 5 diresampled dari distribusi normal yang sama. Jadi, untuk fungsi umum, saya bisa melakukan hal berikut. rnorm(n=100, mean=3, …

10 r  normal-distribution  simulation  truncation 

Adakah yang tahu fungsi atau paket R yang dapat membantu saya mengubah skor z menjadi skor persentil? Tujuan akhir adalah untuk mengklasifikasikan atau memberi peringkat sekelompok responden ke dalam empat kategori berdasarkan ketinggian skor z mereka (skor terendah 20%, 30%, 30%, 20% skor tertinggi). Ada ide? Banyak terima kasih sebelumnya!

10 r  quantiles  change-scores