jags

0열

1답변

%dopar% 및 foreach을 사용할 예정이며 출력을 결합해야합니다. 병렬로 호출 할 함수는 각 호출에 대해 일정한 길이를 갖는 목록을 출력합니다. 그러나이 목록의 요소 길이는 항상 일정하지는 않습니다. 결합 후 결과를 최대한 단순화하고 각 요소의 원본 목록 (반복)을 식별 할 수있게하고 싶습니다. B/C 길이가 이보다 깊은 곳에서 다르다면, answe

4열

2답변

R에서 하나의 체인으로 rjags에서 여러 체인을 결합하는 방법?

진단 목적 (예 : 4 사슬)을 위해 여러 체인이 포함 된 rjags에서 일반적으로 JAGS를 호출합니다. 그 후 나는 종종 사후 매개 변수 추정치 (예 : 예측 된 값 사용, 추가 통계 계산)에 대한 사후 처리를 원합니다. 그러나이 시점에서 사슬을 목록에 저장하는 것은 불편합니다. 체인을 단일 매개 변수 목록으로 결합하는 좋은 방법은 무엇입니까?

1열

1답변

JAGS를 사용하는 베이지안 일반화 된 포아송에 대한 제로 - 하나의 트릭

제로스 트릭을 사용하여 JAGS의 일반 포아송 분포에 카운트 데이터를 맞추고 있습니다. 나는 "베이지안 모델링의 Winbugs 사용"책에 코드를 다음입니다 (286 페이지) 문서라도 내 코드입니다 : GPoisson.model <- function(){ C <- 10000 for(i in 1:N){ zeros[i] <- 0 zeros[

5열

3답변

BUGS/JAGS/STAN 언어의 새로운 함수와 배포판 정의하기

저는 매우 새로운 통계 분석 세계이며 최근에 BUGS/JAGS/STAN 모델 언어에 관심을 보였습니다. 정말로 놀라운 점은 코드 중복을 피하기 위해 정의 된 새로운 기능이나 배포판의 예는 보지 못했다는 것입니다. 예를 들어, 내가 자주 포아송 분포의 제곱을 사용한다고 가정하면, 어쨌든 다음을 수행 할 수 있습니까? dsqpo <- function(lambd

1열

2답변

JAGS : 관측 불가능한 변수와 관찰 가능한 변수 모두를 가진 다 변수 일반?

나를 미치게 만드는 JAGS에 "간단한"문제가 있습니다. u는 잠재 변수 동안 여기 x2[i] ~ dnorm(mu[i,1], tau1); u[i] ~ dnorm(mu[i,2], tau2); , X2는 관찰 변수 (즉, 데이터)이다 본질적으로, 다음 예작동을 고려한다. 이 예에서는 두 가지가 서로 다른 두 개의 정규 분포에서 독립적으로 그려집니다. 그러

1열

1답변

for 루프를 사용하여 여러 매개 변수에 대한 사전 분포 할당

JAGS 설명서를 읽었으나 JAGS/R2JAGS 모델의 여러 매개 변수에 동일한 이전 배포를 할당하는 방법을 찾지 못했습니다. 예를 들어, 현재 나는 다음과 같은 코드를 많이 반복해야 : reg.model <- function() { # Model structure for(i in 1:N){ Y[i] ~ dnorm(mu[i], p

0열

1답변

runjags 요약에서 psrf 값이 너무 낮습니까?

Runjags 분명히 컨버전스 문제가 체인에 대한 psrf = 1.0047 매우 낮은보고 : > print(o, vars = "q_date2") JAGS model summary statistics from 3000 samples (chains = 3; adapt+burnin = 1000): Lower95 Median Upper95 Mea

0열

1답변

JAGS에 대한 변경 내역이 있습니까?

JAGS에 대한 변경 내역이 있습니까? 나는 JAGS 홈페이지와 소스 포지를 둘러 보았고 아무것도 발견하지 못했다. 웹 또는 릴리스에는 변경 로그가 없습니다

0열

1답변

mcmc RJAGS에서 데이터를 추출하는 방법

각 폴링 주에 대해 측정 한 다음 모델을 실행 중입니다. 각 파티마다 100 개 이상의 알파 벡터를 생성합니다. 질문 할 수 있습니다. 선 그래프를 그릴 수있는 방법으로 어떻게 비축 할 수 있습니까? model{ ## measurement for(i in 1:NPOLLS){ p1[i] ~ dnorm(alpha1[WoY[i]] + pol

1열

2답변

베타 배포판을 사용하는 그룹 평균을 산정하기위한 JAGS 코드

JAGS를 사용하여 13 개 장소 (새들은 9 개, 잠재적 서식지는 4 개)에 대한 캐노피 비율의 평균과 SD를 추정하고 싶습니다. 데이터가 0과 1로 묶여 있다는 사실을 설명하기 위해 베타 배포판을 사용하고 있습니다. 다른 배포판 (포아송 및 로그 - 표준)에서 완벽하게 작동하는 모델 문에 대한 코드가 있으며 시도 중입니다. 그 코드를 적응 시키는데 나는