我正在尝试借助 fixst 包对具有固定效应的泊松最大似然模型进行多重估计。
我已经分别运行了所有不同的估计，它们工作得非常好！还包括对因变量 hot、temp 和 heat_index 的逐步函数 sw() 工作得非常好。就在我尝试在方程的固定效应部分使用 sw() 时，会发生以下错误。我不明白发生了什么，尤其是单独运行所有这些固定效果（没有 sw() 函数）也可以很好地工作。根据 fixst 文档，应该可以将 sw() 包含到代码的固定效果部分。
产生错误的代码行： m <- femlm(food_drinks ~ sw(hot,temp,heat_index) |county_fips + sw0(month,week,statemonth,stateweek), data = data,family="poisson")
错误：femlm 中的错误（food_drinks ~ sw(hot, temp, heat_index) | sw(countymonth, : unlist(fixef_id) 中的错误 - 1 : 二元运算符的非数字参数此错误是函数 femlm 的作者未预见到的。如果你认为你对该函数的调用是合法的，你能报告吗？
该错误使我认为我必须将用作固定效果的数据框列转换为分类/因子变量，但它们已经是类型字符，正如我所说，估计在单独运行时有效。
我会很感激任何提示！谢谢！