95％HPD Region R TeachingDemosパッケージを使用

RのTeachingDemosパッケージを使用して95％HPD領域を見つける必要があります。ガンマ分布に従う事後分布があります。95％HPD Region R TeachingDemosパッケージを使用

、ライブラリ（TeachingDemos）パッケージとタイピングをインストールした後、私はなかった：

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を使ったときには初心者ですから、どうしたらいいですか？

正確には、その意味です。ヘルプには '？hpd'をRコンソールに入力してください。これはプログラミング上の質問であり、統計的な質問ではありません。StackOverflowに 'r'タグを使ってこれらを投稿する方が良いでしょう。 – shadowtalker

@ssdecontrolエラーメッセージが何を言おうとしているのか分かりません。？hpdは私に何か新しいことを教えていません。 qgammaをqbetaに変更すると、私のコードは動作しますが、私はqgammaで動作するように自分のコードが必要です。 – sucksatmath

再現可能な例なしでRを使用する方法に関するので、この質問を議論の対象外とすることにしました。 – gung

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)

2016-11-01 03:01:51 shadowtalker

答えて