我的代码旨在生成时间序列分解的多个输出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个图表,即第一个小时级别,不再有。 我的问题:如何更改代码以使其生成一天中的小时(小时变量)的输出图? 我几个小时前发布了一个类似的问题,代码完全不同,我无法使代码工作,因为加载ggplot2的问题 谢谢
答案 0 :(得分:0)
请注意,您在函数调用中使用unique(Dataset$Hour)
而不是d
。