«r» 태그된 질문

lmer () 모델에서 예측 주위의 예측 간격을 얻고 싶습니다. 이에 대한 토론을 찾았습니다. http://rstudio-pubs-static.s3.amazonaws.com/24365_2803ab8299934e888a60e7b16113f619.html http://glmm.wikidot.com/faq 그러나 그들은 무작위 효과의 불확실성을 고려하지 않은 것으로 보인다. 구체적인 예는 다음과 같습니다. 나는 금붕어를 경주하고 있습니다. 지난 100 개 인종에 대한 데이터가 있습니다. RE 추정치와 FE 추정치의 불확실성을 고려하여 101st를 예측하고 싶습니다. 나는 물고기에 …

37 r  mixed-model  prediction  prediction-interval  lme4-nlme 

발트해의 해수면 온도가 해마다 동일하다고 가정하고 함수 / 선형 모델로 설명합니다. 내가 가진 아이디어는 연도를 10 진수 (또는 num_months / 12)로 입력하고 그 시간에 대한 온도를 알아내는 것이 었습니다. R의 lm () 함수에 던지면 정현파 데이터를 인식하지 않으므로 직선을 생성합니다. 그래서 sin () 함수를 I () 괄호 안에 넣고 함수에 …

37 r  regression  time-series  lm 

처리 후 ( Curevs No Cure) 결과 변수를 치료하는 로지스틱 회귀 분석을 작성했습니다 . 이 연구의 모든 환자는 치료를 받았다. 당뇨병이이 결과와 관련이 있는지 확인하고 싶습니다. R에서 로지스틱 회귀 출력은 다음과 같습니다. Call: glm(formula = Cure ~ Diabetes, family = binomial(link = "logit"), data = All_patients) ... Coefficients: Estimate Std. …

37 r  hypothesis-testing  logistic  generalized-linear-model  odds-ratio 

이 질문에는 이미 답변이 있습니다. 난수 생성기의 씨앗은 정확히 무엇입니까? 답변 3 개 나는 set.seed()의사 난수 생성을 위해 R에서 사용한다는 것을 알고 있습니다. 또한 같은 숫자를 사용하면 set.seed(123)결과를 재현 할 수 있습니다. 그러나 내가 얻지 못하는 것은 값 자체가 의미하는 것입니다. 나는 몇 가지 기능을 재생, 일부 사용하고 set.seed(1)또는 set.seed(300)나 …

37 r  random-generation 

나는 다양한 모델로 내 데이터에 적합하고 있음을 파악하려고했다 fitdistr라이브러리 함수 MASS의이 R저를주는 Negative Binomial가장 적합한다. 이제 위키 페이지에서 정의는 다음과 같습니다. NegBin (r, p) 분포는 마지막 시험에서 성공한 k + r Bernoulli (p) 시험에서 k 실패 및 r 성공 확률을 설명합니다. R모델 피팅을 수행하는 데 사용하면 두 개의 매개 변수 …

37 r  distributions  modeling  negative-binomial 

R 의 캐럿 패키지 를 사용하여 분류 및 회귀에 대한 예측 모델을 작성했습니다. Caret는 교차 검증 또는 부트 스트래핑을 통해 모델 하이퍼 파라미터를 조정할 수있는 통합 인터페이스를 제공합니다. 예를 들어 분류를 위해 간단한 '가장 가까운 이웃'모델을 구축하는 경우 몇 개의 이웃을 사용해야합니까? 2? 10? 100? Caret은 데이터를 다시 샘플링하고 다른 …

37 r  time-series  cross-validation 

현재 R 패키지 lme4를 사용하고 있습니다. 임의의 효과가있는 선형 혼합 효과 모델을 사용하고 있습니다. library(lme4) mod1 <- lmer(r1 ~ (1 | site), data = sample_set) #Only random effects mod2 <- lmer(r1 ~ p1 + (1 | site), data = sample_set) #One fixed effect + # random effects mod3 <- lmer(r1 …

37 r  regression  mixed-model  random-effects-model 

Effectspackage 는 package 를 통해 얻은 선형 혼합 효과 모델 결과 를 플로팅 하기위한 매우 빠르고 편리한 방법을 제공 lme4합니다 . 이 effect함수는 신뢰 구간 (CI)을 매우 빠르게 계산하지만 이러한 신뢰 구간은 얼마나 신뢰할 수 있습니까? 예를 들면 다음과 같습니다. library(lme4) library(effects) library(ggplot) data(Pastes) fm1 <- lmer(strength ~ batch + …

36 r  mixed-model  confidence-interval  effects  lme4-nlme 

내 상황에 가장 적합한 교차 유효성 검사 방법을 찾으려고합니다. 다음 데이터는 문제를 해결하기위한 예제 (R)이지만 실제 X데이터 ( xmat)는 서로 상관 관계가 있으며 y변수 ( ymat)를 사용하여 다른 정도와 상관 관계가 있습니다. R 코드를 제공했지만 R에 대한 질문이 아니라 메서드에 대한 질문입니다. XmatX 변수 V1-V100을 ymat포함하고 단일 y 변수 를 …

36 r  regression  cross-validation  linear-model 

내가 가진 일부 데이터에 2 차 다항식 적합을 만들려고합니다. 이 적합도를 플로팅한다고 가정 해 봅시다 ggplot(). ggplot(data, aes(foo, bar)) + geom_point() + geom_smooth(method="lm", formula=y~poly(x, 2)) 나는 얻다: 따라서 두 번째 주문 적합은 아주 잘 작동합니다. R로 계산합니다. summary(lm(data$bar ~ poly(data$foo, 2))) 그리고 나는 얻는다 : lm(formula = data$bar ~ poly(data$foo, …

36 r  regression  interpretation  regression-coefficients 

두 가지 다른 함수 ( prcomp및 princomp)를 사용하여 R을 사용하여 주성분 분석 (PCA)을 수행 했으며 PCA 점수가 다른 것으로 나타났습니다. 어떻게 할 수 있습니까? 이걸 고려하세요: set.seed(999) prcomp(data.frame(1:10,rnorm(10)))$x PC1 PC2 [1,] -4.508620 -0.2567655 [2,] -3.373772 -1.1369417 [3,] -2.679669 1.0903445 [4,] -1.615837 0.7108631 [5,] -0.548879 0.3093389 [6,] 0.481756 0.1639112 [7,] 1.656178 …

36 r  pca  factor-analysis 

샘플 크기가 다른 ( , , ) 3 개의 조류 바이오 매스 그룹 ( , , ) 에서 얻은 데이터를 가지고 그룹이 같은 집단에 속하는지 비교하고 싶습니다.AAAC n A = 15 n B = 13 n C = 12BBBCCCnA=15nA=15n_A=15nB=13nB=13n_B=13nC=12nC=12n_C=12 일원 분산 분석은 확실히 갈 길이지만 내 데이터에 대한 정규성 테스트를 …

36 r  anova  data-transformation  heteroscedasticity 

plot(density(rexp(100)) 분명히 0 왼쪽의 모든 밀도는 바이어스를 나타냅니다. 통계가 아닌 사람들을 위해 일부 데이터를 요약하려고하는데, 음이 아닌 데이터가 왜 밀도가 0의 왼쪽인지에 대한 질문을 피하고 싶습니다. 플롯은 무작위 검사 용입니다. 처리 및 제어 그룹별로 변수 분포를 보여주고 싶습니다. 분포는 종종 지수 적입니다. 히스토그램은 여러 가지 이유로 까다 롭습니다. 빠른 Google …

36 r  pdf  gamma-distribution  kernel-smoothing 

변수가 50,000 개 이상인 모델에 올가미 또는 능선 회귀를 사용하고 싶습니다. R의 소프트웨어 패키지를 사용하고 싶습니다. 수축 매개 변수 ( )를 어떻게 추정 할 수 있습니까?λλ\lambda 편집 : 여기에 내가 지적한 요점이있다 : set.seed (123) Y <- runif (1000) Xv <- sample(c(1,0), size= 1000*1000, replace = T) X <- matrix(Xv, …

36 r  lasso  ridge-regression  high-dimensional 

Amazon EC2를 기반으로 구축 된 클러스터에서 R, python 또는 octave 스크립트를 실행하기위한 플랫폼을 제공하는 작은 회사 목록이 있습니다. 추가해야 할 다른 이름이 있습니까? 구름 번호 오 파니 crdata

36 r  machine-learning