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

忽略inf值并运行lm回归

Laura · 技术社区 · 5 年前

这些是我的变量:

> dput(y)
c(-22.0713165394207, 14.0880914427811, 10.9650636244176, -1.96648890706268, 
-5.30593850426708, -7.54651916037787, 3.84914747321197, 4.4986386904214, 
1.73067625014435, 2.5585960595839, -2.72766183793304, -3.10167452216202, 
2.68853838208521, 1.12662203717498, 1.24951279248057, 3.70075666289518, 
-6.11243972144607, -6.91019769671849, 6.46767794752582, 8.84874735514293, 
2.95606352319898, 3.23883851668917, -2.61692776879569)
> dput(x)
c(`1` = 0.0520523266234464, `2` = Inf, `3` = 0.0520523266234462, 
`4` = 0.0520523266234463, `5` = 0.0520523266234463, `6` = 0.0520523266234461, 
`7` = 0.0520523266234463, `8` = 0.0520523266234466, `9` = 0.0520523266234465, 
`10` = 0.0520523266234465, `11` = 0.0520523266234465, `12` = 0.0520523266234466, 
`13` = 0.0520523266234468, `14` = 0.0520523266234466, `15` = 0.0520523266234467, 
`16` = 0.0520523266234464, `17` = 0.0520523266234463, `18` = 0.0520523266234465, 
`19` = 0.0520523266234466, `20` = 0.0520523266234463, `21` = 0.0520523266234464, 
`22` = 0.0520523266234465, `23` = 0.0520523266234464)

当我运行回归时: summary(lm(ex.return ~ ex.return.skew))

我收到以下错误消息:

Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : 
  NA/NaN/Inf in 'x'

我试着删除position2,然后运行回归,但我不能这样做,因为数据要大得多。所以我正在寻找一种忽略inf/na值并运行回归的方法。

我该怎么做?

有什么帮助吗?

2 回复 | 直到 5 年前

akrun 5 年前

我们可以把无限值转换成 NA 它应该会起作用的

x[is.infinite(x)] <- NA
summary(lm(y ~ x))

rookie 5 年前

您需要删除x和y中的无限项,如下所示:

summary(lm(ex.return[is.finite(df$ex.return)] ~ ex.return.skew[is.finite(df$ex.return)]))

但是,更好的方法是将它们放入data.frame中,并将该data.frame添加到lm的data参数中,该参数过滤掉data.frame的行。

df <- data.frame(ex.return, ex.return.skew)
summary(lm(ex.return ~ ex.return.skew, df[is.finite(df$ex.return),]))

注意,is.finite()适用于na值和-inf/inf

is.finite(c(NA, Inf, 10))
[1] FALSE FALSE  TRUE

如果在data.frame中的任何列的任何行中(即不只是在ex.return中)都可能有inf/-inf和na,则可能需要执行以下操作:

summary(lm(ex.return ~ ex.return.skew, df[is.finite(rowSums(df)),]))

推荐文章

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)”?

1 年前

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

1 年前

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

2 年前

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

2 年前