1

1 年のさまざまな季節 (寒い季節または暖かい季節にグループ化) の時系列データを同じプロットにプロットし、2 つの期間を識別したいと考えています。どうやってやるの？

再現可能なデータは次のとおりです。

library(gamair) 
data(chicago) 
chicago$date<-seq(from=as.Date("1987-01-01"), to=as.Date("2000-12-31"),length=5114)
data.cold <- subset(chicago, quarters(date) %in% c("Q1", "Q4"))
data.warm <- subset(chicago, quarters(date) %in% c("Q2", "Q3"))

以下のコードは、2 つの期間に対して 2 つの別々のプロットを作成します。

with(data.cold,plot(date,death,pch=".", ylab= expression("Mortality Count"), main = "Daily Mortality in cold season"))
with(data.warm ,plot(date,death,pch=".", ylab= expression("Mortality Count"), main = "Daily Mortality in warm season"))

次のコードは 1 つのプロットを作成しますが、2 つの期間の明確な境界はありません。私の願いは、1 つのプロットで 2 つの期間の識別可能なプロットを作成することです。

with(chicago,plot(date,death,pch=".", ylab= expression("Mortality Count"), main = "Daily Mortality from  1987-2000"))
4

1 に答える 1

1

見てください?points：

with(data.cold,plot(date,death,pch=".", ylab= expression("Mortality Count"), main = "Daily Mortality in cold season"))
with(data.warm ,points(date,death,pch=".", ylab= expression("Mortality Count"), main = "Daily Mortality in warm season", col="red"))
legend("topright", c("cold", "warm"), fill=c("black", "red"))

それはあなたのために働きますか？

于 2013-05-22T10:25:38.247 に答える