«r» 태그된 질문

-3, -1, 1, 3 대비로 개체 간 분산 분석을 수행 할 수있는 R 코드를 제공하십시오. 그러한 분석에 적절한 SS (Sum of Squares) 유형에 관한 논쟁이 있음을 이해합니다. 그러나 SAS 및 SPSS (Type III)에 사용되는 기본 SS 유형은 필자의 지역 표준으로 간주됩니다. 따라서이 분석 결과가 해당 통계 프로그램에 의해 생성 된 …

26 r  anova  contrasts  sums-of-squares 

이 사이트의 다른 곳에서 한 질문 에 따르면, AIC는 LOO (Leave-One-Out) 교차 검증과 동일하고 BIC는 K- 폴드 교차 검증과 동일합니다. LOO 및 K-fold에 관련된 기술이 명확하고 AIC 및 BIC 값과 동등한 것으로 입증되도록 R에서 이것을 경험적으로 입증하는 방법이 있습니까? 주석이 달린 코드는 이와 관련하여 도움이 될 것입니다. 또한 BIC를 시연 …

26 r  aic  cross-validation  bic 

결정 요인이 0 인 보안 수익의 상관 관계 행렬이 있습니다. (이는 샘플 상관 행렬과 해당 공분산 행렬이 이론적으로 양의 명확한 값이어야하기 때문에 약간 놀라운 일입니다.) 내 가설은 적어도 하나의 보안이 다른 유가 증권에 선형 적으로 의존한다는 것입니다. R에 선형 의존성을 위해 각 열의 행렬을 순차적으로 테스트하는 함수가 있습니까? 예를 들어, …

26 r  correlation  pca  linear-model  svd 

나는 R을 처음 사용하고 로지스틱 회귀를 주문했다 polr. polr 에 대한 도움말 페이지 하단의 "예제"섹션 (로지스틱 또는 프로 빗 회귀 모델을 정렬 된 요인 반응에 적합) options(contrasts = c("contr.treatment", "contr.poly")) house.plr <- polr(Sat ~ Infl + Type + Cont, weights = Freq, data = housing) pr <- profile(house.plr) plot(pr) pairs(pr) …

26 r  logistic 

R에서 LOESS 회귀 모델을 실행 중이며 12 가지 모델의 출력을 다양한 샘플 크기와 비교하려고합니다. 질문에 대답하는 데 도움이되는 경우 실제 모델을 더 자세히 설명 할 수 있습니다. 샘플 크기는 다음과 같습니다. Fastballs vs RHH 2008-09: 2002 Fastballs vs LHH 2008-09: 2209 Fastballs vs RHH 2010: 527 Fastballs vs LHH 2010: …

26 r  regression  loess 

R의 피어슨 상관 관계에서 p- 값을 찾을 수 있습니까? 피어슨 상관 관계를 찾기 위해 나는 보통 이렇게합니다 col1 = c(1,2,3,4) col2 = c(1,4,3,5) cor(col1,col2) # [1] 0.8315218 그러나 어떻게 이것의 p- 값을 찾을 수 있습니까?

26 r  correlation  p-value  pearson-r 

많은 도움말 사이트를 and이 뒤져서 혼합 모델에서 더 복잡한 중첩 용어를 지정하는 방법에 대해 여전히 혼란 스럽습니다. 또한의 사용으로 혼란 스러워요 :과 /와 |상호 작용을 지정하고 사용하여 임의 요소와 중첩에 lmer()에서 lme4의 패키지 R. 이 질문의 목적을 위해이 표준 통계 모델로 내 데이터를 정확하게 묘사했다고 가정합니다. 고정되고, 및 무작위입니다. 안에 …

26 r  mixed-model  lme4-nlme 

이것은 질문 캐럿 리샘플링 방법 과 유사 하지만, 실제로이 부분에 대해서는 합의 된 방식으로 답변 한 적이 없습니다. 캐럿의 열차 기능 제공 cv및 repeatedcv. 말하는 것과의 차이점은 무엇입니까? MyTrainControl=trainControl( method = "cv", number=5, repeats=5 ) vs MyTrainControl=trainControl( method = "repeatedcv", number=5, repeats=5 ) cv세트를 k- 폴드 (parameter number) 로 나누고 …

26 r  machine-learning  caret 

이전에 예측 프로 를 사용 하여 일 변량 시계열을 예측했지만 워크 플로를 R로 전환하고 있습니다. R에 대한 예측 패키지에는 유용한 기능이 많이 포함되어 있지만 자동으로 실행하기 전에 데이터 변환이 필요하지 않습니다. .arima (). 경우에 따라 예측 전문가는 예측을 수행하기 전에 변환 데이터를 로그하기로 결정하지만 아직 이유를 찾지 못했습니다. 그래서 내 …

26 r  time-series  data-transformation  forecasting  arima 

TED 에서 Peter Donnelly의 이야기에서 영감을 받아 일련의 동전 던지기에 특정 패턴이 나타나는 데 걸리는 시간에 대해 이야기합니다. 나는 R에서 다음 스크립트를 만들었습니다. 이 패턴 중 하나에 도달하기까지 평균 시간 (예 : 동전 던지기 횟수)을 계산합니다. coin <- c('h','t') hit <- function(seq) { miss <- TRUE fail <- 3 trp …

26 r  probability  stochastic-processes 

혼합 효과 모델링 분야의 리더로부터 아래 블록 인용문은 랜덤 효과 ( 'ZCP'모델) 간의 상관 관계가 0 인 모델의 좌표 이동이 모델 예측을 변경한다고 주장합니다. 그러나 누군가가 자신의 주장을 자세히 설명하거나 정당화 할 수 있습니까? 문제의 진술은 7 페이지, 두 번째 단락 ( 다운로드 링크 ) 에 대한 Bates et al의 …

25 r  mixed-model  lme4-nlme 

y정규 예측 변수의 카운트 데이터 를 예측하는 장애물 모델을 고려하십시오 x. set.seed(1839) # simulate poisson with many zeros x <- rnorm(100) e <- rnorm(100) y <- rpois(100, exp(-1.5 + x + e)) # how many zeroes? table(y == 0) FALSE TRUE 31 69 이 경우 69의 0과 31의 양의 카운트로 …

25 r  count-data  zero-inflation 

특정 변수의 값 (~ 15)이 시간이 지남에 따라 어떻게 변하는 지 보여주고 싶지만, 매년 변수가 어떻게 다른지 보여주고 싶습니다. 그래서 나는이 줄거리를 만들었습니다. 그러나 색 구성표를 변경하거나 다른 선 / 모양 유형을 추가 할 때도 지저분 해 보입니다. 이런 종류의 데이터를 시각화하는 더 좋은 방법이 있습니까? R 코드를 사용한 테스트 …

25 r  data-visualization 

이 lmerTest패키지는 anova()선택적으로 Satterthwaite (기본값) 또는 Kenward-Roger의 자유도 (df) 근사값을 갖는 선형 혼합 모형에 대한 기능을 제공합니다 . 이 두 접근법의 차이점은 무엇입니까? 언제 선택해야합니까?

25 r  anova  mixed-model  lme4-nlme  degrees-of-freedom 

현재 혼합 효과 선형 모델을 사용하고 있습니다. R에서 "lme4"패키지를 사용하고 있습니다. 내 모델은 다음과 같은 형식을 취합니다. model <- lmer(response ~ predictor1 + predictor2 + (1 | random effect)) 모형을 실행하기 전에 예측 변수 간의 가능한 다중 공선 성을 확인했습니다. 나는 이것을함으로써 이것을했다 : 예측 변수의 데이터 프레임 만들기 dummy_df …

25 r  correlation  mixed-model  lme4-nlme  multicollinearity