代码之家 › 专栏 › 技术社区 › DanY

为什么我要和ifelse一起得到“警告南生产”?

DanY · 技术社区 · 6 年前

我有一个非常简单的分段函数,它接受向量参数 omega

btilde <- function(omega) {
        if(min(omega)<=0) stop("ALL ELEMENTS OF OMEGA MUST BE POSITIVE")
        ifelse( 0 < omega & omega <=  1, 0.233*omega^2,
        ifelse( 1 < omega & omega <=  3, 0.005*omega^4 - 0.06906*omega^3 + 0.3167*omega^2 - 0.02326*omega,
        ifelse( 3 < omega & omega <= 40, 0.705*sqrt(omega)*log(omega),
        ifelse(40 < omega,               sqrt(3*omega*log(omega) - log(8*pi) - 2/log(omega) - 170),
        NA))))
    }

我得到一个“警告”如果 omega=c(1,41) 但如果 omega=1 或 omega=41 .

> btilde(1)
[1] 0.233
> btilde(41)
[1] 16.8228
> btilde(c(1,41))
[1]  0.2330 16.8228
Warning message:
In sqrt(3 * omega * log(omega) - log(8 * pi) - 2/log(omega) - 170) :
  NaNs produced

NaN ifelse()

2 回复 | 直到 6 年前

mickey 6 年前

从 ifelse 帮助页,

当且仅当测试的任何元素为真时,将评估“是”,类似地,将评估“否”。

FALSE ,如果任何其他元素 TRUE

x = c(0, -1) 
ifelse(x == 0, sqrt(x), x)

[1]  0 -1
Warning message:
In sqrt(x) : NaNs produced

它会产生与代码类似的警告。

Mohammed Ali 6 年前

我建议你关注算法中的一个元素,然后 map 你的函数被添加到一个列表中。

library(purrr)

btilde <- function(omega) {
  if (omega <= 0) stop("ALL ELEMENTS OF OMEGA MUST BE POSITIVE")
  if (omega <= 1) return (0.233*omega^2)
  if (omega <= 3) return (0.005*omega^4 - 0.06906*omega^3 + 0.3167*omega^2 - 0.02326*omega)
  if (omega <= 40) return(0.705*sqrt(omega)*log(omega))
  return(sqrt(3*omega*log(omega) - log(8*pi) - 2/log(omega) - 170))
}

输出如下:

> map_dbl(c(1,41), ~ btilde(.x))
[1]  0.2330 16.8228

推荐文章

Marc B. · 使用ggplot2创建条形图时“缺少值”

1 年前

deschen · tidyverse与外部向量发生突变,该外部向量的元素是数据帧中的列值

1 年前

Laura · 在Shiny中使用可排序的包拖放名称,这些名称将成为图表

1 年前

Mallikarjun M · 如何使用随机森林进行时间序列预测?

1 年前

ly li · 模型摘要:当表格形状改变时,拟合优度消失

1 年前

C.Robin · 将marginaffects::predictions()的结果连接回main df?

1 年前

monotonic · 如何将格式为“col1+col3+col4”的数据帧的行名转换为一列数字向量“c(1,3,4)”?

2 年前

Shawn Hemelstrand · 为什么我的自定义errorbar函数不能在R中工作?

2 年前

RoyBatty · 统计每个字符在整个数据集中出现的次数

2 年前

stats_noob · R: 记录某个“行为”发生的循环的索引?

2 年前