«r» 태그된 질문

MCMCglmm혼합 효과 모델을 실행하기 위해 R 의 패키지를 사용하는 성능 문제가 있습니다. 코드는 다음과 같습니다. MC1<-MCMCglmm(bull~1,random=~school,data=dt,family="categorical" , prior=list(R=list(V=1,fix=1), G=list(G1=list(V=1, nu=0))) , slice=T, nitt=iter, ,burnin=burn, verbose=F) 데이터에는 약 20,000 개의 관측치가 있으며 약 200 개의 학교에 모여 있습니다. 실행하기 전에 데이터 프레임에서 사용하지 않는 모든 변수를 삭제하고 메모리에서 다른 모든 개체를 …

12 r  mixed-model  mcmc  computational-statistics 

다양한 기능적 데이터 분석 방법을 테스트하려고합니다. 이상적으로는 시뮬레이션 된 기능 데이터에 대한 접근 방식의 패널을 테스트하고 싶습니다. 합산 가우시안 노이즈 (아래 코드)를 기반으로 한 접근법을 사용하여 시뮬레이션 FD를 생성하려고 시도했지만 결과 곡선이 실제 와 비교하여 너무 견고 해 보입니다 . 누군가가보다 사실적으로 보이는 시뮬레이션 기능 데이터를 생성하기 위해 함수 / …

12 r  simulation  functional-data-analysis 

주어진 데이터 셋 : x <- c(4.9958942,5.9730174,9.8642732,11.5609671,10.1178216,6.6279774,9.2441754,9.9419299,13.4710469,6.0601435,8.2095239,7.9456672,12.7039825,7.4197810,9.5928275,8.2267352,2.8314614,11.5653497,6.0828073,11.3926117,10.5403929,14.9751607,11.7647580,8.2867261,10.0291522,7.7132033,6.3337642,14.6066222,11.3436587,11.2717791,10.8818323,8.0320657,6.7354041,9.1871676,13.4381778,7.4353197,8.9210043,10.2010750,11.9442048,11.0081195,4.3369520,13.2562675,15.9945674,8.7528248,14.4948086,14.3577443,6.7438382,9.1434984,15.4599419,13.1424011,7.0481925,7.4823108,10.5743730,6.4166006,11.8225244,8.9388744,10.3698150,10.3965596,13.5226492,16.0069239,6.1139247,11.0838351,9.1659242,7.9896031,10.7282936,14.2666492,13.6478802,10.6248561,15.3834373,11.5096033,14.5806570,10.7648690,5.3407430,7.7535042,7.1942866,9.8867927,12.7413156,10.8127809,8.1726772,8.3965665) .. 모수를 추정하여 가장 적합한 확률 분포 (감마, 베타, 정규, 지수, 포아송, 카이 제곱 등)를 결정하고 싶습니다. R을 사용하여 솔루션을 제공하는 다음 링크에 대한 질문을 이미 알고 있습니다 : https : //.com/questions/2661402/given-a-set-of-random-numbers-drawn-from-a- 연속 일 변량 분포- 제안 된 최상의 솔루션은 다음과 같습니다. > …

12 r  matlab  distributions 

JAGS 언어로 지정된 모델에서 rjags를 사용하여 MCMC를 실행했습니다. 해당 모형을 추출하고 모형으로 예측을 수행하는 좋은 방법이 있습니까 (내 매개 변수의 사후 분포를 사용하여)? R에서 모델을 다시 지정하고 매개 변수 후부의 모드를 연결할 수 있습니다. 이 작업을 덜 중복하는 방법이 있는지 궁금합니다. http://sourceforge.net/p/mcmc-jags/discussion/610037/thread/0ecab41c 가 같은 질문을하고 있다고 생각 합니다.

12 r  jags 

참고 :이 질문은 법적 이유로 인해 이전 질문을 삭제해야했기 때문에 다시 게시되었습니다. SAS의 PROC MIXED를 R lme의 nlme패키지 기능과 비교하는 동안 다소 혼란스러운 차이점을 발견했습니다. 구체적으로는, 다른 시험에서 자유도간에 상이 PROC MIXED하고 lme, 그리고 왜 생각해. 다음 데이터 세트에서 시작하십시오 (아래 제공된 R 코드). ind : 측정 대상을 나타내는 계수 …

12 r  mixed-model  sas  degrees-of-freedom  pdf  unbiased-estimator  distance-functions  functional-data-analysis  hellinger  time-series  outliers  c++  relative-risk  absolute-risk  rare-events  regression  t-test  multiple-regression  survival  teaching  multiple-regression  regression  self-study  t-distribution  machine-learning  recommender-system  self-study  binomial  standard-deviation  data-visualization  r  predictive-models  pearson-r  spearman-rho  r  regression  modeling  r  categorical-data  data-visualization  ggplot2  many-categories  machine-learning  cross-validation  weka  microarray  variance  sampling  monte-carlo  regression  cross-validation  model-selection  feature-selection  elastic-net  distance-functions  information-theory  r  regression  mixed-model  random-effects-model  fixed-effects-model  dataset  data-mining 

이전 글 에서 EQ-5D 점수 를 다루는 방법에 대해 궁금했습니다 . 최근에 Bottai와 McKeown 이 제안한 로지스틱 Quantile 회귀 분석을 통해 우연한 결과를 다루는 우아한 방법을 소개했습니다. 공식은 간단합니다. logit(y)=log(y−yminymax−y)logit(y)=log(y−yminymax−y)logit(y)=log(\frac{y-y_{min}}{y_{max}-y}) log (0) 및 0으로 나누지 않도록 범위를 작은 값인 확장하십시오 . 이것은 점수의 경계를 존중하는 환경을 제공합니다.ϵϵ\epsilon 문제는 모든 가 …

12 r  logistic  data-visualization  logit  quantile-regression 

에 R의해 생성 된 데이터 샘플이 rnorm(50,0,1)있으므로 데이터가 분명히 정규 분포를 취합니다. 그러나 R데이터에 대한이 배포 정보를 "알지"않습니다. R샘플이 어떤 종류의 분포에서 나오는지 추정 할 수 있는 방법 이 있습니까? 그렇지 않다면이 shapiro.test기능 을 사용 하고 계속 진행할 것입니다.

12 r  distributions 

나는 부분 최소 제곱 (PLS)에서 아주 새로운 오전 나는 R 함수의 출력을 이해하려고 노력 plsr()에서 pls패키지를. 데이터를 시뮬레이션하고 PLS를 실행하겠습니다 : library(pls) n <- 50 x1 <- rnorm(n); xx1 <- scale(x1) x2 <- rnorm(n); xx2 <- scale(x2) y <- x1 + x2 + rnorm(n,0,0.1); yy <- scale(y) p <- plsr(yy …

12 r  regression  partial-least-squares 

R에서 최근 성, 빈도 및 금전적 가치 (RFM) 모델링 및 고객 가치 모델링을 수행하는 방법을 아는 사람이 있습니까? 또한 누군가 나에게 그것에 대한 몇 가지 문헌을 언급 할 수 있습니까?

12 r  references  marketing 

최근에는 다단계 모델링에 대한 입문 과정을 수강했습니다. 우리가 사용한 대부분의 데이터 세트와 예제는 사회 과학에서 얻은 것입니다. 나는 생물 통계학과에서 2 주 인턴쉽을 얻었습니다. 병원 간 및 5 년 이상 사망률이 높은 응급 상황에서 환자 결과의 병원 수준 변화에 관한 프로젝트를 시작하기를 원합니다. 시간 범위. 나는 다음 주에 인턴쉽을 시작하고 …

12 r  stata  multilevel-analysis  epidemiology 

한 그룹이 다른 그룹보다 더 많은 훈련 데이터를 사용하여 SVM을 구축하려고합니다. 그러나 그룹은 최종 테스트 데이터에서 동일하게 표시됩니다. 따라서 R 패키지 인터페이스 의 class.weights매개 변수 를 사용 하여 훈련 데이터에서 두 그룹의 영향을 균형있게 조정 하고 싶습니다 .e1071libsvm 이 가중치를 정확히 어떻게 지정해야하는지 잘 모르기 때문에 약간의 테스트를 설정했습니다. 일부 …

12 r  machine-learning  svm 

MC 시뮬레이션을 사용하여 다음 적분을 어떻게 근사합니까? ∫1− 1∫1− 1| x-y|d xD Y∫−11∫−11|x−y|dxdy \int_{-1}^{1} \int_{-1}^{1} |x-y| \,\mathrm{d}x \,\mathrm{d}y 감사! 편집 (일부 상황) : 시뮬레이션을 사용하여 적분을 근사화하는 방법을 배우려고 노력하고 있으며 어려움을 겪을 때 연습을하고 있습니다. 편집 2 + 3 : 어떻게 든 혼란스러워서 적분을 별도의 부분으로 나눌 필요가 있다고 …

12 r  self-study  monte-carlo 

잠김 . 이 질문과 주제는 주제가 다르지만 역사적으로 중요하기 때문에이 질문과 답변은 잠겨 있습니다. 현재 새로운 답변이나 상호 작용을받지 않습니다. 세 개의 열로 구성된 ascii 데이터 세트가 있지만 마지막 두 개만 실제 데이터입니다. 이제를 사용하여 데이터의 도트 차트를 만들고 싶습니다 read.csv(file = "result1", sep= " "). R은 세 열을 모두 …

12 r 

자동차가 움직이지 않을 때 1과 자동차가 움직일 때 0으로 이진 시계열이 있습니다. 최대 36 시간 전과 매 시간마다 수평선을 예측하고 싶습니다. 첫 번째 접근 방식은 t-24 (일별 계절), t-48 (주간 계절), 하루 중 시간을 사용하여 Naive Bayes를 사용하는 것입니다. 그러나 결과는 그리 좋지 않습니다. 이 문제에 대해 어떤 기사 나 …

12 r  time-series  forecasting  binary-data 

R에서 Anova를 반복적으로 측정하고 그 데이터 세트에 대한 특정 대비를 수행하려고합니다. 올바른 방법은 Anova()자동차 패키지에서 사용하는 것이라고 생각합니다 . 데이터 를 ?Anova사용하여 얻은 예제로 내 질문을 설명 해 보겠습니다 OBrienKaiser(참고 : 예에서 성별 요인을 생략했습니다) : 우리는 주제 요인, 치료 (3 수준 : 통제, A, B) 및 2 반복 중 …

12 r  anova  repeated-measures  contrasts  sums-of-squares