«r» 태그된 질문

시간 의존적 공변량을 포함하는 Cox 비례 위험 모델에서 생존 시간을 생성하고 싶습니다. 모델은 h ( t |엑스나는) =h0( t ) 특급( γ엑스나는+ α미디엄나는( t ) )h(티|엑스나는)=h0(티)특급⁡(γ엑스나는+α미디엄나는(티))h(t|X_i) =h_0(t) \exp(\gamma X_i + \alpha m_{i}(t)) 여기서 는 Binomial (1,0.5) 및 됩니다.엑스나는엑스나는X_i미디엄나는( t ) =β0+β1엑스나는+β2엑스나는티미디엄나는(티)=β0+β1엑스나는+β2엑스나는티m_{i}(t)=\beta_0 + \beta_1 X_{i} + \beta_2 X_{i} t 실제 매개 …

9 r  survival  cox-model  time-varying-covariate 

내 질문은 일반적으로 SVD (Singular Value Decomposition), 특히 LSI (Latent Semantic Indexing)에 관한 것입니다. 예 를 들어, 7 개의 문서에 대해 5 개의 단어 빈도를 포함하는 가 있습니다.ㅏw o r d× do c u m e n tAword×document A_{word \times document} A = matrix(data=c(2,0,8,6,0,3,1, 1,6,0,1,7,0,1, 5,0,7,4,0,5,6, 7,0,8,5,0,8,5, 0,10,0,0,7,0,0), ncol=7, byrow=TRUE) …

9 r  svd  natural-language  latent-semantic-indexing 

로지스틱 회귀를 사용하여 조정 된 승산 비를 추정하는 방법과 유사한 조정 된 위험률을 추정하는 데 관심이 있습니다. 일부 문헌 (예 : this )은 Huber-White 표준 오류와 함께 포아송 회귀 분석을 사용하는 것이 모델 기반 방법임을 나타냅니다. 연속 공변량을 조정하는 것이 어떻게 영향을 미치는지에 대한 문헌을 찾지 못했습니다. 다음의 간단한 시뮬레이션은이 …

9 r  references  biostatistics  poisson-regression  relative-risk 

200 과목과 1000 변수로 내 데이터 예측의 교차 유효성 검사를 진행하고 있습니다. 변수 수 (사용하고 싶습니다)가 샘플 수보다 많기 때문에 능선 회귀에 관심이 있습니다. 그래서 수축 견적 도구를 사용하고 싶습니다. 다음은 예제 데이터로 구성됩니다. #random population of 200 subjects with 1000 variables M <- matrix(rep(0,200*100),200,1000) for (i in 1:200) { …

9 r  cross-validation  prediction  ridge-regression 

R 패키지에 glmer모델을 장착하고 있습니다 lme4. p- 값이 표시된 anova 테이블을 찾고 있지만 적합한 패키지를 찾을 수 없습니다. R에서 할 수 있습니까? 내가 맞는 모델은 다음과 같은 형식입니다. model1<-glmer(dmn~period*teethTreated+(1|fullName), family="poisson", data=subset(dataset, group=='Four times a year'), control=glmerControl(optimizer="bobyqa"))

9 r  anova  lme4-nlme  mixed-model 

약 5000 개의 기능으로 구성된 데이터 세트가 있습니다. 해당 데이터에 대해 먼저 피처 선택에 Chi Square 테스트를 사용했습니다. 그 후 응답 변수와의 유의미한 관계를 나타내는 약 1500 개의 변수를 얻었습니다. 이제 로지스틱 회귀 분석을 적용해야합니다. R에 glmulti 패키지를 사용하고 있습니다 (glmulti 패키지는 vlm에 효율적인 하위 세트 선택을 제공합니다). 한 번에 …

9 r  logistic  generalized-linear-model  modeling  regression-strategies 

glmR에서 함수를 사용할 때 기본적 으로이 함수는 반복적으로 가중 된 최소 제곱 (IWLS) 방법을 사용하여 모수를 추정 할 수있는 최대 가능성을 찾습니다. 이제 두 가지 질문이 있습니다. IWLS 추정은 가능성 함수의 전체 최대 값을 보장합니까? 이 프레젠테이션 의 마지막 슬라이드를 기반으로 , 나는 그렇지 않다고 생각합니다! 나는 단지 그것을 확인하고 …

9 r  estimation  generalized-linear-model  maximum-likelihood  optimization 

누군가 구매할 품목 그룹을 예측하려고합니다. 즉, 여러 개의 선형 선형 변수가 있습니다. 7 개 항목을 구매할 확률을 예측하기 위해 7 개 정도의 독립적 인 모델을 구축 한 다음 결과를 결합하기보다는 7 개의 관련 종속 변수 ( 그들이 구입할 수있는 것). R을 프로그래밍 언어로 사용하고 있으므로 R 관련 조언을 부탁드립니다.

9 r  machine-learning  multivariate-analysis 

Kruschke의 Doing Bayesian Data Analysis (특히 Poisson 지수 분산 분석) 의 예제를 진행 하고 있습니다. 도 22는 우발 사태 테이블에 대한 빈번한 카이 제곱 독립성 테스트의 대안으로 제시하고있다. 변수가 독립적 인 경우 (예 : HDI가 0을 제외 할 때) 예상보다 많거나 적은 상호 작용에 대한 정보를 얻는 방법을 알 수 …

9 r  bayesian  effect-size  contingency-tables 

일, 주 또는 월과 같은 기간이 많은 시계열이 있습니다. stl()기능을 사용하거나 사용하면 loess(x ~ y)특정 시계열의 추세가 어떻게 보이는지 볼 수 있습니다. 시계열 추세가 증가 또는 감소하는지 감지해야합니다. 어떻게 관리 할 수 ​​있습니까? 나는 선형 회귀 계수를 계산 lm(x ~ y)하고 경사 계수로 재생 하려고했습니다 . ( If |slope|>2 and …

9 r  time-series  trend 

시변 계수, 즉 일반적인 선형 회귀에 대한 확장으로 DLM을 맞추고 싶습니다. yt=θ1+θ2x2yt=θ1+θ2x2y_t = \theta_1 + \theta_2x_2 입니다. 1950 ~ 2011 년에 각각 예측 변수 ( )와 반응 변수 ( ), 해양 및 내륙 연간 어획량이 있습니다. DLM 회귀 모델을 따르고 싶습니다.x2x2x_2ytyty_t yt=θt,1+θt,2xtyt=θt,1+θt,2xty_t = \theta_{t,1} + \theta_{t,2}x_t 시스템 진화 방정식은 θt=Gtθt−1θt=Gtθt−1\theta_t …

9 r  regression  time-series  dlm  dynamic-regression 

계산한다고 주장하는 Excel 시트를보고 있지만이 방법을 인식하지 못하고 뭔가 빠졌는지 궁금합니다.χ2χ2\chi^2 분석하는 데이터는 다음과 같습니다. +------------------+----------+----------+ | Total Population | Observed | Expected | +------------------+----------+----------+ | 2000 | 42 | 32.5 | | 2000 | 42 | 32.5 | | 2000 | 25 | 32.5 | | 2000 | 21 …

9 r  chi-squared  excel 

내 동료가 나에게이 문제를 보냈고 분명히 인터넷에서 라운드를 진행했다. If $3 = 18, 4 = 32, 5 = 50, 6 = 72, 7 = 98$, Then, $10 =$ ? 답은 200 인 것 같습니다. 3*6 4*8 5*10 6*12 7*14 8*16 9*18 10*20=200 R에서 선형 회귀를 할 때 : data <- …

9 r  regression  lm 

통계에 익숙하지 않으며 분산과 선형 회귀의 차이점을 이해하려고합니다. 나는 이것을 탐구하기 위해 R을 사용하고 있습니다. 나는 왜 분산과 회귀가 다르지만 여전히 똑같고 어떻게 시각화 될 수 있는지에 대한 다양한 기사를 읽었습니다. ANOVA는 그룹 내 분산과 그룹 간 분산을 비교하여 테스트 된 그룹간에 차이가 있는지 여부를 확인합니다. ( https://controls.engin.umich.edu/wiki/index.php/Factor_analysis_and_ANOVA ) 선형 …

9 r  regression  anova 

R에서 다른 표본 크기에 대한 정규성 검정의 정확성을 평가하고 싶습니다 (정상 성 테스트가 오도 할 수 있음을 알고 있습니다 ). 예를 들어 Shapiro-Wilk 테스트를 살펴보기 위해 다음과 같은 시뮬레이션을 수행하고 결과를 플로팅하고 샘플 크기가 증가함에 따라 null 감소를 거부 할 확률이 예상됩니다. n <- 1000 pvalue_mat <- matrix(NA, ncol = …

9 r  simulation  power-analysis  normality-assumption