如何在R中繪製覆蓋時間序列？

例如，具有兩個數據集：

數據集1（日期，值）和 dataset2（日期，值）

爲了繪製剛第一，我們可以執行以下操作：

x.Date <- as.Date(dataset1$Date) 
    x <- zoo(dataset1$Value, x.Date) 
    plot(x)

我想添加到同一個窗口（dataset2 $ value，dataset2 $ Date），並偶然設置不同的顏色。

values dataset1 $ Date和dataset2 $ Date不一定相同（有些日子可能重疊，有些日子可能重疊），例如dataset1 $ Date可能包含（dec01，dec02，dec03，dec05）和dataset2 $ Date（dec02 ，dec03，dec06）。

有誰知道如何在同一窗口中繪製兩個（或幾個）時間圖嗎？

來源

2012-09-13 Niko Gamulin

有幾種選擇。以下是與zoo對象配合使用的三個選項。

set.seed(1) 
xz = zoo(ts(rnorm(20), frequency = 4, start = c(1959, 2))) 
yz = zoo(ts(rnorm(20), frequency = 4, start = c(1959, 2))) 
# Basic approach 
plot(xz) 
lines(yz, col = "red") 
# Panels 
plot.zoo(cbind(xz, yz)) 
# Overplotted 
plot.zoo(cbind(xz, yz), 
     plot.type = "single", 
     col = c("red", "blue"))

如果你正在策劃定期ts對象，你還可以遊覽ts.plot：

有關使用GGPLOT2

set.seed(1) 
x = ts(rnorm(20), frequency = 4, start = c(1959, 2)) 
y = ts(rnorm(20), frequency = 4, start = c(1959, 2)) 
ts.plot(x, y, gpars = list(col = c("black", "red")))

來源

2012-09-13 09:05:31 A5C1D2H2I1M1N2O1R2T1

是的！這是我想要的。謝謝！ –

@NikoGamulin，我用一些更好的選項更新了我的答案。 – A5C1D2H2I1M1N2O1R2T1

請注意，如果'x'是一個動物園對象，'plot（x）'和'plot.zoo（x）'具有相同的效果，因爲這將導致'plot'派發'plot.zoo'。請注意，'xyplot'也有一個動物園方法。 –

什麼？

library(ggplot2) 

dta <- data.frame(year = 2012, 
    month = rep(seq(1,9),each=10), 
    day = sample(seq(1,20),90,replace=T), 
    Group = sample(c('A','B'),90,replace=T), 
    Value = seq(1,90) 
) 

dta$Date <- apply(dta[,c('year','month','day')],1,paste,sep='',collapse='-') 
dta$Date <- as.Date(dta$Date) 

qplot(Date,Value,data=dta,geom='line',color=Group)

在你所描述的情況

，我會以下列方式

dataset1$Group <- 'dataset1' 
dataset2$Group <- 'dataset2' 
dta <- rbind(dataset1,dataset2) 
dta$Date <- as.Date(dta$Date)

來源

2012-09-13 09:21:42

我手頭相同的任務組裝dta和一些研究之後，我遇到ts.plot來到{}統計功能這非常有幫助。

函數的用法如下：

ts.plot(..., gpars = list())

gpars是圖形參數，你可以指定情節的圖形組件。

我不得不與此類似，並存儲在一個變量調用時數據：

  [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] 
    V3  1951 1100 433 5638 1760 2385 2602 11007 2490 421 
    V5  433 880 216 4988 220 8241 13229 18704 6289 421 
    V7  4001 440 433 3686 880 9976 12795 21036 13229 1263 
    V9  2385 1320 650 8241 440 12795 13229 19518 11711 1474 
    V11 4771 880 1084 6723 0 17783 17566 27326 11060 210 
    V13 6940 880 2168 2602 1320 21036 16265 10843 15831 1474 
    V15 3903 1760 1951 3470 0 18217 14964  0 13663 2465 
    V17 4771 440 2819 8458 880 25591 24940 1518 17783 1895 
    V19 7807 1760 5205 2385 0 14096 22771 13880 12578 1263 
    V21 5205 880 5205 6506 880 28410 18217 13229 19952 1474 
    V23 6506 1760 5638 7590 880 14747 26675 11928 12795 1474 
    V25 7373 440 5855 10626 0 19301 21470 15398 19952 1895 
    V27 5638 2640 6289  0 880 16482 20603 30796 14313 2316 
    V29 8241 440 6506 6723 880 11277 35784 25157 23205 4423 
    V31 7373 2640 6072 8891 220 17133 27109 31013 27287 4001 
    V33 6723 660 5855 14313 660 6940 26892 17566 24111 4844 
    V35 9325 2420 9325 12578 0 6506 30796 34483 23422 5476 
    V37 4771 440 6872 12361 880 9325 36218 25808 30362 4844 
    V39 9976 2640 7658 12361 440 11277 36001 31013 40555 4633 
    V41 10410 880 6506 12795 440 26241 33398 27976 24940 5686 
    V43 5638 2200 7590 14313 0 9976 34483 29928 33832 6108 
    V45 10843 440 8675 11711 440 7807 29278 24940 43375 4633 
    V47 8675 1760 8891 13663 0 9108 38386 31230 33398 4633 
    V49 10410 1760 9542 13880 440 8675 39051 31446 42507 5476 
     . . . . . . . . .

而且我不得不對同一地塊的各列的時間序列圖。代碼如下：

ts.plot(time,gpars= list(col=rainbow(10)))

該圖的結果如下所示。 Multiple time series in a single plot

來源

2015-07-16 13:43:30 kirancodify

如何在R中繪製覆蓋時間序列？

回答

相關問題