0

我的代码旨在生成时间序列分解的多个输出 jpeg 图（一天中的每个小时）。我运行代码，但没有收到错误消息。事实上，我收到一些迹象表明产出已经产生。

这是创建的函数（stackoverflow 专家给我的）：

HR.plot <- function(y,z){
    ARMAXpath<-file.path("C:", "Desktop","RPLOTS", paste0("HRplot_", z, ".jpg"))
    jpeg(file = ARMAXpath)
    plot(stl(y,s.window="periodic"))
    dev.off()
}

我在 by() 语句中使用此代码：

by(Dataset, Dataset$Hour, function(d) HR.plot(Dataset$UQ,unique(Dataset$Hour)))

它返回 1 个图表，第一个 Hour 级别，仅此而已。我的问题：如何更改代码以使其生成一天中每个小时的输出图（小时变量）？几个小时前，我用完全不同的代码发布了一个类似的问题，但由于加载 ggplot2 的问题，我无法让该代码工作，谢谢

4

1 回答 1

0

请注意，您unique(Dataset$Hour)在函数调用中使用而不是使用d.

于 2013-09-08T19:25:39.690 回答