R：为什么这不起作用？矩阵、舍入误差？

Question

为什么这有效：

ncota <- 1
nslope <- 29
resul <- matrix(rep(0,ncota*nslope*4),ncota*nslope,4)

但这不起作用？

ncota <- 1
sini <- 0.1; sfin <- 1.5; spaso <- 0.05; nslope <- 1+((sfin-sini)/spaso)
resul <- matrix(rep(0,ncota*nslope*4),ncota*nslope,4)

我猜问题是除法给出了一个非整数。 我怎样才能让第二个工作？ 我需要创建一个零矩阵，其大小是通过方程计算计算出来的。

干杯

Answer 1

如果您所要做的就是创建一个零矩阵，则不需要提供正确数量的零，只需提供一个并让 R 将其回收到所需的长度：

matrix(0, ncota*nslope, 4)

第二个失败的原因是 ncota * nslope * 4 不完全是 116：

> (ncota * nslope * 4) == 116
[1] FALSE
> all.equal(ncota * nslope * 4, 116)
[1] TRUE

all.equal 显示，如果您允许浮点误差，则这些是相等的。

?rep 包括以下内容：

 Non-integer values of ‘times’ will be truncated towards zero.  If
 ‘times’ is a computed quantity it is prudent to add a small fuzz.

如果我们按照它所说的去做并添加一个小绒毛，rep 确实会给出所需的 0 数量：

> length(rep(0, times = ncota*nslope*4 + 0.00000001))
[1] 116

正如 Hadley 所指出的（在评论中） ），可以使用 zapsmall 函数轻松添加此模糊：

> length(rep(0, times = zapsmall(ncota*nslope*4)))
[1] 116

Answer 2

您不需要使用rep。这工作得很好：

resul <- matrix(0,ncota*nslope,4)