另一个我无法解决的棘手问题。
我有两个数据帧(“阳光”和“霜”)。一个是年累积日照时数,另一个是冬季极端温度(12月-2月)。我现在要做的是,检查这些data.frames中是否有Frost。在这种情况下,如果温度低于零下19度,我们就会说到Frost。然而,如果一年中日照时数超过300小时,我们就会说是零下21度的霜。然而,另一个前提条件是,必须检查数据帧Frost中的年份,以确定当前年份的结束和下一年的开始(12月至2月)。这意味着,如果1962年日照时数超过300hs,我需要检查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年日照超过300hs,所以气温上升到-21°。还有1964年的霜,因为1963年日照超过300小时。
1963-01-17 -21.47
1964-01-31 -20.54有谁有办法解决这个问题吗?
发布于 2012-12-04 04:06:15
采纳了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.47https://stackoverflow.com/questions/13689356
复制相似问题