我正在尝试使用效果包来提取 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
这可能与这个问题有关,但错误是不同的。