各要素のガンマ関数を含む R の合計を計算する必要があります。ガンマの引数が増加すると、結果として NaN が得られます。問題はガンマの評価に関する数値的なものであると思われます。Rmpfr と gmp のドキュメントは既に読みましたが、整数の階乗しか見つかりませんでした。また、ここにコードを投稿します。エラーの原因についてより良い考えがあるかもしれません。
V1<-w1*(b1^a1)/gamma(a1)
VV1<-outer(V1,V1)
a1mat<-outer(a1, a1-1, FUN="+")
b1mat<-outer(b1, b1, FUN="+")
A<-sum(VV1*gamma(a1mat)/(b1mat^a1mat))
w1 は合計が 1 になる正の実数の配列で、a1 と b1 は正の値のベクトルです。a1 (および a1mat) が長くなり、高い値 (~150) になると、A は NaN になります。