我想使用调查包在数据框中将另外两个变量分组后找到变量的中值。
如果不使用调查包,这是我要使用的代码:
library(dplyr)
df %>%
group_by(a,b) %>%
summarise(MedianC=(median(c)))
我只想使用调查方法得出相同的结果
这是我到目前为止尝试过的:
design_df <- svydesign(ids = ~id, strata = ~strata, weights = ~weight, data = df)
medians_c <- svyby(~c, ~a+b, design_df, svyquantile, quantiles = 0.5, ci = TRUE, vartype = "ci") %>% as_tibble()
如果我在svyby语句中仅包含a而不是〜a + b,则获得a的每个值的c的中值。但是我需要a和b的每个组合的c的中值。