我正在使用可以找到here的Mastering Metrics数据的第5.2和5.3章中的deaths文件。我正在尝试按状态和年份对pop进行分组，并返回平均值pop。我尝试了几种方法但是当我合并组时，结果看起来不合适。
mlda <- read.dta("C:/Users/perdue/Desktop/Adv.MicroEconometrics/HA 8/deaths.dta") ## reading in my data
mlda$avg.pop <- aggregate(mlda[, 7], list(mlda$state), mean)
Error in `$<-.data.frame`(`*tmp*`, avg.pop, value = list(Group.1 = c(1L, :
replacement has 51 rows, data has 24786