首页 > 解决方案 > 如何计算 r 的滞后投资

问题描述

我创造了df1

year gvkey  capex   ppent
2004  1004 13.033 139.137
2005  1004 16.296 213.380
2006  1004 29.891 260.167
2007  1004 30.334 310.393
2008  1004 27.535 245.586
2009  1004 28.855 334.430
...

我创造了df2

year gvkey    ROA
2005  1004 0.02796478
2006  1004 0.04665171
2007  1004 0.05976127
2008  1004 0.06255035
2009  1004 0.03549220
2005  1013 0.06882688
...

我想创造df3

year gvkey    ROA               lag_investment
2005  1004 0.02796478  capex from 2004 / ppent from 2004
2006  1004 0.04665171  capex from 2005 / ppent from 2005
2007  1004 0.05976127  capex from 2006 / ppent from 2006
2008  1004 0.06255035  capex from 2007 / ppent from 2007
2009  1004 0.03549220  capex from 2008 / ppent from 2008
2005  1013 0.06882688  capex from 2004 / ppent from 2004
...

我有超过 2,000 家公司年。gvkey = firm id

我基本上想做的是以下几点:

1) 计算上一年的投资df1

2)创建一个名为“lag_investment”的列df2

2) 将步骤 1) 中的值插入到当前年份行中df2

附加问题:

如果我想执行以下操作,代码会是什么样子?

我创造了 df1

  year gvkey        ROA   ppent  capex
1 2004  1004 0.01320911 139.137 13.033
2 2005  1004 0.03005708 213.380 16.296
3 2006  1004 0.05014214 260.167 29.891
4 2007  1004 0.06423255 310.393 30.334
5 2008  1004 0.06723031 245.586 27.535
6 2009  1004 0.03814769 334.430 28.855
...

我想添加一个变量df1

  year gvkey        ROA   ppent  capex         lag_investment
1 2004  1004 0.01320911 139.137 13.033
2 2005  1004 0.03005708 213.380 16.296  capex from 2004 / ppent from 2004
3 2006  1004 0.05014214 260.167 29.891  capex from 2005 / ppent from 2005
4 2007  1004 0.06423255 310.393 30.334  capex from 2006 / ppent from 2006
5 2008  1004 0.06723031 245.586 27.535  capex from 2007 / ppent from 2007
6 2009  1004 0.03814769 334.430 28.855  capex from 2008 / ppent from 2008
...

我想计算除 2004 年以外所有年份的 lag_investment。

非常感谢!!!

标签: rtransformlag

解决方案


我想你可以lag使用dplyr

library(dplyr)
df1 %>% mutate(lag_investment = lag(capex)/lag(ppent))

#  year gvkey    ROA ppent capex lag_investment
#1 2004  1004 0.0132   139  13.0             NA
#2 2005  1004 0.0301   213  16.3         0.0937
#3 2006  1004 0.0501   260  29.9         0.0764
#4 2007  1004 0.0642   310  30.3         0.1149
#5 2008  1004 0.0672   246  27.5         0.0977
#6 2009  1004 0.0381   334  28.9         0.1121

如果数据框未排序,arrange请先使用按年排序。

df1 %>% arrange(year) %>% mutate(lag_investment = lag(capex)/lag(ppent))

shiftdata.table

library(data.table)
setDT(df1)[, lag_investment := shift(capex)/shift(ppent)]

推荐阅读