R：数学函数的非数字参数

Question

我写了一个用户定义的函数：

epf <- function(z,x,noise=std_noise){
    z_dims <- length(z)
    std_noise <- 0.5*matrix(1,1,z_dims)
    std_noise <- as.data.frame(std_noise)

    obs_prob <- dnorm(z,x[1:z_dims],noise)

    error <- prod(cbind(1,obs_prob))
    return(error)
}

这个函数在另一个函数的for循环中被调用：

w <- matrix(0,N,1)
        for (i in 1:N){
            w[i] <- epf(z,p[i,],R_noise) 
        }

其中z 是一个二维向量，N=1000、p 是一个包含 1000 个观察值和 4 个变量的数据框，R_noise 是一个包含 1 个观察值和 4 个变量的数据框。

这里我得到错误：“数学函数的非数字参数”，对于obs_prob <- dnorm(z,x[1:z_dims],noise)这一行

谁能帮我找出错误？

我已经查看了与我类似的问题，但我仍然无法在我的代码中找到错误。

编辑：
增加了N的定义

Answer 1

dnorm(as.matrix(z), x[1:a_dims], noise) 可能会更好。

更广泛地说，一行两列的数据框可能更好地表示为向量。数据框看起来像矩阵，正如您所说的“二维向量”，但它们在重要方面有所不同。

可能会发生相同的错误，因为您通过传递 R_noise 在其最后一个参数 noise 中提供了第二个数据帧 dnorm。

另外，考虑p[i, ] 有四个值。它被obs_prob 和x[1:z_dims] 子集化。在这种情况下，z_dims 将等于 2，因为 length(z) 是 2。所以您正在评估 dnorm(data.frame(z), p[1, ][1:2], data.frame(R_noise))。