«r» etiketlenmiş sorular

MCMCglmmKarışık efektler modeli çalıştırmak için R paketini kullanarak performans sorunları yaşıyorum . Kod şöyle görünür: MC1<-MCMCglmm(bull~1,random=~school,data=dt,family="categorical" , prior=list(R=list(V=1,fix=1), G=list(G1=list(V=1, nu=0))) , slice=T, nitt=iter, ,burnin=burn, verbose=F) Verilerde yaklaşık 20.000 gözlem vardır ve bunlar yaklaşık 200 okulda kümelenmiştir. Kullanılmamış tüm değişkenleri veri çerçevesinden bıraktım ve çalıştırmadan önce diğer tüm nesneleri bellekten kaldırdım. …

12 r  mixed-model  mcmc  computational-statistics 

Aşağıdaki olasılık işlevim var: Prob = 11 + e- zprob=11+e-z\text{Prob} = \frac{1}{1 + e^{-z}} nerede z= B0+ B1X1+ ⋯ + BnXn.z=B0+B1X1+⋯+BnXn.z = B_0 + B_1X_1 + \dots + B_nX_n. Modelim şöyle görünüyor Pr ( Y= 1 ) = 11 + tecrübe( - [ - 3.92 + 0.014 × ( teklif …

12 r  probability  data-visualization  logistic 

Çeşitli fonksiyonel veri analizi yaklaşımlarını test etmeye çalışıyorum. İdeal olarak, simüle edilmiş fonksiyonel veriler üzerinde sahip olduğum yaklaşımlar panelini test etmek istiyorum. Bir Gauss gürültüsü (aşağıdaki kod) toplayan bir yaklaşım kullanarak simüle edilmiş FD oluşturmaya çalıştım, ancak ortaya çıkan eğriler gerçek şeye kıyasla çok sağlam görünüyor . Birisinin daha gerçekçi …

12 r  simulation  functional-data-analysis 

JJS dilinde belirtilen bir modelde MCMC çalıştırmak için rjags kullandım. Bu modeli çıkartmanın ve onunla tahminler yapmanın iyi bir yolu var mı (parametrelerimin posterior dağılımlarını kullanarak)? Modeli R'de yeniden belirtebilir ve parametre posterlerimin modlarını takabilirim; Sadece bunu yapmanın daha az gereksiz bir yolu olup olmadığını merak ediyorum. Http://sourceforge.net/p/mcmc-jags/discussion/610037/thread/0ecab41c aynı soruyu …

12 r  jags 

Bir veri kümesi verildi: x <- c(4.9958942,5.9730174,9.8642732,11.5609671,10.1178216,6.6279774,9.2441754,9.9419299,13.4710469,6.0601435,8.2095239,7.9456672,12.7039825,7.4197810,9.5928275,8.2267352,2.8314614,11.5653497,6.0828073,11.3926117,10.5403929,14.9751607,11.7647580,8.2867261,10.0291522,7.7132033,6.3337642,14.6066222,11.3436587,11.2717791,10.8818323,8.0320657,6.7354041,9.1871676,13.4381778,7.4353197,8.9210043,10.2010750,11.9442048,11.0081195,4.3369520,13.2562675,15.9945674,8.7528248,14.4948086,14.3577443,6.7438382,9.1434984,15.4599419,13.1424011,7.0481925,7.4823108,10.5743730,6.4166006,11.8225244,8.9388744,10.3698150,10.3965596,13.5226492,16.0069239,6.1139247,11.0838351,9.1659242,7.9896031,10.7282936,14.2666492,13.6478802,10.6248561,15.3834373,11.5096033,14.5806570,10.7648690,5.3407430,7.7535042,7.1942866,9.8867927,12.7413156,10.8127809,8.1726772,8.3965665) .. Parametrelerin bir tahminiyle en uygun olasılık dağılımını (gama, beta, normal, üstel, poisson, ki-kare, vb.) Belirlemek istiyorum. R kullanarak bir çözümün sağlandığı aşağıdaki bağlantıdaki sorunun zaten farkındayım: /programming/2661402/given-a-set-of-random-numbers-drawn-from-a- sürekli-tek-değişken-dağılım-f en iyi önerilen çözüm şudur: > library(MASS) > fitdistr(x, 't')$loglik #$ > fitdistr(x, 'normal')$loglik …

12 r  matlab  distributions 

Not: önceki sorumun yasal nedenlerle silinmesi gerektiğinden, bu soru bir gönderidir. Fonksiyonlu SAS PROC MIXED karşılaştırarak birlikte lmegelen nlmeR paketin, bazı çok kafa farklılıklar tökezledi. Daha spesifik olarak, farklı testlerdeki özgürlük dereceleri ve arasında farklılık gösterir PROC MIXEDve lmenedenini merak ettim. Aşağıdaki veri kümesinden başlayın (R kodu aşağıda verilmiştir): ind: …

12 r  mixed-model  sas  degrees-of-freedom  pdf  unbiased-estimator  distance-functions  functional-data-analysis  hellinger  time-series  outliers  c++  relative-risk  absolute-risk  rare-events  regression  t-test  multiple-regression  survival  teaching  multiple-regression  regression  self-study  t-distribution  machine-learning  recommender-system  self-study  binomial  standard-deviation  data-visualization  r  predictive-models  pearson-r  spearman-rho  r  regression  modeling  r  categorical-data  data-visualization  ggplot2  many-categories  machine-learning  cross-validation  weka  microarray  variance  sampling  monte-carlo  regression  cross-validation  model-selection  feature-selection  elastic-net  distance-functions  information-theory  r  regression  mixed-model  random-effects-model  fixed-effects-model  dataset  data-mining 

RTarafından oluşturulan veri bir örnek var rnorm(50,0,1), bu yüzden veri açıkça normal bir dağılım alır. Ancak Rverilerle ilgili bu dağıtım bilgisini "bilmez". RÖrneklemimin ne tür bir dağıtımdan geldiğini tahmin edebilecek bir yöntem var mı ? Değilse, shapiro.testişlevi kullanacağım ve bu şekilde ilerleyeceğim.

12 r  distributions 

Önceki bir yazıda EQ-5D puanlarıyla nasıl başa çıkacağımı merak ettim . Son zamanlarda sınırlı sonuçlarla başa çıkmak için zarif bir yol sunan Bottai ve McKeown tarafından önerilen lojistik kantil regresyon üzerine tökezledim . Formül basit: l o gben t ( y) = l o g( y- ym i nym bir …

12 r  logistic  data-visualization  logit  quantile-regression 

Ben kısmi en küçük kareler (PLS) çok yeni ve R işlevin çıktısını anlamaya çalışma plsr()içinde plspaketin. Verileri simüle edip PLS'yi çalıştıralım: library(pls) n <- 50 x1 <- rnorm(n); xx1 <- scale(x1) x2 <- rnorm(n); xx2 <- scale(x2) y <- x1 + x2 + rnorm(n,0,0.1); yy <- scale(y) p <- plsr(yy …

12 r  regression  partial-least-squares 

Biri bana R'de yenilik, frekans ve parasal değer (RFM) modellemesi ve müşteri değeri modellemesi yapmayı söyleyebilir mi? Ayrıca, birisi bana bazı literatürü önerebilir mi?

12 r  references  marketing 

Son zamanlarda çok seviyeli modelleme üzerine bir giriş dersi aldım. Kullandığımız veri kümelerinin ve örneklerin çoğu sosyal bilimlerdendi. Bir biyoistatistik departmanında 2 haftalık bir staj aldım, burada hem hastaneler arasında hem de 5 yıl boyunca yüksek mortalite oranına sahip bir acil durum için hasta düzeyindeki hastane sonuçlarındaki değişiklikle ilgili bir …

12 r  stata  multilevel-analysis  epidemiology 

Bir grubun diğerinden daha fazla temsil edildiği eğitim verilerinden bir SVM oluşturmaya çalışıyorum. Ancak, gruplar nihai test verilerinde eşit olarak temsil edilecektir. Bu nedenle, iki paketin eğitim verilerindeki etkisini dengelemek için R paketi arabiriminin class.weightsparametresini kullanmak istiyorum .e1071libsvm Bu ağırlıkların tam olarak nasıl belirlenmesi gerektiğinden emin olmadığım için küçük bir …

12 r  machine-learning  svm 

MC simülasyonunu kullanarak aşağıdaki integrali nasıl yaklaştırabilirim? ∫1−1∫1−1|x−y|dxdy∫−11∫−11|x−y|dxdy \int_{-1}^{1} \int_{-1}^{1} |x-y| \,\mathrm{d}x \,\mathrm{d}y Teşekkürler! Düzenleme (Bazı bağlamlar): Simülasyonun integralleri yaklaşık olarak nasıl kullanacağımı öğrenmeye çalışıyorum ve bazı zorluklarla karşılaştığımda bazı alıştırmalar yapıyorum. Edit 2 + 3 : Bir şekilde kafam karıştı ve integrali ayrı parçalara bölmem gerektiğini düşündüm. Yani, aslında …

12 r  self-study  monte-carlo 

Kilitli . Bu soru ve cevapları kilitlidir çünkü soru konu dışıdır, ancak tarihsel önemi vardır. Şu anda yeni yanıtları veya etkileşimleri kabul etmiyor. Üç sütundan oluşan bir ascii veri kümem var, ancak yalnızca son ikisi gerçek verilerdir. Şimdi kullanarak bir nokta grafiği oluşturmak istiyorum read.csv(file = "result1", sep= " "). …

12 r 

Araba hareket etmiyorken 1, araba hareket ettiğinde 0 olan bir ikili zaman serim var. 36 saat ileri ve her saat için bir zaman ufku tahmin etmek istiyorum. İlk yaklaşımım şu girişleri kullanarak bir Naive Bayes kullanmaktı: t-24 (günlük mevsimsel), t-48 (haftalık mevsimsel), günün saati. Ancak, sonuçlar çok iyi değil. Bu …

12 r  time-series  forecasting  binary-data