我正在尝试使用R.中的单个或最小命令创建多个图[6个图或多个图]
让我们先看一下每小时的数据:-
str(ZZZ)
'data.frame': 291960 obs. of 9 variables:
$ TRADE_DT : POSIXct, format: "2014-11-01" "2014-11-01" "2014-11-01" "2014-11-01" ...
$ YEAR : int 2014 2014 2014 2014 2014 2014 2014 2014 2014 2014 ...
$ MONTH : int 11 11 11 11 11 11 11 11 11 11 ...
$ hour_num : int 1 1 1 1 1 1 1 1 1 1 ...
$ source : Factor w/ 5 levels "AB","EF","EI",..: 1 1 1 1 1 1 1 1 1 1 ...
$ LSE_CD : int 116 116 116 116 116 116 135 135 135 135 ...
$ utility_cd: Factor w/ 6 levels "CPL","SHARY",..: 1 2 3 4 5 6 1 4 5 6 ...
$ load : num 12.834 0.502 31.436 13.948 31.314 ...
$ total_load: num 13.929 0.524 35.864 14.77 33.161 ...
dput(头(ZZZ))
structure(list(TRADE_DT = structure(c(1414818000, 1414818000,
1414818000, 1414818000, 1414818000, 1414818000), class = c("POSIXct",
"POSIXt"), tzone = ""), YEAR = c(2014L, 2014L, 2014L, 2014L,
2014L, 2014L), MONTH = c(11L, 11L, 11L, 11L, 11L, 11L), hour_num = c(1L,
1L, 1L, 1L, 1L, 1L), source = structure(c(1L, 1L, 1L, 1L, 1L,
1L), .Label = c("AB", "EF", "EI", "IB", "ST"), class = "factor"),
LSE_CD = c(116L, 116L, 116L, 116L, 116L, 116L), utility_cd = structure(1:6, .Label = c("CPL",
"SHARY", "TNMP", "TXRL", "TXTU", "WTU"), class = "factor"),
load = c(12.83423, 0.501589, 31.435567, 13.947688, 31.314148,
2.237439), total_load = c(13.928702, 0.524432, 35.864181,
14.770245, 33.161105, 2.417721)), .Names = c("TRADE_DT",
"YEAR", "MONTH", "hour_num", "source", "LSE_CD", "utility_cd",
"load", "total_load"), row.names = c(NA, 6L), class = "data.frame")
我有兴趣根据每个实用程序覆盖我的来源(AB、EI、EF等)。对于6个实用程序,它应该生成6个图,其中每个图将有5行(或者根据需要有2或3行)每个效用有一个图&每个图都应该有多条基于源的线听起来很简单,但当数据是这种格式时,我还没能做到这一点。
我能够在图表中叠加多条线
然而,当我把我的来源(5个因素)变成5个不同的列,并从图片中删除小时时,我能够做到这一点;每天汇总。
str(YYY)
'data.frame': 102 obs. of 5 variables:
$ TRADE_DT: POSIXct, format: "2014-01-01" "2014-01-02" "2014-01-03" ...
$ AB : num 289 336 356 258 316 ...
$ EI : num 306 347 370 282 335 ...
$ IB : num 282 325 299 250 307 ...
$ EF : num 304 348 367 281 335 ...
ggplot(YYY, aes(TRADE_DT)) +
geom_line(aes(y = AB, colour = "AB")) +
geom_line(aes(y = EI, colour = "EI")) +
geom_line(aes(y = IB, colour = "IB")) +
geom_line(aes(y = EF, colour = "EF"))
但是,上面的方法并没有像我想要的那样通过utility_cd或LSE_cd来分离图,而且我不得不去掉hour。我见过人们使用";通过命令";在SAS中一次创建这些多个图。
有一个神奇的";命令";对于这种类型的交易,在R中?我将把我所有的图形输出到一个大pdf中,我可以自己处理。
如果有人能分享用这些标准生成多个图的秘密,我会非常感激。此外,当我有24小时的数据时,这条线看起来不像线,它们看起来像是通过一条倾斜的水平线连接在一起的。
再次感谢!
最佳,Gyve
请提供dput(head(YOUR DATA SET))
而不是str
,因为str
不太适合复制您的数据。如何制作一个伟大的R可复制的例子?
希望这能有所帮助:
1.绘制数据对于ggplot
,您需要重新整形2包装的措辞中的熔融数据集
获取第二个数据集:
YYY <- data.frame(TRADE_DT = seq(as.Date("2014-01-01"),as.Date("2014-01-05"), length.out = 5),
AB = c(289,336,356,258,316),
EI = c(306,347,370,282,335),
IB = c(282,325,299,250,307),
EF = c(304,348,367,281,335))
现在我们使用melt
来根据我们的需求进行调整:
require(reshape2)
YYY_molten <- melt(YYY,"TRADE_DT")
> head(YYY_molten)
TRADE_DT variable value
1 2014-01-01 AB 289
2 2014-01-02 AB 336
3 2014-01-03 AB 356
4 2014-01-04 AB 258
5 2014-01-05 AB 316
6 2014-01-01 EI 306
现在您可以使用ggplot
require(ggplot2)
ggplot(YYY_molten, aes(x = TRADE_DT, y = value, col = variable)) + geom_line()
2.按实用程序绘图
假设_utility_cd_是包含实用程序数据的列,您可以执行以下操作:
ZZZ_split <- split(ZZZ, f = ZZZ$utility_cd)
lapply(ZZZ_split, function(subset){
# function that melts and plots your subset/utility
})
如果我正确解读你的str,它应该是:
lapply(ZZZ_split, function(subset){
print(ggplot(subset, aes(x=TRADE_DT, y=LSE_CD, col = source)) + geom_line())
})