JAGS에서 검열 문제를 해결하는 방법에 대한 질문이 있습니다.
X 값에 측정 오류가있는 이변 량 혼합 법선이 관찰됩니다. 관측 된 검열 된 값의 진정한 기본 '평균'을 모델링하고 싶습니다.
여기 내가 지금 가진 것입니다 :
for (i in 1:n){
x[i,1:2]~dmnorm(mu[z[i],1:2], tau[z[i],1:2,1:2])
z[i]~dcat(prob[ ])
}
Y에도 측정 오류가 있습니다. 내가하고 싶은 것은 다음과 같습니다.
for (i in 1:n){
x_obs[i] ~ dnorm(x_true[i],prec_x)I(x_true[i],)
y_obs[i] ~ dnorm(y_true[i],prec_y)
c(x_true[i]:y_true[i])~dmnorm(mu[ z [ i ],1:2], tau[z[i],1:2,1:2])
z[i]~dcat(prob[ ])
}
#priors for measurement error
e_x~dunif(.1,.9)
prec_x<-1/pow(e_x,2)
e_y~dunif(2,4)
prec_y<-1/pow(e_y,2)
분명히 c 명령은 JAGS에서 유효하지 않습니다.
미리 감사드립니다.
3
자르려면 T (-,-)를 사용하되, censuring 및 truncationq에 대한 정보는 사용자 매뉴얼을 읽으십시오
—
David LeBauer