Betareg导致内存分配问题。

Question

我试图使用betareg包在R中运行一个beta回归模型。然而，我意外地遇到了内存大小问题。让我解释一下为什么这让我吃惊。

我使用的是Windows 7，64位，R-64，内存中有32 GB .

我正在运行的betareg命令是：

br1 <- betareg(dfp ~ ago + evl + spe + low + poly(fpt, 4, raw = T), data = tt[tt$zero_one_ind == 1, ], model = T, y = F, x = F)

betareg模型的对象大小是：

打印(object.size(Br1)，单位= " Mb ") 46 Mb

这是我正在收到的错误消息：

Error: cannot allocate vector of size 344.1 Gb
In addition: Warning messages:
1: In diag(x %*% xwx1 %*% t(x)) :
Reached total allocation of 32693Mb: see help(memory.size)
2: In diag(x %*% xwx1 %*% t(x)) :
 Reached total allocation of 32693Mb: see help(memory.size)
3: In diag(x %*% xwx1 %*% t(x)) :
 Reached total allocation of 32693Mb: see help(memory.size)
4: In diag(x %*% xwx1 %*% t(x)) :
Reached total allocation of 32693Mb: see help(memory.size)

betareg模型成功地在R中运行并估计了系数，据我所知，所有的时隙都被填充了，但看起来R无法构造方差协方差矩阵。这里有什么不对劲的地方吗？

Answer 1

我也有同样的问题，解决办法很简单。

从手册中：

请注意，在大样本中，默认的残差"sweighted2“可能需要计算，因此可能需要在此类应用程序中进行修改。

例如，您可以在摘要中使用其他选项之一：

类型=c(“皮尔森”、“越轨”、“响应”、“加权”、“抽奖”、“sweighted2”)

Answer 2

我在使用betareg时也遇到了同样的问题。我只对系数意义的p值感兴趣，所以这对我来说是一种有用的方法：

fit_frst_spnd_model <- betareg(formula = frst_spnd_util_pc2 ~ .
                        ,data = train_data_frst_txn2_2
                        )
library(lmtest)
coeftest(fit_frst_spnd_model)