R: Using aggregate within a function, not working

你。 提交于 2021-02-11 12:28:17

问题


I am fairly new to R so apologies if an answer to this already exists that I am unable to find.

I cannot replicate the exact error I have with my own dataset, but since an error is produced nonetheless here we go. What I am trying to do is to create a function to calculate the average of several columns conditionnal on the values of other ones. Let's say that

d1 <- c(1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4)
d2 <- c(1:12)
d3 <- c(1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2)
df <- cbind(d1, d2, d3)

indicators_plotting1 <- function(a, b = d2, c = 1) {
  temp_data <- df %>%
    filter(d3 = c)
    for_plot <- aggregate(b ~ a, df, mean)
    plot(for_plot$b, for_plot$a)
}
indicators_plotting1(a = d1)

It works well if I do not include the %>% filter(d3 = c) bit. When I include it I get Error in UseMethod("filter_") : no applicable method for 'filter_' applied to an object of class "c('matrix', 'double', 'numeric')".

That when I do it with the dataset I propose here. When I do it with my own larger one and with four filters, what I get is Error in eval(predvars, data, env) : object 'd1' not found.

Any thoughts? Thanks in advance


回答1:


Consider following base R adjustment using dynamic formula with reformulate in aggregate() and plot(). For below solution, you need to pass column reference as strings. And below works on named matrices or data frames:

indicators_plotting1 <- function(a, b = "d2", c = 1) {
  temp_data <- subset(df, d3 == c)
  for_plot <- aggregate(reformulate(a, b), temp_data, mean)
  plot(reformulate(b, a), for_plot)
}

indicators_plotting1(a = "d1")

Alternatively a one-line function:

indicators_plotting1 <- function(a, b = "d2", c = 1) {
  plot(reformulate(b, a), 
       aggregate(reformulate(a, b), subset(df, d3 == c), mean))
}

Online Demo (click Run it (F8) at bottom for plot)



来源:https://stackoverflow.com/questions/61213456/r-using-aggregate-within-a-function-not-working

标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!