Speedlm更新“需要呼叫组件的对象”

Question

我试图使用speedlm(speedglm package)的更新选项，因为我没有足够的RAM一次计算整个模型，而biglm只使用一个CPU。
下面的代码是出错的可重现示例。

library(speedglm) 
formula <- Sepal.Length ~ Sepal.Width 
chunk1 <- iris[1:10,] 
chunk2 <- iris[11:20,] 
chunk3 <- iris[21:30,] 
lmfit <- speedlm(formula, chunk1) 
lmfit <- update(lmfit, chunk2) 
lmfit <- update(lmfit, chunk3) 

，我发现了以下错误：

> lmfit <- speedlm(formula, chunk1) 
> lmfit <- update(lmfit, chunk2) 
> lmfit <- update(lmfit, chunk3) 
Error in update.default(lmfit, chunk3) : 
    need an object with call component 
> 

如果是由于有update代替updateWithMoreData我已经期待着与chunk2更新后的错误。

想知道解决这个问题的方法，或者如果我必须使用替代方法。
在此先感谢！

---

使用updateWithMoreData得到以下错误：

> lmfit <- speedlm(formula, chunk1) 
> lmfit <- updateWithMoreData(lmfit, chunk2) 
Error: object of type 'symbol' is not subsettable 
> lmfit <- updateWithMoreData(lmfit, chunk3) 
Error: object of type 'symbol' is not subsettable 
> 

---

下面的代码工作，道具@LyzandeR

> library(speedglm) 
> chunk1 <- iris[1:10,] 
> chunk2 <- iris[11:20,] 
> chunk3 <- iris[21:30,] 
> lmfit <- speedlm(Sepal.Length ~ Sepal.Width, chunk1) 
> 
> for (i in list(11,20, 21:30)){ 
+ lmfit2 <- updateWithMoreData(lmfit, iris[i,]) 
+ } 
> lmfit2 
Linear Regression Model of class 'speedlm': 

Call: speedlm(formula = Sepal.Length ~ Sepal.Width, data = chunk1) 

Coefficients: 
(Intercept) Sepal.Width 
    2.9876  0.5813 

> 

Answer 1

为了使用update来更新你的你的模型根据@ Roland的评论，需要使用updateWithMoreData。虽然有一个问题。

如果使用这样的代码，你会得到：

library(speedglm) 
formula <- Sepal.Length ~ Sepal.Width 
chunk1 <- iris[1:10,] 
chunk2 <- iris[11:20,] 
chunk3 <- iris[21:30,] 
lmfit <- speedlm(formula, chunk1) 
#runs fine up to here 

#but this gives an error 
lmfit2 <- updateWithMoreData(lmfit, chunk2) 
Error: object of type 'symbol' is not subsettable 

显然，这是造成的，因为它试图把这个（从回溯）：

as.formula(object$call[[2]])

其失败的原因lmfit$call[[2]]返回formula。但是，如果你把它改成这样它的工作原理：

library(speedglm) 
chunk1 <- iris[1:10,] 
chunk2 <- iris[11:20,] 
chunk3 <- iris[21:30,] 
#use the actual formula below so lmfit$call[[2]] will return it 
#and now all the three lines below work fine 
lmfit <- speedlm(Sepal.Length ~ Sepal.Width, chunk1) 
lmfit2 <- updateWithMoreData(lmfit, chunk2) 
lmfit3 <- updateWithMoreData(lmfit, chunk3) 

注意，当你打印出来，既会说他们是因为呼叫保持不变，但结果是正确的块1：

输出：

> lmfit2 
Linear Regression Model of class 'speedlm': 

Call: speedlm(formula = Sepal.Length ~ Sepal.Width, data = chunk1) 

Coefficients: 
(Intercept) Sepal.Width 
    1.8398  0.9181 

> lmfit 
Linear Regression Model of class 'speedlm': 

Call: speedlm(formula = Sepal.Length ~ Sepal.Width, data = chunk1) 

Coefficients: 
(Intercept) Sepal.Width 
    2.3882  0.7468