«r» 태그된 질문

R에서 glmnet을 사용하여 LASSO 모델을 올바르게 훈련하고 테스트하는 방법에 대해 궁금합니다. 특히 외부 테스트 데이터 세트가 부족한 경우 교차 검증 (또는 다른 유사한 방법)을 사용하여 LASSO 모델을 테스트 하는 방법을 궁금합니다 . 시나리오를 정리하겠습니다 : 내 glmnet 모델을 알리고 훈련 할 데이터 세트가 하나뿐입니다. 결과적으로 교차 검증을 사용하여 데이터를 분할하여 …

10 r  cross-validation  glmnet  lasso 

제목에서 알 수 있듯이 라이브러리의 LBFGS 옵티 마이저를 사용하여 glmnet linear의 결과를 복제하려고합니다 lbfgs. 이 옵티마이 저는 목적 함수 (L1 정규화 용어가없는)가 볼록한 한, 미분에 대해 걱정할 필요없이 L1 정규화 용어를 추가 할 수 있습니다. 의 탄성 순 선형 회귀 문제 glmnet 용지가 주어진다 여기서 X \ in \ mathbb …

10 r  regression  regularization  glmnet  elastic-net 

이 질문에는 하나 이상의 심각한 오해가있을 수 있지만 계산을 올바르게하는 것이 아니라 시계열을 배우는 데 중점을 둡니다. 시계열의 적용을 이해하려고 할 때, 데이터의 비추 세화로 미래의 가치를 예측할 수없는 것처럼 보입니다. 예를 들어 패키지의 gtemp시계열은 astsa다음과 같습니다. 예측 된 미래 가치를 표시 할 때 지난 수십 년간 상승 추세를 고려해야합니다. …

10 r  time-series  data-visualization 

우선이 포럼에서 몇 가지 조사를 해본 결과 매우 비슷한 질문에 대한 답변을 받았지만 일반적으로 제대로 답변되지 않았거나 때로는 답변이 상세하지 않은 경우가 있습니다. 그래서 이번에는 내 질문은 : 나는 두 개의 데이터 세트를 가지고 있는데, 각각에 대해 다항식 회귀를 수행합니다. Ratio<-(mydata2[,c(2)]) Time_in_days<-(mydata2[,c(1)]) fit3IRC <- lm( Ratio~(poly(Time_in_days,2)) ) 다항식 회귀 도표는 …

10 r  regression  statistical-significance  regression-coefficients  group-differences 

PCA가 제공하는 좌표 공간이 아니라 약간 다른 (회전) 벡터 세트에 대해 데이터 세트의 분산 백분율을 얻는 방법을 이해하고 싶습니다. set.seed(1234) xx <- rnorm(1000) yy <- xx * 0.5 + rnorm(1000, sd = 0.6) vecs <- cbind(xx, yy) plot(vecs, xlim = c(-4, 4), ylim = c(-4, 4)) vv <- eigen(cov(vecs))$vectors ee …

10 r  variance  pca  linear-algebra 

482 개의 관측치 데이터 세트가 있습니다. data=Populationfull 3 개의 SNP에 대한 유전자형 연관 분석을 수행하려고합니다. 내 분석을 위해 모델을 작성하려고 시도하고 aov (y ~ x, data = ...)를 사용하여 Im. 하나의 특성에 대해 다음과 같이 모델에 포함 된 여러 고정 효과 및 공변량이 있습니다. Starts <- aov(Starts~Sex+DMRT3+Birthyear+Country+Earnings+Voltsec+Autosec, data=Populationfull) summary(Starts) Df …

10 r  anova 

이 유형의 줄거리를 무엇이라고 부르며 R로 만들 수 있습니까? 편집 : 많은 감사합니다-매우 도움이됩니다. 지금까지 최고의 타이틀 : 정량화 된 바이올린 음모!

10 r  data-visualization 

교육 게임에 대한 연구를하고 있으며 현재 진행중인 일부 프로젝트에는 BoardGameGeek (BGG) 및 VideoGameGeek (VGG)의 데이터를 사용하여 게임의 디자인 요소 (예 : "제 2 차 세계 대전에서 설정 됨", "롤링 주사위 포함") 간의 관계를 조사합니다. ) 및 해당 게임의 플레이어 등급 (예 : 10 점 만점). 이러한 각 디자인 요소는 BGG …

10 r  categorical-data  factor-analysis  binary-data 

에 의해 추정 된 모델의 p,d and q값을 어떻게 얻을 수 있습니까?ARIMA(p,d,q)auto.arima(mytimeseries) arima_model <-auto.arima (mytimeseries, ic = 'bic') 우리가 출력을 보면 arima_model $ arma 우리는 얻는다 [1] 1000 0 2 0 위 순서대로 숫자의 의미는 무엇입니까?

10 r  arima 

비선형 경계의 분류 문제는 단순한 퍼셉트론 으로 해결할 수 없습니다 . 다음 R 코드는 설명을위한 것이며 Python 의이 예제 를 기반으로합니다 . nonlin <- function(x, deriv = F) { if (deriv) x*(1-x) else 1/(1+exp(-x)) } X <- matrix(c(-3,1, -2,1, -1,1, 0,1, 1,1, 2,1, 3,1), ncol=2, byrow=T) y <- c(0,0,1,1,1,0,0) syn0 …

10 r  neural-networks  svm  kernel-trick  perceptron 

p- 값을 결합 하기 위해 Fisher 또는 Stouffer의 방법 을 구현하는 R 패키지 (또는 기본 R 함수)가 있습니까? 이것을 코딩하는 것은 거의 사소한 것이지만 패키지를 사용하고 인용하고 싶습니다. 이 질문의 예제 코드 : p- 값을 결합하는 Fisher의 방법-아래쪽 꼬리는 어떻습니까?

10 r  combining-p-values 

1 분 재고 반품에 대해 2 개의 표본 평균을 비교하고 싶습니다. 나는 그들이 Laplace 분산되어 있다고 가정하고 (이미 확인) 반환을 2 그룹으로 나눕니다. 그것들이 크게 다른지 어떻게 확인할 수 있습니까? 값이 300 개가 넘더라도 QQ 플롯은 정규 분포와 큰 차이가 있기 때문에 정규 분포처럼 취급 할 수 없다고 생각합니다.

10 r  confidence-interval  laplace-distribution 

에서 lmer내 기능 lme4에 R임의 효과의 모델 매트릭스 구성에 대한 요구가 같은 설명, 여기 , 페이지 7-9.ZZZ 계산 하는 것은 J_i 및 X_i 의 두 행렬의 KhatriRao 및 / 또는 Kronecker 곱을 수반 합니다. J i X iZZZJiJiJ_iXiXiX_i 행렬 JiJiJ_i 는 한 입 가득입니다. "그룹화 인자 지수의 지표 행렬"이지만 더 …

10 r  mixed-model  lme4-nlme  matrix 

다음과 같은 간단한 X 및 Y 벡터가 있습니다. > X [1] 1.000 0.063 0.031 0.012 0.005 0.000 > Y [1] 1.000 1.000 1.000 0.961 0.884 0.000 > > plot(X,Y) log of X를 사용하여 회귀를 원합니다. log (0)을 얻지 않으려면 +1 또는 +0.1 또는 +0.00001 또는 +0.000000000000001을 넣으십시오. > summary(lm(Y~log(X))) Error …

10 r  regression  lognormal 

상황을 상상해보십시오. 우리는 3 개의 광산에 대한 역사적 기록 (20 년)을 가지고 있습니다. 은이 있으면 내년에 금을 찾을 확률이 높아 집니까? 그러한 질문을 테스트하는 방법? 예제 데이터는 다음과 같습니다. mine_A <- c("silver","rock","gold","gold","gold","gold","gold", "rock","rock","rock","rock","silver","rock","rock", "rock","rock","rock","silver","rock","rock") mine_B <- c("rock","rock","rock","rock","silver","rock","rock", "silver","gold","gold","gold","gold","gold","rock", "silver","rock","rock","rock","rock","rock") mine_C <- c("rock","rock","silver","rock","rock","rock","rock", "rock","silver","rock","rock","rock","rock","silver", "gold","gold","gold","gold","gold","gold") time <- seq(from = 1, to …

10 r  time-series  hypothesis-testing  stochastic-processes