我正在尝试使用效果包来提取PLM模型的预测值。不幸的是,只要我包含一个因子变量,效果包就会出现错误消息
Error in Summary.factor(c(2L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 1L,
‘range’ not meaningful for factors
有没有人找到解决此问题的方法?
以下是复制问题的最低工作示例
data("Produc", package = "plm")
Produc$dummy1<-(rbinom(nrow(Produc), 1 , 0.4)) #Generate Random Dummy
Produc$dummy2<-as.factor(Produc$dummy1)
mod1 <- plm(gsp ~ pcap * dummy1,data = Produc, index = c("state","year"), model="fd")
mod2 <- plm(gsp ~ pcap * dummy2,data = Produc, index = c("state","year"), model="fd")
Effect(c("pcap","dummy1"),mod1) #This works
Effect(c("pcap","dummy2"),mod2) #This does not
这可能与此question有关,但错误不同。