저는 매우 새로운 통계 분석 세계이며 최근에 BUGS/JAGS/STAN 모델 언어에 관심을 보였습니다. 정말로 놀라운 점은 코드 중복을 피하기 위해 정의 된 새로운 기능이나 배포판의 예는 보지 못했다는 것입니다. 예를 들어, 내가 자주 포아송 분포의 제곱을 사용한다고 가정하면, 어쨌든 다음을 수행 할 수 있습니까? BUGS/JAGS/STAN 언어의 새로운 함수와 배포판 정의하기
dsqpo <- function(lambda) {
tmp ~ dpois(lambda)
tmp2 <- tmp * tmp
return(tmp2)
}
하고 나중에 새로운 임시 변수 각 시간을 정의하지 않고
model{
...
x ~ dsqpo(alpha)
y ~ dsqpo(beta)
}
에
.
저는 학계가 아닌 사람들을 위해서이 튜토리얼을 '사야'합니다. –
저자에게 이메일을 보내거나 scihub에서 검색 할 수 있습니다. –