文章/答案/技术大牛

发布

社区首页 >问答首页 >使用lapply来标准化ts()函数

问使用lapply来标准化ts()函数
EN

Stack Overflow用户

提问于 2018-06-21 06:11:20

回答 1查看 222关注 0票数 0

我是R和编写函数的新手。我花了几个小时试图弄清楚这一点，并在谷歌上搜索，但似乎什么也找不到。希望你能帮上忙？我想使用lapply()通过ts()函数分析下面的数据。

我的代码如下所示：

library(dplyr)

#group out different sites
 mylist <- data %>%
   group_by(Site)
 mylist

#Write ts() function
 alpha_function = function(x) {
    ts_alpha = ts(x$Temperature, frequency=12, start=c(0017, 7, 20))
     return(data.frame(ts_alpha))    
 }

 #Run list through lapply()
 results = lapply(mylist, alpha_function())

但是我得到了这个错误:参数"x“没有缺省。

我有一个数据集，如下所示：

Site(factor)    Date(POSIXct)         Temperature(num)
1               0017-03-04            2.73
2               0017-03-04            3.73
3               0017-03-04            2.71
4               0017-03-04            2.22
5               0017-03-04            2.89
etc.

我有5个不同地点不同日期的3000多个温度读数。

提前感谢！

lapply

回答 1

Stack Overflow用户

回答已采纳

发布于 2018-06-21 21:54:46

使用dplyr和tidyverse时，建议的方法是将内容保存在数据帧中：

library(tidyverse)
library(zoo)

dat %>% 
  nest(-Site) %>% 
  mutate(data = map(data, ~ zoo(.x$Temperature, .x$Date)))

# # A tibble: 5 x 2
#   Site  data     
#   <fct> <list>   
# 1 a     <S3: zoo>
# 2 b     <S3: zoo>
# 3 c     <S3: zoo>
# 4 d     <S3: zoo>
# 5 e     <S3: zoo>

或者，如果我们必须使用ts而不是zoo对象，我们可以使用as.ts(zoo(...))。

如果我们仍然喜欢使用常规列表，我们可以使用基本split()和lapply()

dat %>% 
  split(.$Site) %>% 
  lapply(function(.x) zoo(.x$Temperature, .x$Date))

# List of 5
#  $ a:‘zoo’ series from 2017-03-04 12:00:00 to 2017-05-06 00:30:00
#   Data: num [1:3000] 5.37 5.49 5.32 5.44 5.43 ...
#   Index:  POSIXct[1:3000], format: "2017-03-04 12:00:00" ...
#  $ b:‘zoo’ series from 2017-03-04 12:00:00 to 2017-05-06 00:30:00
#   Data: num [1:3000] 5.36 5.22 5.15 5.41 5.41 ...
#   Index:  POSIXct[1:3000], format: "2017-03-04 12:00:00" ...
#  $ c:‘zoo’ series from 2017-03-04 12:00:00 to 2017-05-06 00:30:00
#   Data: num [1:3000] 6.08 6.11 6.22 6.13 6.03 ...
#   Index:  POSIXct[1:3000], format: "2017-03-04 12:00:00" ...
#  $ d:‘zoo’ series from 2017-03-04 12:00:00 to 2017-05-06 00:30:00
#   Data: num [1:3000] 5.06 4.96 5.23 5.16 5.29 ...
#   Index:  POSIXct[1:3000], format: "2017-03-04 12:00:00" ...
#  $ e:‘zoo’ series from 2017-03-04 12:00:00 to 2017-05-06 00:30:00
#   Data: num [1:3000] 5.1 5.08 5.14 5.13 5.22 ...
#   Index:  POSIXct[1:3000], format: "2017-03-04 12:00:00" ...

(其中，dat的生成方式如下：

n_sites <- 5
n_dates <- 3000

set.seed(123) ; dat <- tibble(
  Site = factor(rep(letters[1:n_sites], each = n_dates)), 
  Date = rep(seq.POSIXt(as.POSIXct("2017-03-04 12:00:00"), by = "30 min", length.out = n_dates), times = n_sites),
  Temperature = as.vector(replicate(n_sites, runif(1, 5, 6) + cumsum(rnorm(n_dates, 0, 0.1))))
)

票数 1

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/50957786

复制

相似问题

问使用lapply来标准化ts()函数
EN

回答 1

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问使用lapply来标准化ts()函数EN

回答 1

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问使用lapply来标准化ts()函数
EN