«r» etiketlenmiş sorular

2 ilacın ( drug1, drug2) hastanın düşme olasılığı üzerindeki etkisini tahmin etmeye çalışıyorum ( event). Hastalar bir kereden fazla düşebilir ve herhangi bir noktada ilaçlara konabilir veya çıkarılabilir. Benim sorum, verinin zaman dilimi (günler), özellikle de günler arasında çakışma olması gerekip gerekmediğine göre nasıl yapılandırılması gerektiğidir. Yapımın yanlış olduğunu düşündüğüm …

9 r  survival  cox-model 

Dixon, Coles ( 1997 ) 'de, futboldaki skorları modellemek için (4.3)' te değiştirilmiş iki bağımsız Poisson modeli için maksimum olabilirlik tahminini kullanmışlardır. Alfa ve beta yanı sıra ev etkisi parametreleri (s. 274, Tablo 4) herhangi bir paket kullanmadan (her zamanki gibi bağımsız Poisson modelleri kullanarak çok iyi) "çoğaltmak" için R …

9 r  modeling  maximum-likelihood  games 

Bir lojistik regresyon modelim var (elastik net regülasyonlu R'de glmnet ile uyumlu) ve gerçek pozitifler ve yanlış pozitifler arasındaki farkı en üst düzeye çıkarmak istiyorum. Bunu yapmak için aşağıdaki prosedür akla geldi: Standart lojistik regresyon modeline uygun Tahmin eşiğini 0,5 olarak kullanarak, tüm pozitif tahminleri belirleyin Olumlu tahmin edilen gözlemler …

9 r  regression  logistic  classification  glmnet 

77 çeyrek boyunca 33 kişiyle panel verilerine dayalı vektör otomatik regresyon (VAR) ve dürtü yanıt işlevi (IRF) tahmini üzerinde çalışıyorum. Bu tür bir durum nasıl analiz edilmelidir? Bu amaçla hangi algoritmalar var? Bu analizleri R'de yapmayı tercih ederim, bu yüzden herhangi biri R kodu veya bu amaç için tasarlanmış bir …

9 r  econometrics  panel-data  var  impulse-response 

Belirli bir uyumlulaştırma modeliyle tutarlı verileri simüle etmek için bir prosedür bulmakla ilgileniyorum. Barron ve Kenny'nin (1986) ilk önce özetlediği ve Judd, Yzerbyt ve Muller (2013) gibi başka yerlerde açıklanan arabuluculuk modellerinin test edilmesine yönelik genel doğrusal yapısal denklem modeli çerçevesine göre , sonucu için arabuluculuk modelleri , arabulucu \ …

9 r  regression  simulation  random-generation  mediation 

İntertidal topluluklarla ilgili bir veri setini analiz ediyorum. Veriler kuadratlarda yüzde örtüsü (deniz yosunu, kıskaç, midye, vb.) Tür sayısı açısından yazışma analizi (CA) ve doğrusal çevresel (tür değil) eğilimler için daha yararlı bir şey olarak temel bileşen analizi (PCA) hakkında düşünmeye alışkınım . PCA veya CA'nın yüzde kapak için daha …

9 r  pca  multivariate-analysis  correspondence-analysis 

Burada aşağıdaki veriler var . Hidrokarbon yüzdesi 1.0 olduğunda ortalama saflıkta% 95 güven aralığını hesaplamaya çalışıyorum. R harfine şunu giriyorum. > predict(purity.lm, newdata=list(hydro=1.0), interval="confidence", level=.95) fit lwr upr 1 89.66431 87.51017 91.81845 Ancak, bu sonucu kendim nasıl elde edebilirim? Aşağıdaki denklemi kullanmaya çalıştım. sn e w=s2( 1 +1N-+(xn e w-x¯)2∑ …

9 r  regression  confidence-interval  prediction-interval 

Aşağıdaki gibi bir örnek veri kümesi var: Volume <- seq(1,20,0.1) var1 <- 100 x2 <- 1000000 x3 <- 30 x4 = sqrt(x2/pi) H = x3 - Volume r = (x4*H)/(H + Volume) Power = (var1*x2)/(100*(pi*Volume/3)*(x4*x4 + x4*r + r*r)) Power <- jitter(Power, factor = 1, amount = 0.1) plot(Volume,Power) Şekilde, …

9 r  hypothesis-testing  nonlinear-regression  goodness-of-fit  heteroscedasticity 

Bu belki de aptalca bir soru olabilir, ancak caret ile bir model oluştururken LOOCVveya veya (hatta daha da fazlası) gibi bir şey kullanırken LGOCV, aslında çapraz doğrulama adımının bu olması durumunda verileri tren ve test setlerine bölmenin yararı nedir? neyse? İlgili soruların bazılarını okudum ve çapraz onaylama yöntemlerinden bazılarının (örneğin …

9 r  machine-learning  cross-validation  caret 

Metafor paketini kullanarak R etki boyutlarının d meta-analizini yapıyorum . d , hastalar ile sağlıklı arasındaki bellek puanlarındaki farklılıkları temsil eder. Bununla birlikte, bazı çalışmalar sadece ilgili d ölçümünün alt skorlarını rapor etmektedir (örneğin, birkaç farklı bellek skoru veya üç ayrı bellek testi bloğundan alınan skorlar). Lütfen çalışmaların etki boyutlarını …

9 r  meta-analysis  effect-size  meta-regression 

Monte Carlo simülasyonunun chisq.test()R'deki fonksiyonda kullanımını anlamak istiyorum . 128 seviye / sınıfa sahip nitel bir değişkenim var. Örneklem boyutum 26 (daha fazla "birey" örnekleyemedim). Açıkçası, ben 0 "bireyler" ile bazı seviyeleri olacak. Ancak gerçek şu ki, mümkün olan 127 ders arasından sadece çok az sayıda sınıfım var. Ki-kare testi …

9 r  chi-squared  monte-carlo 

Arka fon Epstein'dan (1991) bir makalede: Aylık ortalamalardan günlük klimatolojik değerler elde edildiğinde, periyodik ve eşit aralıklı değerler için Fourier enterpolasyonunun hesaplanması için formülasyon ve bir algoritma verilmektedir. Bu makalede amaç, enterpolasyon yoluyla aylık ortalama değerlerden günlük değerler elde etmektir . Kısacası, bilinmeyen günlük değerlerin harmonik bileşenlerin toplamı ile temsil …

9 r  interpolation  fourier-transform 

randomForestYeni bir vakanın sınıfını öngören bir uygulamada kullanmak istediğim iyi bir sınıflandırma modelim var. Yeni vaka kaçınılmaz olarak eksik değerlere sahip. Tahmin NA'lar için böyle çalışmaz. Bunu nasıl yapmalıyım? data(iris) # create first the new case with missing values na.row<-45 na.col<-c(3,5) case.na<-iris[na.row,] case.na[,na.col]<-NA iris.rf <- randomForest(Species ~ ., data=iris[-na.row,]) # …

9 r  random-forest  prediction  missing-data 

Durum Bir bağımlı ve bir bağımsız değişken ile bir veri kümesi var . Birlikte regresyon doğrusal sürekli parçalı uygun istediğiniz meydana gelen bilinen sabit / kesme noktaları . Kesme noktaları belirsizlik olmadan bilinir, bu yüzden onları tahmin etmek istemiyorum. Sonra İşte bir örnekyyyxxxkkk(a1,a2,…,ak)(a1,a2,…,ak)(a_{1}, a_{2}, \ldots, a_{k})yi=β0+β1xi+β2max(xi−a1,0)+β3max(xi−a2,0)+…+βk+1max(xi−ak,0)+ϵiyi=β0+β1xi+β2max⁡(xi−a1,0)+β3max⁡(xi−a2,0)+…+βk+1max⁡(xi−ak,0)+ϵi y_{i} = \beta_{0} + …

9 r  regression  standard-error  piecewise-linear 

EDIT: Bu yazı yaptığımdan beri, burada ek bir yazı ile takip ettim . Aşağıdaki metnin özeti: Bir model üzerinde çalışıyorum ve doğrusal regresyon, Box Cox dönüşümleri ve GAM'ı denedim, ancak fazla ilerleme kaydetmedim Kullanarak R, şu anda büyük lig (MLB) düzeyinde küçük lig beyzbol oyuncularının başarısını tahmin etmek için bir …

9 r  multiple-regression  data-transformation  nonlinear-regression  gam