问题描述
我正在尝试使用 R 中的 dgamma 函数来获取值矩阵中的 gamma 密度函数值。因为我的值通常很小,所以我也需要在日志空间中做这个过程,以避免函数内的下溢。
例如:
values = matrix(data = seq(1:9),nrow = 3,ncol = 3)
dgamma(values,shape = 10,scale = 1/10,log = T)
# [,1] [,2] [,3]
#[1,] 0.2240234 -17.29933 -42.26279
#[2,] -3.5376519 -25.29104 -51.06100
#[3,] -9.8884660 -33.65014 -60.00096
但是,因为我在一个具有相当大维度(通常约为 1000x10000)的矩阵上执行此操作,并且对不同值的矩阵重复执行此操作,所以此过程是我代码中的瓶颈。
我想知道是否有人有任何聪明的技巧来加速 R 中的这个过程?
谢谢。
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)