R 线性概率模型

我有一个代表一些国际奥林匹克运动会参与者的数据集，它有列性别（虚拟变量）、奥林匹克年、一个人的名字和这个人获得的奖品。我需要估计以下线性概率模型：1{gold_ind_I=1}=α+β×1{d_female_I=1}+∑_（j来自1968年至2018年）μ_j×1{Year j=Year i}+εi
，也就是说，我的因变量是虚拟变量，用于描述某个人是否获得了金牌，预测因子是虚拟性别（如果一个人是女性，则等于1）和每年的虚拟变量。例如，在1999年，该年的dummy等于1，所有其他年份的Dummie为0。

---

我可以对变量性别进行回归。问题在于将假人包含在与年份相关的模型中，因为我有这些年份的大量样本（从1960年代到2000年代）

找到了一个解决方案：只需使用“fastDummies”软件包。

如果您包含一个简单的样本输入和所需的输出，可以用来测试和验证可能的解决方案，那么就更容易帮助您。