我想改变的精度在R的计算.例如,我想计算x^6与x = c(-2.5e+59, -5.6e+60).为了计算它我应该改变R的精度,否则结果是Inf,我不知道该怎么做.

我保证这不仅仅是另一个掷骰子的家庭作业问题。我实现了一个函数来计算在掷骰子s时获得小于总和的概率n m。我的函数适用于 的小值，n但我发现大值的n. 见附图。任何人都了解发生了什么？

我的概率函数

从此数学堆栈交换中实现

probability <- function(s, m, n) {

  i <- 0:((s-1-n) / m)
  m^(-n) * sum((-1)^i * choose(n, i) * choose(s - 1 - i * m, n))

}

开始打破 ~ n > 80

n <- 1:90 # number of dice
m <- 6 # number of sides
s <- floor(mean(1:m)*n) # sum of faces
p <- mapply(probability, s = s, m = m, n …

我正在阅读,我自己(不是HW)关于编程,一个练习涉及在R中编写Pascal的三角形.我的第一个想法是制作一个列表,然后将其添加到其中,但这并不能很好地工作.然后我想到了从一个向量开始,并在最后制作一个列表.然后我想到制作一个矩阵,并在最后制作一个列表.

不知道哪种方式可以解决这个问题.

任何提示？

谢谢