1

函数 my.ccf 会生成一个绘图等。

my.ccf <- function(dat)
{
require(forecast)
modx <- auto.arima(dat[,1])
modx$series <- colnames(dat)[1]
mody <- Arima(dat[,2], model=modx)
mody$series <- colnames(dat)[2]
ccf1 <- ccf(residuals(modx), residuals(mody), type="correlation", ylab="CCF", 
            main=paste(colnames(dat)[1], "&", colnames(dat)[2], sep=" "))
return(list(modx=modx, mody=mody, ccf=ccf1))    
}
dat1 <- data.frame(p=rnorm(50), q=rnorm(50)); my.ccf(dat1)

生成的图形以“p & x”等为标题

我正在使用函数year.ccf 中的do.call 将这个函数应用于另一个data.frame (dat2)。首先，它将按年对传递的 data.frame 进行子集化（我有单独的函数来执行此操作。出于说明目的，我在此处使用 for 循环），然后将 my.ccf 应用于名为“x”的列与给定列的配对组合在 var.do. var.do 的长度从 3 到 9 不等。我必须在不同的数据帧上重复多次（超过 500 次）并尝试尽可能地自动执行。

year.ccf1 <- function(dat, var.do=c('a','b'))
{
year.form <- format(dat2$date, '%Y')
dat$year <- factor(year.form, levels=unique(year.form), ordered=TRUE)
yl <- levels(dat$year)
ydat <- list()
for(i in 1:length(yl))
ydat[[i]] <- subset(dat,year==yl[i])
wdat <- list()
for(i in 1:length(yl)){  
a <- list()
for(j in 1:length(var.do)){
#par(mfrow=c(2,2))
a[[j]] <- do.call(my.ccf, list(ydat[[i]][,eval(c(var.do[j],'x'))]))
}
names(a) <- var.do
wdat[[i]] <- a         
}
names(wdat) <- yl
return(wdat)
}
date <- seq(as.Date("2011/1/1"), as.Date("2013/8/1"),  by = "months")
dat2 <- data.frame(date=date,x =rnorm(length(date)), y=rnorm(length(date)),        
a=rnorm(length(date)),b=rnorm(length(date)))
pdf("./temp/dat2.pdf")
A <- year.ccf1(dat2, var.do=c('a', 'b'))
dev.off()

我正在尝试将 dat2.pdf 中每个图的标题更改为带有列名的年份标记，例如“2011 年：a 和 x 的 CCF”。不确定我是否遗漏了一些明显的东西了解如何在一个 do.call 中进行辩论，如下所示。

do.call(my.ccf, main=paste("Year ", yl[i], var.do[j], " & x", sep=" "), list(ydat[[i]
[,eval(c(var.do[j],'x'))]))

还尝试将多个绘图放在一页中以使 pdf 文件更小。我尝试使用 par(mfrow) （未注释）来处理第二个问题，但失败了。

任何帮助深表感谢。

4

1 回答 1

1

我认为你应该首先给my.ccf. 例如：

my.ccf <- function(dat,year=NA)
{
   ...
  title=paste(colnames(dat)[1], "&", colnames(dat)[2], sep=" ")
  if(!is.na(year))
    title <- paste('Year', year,': CCF',title)
  ccf1 <- ccf(residuals(modx), ..., main=title)
  .....
}

然后你改变你do.call的东西是这样的：

 do.call(my.ccf, list(dat=list(ydat[[i]][,eval(c(var.do[j],'x'))]),year=i))

就是说认为您的year.ccf1功能过于复杂，而且效率极低。尝试使用outer和mapply例如重写它。

year.ccf1 <- function(dat, var.do=c('a','b'))
{
  dat$year <- factor(format(dat$date, '%Y'))
  ydat <- split(dat,dat$year)
  outer(names(ydat),var.do,function(x,y)
    Map(function(i,j)
        do.call(my.ccf, list(dat=ydat[[i]][,c(j,'x')],year=i))
           ,x,y)
  )
}
于 2013-08-02T09:32:03.033 回答