R 线性回归中的聚类标准误差

新手问题：我目前正在写一篇关于分析工具对公司粉丝数量影响的论文。我创建了以下回归模型：

Fans~eventWeek*analysistool

---

这个函数的问题在于它假设数据中的每一行都是不相关的，而事实并非如此。因此，我的论文指导建议我在R中使用聚集标准错误，并将其聚集在公司编号上。有人知道怎么做吗

我相信可以用GLMs回答您的问题。

链接已断开。