r - 简单函数给出的结果为 NaN？

Question

我写了一个简单的函数，但我想知道为什么我会得到 NaN 作为结果。这是代码：

数据

      en1=c(1,9,5,9,8,7)
      en=c(1,5,4,9,3,7)
      en2=c(1,4,5,4,8,7)

功能：

跑：

        err1(en,en1,en2)

结果：

         [1] NaN
           Warning message:
       In sqrt(mean(((pr - ob) * (ob - ref)), na.rm = na.rm)) : NaNs produced

请对此有任何想法

score 5 · Accepted Answer

NaN由于负数的平方根未定义（对于实数），因此您获得。

您可能希望abs在您的函数中使用：

err1 <- function(pr,ob,ref,na.rm=TRUE){
    sqrt(mean(abs((pr-ob)*(ob-ref)), na.rm = na.rm))}     

err1(en,en1,en2)
[1] 1.825742

1 回答 1