R语言带Newey-West标准误的线性回归模型缺失洪水交互项系数求助

我正在做自然灾害对收入影响的回归分析，完成自相关（AC）和异方差检验后，为模型应用了Newey-West（NW）标准误。但用stargazer输出模型结果时，洪水虚拟变量Lagged1y_Disaster_Dummy_Flood对应的Interaction_Flood系列交互项完全没有显示系数。奇怪的是，单独拟合仅包含洪水交互项的模型时，这些项均显著；不构建交互项时，相关变量也显著。以下是我的代码、回归输出关键部分及相关系数矩阵，恳请帮忙解决：

模型代码 #

model_pooled_ols_5e <- lm(log_Average_Wage_and_Salary_Income ~ Lagged1y_Disaster_Dummy_Storm + Lagged1y_Disaster_Dummy_Tropical_Cyclone_Storm + Lagged1y_Disaster_Dummy_Wildfire + Lagged1y_Disaster_Dummy_Flood + Intensity_Score + Duration_Under_1_Week_Dummy + Duration_1_to_3_Weeks_Dummy + Duration_1_Month_Dummy + Duration_Over_1_Month_Dummy + Inner_Regional_Dummy + Major_Cities_Dummy + Outer_Regional_Dummy + Remote_Very_Remote_Dummy + Interaction_Storm_Inner_Regional + Interaction_Storm_Major_Cities + Interaction_Storm_Outer_Regional + Interaction_Storm_Remote_Very_Remote + Interaction_Tropical_Cyclone_Inner_Regional + Interaction_Tropical_Cyclone_Major_Cities + Interaction_Tropical_Cyclone_Outer_Regional + Interaction_Tropical_Cyclone_Remote_Very_Remote + Interaction_Wildfire_Inner_Regional + Interaction_Wildfire_Major_Cities + Interaction_Wildfire_Outer_Regional + Interaction_Wildfire_Remote_Very_Remote + Interaction_Flood_Inner_Regional + Interaction_Flood_Major_Cities + Interaction_Flood_Outer_Regional + Interaction_Flood_Remote_Very_Remote + Year, data=Merged_ABS_EMDAT_simple_duplicates_cleaned_v5)

回归输出（关键部分） #

因变量：log_Average_Wage_and_Salary_Income
...
Interaction_Flood_Inner_Regional
Interaction_Flood_Major_Cities
Interaction_Flood_Outer_Regional
Interaction_Flood_Remote_Very_Remote
...
注：已应用Newey-West标准误

相关系数矩阵（洪水交互项部分） #

log_Average_Wage_and_Salary_Income Lagged1y_Disaster_Dummy_Storm
Interaction_Flood_Inner_Regional                            -0.18158278                   -0.03057133
Interaction_Flood_Major_Cities                               0.09542833                   -0.05145906
Interaction_Flood_Outer_Regional                            -0.18727074                   -0.02185281
Interaction_Flood_Remote_Very_Remote                        -0.08479981                   -0.01270038
...

• 多重共线性排查：你同时放入了洪水主效应、四个区域虚拟变量以及对应的四个洪水-区域交互项，这会导致完全共线性（主效应+区域虚拟变量+交互项存在线性相关关系）。建议去掉一个区域基准项（比如Remote_Very_Remote_Dummy），同时移除对应基准组的洪水交互项（Interaction_Flood_Remote_Very_Remote），重新拟合模型。
• stargazer参数检查：确认是否设置了keep/omit参数误过滤了这些交互项；或者因系数过小被自动省略，可尝试增加digits参数的小数位数，或调整keep.stat参数确保输出完整。
• Newey-West标准误传递验证：确保你是通过sandwich包计算NW标准误后正确传递给stargazer，比如用coeftest(model_pooled_ols_5e, vcov = NeweyWest(model_pooled_ols_5e))生成带NW标准误的结果，再传入stargazer，避免计算过程中丢失变量信息。
• 变量变异检查：用table(Interaction_Flood_Inner_Regional, data=Merged_ABS_EMDAT_simple_duplicates_cleaned_v5)这类命令，查看洪水交互项是否存在全0/全1的情况，若变量无变异，模型无法估计系数。

内容的提问来源于stack exchange，提问作者Lucy