0

我正在尝试使用效果包来提取 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

这可能与这个问题有关,但错误是不同的。

4

0 回答 0