关于R语言滞后的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请教一个问题。
我想对某一列求滞后，但是是对不同企业分别滞后。我用了lag（x，n=1L），但是这一列全滞后了，没有按企业分开。怎么解决？求帮忙。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：R语言 lag 求帮忙

建议你把数据结构补充一下方便别人回答你的问题。我这里假定你是一个data frame叫DF，你有一个包含企业名字的列，假定叫“companies", 然后你要求lag的列叫"V1",那么你可以用dplyr，先按照"companies"列分组，再进行lag，那么结果就是会对不同的企业分别lag。

1. DF %>%
   
2.   group_by(companies) %>%
   
3.   mutate(V2 = lag(V1, n=1L) ->
   
4.   DF

复制代码

suimong 发表于 2015-12-10 14:43
建议你把数据结构补充一下方便别人回答你的问题。我这里假定你是一个data frame叫DF，你有一个包含企业名字 ...

你好，非常感谢你，按照你的方法可以了。但是我是想按时间滞后的，按季度滞后，公司是可以了，但是季度滞后没有解决，出现了这种情况。我想让如果没有3月份，6月份的滞后一期就是NA,怎么同时做到这两个？？

suimong 发表于 2015-12-10 14:43
建议你把数据结构补充一下方便别人回答你的问题。我这里假定你是一个data frame叫DF，你有一个包含企业名字 ...

Stkcd          Accper             size                            V1
4        2014/3/31          253126108.2       
4        2014/6/30                275964182.4     253126108.2
4        2014/9/30                283038604.3         275964182.4
4        2014/12/31        338282581.1        283038604.3
5        2005/3/31               1800641798                     NA
5        2005/6/30          1773051158        1800641798
5        2005/9/30          1781712610        1773051158
5        2005/12/31        1643873315        1781712610
5        2006/6/30               1559309892         1643873315
5        2006/12/31        1528293792         1559309892

不知道是不是图片看不到。比如这样，2006年6月份的不应该有滞后一期的，因为没有3月份的数据。该怎么同时和按企业滞后一起做到？

suimong 发表于 2015-12-10 14:43
建议你把数据结构补充一下方便别人回答你的问题。我这里假定你是一个data frame叫DF，你有一个包含企业名字 ...

Stkcd          Accper             size                            V1
4        2014/3/31          253126108.2       
4        2014/6/30                275964182.4     253126108.2
4        2014/9/30                283038604.3         275964182.4
4        2014/12/31        338282581.1        283038604.3
5        2005/3/31               1800641798                     NA
5        2005/6/30          1773051158        1800641798
5        2005/9/30          1781712610        1773051158
5        2005/12/31        1643873315        1781712610
5        2006/6/30               1559309892         1643873315
5        2006/12/31        1528293792         1559309892

不知道是不是图片看不到。比如这样，2006年6月份的不应该有滞后一期的，因为没有3月份的数据。该怎么同时和按企业滞后一起做到？

suimong 发表于 2015-12-10 14:43
建议你把数据结构补充一下方便别人回答你的问题。我这里假定你是一个data frame叫DF，你有一个包含企业名字 ...

Stkcd          Accper             size                            V1
4        2014/3/31          253126108.2       
4        2014/6/30                275964182.4     253126108.2
4        2014/9/30                283038604.3         275964182.4
4        2014/12/31        338282581.1        283038604.3
5        2005/3/31               1800641798                     NA
5        2005/6/30          1773051158        1800641798
5        2005/9/30          1781712610        1773051158
5        2005/12/31        1643873315        1781712610
5        2006/6/30               1559309892         1643873315
5        2006/12/31        1528293792         1559309892

不知道是不是图片看不到。比如这样，2006年6月份的不应该有滞后一期的，因为没有3月份的数据。该怎么同时和按企业滞后一起做到？

这样的话你只要把日期列给补齐就可以了。假设日期从2005/3/31开始，到2014/12/31结束，那么你先生成一个日期序列：

1. # 在起始日期上先加一天是为了避免碰到小月时溢出至下个月
   
2. dt<-seq(as.Date("2005/3/31")+1,
   
3. as.Date("2014/12/31")+1, by = "1 quarter")-1
   
4. dt<-data.frame(DATE = dt)
   
5. 
   
6. # 假设你的数据为DF，日期列是Accper?
   
7. DF %>%
   
8.   right_join(dt, by = c("Accper" = "DATE")) ->
   
9. DF

复制代码

right_join保证了所有dt内的值都会出现在DF$Accper中，同时对于原来不存在的Accper中的值（例如2006/3/31)，该行的其他列就会被自动填充NA。这时候你再lag的话就会有你想要的效果了。

suimong 发表于 2015-12-11 02:52
这样的话你只要把日期列给补齐就可以了。假设日期从2005/3/31开始，到2014/12/31结束，那么你先生成一个日期 ...

谢谢。我试了一下，但是我发现并没有补齐。因为dt里的所有四十个日期在Accper里的日期都有，所有合并后没有什么变化。它并不是按照每家企业分别合并，所以原来少日期的企业依然少日期。我按企业分组了(group_by),仍然是这样。怎么解决？

我明白了，你是针对每个企业都有一组日期，范围在2005/3/31-2014/12/31，但是可能企业A缺了2014/06/30，企业B缺了2013/09/30，是不是这样？如果是这种情况的话，那就是要按照分组inner_join了，这个任务dplyr就没法胜任了。我们可以用plyr里的ddply：

1. library(plyr)
   
2. ddply(.data = DF,
   
3.       .variables = "companies",
   
4.       .fun = function(x) right_join(x, dt, by = c("Accper" = "DATE")))

复制代码

suimong 发表于 2015-12-11 16:18
我明白了，你是针对每个企业都有一组日期，范围在2005/3/31-2014/12/31，但是可能企业A缺了2014/06/30，企业 ...

您好，我生成了日期序列DATE，然后运行了下面程序：
library(plyr)
A13 <- ddply(.data = A12,
      .variables = "Stkcd",
      .fun = function(A12) right_join(A12, dt, by = c("Accper" = "DATE")))
但是出现了这样的结果：
Error in .fun(piece, ...) : could not find function "right_join"
而且把right_join改成inner_join也是错。请问怎么回事？