R语言数据处理的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

有一个表格df1，包括id, t, dm, s这4列。其中

* id有重复。这里只列出id = 1,2的情况，但是真实数据中还有很多，id=1,2,...,100.

* t可以理解为时间，1-7天，或者1-7月等等。

* dm是另一个参数，要根据dm做一些lag。

* s是一些数值。

df1

根据这些来计算df2中x的值。具体算法如下：

当id = 1时，dm = 1。x的值比s的值滞后(lag)一行。

当id = 2时，dm有两个值(2和3)。

当dm = 2时， x的值比s的值滞后(lag)两行。

当dm = 3时， x的值比s的值滞后(lag)三行。

其中第12行的x值为0

df2

这个如何实现呢？

我已经试过left join,shift等办法，都没有搞定。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：数据处理 R语言 left join shift left

你的描述并不清楚，为什么dm = 1，你是所有的数据都往后移一格，但是dm = 2的时候只有前两个数据往后移两格，而不是所有数据都往后移两格，如果id = 2, dm = 1，3或者1，2的话又会是怎么样？

谢谢关注。你说的这两种情况我在下面列出来了。

当id = 2, dm = 1，3时，x的值为0 4 8 0 0 3 5.

当id = 2, dm = 1，2时，x的值为0 4 8 0 3 5 45

wx2123 发表于 2022-2-25 10:17
谢谢关注。你说的这两种情况我在下面列出来了。

当id = 2, dm = 1，3时，x的值为0 4 8 0 0 3 5.

实际上你并没有回答我的问题，只是把其中的一种结果列出来了而已，关键是

dm = 2的时候只有前两个数据往后移两格，而不是所有数据都往后移两格

这个问题本来是要了解x移动的规则，但你现在只列出了两种情况，dm = {1,1,2,2,2,2,2} 和 dm = {2,2,3,3,3,3,3}
那如果是dm = {1,1,1,1,1,2,2}，dm = {2,2,2,2,2,1,1}, dm = {1,1,1,1,1,3,3}, dm = {3,3,3,3,1,1,1}，dm = {1,2,3,3,3,3,3,3}, dm = {3,3,3,2,2,2,1}这些情况呢，dm的取值和组合有什么限制吗？
请把你的移动规则描述清楚（或者设计清楚），否则确实不知道你想干什么。

根据你目前的描述，可以这样处理，当然，随着描述的进一步准确，潜在的bug可以修复

1. library(dplyr)
   
2. df1 <- data.frame(id = rep(1:2, each = 7),
   
3.                   t = rep(1:7, times = 2),
   
4.                   dm = c(rep(1,7), rep(1,2), rep(2,5)),
   
5.                   s = c(15,66,77,80,97,26,15,4,8,3,5,45,85,35))
   
6. 
   
7. dfs <- split(df1, df1$id)
   
8. 
   
9. x_val <- function(dat){
   
10.   
    
11.   rle_dm <- rle(dat$dm)
    
12.   l <- rle_dm$lengths[1]
    
13.   n <- rle_dm$values
    
14.    
    
15.   if(length(n) == 1)
    
16.     dat <- dat %>% mutate(x = lag(s, default = 0))
    
17.   
    
18.   else{
    
19.     s <- dat$s
    
20.     x1 <- lag(s, n = n[1], default = 0)
    
21.     x2 <- lag(s[-c(1:l)], n = n[2], default = 0)
    
22.     x <- c(head(x1, n[1] + l), tail(x2, length(s) - n[1] - l))
    
23.     dat <- dat %>% mutate(x)
    
24.   }
    
25. }
    
26. 
    
27. dfx <- lapply(dfs, x_val)
    
28. (df <- bind_rows(dfx))

复制代码

感谢回复，我先试试看。然后在想想lag的规则到底是什么。我也在摸索。

owenqi 发表于 2022-2-26 03:36
根据你目前的描述，可以这样处理，当然，随着描述的进一步准确，潜在的bug可以修复

你说的对，我也在探索x移动的规则。而且我可能在误导你，因为我自己说的也不一定对。
这个问题可能没有那么复杂。

不知道你SAS熟不熟，我就是想产生下面SAS程序的结果：

1. data One;
   
2.   input  id t dm e s f;
   
3.   cards;
   
4. 1 1  1  10 15 0                                                                                                               
   
5. 1 2  1  20 66 0                                                                                                         
   
6. 1 3  1  30 77 0                                                                                                            
   
7. 1 4  1  40 80 0
   
8. 1 5  1  40 97 0   
   
9. 1 6  1  40 26 0
   
10. 1 7  1  10 15 0                                                                                          
    
11. 2 1  1  90 4  2
    
12. 2 2  1  80 8  1
    
13. 2 3  2  66 3  2
    
14. 2 4  2  44 5  1
    
15. 2 5  2  90 45  2
    
16. 2 6  2  80 85  1
    
17. 2 7  2  66 35  2                                                                                                            
    
18. ;
    
19. run;
    
20. 
    
21. data Two;
    
22.   set One;
    
23.   by id;
    
24. 
    
25.   array o{0:150}  _temporary_;
    
26.   array dc{0:150} _temporary_;
    
27.        
    
28.   if first.id then do;
    
29.     do i = 0 to 10;
    
30.         dc[i] =0;
    
31.         o[i] = 0;
    
32.     end;
    
33.   end;
    
34.        
    
35.         /* 主要问题在这里 */

复制代码

输出结果如下。注意每个id的最后一行并非是移动，而是要求和。谢谢！

我发现有一部分程序没有贴上去，而且是重要的部分。我把Data Two再贴一次。

1. data Two;
   
2.   set One;
   
3.   by id;
   
4.   array o{0:150}  _temporary_;
   
5.   array dc{0:150} _temporary_;
   
6.   if first.id then do;
   
7.     do i = 0 to 10;
   
8.       dc[i] =0;
   
9.       o[i] = 0;
   
10.     end;
    
11.   end;
    
12.        
    
13.         /* 主要问题在这里 */
    
14.   dc[t+dm] = dc[t+dm] + s;
    
15.   if last.id then do;
    
16.      x    = 0;
    
17.      do i = 0 to 3;
    
18.          x = x + dc[t+i];
    
19.          d3 = dc[t+i];
    
20.         put '------------------';
    
21.         put 'i t d3 x';
    
22.         put  i t d3 x;
    
23.      end;
    
24.   end;
    
25.   else do;
    
26.       x = dc[t];
    
27.       put '********************';
    
28.       put x;
    
29.    end;
    
30. run;

复制代码

wx2123 发表于 2022-2-26 13:33
你说的对，我也在探索x移动的规则。而且我可能在误导你，因为我自己说的也不一定对。
这个问题可能没有那 ...

你的代码没办法输出你给的图，你的没有定义x.

owenqi 发表于 2022-2-27 05:20
你的代码没办法输出你给的图，你的没有定义x.

x是在data two里面算出来的