Question

我正在研究营销组合建模，并且正在关注文章

https://analyticsartist.wordpress.com/2014/01/31/adstock-rate-deriving-with-analytical-methods/

本文定义了adstock函数，如下所示：

adstock <- function(x, rate=0){
  return(as.numeric(filter(x=x, filter=rate, method="recursive")))
}

并进一步使用R中nlsm包中的minpack.lm来计算比率和系数。

model1 <- nlsLM(Applications~b0 + b1 * adstock(Media1, r1) + b2 * adstock(Media2, r2) +
                  b3 * adstock(Media3, r3) + b4 * adstock(Media4, r4) + b5 * adstock(Media5, r5) +
                  b6 * adstock(Media6, r6) + b7 * adstock(Media7, r7),
                algorithm = "LM",
                start     = c(b0=   1, b1=   1, b2=   1, b3 = 1, b4 = 1, b5 =1, b6= 1, b7= 1, r1=0, r2=0, r3=0, r4=0, r5=0, r6=0, r7=0),
                lower     = c(b0=-Inf, b1=-Inf, b2=-Inf, b3 = -Inf, b4 = -Inf, b5 =-Inf, b6= -Inf, b7= -Inf, r1=0, r2=0, r3=0, r4=0, r5=0, r6=0,     r7=0),
                upper     = c(b0= Inf, b1= Inf, b2= Inf, b3 = Inf, b4 = Inf, b5 =Inf, b6= Inf, b7= Inf, r1=0.5, r2=0.5, r3=0.5, r4=0.5, r5=0.5, r6=0.5, r7=0.5))

但是，模型不断失败，并显示以下错误

Error in filter_(.data, .dots = compat_as_lazy_dots(...)) : 
argument ".data" is missing, with no default

似乎错误来自adstock函数，但我不确定如何解决。

我真的很希望有人能帮助解决这个问题。

非常感谢！

Answer 1

（这是一个常见的问题，但是由于我找不到重复项，因此我现在提供一个答案。）

您在这里看到的错误来自dplyr::filter，而不是您期望使用的错误：stats::filter。加载dplyr时，您应该已经看到以下内容：

library(dplyr)
# Attaching package: 'dplyr'
# The following objects are masked from 'package:stats':
#     filter, lag
# The following objects are masked from 'package:base':
#     intersect, setdiff, setequal, union

使用非基本函数时，要避免这种情况（在将程序包发布到CRAN时鼓励/强制使用）是明确的。我通常以为stats::可以不受此影响，但是使用dplyr无疑是必须的。

因此，在filter附近的任何地方使用dplyr时，只需简单地明确代码的正确性：

adstock <- function(x, rate=0){
  return(as.numeric(stats::filter(x=x, filter=rate, method="recursive")))
}

FWIW，R的名称空间管理和与python更明确的方法的大致等效性：

R                          Python
----------------------     ----------------------
                           import pkgname         | explicit namespace use
pkgname::function(...)     pkgname.function(...)  |

                           import pkgname as p    | no R equivalent?
                           p.function(...)        |

library(pkgname)           import * from pkgname  | permissive namespace,
function(...)              function(...)          |   enables masking

过滤器功能不断失败，缺少数据参数

1 个答案: