95 % HPD 지역 R TeachingDemos 패키지 사용

Question

R에서 TeachingDemos 패키지를 사용하여 95 % HPD 지역을 찾아야합니다. 나는 감마 분포를 따르는 사후 분포를 가지고 있습니다.

라이브러리 (TeachingDemos)에서 압축하고 입력을 설치 한 후에 I는 않았다

a = 200 
b = 20 
hpd(qgamma,shape1=a,shape2=b, conf=0.95) 

(a 및 b는 감마 분포의 알파 및 베타 값이다)

I을 코드를 실행할 때 계속해서 다음과 같은 오류 메시지가 나타납니다.

Error in posterior.icdf(1 - conf + x, ...) : 
    unused arguments (shape1 = 200, shape2 = 20) 

나는 R을 사용하는 데 초보자입니다. 그래서 내가 뭘 잘못하고 있니?

Answer 1

hpd의 서명은

hpd(posterior.icdf, conf=0.95, tol=0.00000001,...) 

입니다 그리고 문서 ...이라고 주장한다 "posterior.icdf에 전달 된 추가 인수." 그 어딘가에 hpd 안에 인수가 hpd`` apart from conf의 and tol`에 전달 된 어떤 ... 단지입니다

posterior.icdf(x, ...) 

같이 보입니다 코드의 라인이, 있다는 것을 의미한다.

qbeta의 함수 시그니처를 보면 shape1 및 shape2이라는 인수가 있음을 알 수 있습니다. 반면에 qgamma을 보면 이 아니고 유효한 인수 이름이 인 것을 알 수 있습니다.

qgamma과 같은 인수를 직접 호출하면 동일한 오류 메시지가 표시됩니다.

qgamma(0.5, shape1=200, shape2=20)