«r» etiketlenmiş sorular

R'de 2 taraflı Kolmogorov Smirnov testi için güç analizi yapmak mümkün müdür? İki ampirik dağılımın ks.test () kullanarak farklılık gösterip göstermediğini test ediyorum ve bir güç analizi eklemek istiyorum. R'deki KS testleri için yerleşik güç analizleri bulamadım. Herhangi bir öneriniz var mı? Düzenleme : Bunlar rasgele oluşturulan dağıtım verilerimi yakından …

10 r  power-analysis  kolmogorov-smirnov 

Bayes metodolojisi kullanarak AB testi sonuçlarını hesaplamak için Bayesian ab test formülünü kullanıyorum. Pr(pB>pA)=∑i=0αB−1B(αA+i,βB+βA)(βB+i)B(1+i,βB)B(αA,βA)Pr(pB>pA)=∑i=0αB−1B(αA+i,βB+βA)(βB+i)B(1+i,βB)B(αA,βA) \Pr(p_B > p_A) = \sum^{\alpha_B-1}_{i=0} \frac{B(\alpha_A+i,\beta_B+\beta_A)}{(\beta_B+i)B(1+i,\beta_B)B(\alpha_A, \beta_A)} nerede αAαA\alpha_A bir arada artı A için başarı sayısı βAβA\beta_A bir arada artı A için başarısızlık sayısı αBαB\alpha_B bir artı B için başarı sayısı βBβB\beta_B bir artı B için …

10 r  bayesian  ab-test 

Ben kullanarak bazı davranışsal verilerin analiz ediyorum lme4yılında Rçoğunlukla ardından Bodo Kış mükemmel öğreticiler , ama düzgün etkileşimleri hallediyorum ben anlamıyorum. Daha kötüsü, bu araştırmaya dahil olan hiç kimse karışık modeller kullanmaz, bu yüzden işlerin doğru olduğundan emin olmak konusunda biraz şaşırıyorum. Sadece yardım için bir çığlık atmak yerine, sorunu …

10 r  mixed-model  p-value  lme4-nlme 

R kullanarak zaman serisi analizi yapıyorum. Verilerimi trend, mevsimsel ve rastgele bileşenlere ayırmalıyım. 3 yıllık haftalık verilerim var. R - stl()ve 'de iki fonksiyon buldum decompose(). Bunun stl()çarpımsal ayrışma için iyi olmadığını okudum . Biri bana bu fonksiyonların hangi senaryoda kullanılabileceğini söyleyebilir mi?

10 r  time-series 

Robert Kabacoff en düşük Hızlı-R Ben # Bootstrap 95% CI for regression coefficients library(boot) # function to obtain regression weights bs <- function(formula, data, indices) { d <- data[indices,] # allows boot to select sample fit <- lm(formula, data=d) return(coef(fit)) } # bootstrapping with 1000 replications results <- boot(data=mtcars, statistic=bs, …

10 r  regression  p-value  bootstrap 

İki örnek oranlarının karşılaştırılması, örnek büyüklüğü tahmini: R ve Stata Örnek boyutları için aşağıdaki gibi farklı sonuçlar aldım: içinde R power.prop.test(p1 = 0.70, p2 = 0.85, power = 0.90, sig.level = 0.05) Sonuç: n = 160.7777n=160.7777n = 160.7777 (yani 161) her grup için. In Stata sampsi 0.70 0.85, power(0.90) alpha(0.05) …

10 r  stata  sample-size  power-analysis  jmp 

Ben onları denemek düşündüm bu yüzden vb MNIST, cifar, STL-10, Rastgele Ormanları uygulanması ile ilgili herhangi literatürü bulamadı permütasyon değişmeyen MNIST kendim. İçinde R , I güvenilir: randomForest(train$x, factor(train$y), test$x, factor(test$y), ntree=500) Bu 2 saat sürdü ve% 2.8 test hatası aldı. Ben de denedim scikit-öğrenme ile, RandomForestClassifier(n_estimators=2000, max_features="auto", max_depth=None) 70 …

10 r  machine-learning  classification  random-forest  scikit-learn 

Aşağıdaki kodda glm ve "elle" mle2 kullanarak gruplandırılmış veriler üzerinde lojistik regresyon gerçekleştiriyorum. Neden R'deki logLik işlevi bana bir günlük olasılığı verir logLik (fit.glm) = - 2.336 birinden farklı logLik (fit.ml) = - 5.514 Elle alıyorum? library(bbmle) #successes in first column, failures in second Y <- matrix(c(1,2,4,3,2,0),3,2) #predictor X <- …

10 r  self-study  generalized-linear-model 

SMOTE'u çok sınıflı sınıflandırma sorunumdaki dengesizliği düzeltmek için kullanmaya çalışıyorum. SMOTE, SMOTE yardım belgesine göre iris veri kümesinde mükemmel çalışmasına rağmen, benzer bir veri kümesinde çalışmaz. Verilerim şöyle görünüyor. 1, 2, 3 değerlerine sahip üç sınıfı olduğunu unutmayın. > data looking risk every status 1 0 1 0 1 2 …

10 r  classification  unbalanced-classes  oversampling 

Ben böyle modeller var: require(nlme) set.seed(123) n <- 100 k <- 5 cat <- as.factor(rep(1:k, n)) cat_i <- 1:k # intercept per kategorie x <- rep(1:n, each = k) sigma <- 0.2 alpha <- 0.001 y <- cat_i[cat] + alpha * x + rnorm(n*k, 0, sigma) plot(x, y) m1 <- …

10 r  mixed-model  aic 

R'de ayrık zamanlı bir model yerleştirmeye çalışıyorum, ancak nasıl yapılacağından emin değilim. Bağımlı değişkeni farklı satırlarda, her bir zaman gözlemi için bir tane düzenleyebileceğinizi ve glmbir logit veya cloglog bağlantısıyla işlevi kullanabileceğinizi okudum. Bu anlamda, üç sütun vardır: ID, Event(her zaman atıl 1 ya da 0) ve Time Elapsedek olarak, …

10 r  survival  pca  sas  matlab  neural-networks  r  logistic  spatial  spatial-interaction-model  r  time-series  econometrics  var  statistical-significance  t-test  cross-validation  sample-size  r  regression  optimization  least-squares  constrained-regression  nonparametric  ordinal-data  wilcoxon-signed-rank  references  neural-networks  jags  bugs  hierarchical-bayesian  gaussian-mixture  r  regression  svm  predictive-models  libsvm  scikit-learn  probability  self-study  stata  sample-size  spss  wilcoxon-mann-whitney  survey  ordinal-data  likert  group-differences  r  regression  anova  mathematical-statistics  normal-distribution  random-generation  truncation  repeated-measures  variance  variability  distributions  random-generation  uniform  regression  r  generalized-linear-model  goodness-of-fit  data-visualization  r  time-series  arima  autoregressive  confidence-interval  r  time-series  arima  autocorrelation  seasonality  hypothesis-testing  bayesian  frequentist  uninformative-prior  correlation  matlab  cross-correlation 

R'de iki yönlü tekrarlanan ölçümler ANOVA için hata terimlerini tanımlamakta sorun yaşıyorum Verilerim, bir ağaçtan çıkarılan bir çekirdek boyunca üç radyal konum (iç, orta ve dış) için ahşap yoğunluk tahminlerinden oluşuyor. Toplam 20 ağaç türü, her türden 6 kişi ve her ağaçtan iki çekirdek vardır. Radyal konumun ahşap yoğunluğu üzerindeki …

10 r  anova  repeated-measures 

En iyi modeli veya modelleri seçmek için AIC'yi kullanarak karışık model analizimi nasıl yürüteceğime dair literatürü araştırmaya çalıştıktan sonra bunalmış hissediyorum. Verilerimin bu kadar karmaşık olduğunu düşünmüyorum, ancak yaptığımın doğru olduğunu doğrulamak için arıyorum ve sonra nasıl ilerleyeceğiniz konusunda tavsiyede bulunuyorum. Lme veya lmer kullanmam gerekip gerekmediğinden emin değilim ve …

10 r  lme4-nlme  random-effects-model  aic 

Parametre dağılımı hakkında bir fikriniz yoksa ne olur? Hangi yaklaşımı kullanmalıyız? Çoğu zaman, belirli bir değişkenin belirli bir türün varlığı / yokluğu üzerinde herhangi bir etkisi olup olmadığını ve değişkenin değişken önemine göre kabul edilip edilmediğini anlamayı hedefleriz. Bu, çoğu zaman, bir parametrenin olması gereken hızlı dağıtım üzerinde düşünmediğimiz anlamına …

10 r  bayesian  mcmc  bugs  winbugs 

Sorularım mgcv R paketindeki GAM'lerle ilgilidir . Küçük bir örneklem büyüklüğü nedeniyle, bir kereye mahsus bırakma çapraz doğrulaması kullanarak tahmin hatasını belirlemek istiyorum. Bu makul mi? Bunu nasıl yapabilirim bir paket veya kod var mı? errorest()İşlevi IPRED paketin çalışmaz. Basit bir test veri kümesi: library(mgcv) set.seed(0) dat <- gamSim(1,n=400,dist="normal",scale=2) b<-gam(y~s(x0)+s(x1)+s(x2)+s(x3),data=dat) …

10 r  cross-validation  gam  mgcv