计算每个先前组号的平均值

如何解决计算每个先前组号的平均值

我想计算之前每组的平均含糖量。

library(dplyr)

df <- tibble(grp = c(1,1,2,3,4,5,6,7,8),sog = c(12.00,12.00,11.00,10.00,9.00,2.00,8.80,8.70,7.80,9.54,9.80,10.40,4.00,4.50,3.60,7.00,8.00,0.00),res = c(NA,NA,66.66667,275.0,0)
)

我想要类似的东西

df %>%
    group_by(grp) %>%
    mutate(sog_avg = dplyr::lag(mean(sog),n = 'by every previous group',default = NA))

现在的数据是:

 grp   sog   res

 1     12     NA  
 1     12     NA  
 1     12     NA  
 1     11     NA  
 1     10     NA  
 1      9     NA  
 2      2     NA  
 3     8.8    0  
 3     8.7    0  
 3     7.8    0  
 4     11      0  
 4     11      0  
 4     12      0  
 4     11      0  
 5     9.54   0  
 5     9.8    0  
 5     10.4    0  
 6     4     66.7
 6     4     66.7
 6     4.5   66.7
 6     3.6   66.7
 7     7    275  
 7     8    275  
 7     9    275  
 8     0      0  

我想要的结果是:

 grp   sog   res  sog_avg

 1     12     NA  NA
 1     12     NA  NA
 1     12     NA  NA
 1     11     NA  NA
 1     10     NA  NA
 1      9     NA  NA
 2      2     NA  11
 3     8.8    0    2
 3     8.7    0    2
 3     7.8    0    2
 4     11     0   8.43
 4     11     0   8.43 
 4     12     0   8.43
 4     11     0   8.43
 5     9.54   0   11.25
 5     9.8    0   11.25
 5     10.4   0   11.25
 6     4     66.7 9.91
 6     4     66.7 9.91
 6     4.5   66.7 9.91 
 6     3.6   66.7 9.91
 7     7    275   4.025 
 7     8    275   4.025
 7     9    275   4.025
 8     0      0   8

解决方法

我们可以计算每个组的均值,使用lag将其后退一步,然后加入数据:

library(dplyr)

df %>%
  group_by(grp) %>%
  summarise(sog_avg = mean(sog,na.rm = TRUE)) %>%
  mutate(sog_avg = lag(sog_avg)) %>%
  left_join(df,by = 'grp')

#   grp sog_avg   sog   res
#1    1      NA 12.00    NA
#2    1      NA 12.00    NA
#3    1      NA 12.00    NA
#4    1      NA 11.00    NA
#5    1      NA 10.00    NA
#6    1      NA  9.00    NA
#7    2   11.00  2.00    NA
#8    3    2.00  8.80   0.0
#9    3    2.00  8.70   0.0
#10   3    2.00  7.80   0.0
#11   4    8.43 11.00   0.0
#12   4    8.43 11.00   0.0
#13   4    8.43 12.00   0.0
#14   4    8.43 11.00   0.0
#15   5   11.25  9.54   0.0
#16   5   11.25  9.80   0.0
#17   5   11.25 10.40   0.0
#18   6    9.91  4.00  66.7
#19   6    9.91  4.00  66.7
#20   6    9.91  4.50  66.7
#21   6    9.91  3.60  66.7
#22   7    4.03  7.00 275.0
#23   7    4.03  8.00 275.0
#24   7    4.03  9.00 275.0
#25   8    8.00  0.00   0.0
,

添加了purrr的一个选项可能是:

df %>%
 mutate(sog_avg = map_dbl(.x = grp - 1,~ mean(sog[grp == .x]))) 

     grp   sog   res sog_avg
   <dbl> <dbl> <dbl>   <dbl>
 1     1  12      NA     NaN
 2     1  12      NA     NaN
 3     1  12      NA     NaN
 4     1  11      NA     NaN
 5     1  10      NA     NaN
 6     1   9      NA     NaN
 7     2   2      NA      11
 8     3   8.8     0       2
 9     3   8.7     0       2
10     3   7.8     0       2

或仅使用dplyr

df %>%
 group_by(grp) %>%
 mutate(sog_avg = mean(sog)) %>%
 ungroup() %>%
 mutate(sog_avg = lag(sog_avg)) %>%
 group_by(grp) %>%
 mutate(sog_avg = first(sog_avg)) 

     grp   sog   res sog_avg
   <dbl> <dbl> <dbl>   <dbl>
 1     1  12      NA      NA
 2     1  12      NA      NA
 3     1  12      NA      NA
 4     1  11      NA      NA
 5     1  10      NA      NA
 6     1   9      NA      NA
 7     2   2      NA      11
 8     3   8.8     0       2
 9     3   8.7     0       2
10     3   7.8     0       2

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。

相关推荐


依赖报错 idea导入项目后依赖报错,解决方案:https://blog.csdn.net/weixin_42420249/article/details/81191861 依赖版本报错:更换其他版本 无法下载依赖可参考:https://blog.csdn.net/weixin_42628809/a
错误1:代码生成器依赖和mybatis依赖冲突 启动项目时报错如下 2021-12-03 13:33:33.927 ERROR 7228 [ main] o.s.b.d.LoggingFailureAnalysisReporter : *************************** APPL
错误1:gradle项目控制台输出为乱码 # 解决方案:https://blog.csdn.net/weixin_43501566/article/details/112482302 # 在gradle-wrapper.properties 添加以下内容 org.gradle.jvmargs=-Df
错误还原:在查询的过程中,传入的workType为0时,该条件不起作用 &lt;select id=&quot;xxx&quot;&gt; SELECT di.id, di.name, di.work_type, di.updated... &lt;where&gt; &lt;if test=&qu
报错如下,gcc版本太低 ^ server.c:5346:31: 错误:‘struct redisServer’没有名为‘server_cpulist’的成员 redisSetCpuAffinity(server.server_cpulist); ^ server.c: 在函数‘hasActiveC
解决方案1 1、改项目中.idea/workspace.xml配置文件,增加dynamic.classpath参数 2、搜索PropertiesComponent,添加如下 &lt;property name=&quot;dynamic.classpath&quot; value=&quot;tru
删除根组件app.vue中的默认代码后报错:Module Error (from ./node_modules/eslint-loader/index.js): 解决方案:关闭ESlint代码检测,在项目根目录创建vue.config.js,在文件中添加 module.exports = { lin
查看spark默认的python版本 [root@master day27]# pyspark /home/software/spark-2.3.4-bin-hadoop2.7/conf/spark-env.sh: line 2: /usr/local/hadoop/bin/hadoop: No s
使用本地python环境可以成功执行 import pandas as pd import matplotlib.pyplot as plt # 设置字体 plt.rcParams[&#39;font.sans-serif&#39;] = [&#39;SimHei&#39;] # 能正确显示负号 p
错误1:Request method ‘DELETE‘ not supported 错误还原:controller层有一个接口,访问该接口时报错:Request method ‘DELETE‘ not supported 错误原因:没有接收到前端传入的参数,修改为如下 参考 错误2:cannot r
错误1:启动docker镜像时报错:Error response from daemon: driver failed programming external connectivity on endpoint quirky_allen 解决方法:重启docker -&gt; systemctl r
错误1:private field ‘xxx‘ is never assigned 按Altʾnter快捷键,选择第2项 参考:https://blog.csdn.net/shi_hong_fei_hei/article/details/88814070 错误2:启动时报错,不能找到主启动类 #
报错如下,通过源不能下载,最后警告pip需升级版本 Requirement already satisfied: pip in c:\users\ychen\appdata\local\programs\python\python310\lib\site-packages (22.0.4) Coll
错误1:maven打包报错 错误还原:使用maven打包项目时报错如下 [ERROR] Failed to execute goal org.apache.maven.plugins:maven-resources-plugin:3.2.0:resources (default-resources)
错误1:服务调用时报错 服务消费者模块assess通过openFeign调用服务提供者模块hires 如下为服务提供者模块hires的控制层接口 @RestController @RequestMapping(&quot;/hires&quot;) public class FeignControl
错误1:运行项目后报如下错误 解决方案 报错2:Failed to execute goal org.apache.maven.plugins:maven-compiler-plugin:3.8.1:compile (default-compile) on project sb 解决方案:在pom.
参考 错误原因 过滤器或拦截器在生效时,redisTemplate还没有注入 解决方案:在注入容器时就生效 @Component //项目运行时就注入Spring容器 public class RedisBean { @Resource private RedisTemplate&lt;String
使用vite构建项目报错 C:\Users\ychen\work&gt;npm init @vitejs/app @vitejs/create-app is deprecated, use npm init vite instead C:\Users\ychen\AppData\Local\npm-