r - split-apply-combine 论文中的代码错误 - 如何解决？

Question

为了尝试掌握 R 中的数据操作，我开始阅读 Hadley 关于 split-apply-combine 的论文.

我在第 3 页并尝试通过代码来理解它。不幸的是，代码有误，我的复制是忠实的(我已经完成了 c&p 和手写)。当我试图学习这些东西时，我一开始是对的，我实际上无法说出它有什么问题。我在 R2.5 和 R3.0 上都试过了

library("MASS")
library("plyr")
data(ozone)
one<-ozone[1,1,]
month<-ordered(rep(1:12,length=72))
model<-rlm(one ~ month - 1)
deseas<-resid(model)
deseasf<-function(value) {rlm(value ~ month - 1)}
models<-aaply(ozone,1:2,deseasf)
deseas<-aaply(models,1:2,resid)

模型行错误为 Error: Results must have one or more dimensions.

有人可以告诉我它是否适合他们，或者如果不适合需要修复/修改什么，为什么？





PS - 无法查看 http://plyr.had.co.nz/勘误表，因为我的工作代理目前阻止了该站点!

Answer 1

应该是

models <- alply(ozone, 1:2, deseasf)
deseas <- ldply(models, resid)