我有一个面板数据（美国所有上市公司都有季度数据），如果在投资推荐网站上推荐该公司，则因变量为1，否则为0。自变量是几个股票特征（例如股本回报率，市盈率，销售增长，市值等）。
我使用pglm软件包执行了logit池回归。之后，我尝试使用blorr和pscl软件包来计算伪R平方，但是它不起作用。
所以，我的代码如下：
Panel.set.buy <- pdata.frame(DB_buy, index = c("gvkey", "datacqtr"))
pglm_buy <- pglm(VIC ~ ROA + ROE + LEV + CAPEX + SG + TURN + RET + MKCAP +
PriceEarnings + M_B + NEWS_rs + NANALYSTS, Panel.set.buy, family =
binomial(link = "logit"), model = "pooling")
summary (pglm_buy) <br> </br>
pR2(pglm_buy)
blr_rsq_mcfadden(pglm_buy)
运行最后两行（pR2和blr_rsq_mcfadden）时，出现以下错误：
UseMethod（“ pR2”）中的错误：没有将适用于'pR2'的适用方法应用于类“ c（'maxLik'，'maxim'，'list'）
我在做什么错？我该如何解决问题？
谢谢