多维数组的伽马分布

Question

我有一个CSV，其中包含1000行和5个实数列，它们同时代表5个不同的测量值。

我被要求使用R执行类似

的操作

基于该值创建伽马函数，并计算\ mu和\ gamma。

使用此伽玛系数评估选择新的五个测量值的概率是多少，从而计算出新的测量值与总体之间的差异大于某个误差的概率。

直观地让$ x_ {n + 1} $为新度量$ P（X-x_n）<\ epsilon $。

我发现使用堆栈如何使用R中的gamma函数对单个测量结果进行存档，但是，我不知道如何概括。 在堆栈中的代码下面

x=rgamma(20,shape=1, rate=2) #shape=alpha, #rate=beta. 
#x = (sample.int(101,size=100,replace=TRUE)-1)/100
dloglik=function(alpha){
  return(log(alpha)-digamma(alpha)-log(mean(x))+mean(log(x)))
}
alpha=uniroot(dloglik, lower=0.001, upper=1000)$root
beta=alpha/mean(x)

有人可以帮忙吗？

保罗