«r» 태그된 질문

R을 사용하여 일련의 데이터를 가진 로그 정규 모델을 피팅했습니다. 결과 매개 변수는 다음과 같습니다. meanlog = 4.2991610 sdlog = 0.5511349 이 모델을 Scipy로 이전하고 싶습니다. Scipy를 사용하여 1과 3.1626716539637488e + 90의 모양과 스케일을 얻을 수있었습니다. 매우 다른 숫자입니다. 또한 meanlog 및 sdlog의 exp를 사용하려고했지만 기괴한 그래프를 계속 얻습니다. 나는 scipy에서 …

10 r  python  numpy  scipy 

R을 사용하여 Mann-Whitney U의 임계 값을 찾으려고하면 값이 항상 1+ 임계 값이라는 것을 알았습니다. 예를 들어, 의 경우 (양측) 임계 값은 8이며 의 경우 (양측) ) 임계 값은 22 ( 표 확인 )이지만 다음과 같습니다.α=.05,n=10,m=5α=.05,n=10,m=5\alpha=.05, n = 10, m = 5α=.05,n=12,m=8α=.05,n=12,m=8\alpha=.05, n=12, m=8 > qwilcox(.05/2,10,5) [1] 9 > qwilcox(.05/2,12,8) [1] …

10 r  hypothesis-testing  nonparametric  wilcoxon-mann-whitney 

사용자 정의 대비로 (일종 당) 일원 분산 분석을 수행하고 있습니다. [,1] [,2] [,3] [,4] 0.5 -1 0 0 0 5 1 -1 0 0 12.5 0 1 -1 0 25 0 0 1 -1 50 0 0 0 1 여기서 강도 0.5를 5와 5, 5를 12.5와 비교합니다. 이것들은 내가 작업중 …

10 r  anova  contrasts  generalized-least-squares 

두 가지 다른 시간 변수를 모델링하고 싶습니다. 일부 변수는 내 데이터에서 상당히 공 선적입니다 (나이 + 코호트 = 기간). 이렇게하면에 대한 몇 가지 문제와 lmer상호 작용이 발생 poly()했지만 IIRC lmer와 동일한 결과를 얻었을 것 nlme입니다. 분명히 poly () 함수의 기능에 대한 이해가 부족합니다. 나는 무엇을 이해하고 그것 poly(x,d,raw=T)없이는 raw=T직교 다항식 …

10 r  lme4-nlme  polynomial 

두 주가의 시계열 간의 리드 랙을 분석하려고합니다. 정기적 인 시계열 분석에서는 VECM (Granger Causality) 인 Cross Correlaton을 수행 할 수 있습니다. 그러나 불규칙한 간격의 시계열에서 어떻게 동일하게 처리합니까? 가설은 계측기 중 하나가 다른 쪽을 이끈다는 것입니다. 두 기호에 대한 데이터가 마이크로 초에 있습니다. RTAQ 패키지를 살펴보고 VECM을 적용 해 보았습니다. …

10 r  time-series  cross-correlation  unevenly-spaced-time-series 

동료의 작업을 복제하려고하고 분석을 Stata에서 R로 옮기려고합니다. 그녀가 사용하는 모델은 nbreg 함수 내에서 "cluster"옵션을 호출하여 표준 오류를 군집화합니다. 이 옵션의 의미와 이유에 대한 자세한 설명은 http://repec.org/usug2007/crse.pdf 를 참조하십시오. 내 질문은 R 내에서 음의 이항 회귀에 대해 동일한 옵션을 호출하는 방법입니다. 본 논문의 주요 모델은 Stata에 다음과 같이 명시되어있다 xi: nbreg …

10 r  stata  negative-binomial 

lm(y~a*b)R 구문에서 a이진 변수이고 b숫자 변수 인 데이터와 같이 데이터를 맞추면 교호 a:b작용 항은 y~bat a= 0과 at a= 1 사이의 차이 입니다. 자,하자가 사이의 관계를 말할 y와 b곡선이다. 지금 맞는 경우 lm(y~a*poly(b,2)), 다음 a:poly(b,2)1의 변화의 변화 y~b의 수준을 조건으로 a위와 같이하고 a:poly(b,2)2있는 변화 y~b^2의 수준에 조건 a. 약간의 핸드 …

10 r  multiple-regression  splines 

간단한 확률이 내 문제에 효과가 있는지 또는 로지스틱 회귀와 같은보다 정교한 방법을 사용하고 배우는 것이 더 나은지 결정하려고합니다. 이 문제의 반응 변수는 이항 반응 (0, 1)입니다. 나는 모두 범주적이고 순서가없는 많은 예측 변수를 가지고 있습니다. 예측 변수의 조합이 1의 가장 높은 비율을 산출하는지 확인하려고합니다. 로지스틱 회귀가 필요합니까? 범주 형 예측 …

10 r  probability  logistic 

정규화를 수행하는 방법에는 , L 1 및 L 2 규범 기반 정규화 와 같은 많은 방법이 있습니다 . Friedman Hastie & Tibsharani 에 따르면 , 최고의 정규화 기는 문제, 즉 실제 목표 함수의 특성, 사용 된 특정 기준, 신호 대 잡음비 및 샘플 크기에 따라 달라집니다.엘0엘0L_0엘1엘1L_1엘2엘2L_2 다양한 정규화 방법의 방법과 …

10 r  regression  machine-learning  regularization 

3 회 반복 측정 ANOVA를 수행했습니다. 어떤 사후 분석이 유효합니까? 이것은 개체 내 반복 측정을 갖는 요소 중 하나를 가진 완전 균형 설계 (2x2x2)입니다. R에서 반복 측정 ANOVA에 대한 다변량 접근법을 알고 있지만, 첫 번째 본능은 간단한 aov () 스타일의 ANOVA로 진행하는 것입니다. aov.repeated <- aov(DV ~ IV1 * IV2 …

10 r  anova  repeated-measures  interaction  post-hoc 

두 개의 교차 임의 효과가 있다고 생각되는 일부 데이터를 모델링하고 있습니다. 그러나 데이터 세트의 균형이 맞지 않으므로이를 설명하기 위해 수행해야 할 작업이 확실하지 않습니다. 내 데이터는 일련의 이벤트입니다. 클라이언트가 공급자와 만나 작업을 수행 할 때 이벤트가 발생합니다 (성공 여부). 수천 명의 고객과 제공자가 있으며 각 고객과 제공자는 다양한 이벤트 (대략 …

10 r  mixed-model  random-effects-model  logistic 

SEM을 적용하려는 하나의 데이터 세트에 많은 문제가 있습니다. 우리는 지표가 각각 5 가지 잠재 요인 A, B, C, D, E의 존재를 가정합니다. A1-A5 (순서 계수), B1-B3 (정량적), C1, D1, E1 (E1에 대해 2 단계 만 포함 된 마지막 세 가지 요인 모두) 모든 요소 사이의 공분산에 관심이 있습니다. 나는 OpenMx그렇게 …

10 r  modeling  multiple-regression  sem 

임상 재활 데이터의 데이터 세트를 분석해야합니다. 정량화 된 "입력"(치료량)과 건강 상태의 변화 사이의 가설 중심 관계에 관심이 있습니다. 데이터 세트는 비교적 작지만 (n ~ 70) 두 가지 모두의 시간적 변화를 반영하는 데이터를 반복했습니다. R의 비선형 혼합 효과 모델링에 익숙하지만 여기서 입력과 출력 사이의 "인과적인"관계에 관심이 있으므로 SEM의 반복 측정 응용 …

10 r  repeated-measures  panel-data  sem 

두 가지 "적합도 카이-제곱"테스트를 비교할 솔루션을 찾고 있습니다. 보다 정확하게는 두 개의 독립적 인 실험 결과를 비교하고 싶습니다. 이 실험에서 저자들은 랜덤 추정 (예상 주파수)과 관측 된 주파수를 비교하기 위해 적합도 카이 제곱을 사용했습니다. 두 실험의 참가자 수는 동일하고 실험 절차는 동일하며 자극 만 변경되었습니다. 두 실험 결과는 유의 한 …

10 r  confidence-interval  chi-squared 

약 70 개의 변수가있는 데이터 세트가 있습니다. 내가 찾고있는 것은 CV를 사용하여 다음과 같은 방식으로 가장 유용한 변수를 찾는 것입니다. 1) 20 개의 변수를 임의로 선택합니다. 2) stepwise/ LASSO/ lars/ etc를 사용 하여 가장 중요한 변수를 선택하십시오. 3) ~ 50x를 반복하고 어떤 변수가 가장 자주 선택 (제거되지 않음)되었는지 확인합니다. 이것은 …

10 r  cross-validation  feature-selection  random-forest  stepwise-regression