0

我有一个数据框，其中包含来自 COMPUSTAT 的变量，其中包含各种会计项目的数据，包括来自不同公司的 SG&A 费用。

我想在数据框中创建一个新变量，该变量按时间顺序累积每家公司的 SG&A 费用。我使用 PERMNO 代码作为每个公司的唯一 ID。

我已经尝试过这段代码，但它似乎不起作用：

crsp.comp2$cxsgaq <- crsp.comp2 %>%
  group_by(permno) %>%
  arrange(date) %>%
  mutate_at(vars(xsgaq), cumsum(xsgaq))

（xsgag 是 SG&A 费用的 COMPUSTAT 变量）

非常感谢您的帮助

4

2 回答 2

3

您的示例代码正在尝试将整个数据帧写入crsp.comp2变量crsp.comp2$cxsgaq。

通常vars()函数变量需要“引用”；尽管在您的情况下，请使用标准mutate()函数并在那里分配cxsgaq变量。

crsp.comp2 <- crsp.comp2 %>%
  group_by(permno) %>%
  arrange(date) %>%
  mutate(cxsgaq = cumsum(xsgaq))

鸢尾花数据集的可重现示例：

library(tidyverse)
iris %>% 
  group_by(Species) %>% 
  arrange(Sepal.Length) %>% 
  mutate(C.Sepal.Width = cumsum(Sepal.Width))
于 2021-02-26T14:21:45.177 回答
1

基于@m-viking 的答案，如果使用 WRDS PostgreSQL 服务器，您只需使用window_order(from dplyr) 代替arrange. （我使用 Compustat 公司标识符gvkey代替，permno以便此代码有效，但想法是相同的。）

library(dplyr, warn.conflicts = FALSE)
library(DBI)

pg <- dbConnect(RPostgres::Postgres(), 
                bigint = "integer", sslmode='allow')

fundq <- tbl(pg, sql("SELECT * FROM comp.fundq"))

comp2 <-
  fundq %>%
  filter(indfmt == "INDL", datafmt == "STD",
         consol == "C", popsrc == "D")

comp2 <- 
  comp2 %>%
  group_by(gvkey) %>%
  dbplyr::window_order(datadate) %>%
  mutate(cxsgaq = cumsum(xsgaq))

comp2 %>%
  filter(!is.na(xsgaq)) %>%
  select(gvkey, datadate, xsgaq, cxsgaq)
#> # Source:     lazy query [?? x 4]
#> # Database:   postgres [iangow@wrds-pgdata.wharton.upenn.edu:9737/wrds]
#> # Groups:     gvkey
#> # Ordered by: datadate
#>    gvkey  datadate   xsgaq cxsgaq
#>    <chr>  <date>     <dbl>  <dbl>
#>  1 001000 1966-12-31 0.679  0.679
#>  2 001000 1967-12-31 1.02   1.70 
#>  3 001000 1968-12-31 5.86   7.55 
#>  4 001000 1969-12-31 7.18  14.7  
#>  5 001000 1970-12-31 8.25  23.0  
#>  6 001000 1971-12-31 7.96  30.9  
#>  7 001000 1972-12-31 7.55  38.5  
#>  8 001000 1973-12-31 8.53  47.0  
#>  9 001000 1974-12-31 8.86  55.9  
#> 10 001000 1975-12-31 9.59  65.5  
#> # … with more rows

由reprex 包于 2021-04-05 创建(v1.0.0)

于 2021-04-05T23:44:19.457 回答