split-apply-combine纸张代码出错 - 如何解决？

Question

为了尝试掌握R中的数据操作,我开始阅读Hadley关于分割 - 应用 - 组合的论文.

我在第3页,试图通过代码来理解它.不幸的是,代码是错误的,我的复制是忠实的(我已经完成了c&p和handtyped).因为我正在努力学习这些东西而且我刚开始时我无法确定它有什么问题.我在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)