문제
정보 및 사전 정보를 바탕으로 사후 밀도를 추정하기 위해 베이지안 분석을 수행하는 R 함수를 작성 중입니다. 사용자가 이전을 재고해야 할 경우 경고를 보내는 기능을 원합니다.
이 질문에서 나는 사전 평가 방법을 배우고 싶습니다. 이전의 질문들은 정보화 된 이전의 내용 ( 여기 및 여기 ) 을 다루는 메커니즘을 다루었 다 .
다음과 같은 경우 이전의 재평가가 필요할 수 있습니다.
- 데이터는 이전을 진술 할 때 고려되지 않은 극단적 인 경우를 나타냅니다
- 데이터 오류 (예 : 이전이 kg 일 때 데이터가 g 단위 인 경우)
- 코드의 버그로 인해 사용 가능한 사전 설정에서 잘못된 사전이 선택되었습니다.
첫 번째 경우, 데이터 값이 지원되지 않는 범위 (예 : logN 또는 감마의 경우 <0)에 있지 않는 한, 우선 순위는 일반적으로 데이터가 일반적으로 데이터를 압도 할 정도로 충분히 확산됩니다. 다른 경우는 버그 또는 오류입니다.
질문
- 데이터를 사용하여 이전을 평가하는 유효성 에 관한 문제가 있습니까?
- 이 문제에 가장 적합한 특정 테스트가 있습니까?
예
파란색 데이터는 유효한 사전 + 데이터 조합 일 수 있지만 빨간색 데이터는 음수 값을 지원하는 사전 분포가 필요합니다.
set.seed(1)
x<- seq(0.01,15,by=0.1)
plot(x, dlnorm(x), type = 'l', xlim = c(-15,15),xlab='',ylab='')
points(rnorm(50,0,5),jitter(rep(0,50),factor =0.2), cex = 0.3, col = 'red')
points(rnorm(50,8,0.5),jitter(rep(0,50),factor =0.4), cex = 0.3, col = 'blue')