2

使用 data.table,我正在尝试编写一个函数,该函数将数据表、公式对象和字符串作为参数,并创建和存储多个模型对象。

myData <- data.table(c("A","A","A","B","B","B"),c(1,2,1,4,5,5),c(1,1,2,5,6,4))
## This works.
ModelsbyV1 <- myData[,list(model=list(lm(V2~V3)),by=V1)]

##This does not.
SectRegress <- function (df,eq,sectors) {
  Output <- df[,list(model=list(lm(eq))),
             by=sectors]
  return(Output)
}

Test <- SectRegress(myData,formula(V2~V3),sectors="V1")
##Error in eval(expr, envir, enclos) : object 'X' not found

我已经尝试在函数中附加 df 。但是,这会使按类型分组的能力无效。函数内部的 colnames(df) 包括“X”。我难住了。

4

1 回答 1

2

您必须在环境中对其进行评估.SDlm否则无法“看到”V2 和 V3):

SectRegress <- function (df,eq,sectors) {
    Output <- df[, list(model=list(lm(eq, .SD))), by=sectors]
    return(Output)
}
Test <- SectRegress(myData,formula(V2~V3),sectors="V1")
于 2013-04-26T09:14:15.377 回答