R语言 summarise函数分组计算失败

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

1. library(openxlsx)

复制代码

程序设计希望能将数据按照sp分组，通过自己定义的函数计算不同sp的relativeBA和relativeN，但是并没有成功达到预期，报错如下：

1. > groups.importance_by_sp=summarise(sp_groups,BA_relative=relativeBA(),N_relative=relativeN())
   
2. `summarise()` has grouped output by 'sp'. You can override using the `.groups` argument.
   
3. Warning message:
   
4. Returning more (or less) than 1 row per `summarise()` group was deprecated in dplyr 1.1.0.
   
5. ℹ Please use `reframe()` instead.
   
6. ℹ When switching from `summarise()` to `reframe()`, remember that `reframe()` always returns an ungrouped data frame and adjust
   
7.   accordingly.
   
8. Call `lifecycle::last_lifecycle_warnings()` to see where this warning was generated.
   
9. > groups.importance_by_sp=reframe(sp_groups,BA_relative=relativeBA(),N_relative=relativeN())
   
10. > groups.importance_by_sp
    
11. # A tibble: 21,112 × 3
    
12.    sp    BA_relative N_relative
    
13.    <chr>       <dbl>      <dbl>
    
14. 1 乌桕   0.00000167          1
    
15. 2 乌桕   0.00000211          1
    
16. 3 乌桕   0.00000329          1
    
17. 4 乌桕   0.00000437          1
    
18. 5 乌桕   0.00000147          1

复制代码

提示我改用reframe函数，但修改后，根本没有分组计算，而是对原数据每一行分开进行了处理。请问是不是我的定义函数哪里需要修改呢？求大神帮助~谢谢！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：arise Summa aris rise R语言

1. library(openxlsx)
   
2. comm=read.xlsx("comm_island.xlsx")
   
3. 
   
4. attach(comm)
   
5. comm$branch=as.numeric(branch)
   
6. comm$dbh=as.numeric(dbh)
   
7. comm$h=as.numeric(h)
   
8. 
   
9. unique(comm$sp)
   
10. table(comm$sp)
    
11. 
    
12. comm=comm[comm$branch==0,] #选择主枝
    
13. #前置操作
    
14. 
    
15. BA=function(a){
    
16. a=comm$dbh
    
17. r=a/2
    
18. S=pi*r^2
    
19. return(S)
    
20. }
    
21. #计算胸高断面积的函数
    
22. 
    
23. total_BA=sum(apply(comm,1,BA))#计算总BA
    
24. total_BA=61506864
    
25. total_N=nrow(comm)#计算总个体数
    
26. total_N=728
    
27. #以上正确
    
28. 
    
29. relativeBA=function(a){
    
30. sum_S=0
    
31. a=comm$dbh
    
32. r=a/2
    
33. S=pi*r^2
    
34. sum_S=sum_S+S #计算该物种总胸高断面积（在groupby情况下）
    
35. relativeBA=sum_S/total_BA
    
36. return(relativeBA)
    
37. }
    
38. #计算某物种相对胸高断面积
    
39. 
    
40. relativeN=function(b){
    
41. b=comm$sp
    
42. is.character(b)
    
43. comm_selected=comm[comm$sp=="b"]
    
44. N=nrow(comm_selected)
    
45. relativeN=N/total_N
    
46. return(relativeN)
    
47. }
    
48. #计算某物种相对个体数
    
49. 
    
50. library(dplyr)
    
51. 
    
52. sp_groups=group_by(comm,sp)
    
53. #按照sp分组
    
54. 
    
55. groups.importance_by_sp=summarise(sp_groups,BA_relative=relativeBA(),N_relative=relativeN())

复制代码

海阔天空7890 发表于 2023-10-13 15:30

原代码如上