我在设置R中的面板数据模型（固定效果）时遇到问题。
目前我正在运行以下代码：
fe1 <- summary(lm(qnorm(y) ~ factor(Bank) -1 + factor(Country)*x ,data=PDwideHPI))
每个国家/地区每个银行有6次观察。我希望摘要输出能为我提供特定于银行的拦截，其次是x的国家特定系数。但是，R-console会返回特定于银行的AND特定国家/地区的拦截。
因此我收到的估计数：因子（银行）1因子（银行）2因子（银行）3因子（银行）4因子（国家）1因子（国家）2 x因子（国家）1：x因子（国家） ）2：X
而我的预期如下：因子（银行）1因子（银行）2因子（银行）3因子（银行）4 x因子（国家）1：x因子（国家）2：x
有没有简单的方法来解决这个问题？我对R很新，因此我希望这对于专业人士来说是一个简单的解决方案。提前谢谢
答案 0 :(得分：0)
使用默认的R公式语法，*不仅包括交互术语，还包括各个术语。如果您只想要交互术语，那么您就是:运算符。所以，在你的情况下，你想要
fe1 <- summary(lm(qnorm(y) ~ factor(Bank) -1 + factor(Country):x ,data=PDwideHPI))