【问题标题】:Writing a Function to Lapply Over Large List编写一个函数来覆盖大列表
【发布时间】:2014-02-17 09:32:15
【问题描述】:

基本上我有几天的大型时间序列数据框。我已经编写了一些代码,它在数据框中一次只运行一天,但现在我想对其进行调整,以便它运行一整天。对于我的数据框中的每一天,都有一列包含当天的日出时间和一列包含日落时间。我想利用这些时间将每一天分成白天和晚上。我的日出日落柱子看起来像这样,而且每天都不一样:

Sunrise              Sunset
2010-01-19 08:55:12 2010-01-19 17:26:34

我使用 split 将数据框按日期划分,以给出一个包含 10 个元素(天)的大列表

# Splits data frame by date 
sepdays<- split(df, df$Date)

# Function to split each day into day and night hours
daynight <- function(){
rise <- as.character(df$Sunrise[1])
rise <- substr(rise, 12,19)
set <- as.character(df$Sunset[1])
set <- substr(set, 12,19)
day <- df[df$Time>rise & df$Time<set, ]
df.night1<-df[df$Time<rise,]
df.night2<-df[df$Time>set,]
night <- merge.data.frame(df.night1,df.night2, sort = TRUE, all.x = TRUE, all.y=TRUE)
return(table(day$Activity))
}

# Apply function over list of days
lapply(sepdays,daynight)

当我运行 lapply 时,出现未使用的参数错误:

FUN(X[[1L]], ...) 中的错误:未使用的参数 (X[[1]])

我也不确定这是否是获取每个匹配日的特定日出/日落时间的最佳方式。我意识到我的函数没有参数,但我是 R 新手,所以不太确定自己在做什么。

这是我的数据的样子。

Date       Time      Activity  Sunrise              Sunset
2010-01-19 23:58:00  1         2010-01-19 08:55:12  2010-01-19 17:26:34
2010-01-19 23:59:00  1         2010-01-19 08:55:12  2010-01-19 17:26:34
2010-01-19 00:00:00  0         2010-01-20 08:54:13  2010-01-20 17:28:11
2010-01-19 00:01:00  0         2010-01-20 08:54:13  2010-01-20 17:28:11
2010-01-20 00:02:00  1         2010-01-20 08:54:13  2010-01-20 17:28:11
2010-01-20 00:03:00  0         2010-01-20 08:54:13  2010-01-20 17:28:11
2010-01-20 00:04:00  1         2010-01-20 08:54:13  2010-01-20 17:28:11

我希望我的输出包含每个日期的活动表,例如:

2010-01-19
1  0
2  0

2010-01-20
1  0
2  3

【问题讨论】:

  • 您正在尝试重新发明轮子。使用 R 的日期时间工具和包,如 plyr、dplyr 或 data.table。如果您清楚地说明您的最终目标并在您的问题中添加示例数据(例如,使用dput(head(df))),那么应该会有人向您展示如何轻松做到这一点。
  • 您的示例输出表与示例输入有什么关系?
  • 每天都应该生成一个活动值表。 IE。每天的 1 总数和 0 总数。
  • 错误信息告诉你函数daynight需要有参数。 daynight &lt;- function() 应该是 daynight &lt;- function(df)
  • 我看到两个 1 和两个 0 用于 2010-01-19。与日出/日落时间的关系在哪里?

标签: r dataframe time-series lapply


【解决方案1】:

我不太确定,因为你的问题很模糊,但我认为你可以这样做:

DF <- read.table(text="Date,       Time,      Activity,  Sunrise,              Sunset
2010-01-19, 23:58:00,  1,         2010-01-19 08:55:12,  2010-01-19 17:26:34
2010-01-19, 23:59:00,  1,         2010-01-19 08:55:12,  2010-01-19 17:26:34
2010-01-19, 00:00:00,  0,         2010-01-19 08:55:12,  2010-01-19 17:26:34
2010-01-19, 00:01:00,  0,         2010-01-19 08:55:12,  2010-01-19 17:26:34
2010-01-19, 09:01:00,  0,         2010-01-19 08:55:12,  2010-01-19 17:26:34
2010-01-20, 00:02:00,  1,         2010-01-20 08:54:13,  2010-01-20 17:28:11
2010-01-20, 00:03:00,  0,         2010-01-20 08:54:13,  2010-01-20 17:28:11
2010-01-20, 00:04:00,  1,         2010-01-20 08:54:13,  2010-01-20 17:28:11", header=TRUE, sep=",")

DF$datetime <- as.POSIXct(paste(DF$Date, DF$Time), "%Y-%m-%d %H:%M:%S", tz="GMT")
DF$date <- as.Date(DF$datetime)
DF$Sunrise <- as.POSIXct(DF$Sunrise, "%Y-%m-%d %H:%M:%S", tz="GMT")
DF$Sunset <- as.POSIXct(DF$Sunset, "%Y-%m-%d %H:%M:%S", tz="GMT")


DF$day <- (DF$datetime > DF$Sunrise) & (DF$datetime < DF$Sunset)

#        Date      Time Activity             Sunrise              Sunset            datetime   day       date
#1 2010-01-19  23:58:00        1 2010-01-19 08:55:12 2010-01-19 17:26:34 2010-01-19 23:58:00 FALSE 2010-01-19
#2 2010-01-19  23:59:00        1 2010-01-19 08:55:12 2010-01-19 17:26:34 2010-01-19 23:59:00 FALSE 2010-01-19
#3 2010-01-19  00:00:00        0 2010-01-19 08:55:12 2010-01-19 17:26:34 2010-01-19 00:00:00 FALSE 2010-01-19
#4 2010-01-19  00:01:00        0 2010-01-19 08:55:12 2010-01-19 17:26:34 2010-01-19 00:01:00 FALSE 2010-01-19
#5 2010-01-19  09:01:00        0 2010-01-19 08:55:12 2010-01-19 17:26:34 2010-01-19 09:01:00  TRUE 2010-01-19
#6 2010-01-20  00:02:00        1 2010-01-20 08:54:13 2010-01-20 17:28:11 2010-01-20 00:02:00 FALSE 2010-01-20
#7 2010-01-20  00:03:00        0 2010-01-20 08:54:13 2010-01-20 17:28:11 2010-01-20 00:03:00 FALSE 2010-01-20
#8 2010-01-20  00:04:00        1 2010-01-20 08:54:13 2010-01-20 17:28:11 2010-01-20 00:04:00 FALSE 2010-01-20

table(DF[,c("date", "Activity", "day")])

#, , day = FALSE
#
#            Activity
#date         0 1
#  2010-01-19 2 2
#  2010-01-20 1 2
#
#, , day = TRUE
#
#            Activity
#date         0 1
#  2010-01-19 1 0
#  2010-01-20 0 0

这样更容易阅读,效率更高。

【讨论】:

    猜你喜欢
    • 2022-11-23
    • 2022-11-13
    • 2014-11-16
    • 1970-01-01
    • 1970-01-01
    • 2014-09-08
    • 1970-01-01
    • 1970-01-01
    • 2018-12-29
    相关资源
    最近更新 更多