«r» 태그된 질문

R에서 양면 Kolmogorov Smirnov 테스트를위한 전력 분석을 수행 할 수 있습니까? ks.test ()를 사용하여 두 개의 경험적 분포가 다른지 테스트하고 전력 분석을 추가하려고합니다. R에서 KS 테스트에 대한 기본 제공 전력 분석을 찾을 수 없습니다. 제안 사항이 있습니까? 편집 : 이들은 내 데이터와 거의 비슷한 무작위로 생성 된 분포입니다 (지수 분포에 …

10 r  power-analysis  kolmogorov-smirnov 

Bayesian 방법론을 사용하여 AB 테스트의 결과를 계산하기 위해 Bayesian ab 테스트 의 공식을 사용하고 있습니다. 홍보 (피비>피ㅏ) =∑나는 = 0α비− 1비 (αㅏ+ 나는 ,β비+βㅏ)(β비+ i ) B ( 1 + 나는 ,β비) B (αㅏ,βㅏ)홍보(피비>피ㅏ)=∑나는=0α비−1비(αㅏ+나는,β비+βㅏ)(β비+나는)비(1+나는,β비)비(αㅏ,βㅏ) \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)} 어디 αㅏαㅏ\alpha_A 에 A의 성공 횟수를 더한 것 βㅏβㅏ\beta_A하나의 …

10 r  bayesian  ab-test 

나는 주로 Bodo Winter의 훌륭한 자습서를 따라 lme4in을 사용하여 일부 행동 데이터를 분석하고 있지만 상호 작용을 올바르게 처리하고 있는지 이해할 수 없습니다. 더군다나,이 연구에 참여한 어느 누구도 혼합 모델을 사용하지 않기 때문에 문제가 옳다는 것을 확신 할 때 약간 표류합니다.R 도움을 요청하는 것보다는 문제를 해석하기 위해 최선을 다한 다음에 집단적인 …

10 r  mixed-model  p-value  lme4-nlme 

R을 사용하여 시계열 분석을 수행하고 있습니다. 데이터를 추세, 계절 및 랜덤 구성 요소로 분해해야합니다. 3 년 동안 매주 데이터가 있습니다. 나는 R의 두 가지 기능을 발견했다 - stl()와 decompose(). 나는 stl()곱셈 분해에 좋지 않다는 것을 읽었습니다 . 이 기능을 사용할 수있는 시나리오를 아는 사람이 있습니까?

10 r  time-series 

로버트 Kabacoff의에서 빠른-R 내가 가진 # 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, R=1000, …

10 r  regression  p-value  bootstrap 

비율의 2- 표본 비교, 표본 크기 추정 : R vs Stata 다음과 같이 샘플 크기에 대해 다른 결과를 얻었습니다. 에서 R power.prop.test(p1 = 0.70, p2 = 0.85, power = 0.90, sig.level = 0.05) 결과 : 각 그룹에 대해 (그래서 161)입니다.n = 160.7777엔=160.7777n = 160.7777 에서 STATA sampsi 0.70 0.85, power(0.90) …

10 r  stata  sample-size  power-analysis  jmp 

나는 MNIST, CIFAR, STL-10 등에 Random Forests를 적용하는 것에 대한 문헌을 찾지 못했기 때문에 순열 불변의 MNIST로 직접 시도 할 것이라고 생각했습니다 . 에서 R , 나는 시도했다 : randomForest(train$x, factor(train$y), test$x, factor(test$y), ntree=500) 2 시간 동안 실행되었으며 2.8 % 테스트 오류가 발생했습니다. 나는 또한 시도 scikit가 배울 와 함께, …

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

다음 코드에서는 glm을 사용하고 mle2를 사용하여 "손으로"그룹화 된 데이터에 대해 로지스틱 회귀를 수행합니다. R의 logLik 함수가 왜 logLik (fit.ml) =-5.514와 다른 로그 가능성 logLik (fit.glm) =-2.336을 제공합니까? library(bbmle) #successes in first column, failures in second Y <- matrix(c(1,2,4,3,2,0),3,2) #predictor X <- c(0,1,2) #use glm fit.glm <- glm(Y ~ X,family=binomial (link=logit)) …

10 r  self-study  generalized-linear-model 

SMOTE를 사용하여 멀티 클래스 분류 문제의 불균형을 수정하려고합니다. SMOTE 도움말 문서에 따라 SMOTE가 홍채 데이터 세트에서 완벽하게 작동하지만 유사한 데이터 세트에서는 작동하지 않습니다. 내 데이터 모양은 다음과 같습니다. 값이 1, 2, 3 인 3 개의 클래스가 있습니다. > data looking risk every status 1 0 1 0 1 2 0 …

10 r  classification  unbalanced-classes  oversampling 

나는 이런 모델을 가지고있다 : 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에 이산 시간 모델을 맞추려고하지만 어떻게 해야할지 모르겠습니다. 종속 변수를 각 시간 관찰마다 하나씩 다른 행 glm으로 구성하고 logit 또는 cloglog 링크와 함께 함수를 사용할 수 있다는 것을 읽었습니다. 이런 의미에서, 나는 세 개의 열이 있습니다 : ID, Event(각 시간 경과시 1 또는 0) 및 Time Elapsed(관측 시작부터 ) 그리고 …

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에서 양방향 반복 측정 ANOVA에 대한 오류 항을 정의하는 데 문제가 있습니다. 내 데이터는 나무에서 추출한 코어를 따라 세 개의 방사형 위치 (내부, 중간 및 외부)에 대한 목재 밀도 추정치로 구성됩니다. 총 20 개의 나무 종, 각 종의 6 개 개체 및 각 나무에서 2 개의 코어가 있습니다. 목재 밀도에 …

10 r  anova  repeated-measures 

AIC를 사용하여 최상의 모델을 선택하여 혼합 모델 분석을 실행하는 방법에 대한 문헌을 찾으려고 노력한 후 압도적이라고 생각합니다. 내 데이터가 그렇게 복잡하다고 생각하지는 않지만 내가 한 일이 올바른지 확인한 다음 진행 방법에 대한 조언을 찾고 있습니다. REML 또는 ML을 사용해야하는 경우 lme 또는 lmer를 사용 해야하는지, 그중 하나를 사용 해야하는지 확실하지 …

10 r  lme4-nlme  random-effects-model  aic 

모수 분포에 대한 아이디어가 없으면 어떻게됩니까? 어떤 접근법을 사용해야합니까? 대부분의 경우 특정 변수가 특정 종의 존재 유무에 영향을 미치고 변수의 중요성에 따라 변수가 허용되는지 여부를 미달하려고합니다. 이것은 대부분의 경우 매개 변수가 가지고 있어야하는 배포에 대해 생각하지 않는다는 것을 의미합니다. b1, b2, b3 및 b4가 -2와 2 사이에서 변해야하고 b0이 -5와 …

10 r  bayesian  mcmc  bugs  winbugs 

내 질문은 mgcv R 패키지 의 GAM 을 다루고 있습니다 . 표본 크기가 작기 때문에 leave-one-out 교차 유효성 검사를 사용하여 예측 오류를 결정하고 싶습니다. 이것이 합리적입니까? 이 작업을 수행 할 수있는 패키지 또는 코드가 있습니까? ipred 패키지 의 errorest()기능 이 작동 하지 않습니다. 간단한 테스트 데이터 세트는 다음과 같습니다. library(mgcv) …

10 r  cross-validation  gam  mgcv