带加权因变量的双向方差分析

Question

试图按个人规模比较多年来的 3 个独立人群，我有这样的数据集：

year <- c(rep(2000,5),rep(2001,3),rep(2002,7))
region <- c(1,1,2,3,3,1,2,3,rep(1,3),rep(2,3),3)
size <- c(28,24,26,56,47,85,12,24,68,71,42,59,12,25,33)
count <- c(3,8,9,1,2,4,7,12,4,8,3,2,7,15,4)
df <- data.frame(year, region, size, count)

这给出了：

year region size count
2000      1   28     3
2000      1   24     8
2000      2   26     9
2000      3   56     1
2000      3   47     2
2001      1   85     4
2001      2   12     7
2001      3   24    12
2002      1   68     4
2002      1   71     8
2002      1   42     3
2002      2   59     2
2002      2   12     7
2002      2   25    15
2002      3   33     4

我想做一个 2-Way ANOVA：

model.2way <- lm(size ~ year * region, df)  # example of code
anova(model.2way)

我的问题是变量size 由count 加权：对于每个size，我有count 个人数。我获得了数百万个数据，并且无法轻松地将我的数据转换为具有数百万个 size 值。

您知道用这种加权数据制作 2-Way ANOVA 的方法吗？

提前致谢！

Answer 1

model.2way <- lm(size ~ year * region, df, weights = count)

来自?lm：

...当“权重”的元素是正整数 w_i 时，每个响应 y_i 是 w_i 个单位权重观测值的平均值 ...

换句话说，权重为 2 意味着该案例出现了两次。