大型数据集的面板数据回归（ffdf）

Question

我使用ffdf类来存储相当大的面板数据集。我想知道是否有任何方法可以为此类中的数据运行固定效果回归。

library(ffbase)
size <- 100000
index <- 100
a <- ffrandom(n=index, rfun = runif)
panel_index <- ffseq_len(index)
times_vec <- ffrep.int(size,times=index)
b <- ffrep.int(a, times=times_vec)
panel_index <- ffrep.int(panel_index, times=times_vec)
error <- ffrandom(n=size*index, rfun=rnorm)
intercept <- ffrep.int(10, times=size*index)
x <- ffrandom(n=size*index, rfun=rnorm)
y <- x+intercept+b+error

dataf <- ffdf(y = y, x = x, panel_index = panel_index) ## on disk
dim(dataf)

我想运行像

这样的东西

library(plm)
fixed <- plm(y ~ x, data=dataf, 
    index="panel_index", model="within")

但是我收到以下错误

Error in is.finite(x[[i]]) : 
  default method not implemented for type 'list'

哪些替代软件包可以处理ffdf数据？