这是我想估计加权平均数的列列表:

vlist = c("Public_Welf_Total_Exp", "Welf_Cash_Total_Exp", "Welf_Cash_Cash_Assist", "Welf_Ins_Total_Exp","Total_Educ_Direct_Exp", "Higher_Ed_Total_Exp", "Welf_NEC_Cap_Outlay","Welf_NEC_Direct_Expend", "Welf_NEC_Total_Expend", "Total_Educ_Assist___Sub", "Health_Total_Expend", "Total_Hospital_Total_Exp", "Welf_Vend_Pmts_Medical","Hosp_Other_Total_Exp","Unemp_Comp_Total_Exp", "Unemp_Comp_Cash___Sec", "Total_Unemp_Rev", "Hous___Com_Total_Exp", "Hous___Com_Construct")

这是我一直使用的代码:

df = df %>% group_by(cz) %>% mutate_at(vlist, weighted.mean(., df$Population))

我也尝试过:

df = df %>% group_by(cz) %>% mutate_at(vlist, function(x) weighted.mean(x, df$Population))

df = df %>% group_by(cz) %>% mutate_at(vars(Public_Welf_Total_Exp, Welf_Cash_Total_Exp), weighted.mean(., df$Population))

但是,我尝试的所有操作都会导致以下错误,即使我的任何变量中都没有NAs:

Error in weighted.mean.default(., df$Population) :

'x' and 'w' must have the same length

我知道我可以用lapply来做下面的估计,但是我不知道如何用lapply来对另一个变量进行分组。如果有任何建议,我将不胜感激!