另一个我无法解决的棘手问题。
我有两个数据帧(“阳光”和“霜冻”)。一个是每年累积的日照时数,另一个是冬季极端温度(12月 - 2月)。我现在需要做的是,检查Frost的这些data.frames。 如果温度低于-19度,我们在这种情况下谈论弗罗斯特。但是,如果一年中的日照时数超过300小时,我们会谈到-21度的冰霜。但是,另一个先决条件是,必须检查数据帧Frost中的年份是否在当前和下一年的开始(12月到2月)结束。这意味着,如果1962年的日照时数超过300小时,我需要在1962年12月至1963年2月检查冰霜。 这是data.frames。
"Sunshine"
year hours
1962 396
1963 290
1964 301
"Frost"
Date temperature
1962-12-26 -19.21
1963-01-13 -20.29
1963-01-17 -21.47
1964-01-31 -20.54
现在我遇到的问题是如何将它们与数据帧结合起来,这样它只向我展示了Frost温度的日期,满足条件。
在这种情况下,我们没有在1962年冬季以19,21°进行冰霜事件,因为在一年中阳光超过300小时,触发温度为-21°。 1963年1月我们举行了冰霜活动,因为1962年阳光超过300小时,所以情况升至-21°。同样是弗罗斯特在1964年,因为1963年阳光超过300小时。
1963-01-17 -21.47
1964-01-31 -20.54
有没有人知道如何解决这个问题。
答案 0 :(得分:0)
接受Justin的建议,并假设Frost$Date
有课程Date
:
x <- Frost
x$year <- as.numeric(format(x$Date, "%Y"))
x <- merge(x, Sunshine)
x <- transform(x, threshold = ifelse(hours > 300, -21, -19))
subset(x, temperature < threshold)[,2:3] # Noting that column 1 is year
# Date temperature
#2 1963-01-13 -20.29
#3 1963-01-17 -21.47