我正在使用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")
但我收到以下错误消息：
is.finite(x[[i]]) 中的错误：未为类型“列表”实现默认方法
可以处理数据的替代包有哪些ffdf？